High-Density Genetic Map Construction and QTL Mapping of Leaf and Needling Traits in Ziziphus jujuba Mill

The Chinese jujube (Ziziphus jujuba Mill., 2n = 2x = 24), one of the most popular fruit trees in Asia, is widely cultivated and utilized in China, where it is traditionally consumed as both a fresh and dried food resource. A high-density genetic map can provide the necessary framework for quantitative trait loci (QTL) analyses and map-based gene cloning and molecular breeding. In this study, we constructed a new high-density genetic linkage map via a genotyping-by-sequencing approach. For the consensus linkage map, a total of 3,792 markers spanning 2,167.5 cM were mapped onto 12 linkage groups, with an average marker interval distance of 0.358 cM. The genetic map anchored 301 Mb (85.7%) of scaffolds from the sequenced Z. jujuba "Junzao" genome. Based on this genetic map, 30 potential QTLs were detected, including 27 QTLs for leaf traits and 3 QTLs for needling length. This high-density genetic map and the identified QTLs for relevant agronomic traits lay the groundwork for functional genetic mapping, map-based cloning, and marker-assisted selection in jujube.

Genome Improvement and Genetic Map Construction for Aethionema arabicum, the First Divergent Branch in the Brassicaceae Family

The genus Aethionema is a sister-group to the core-group of the Brassicaceae family that includes Arabidopsis thaliana and the Brassica crops. Thus, Aethionema is phylogenetically well-placed for the investigation and understanding of genome and trait evolution across the family. We aimed to improve the quality of the reference genome draft version of the annual species Aethionema arabicum. Second, we constructed the first Ae. arabicum genetic map. The improved reference genome and genetic map enabled the development of each other. We started with the initially published genome (version 2.5). PacBio and MinION sequencing together with genetic map v2.5 were incorporated to produce the new reference genome v3.0. The improved genome contains 203 MB of sequence, with approximately 94% of the assembly made up of called (non-gap) bases, assembled into 2,883 scaffolds (with only 6% of the genome made up of non-called bases (Ns)). The N₅₀ (10.3 MB) represents an 80-fold increase over the initial genome release. We generated a Recombinant Inbred Line (RIL) population that was derived from two ecotypes: Cyprus and Turkey (the reference genotype. Using a Genotyping by Sequencing (GBS) approach, we generated a high-density genetic map with 749 (v2.5) and then 632 SNPs (v3.0) was generated. The genetic map and reference genome were integrated, thus greatly improving the scaffolding of the reference genome into 11 linkage groups. We show that long-read sequencing data and genetics are complementary, resulting in an improved genome assembly in Ae. arabicum. They will facilitate comparative genetic mapping work for the Brassicaceae family and are also valuable resources to investigate wide range of life history traits in Aethionema.

An Ultra High-Density Arabidopsis thaliana Crossover Map That Refines the Influences of Structural Variation and Epigenetic Features

Many environmental, genetic, and epigenetic factors are known to affect the frequency and positioning of meiotic crossovers (COs). Suppression of COs by large, cytologically visible inversions and translocations has long been recognized, but relatively little is known about how smaller structural variants (SVs) affect COs. To examine fine-scale determinants of the CO landscape, including SVs, we used a rapid, cost-effective method for high-throughput sequencing to generate a precise map of >17,000 COs between the Col-0 and Ler-0 accessions of Arabidopsis thaliana COs were generally suppressed in regions with SVs, but this effect did not depend on the size of the variant region, and was only marginally affected by the variant type. CO suppression did not extend far beyond the SV borders and CO rates were slightly elevated in the flanking regions. Disease resistance gene clusters, which often exist as SVs, exhibited high CO rates at some loci, but there was a tendency toward depressed CO rates at loci where large structural differences exist between the two parents. Our high-density map also revealed in fine detail how CO positioning relates to genetic (DNA motifs) and epigenetic (chromatin structure) features of the genome. We conclude that suppression of COs occurs over a narrow region spanning large- and small-scale SVs, representing an influence on the CO landscape in addition to sequence and epigenetic variation along chromosomes.

An SNP-Based Genetic Map and QTL Mapping for Growth Traits in the Red-Spotted Grouper (Epinephelus akaara)

The red-spotted grouper (Epinephelus akaara) is one of the most commercially important aquatic species in China. However, its seedstock has low larval survival rates, and its stability is confronted with the danger of overexploitation. In this study, a high-density genetic map was constructed using 3435 single nucleotide polymorphisms (SNPs) from 142 first generation (F₁) full-sib offspring and two parents of a red-spotted grouper population. The total genetic length of the map was 2300.12 cM with an average intermarker distance of 0.67 cM. Seventeen genome-wide significant quantitative trait loci (QTLs) for growth-related traits were detected on 24 linkage groups, including 5 QTLs for full length, 7 QTLs for body length, and 5 QTLs for body weight. The contribution values of explained phenotypic variance ranged from 10.7% to 12.9%. Moreover, 13 potential candidate genes for growth-related traits were identified. Collectively, these findings will be useful for conducting marker-assisted selection of the red-spotted grouper in future studies.

Genetic determination of sex and shell color in the Pacific abalone Haliotis discus hannai revealed by an integrated linkage map

Integrated linkage maps for each sex have been constructed for the Pacific abalone Haliotis discus hannai using three F₁ mapping families based on co-dominant markers. A total of 273 markers were placed on the female map, spanning 927.3 cM with an average interval of 3.64 cM, whereas 277 markers were mapped on the male map, covering 727.0 cM with an average spacing of 2.80 cM. Both female and male maps consisted of 18 linkage groups, corresponding well with the number of chromosomes. Furthermore, the sex-determining locus and the green/orange shell color controlling locus were mapped to the linkage group 3 (LG3) and LG9 respectively. A marker completely linked to phenotypic sex was identified, and the sex determination system was further concluded as paternal heterogametic (males XY and females XX). Based on the segregation ratio of the shell color in the progeny, a simple recessive model of epistasis was proposed to explain the distribution of different color morphs (green, orange and blue): the recessive allele determining orange type masks the effect of the locus controlling green and blue types, whereas the dominant allele at the green/orange locus permits the expression of green and blue types controlled by another locus. The current consensus map provides a useful framework for genetic studies in the Pacific abalone. Mapping of the sex-determining locus and the shell color-controlling locus leads to further understanding of the mechanisms underlying these important traits.

A SNP-Based High-Density Genetic Map of Leaf and Fruit Related Quantitative Trait Loci in Wolfberry (Lycium Linn.)

Wolfberry (Lycium Linn. 2n = 24) fruit, Gouqizi, is a perennial shrub, traditional food and medicinal plant resource in China. Leaf and fruit related characteristics are economically important traits that are the focus for genetic improvement, but few studies into the molecular genetics of this crop have been reported to date. Here, an F₁ population (302 individuals) derived from a cross between "NO.1 Ningqi" (Lycium barbarum L.) and "Chinese gouqi" (Lycium chinese Mill.) was constructed. We recorded fruit weight, longitude, diameter and index along with leaf length, width and index for three consecutive years from 2015 to 2017. Based on this population and these phenotypic data, we constructed the first high-density genetic map of Lycium using specific length amplified fragment sequencing (SLAF-seq) and analyzed quantitative trait loci (QTLs). The map contains 6733 single nucleotide polymorphisms and 12 linkage groups (LG) with a total map distance of 1702.45 cM and an average map distance of 0.253 cM. A total of 55 QTLs were mapped for more than 2 years, of which 18 stable QTLs for fruit index on LG 11, spanning an interval of 73.492-90.945 cM, were detected. qFI11-15 for fruit index was an impressive QTL with logarithm of odds (LOD) and proportion of variance explained (PEV) values reaching 11.07 and 19.7%, respectively. The QTLs on LG 11 were gathered tightly, having an average interval of less than 1 cM per QTL, suggesting that there might be a cluster region controlling fruit index. Remarkably, qLI10-2 and qLI11-2 for leaf index were detectable for 3 years. These results give novel insight into the genetic control of leaf and fruit related traits in Lycium and provide robust support for undertaking further positional cloning studies and implementing marker-assisted selection in seedlings.

Development of an integrated 200K SNP genotyping array and application for genetic mapping, genome assembly improvement and genome wide association studies in pear (Pyrus)

Pear (Pyrus; 2n = 34), the third most important temperate fruit crop, has great nutritional and economic value. Despite the availability of many genomic resources in pear, it is challenging to genotype novel germplasm resources and breeding progeny in a timely and cost-effective manner. Genotyping arrays can provide fast, efficient and high-throughput genetic characterization of diverse germplasm, genetic mapping and breeding populations. We present here 200K AXIOM^® PyrSNP, a large-scale single nucleotide polymorphism (SNP) genotyping array to facilitate genotyping of Pyrus species. A diverse panel of 113 re-sequenced pear genotypes was used to discover SNPs to promote increased adoption of the array. A set of 188 diverse accessions and an F₁ population of 98 individuals from 'Cuiguan' × 'Starkrimson' was genotyped with the array to assess its effectiveness. A large majority of SNPs (166 335 or 83%) are of high quality. The high density and uniform distribution of the array SNPs facilitated prediction of centromeric regions on 17 pear chromosomes, and significantly improved the genome assembly from 75.5% to 81.4% based on genetic mapping. Identification of a gene associated with flowering time and candidate genes linked to size of fruit core via genome wide association studies showed the usefulness of the array in pear genetic research. The newly developed high-density SNP array presents an important tool for rapid and high-throughput genotyping in pear for genetic map construction, QTL identification and genomic selection.

Genetic Mapping of Prince Rupprecht's Larch (Larix principis-rupprechtii Mayr) by Specific-Locus Amplified Fragment Sequencing

A high-density genetic linkage map is essential for plant genetics and genomics research. However, due to the deficiency of genomic data and high-quality molecular markers, no genetic map has been published for Prince Rupprecht’s larch (Larix principis-rupprechtii Mayr), a conifer species with high ecological and commercial value in northern China. In this study, 145 F1 progeny individuals from an intraspecific cross between two elite clones of L. principis-rupprechtii and their parents were employed to construct the first genetic map in this important tree species using specific-locus amplified fragment sequencing (SLAF-seq). After preprocessing, the procedure yielded 300.20 Gb of raw data containing 1501.22 M pair-end reads. A total of 324,352 SNP markers were detected and 122,785 of them were polymorphic, with a polymorphism rate of 37.86%. Ultimately, 6099 SNPs were organized into a genetic map containing 12 linkage groups, consistent with the haploid chromosome number of larch and most other species in the Pinaceae family. The linkage map spanned 2415.58 cM and covered 99.6% of the L. principis-rupprechtii genome with an average of 0.4 cM between adjacent markers. To the best of our knowledge, this map is the first reference map for L. principis-rupprechtii, as well as the densest one obtained in larch species thus far. The genome-wide SNPs and the high-resolution genetic map will provide a foundation for future quantitative trait loci mapping, map-based cloning, marker-assisted selection, comparative genomics, and genome sequence assembly for larch trees.

Ecological speciation in sympatric palms: 3. Genetic map reveals genomic islands underlying species divergence in Howea

Although it is now widely accepted that speciation can occur in the face of continuous gene flow, with little or no spatial separation, the mechanisms and genomic architectures that permit such divergence are still debated. Here, we examined speciation in the face of gene flow in the Howea palms of Lord Howe Island, Australia. We built a genetic map using a novel method applicable to long-lived tree species, combining it with double digest restriction site-associated DNA sequencing of multiple individuals. Based upon various metrics, we detected 46 highly differentiated regions throughout the genome, four of which contained genes with functions that are particularly relevant to the speciation scenario for Howea, specifically salt and drought tolerance.

Development and evaluation of a genome-wide Coffee 8.5K SNP array and its application for high-density genetic mapping and for investigating the origin of Coffea arabica L

Coffee species such as Coffea canephora P. (Robusta) and C. arabica L. (Arabica) are important cash crops in tropical regions around the world. C. arabica is an allotetraploid (2n = 4x = 44) originating from a hybridization event of the two diploid species C. canephora and C. eugenioides (2n = 2x = 22). Interestingly, these progenitor species harbour a greater level of genetic variability and are an important source of genes to broaden the narrow Arabica genetic base. Here, we describe the development, evaluation and use of a single-nucleotide polymorphism (SNP) array for coffee trees. A total of 8580 unique and informative SNPs were selected from C. canephora and C. arabica sequencing data, with 40% of the SNP located in annotated genes. In particular, this array contains 227 markers associated to 149 genes and traits of agronomic importance. Among these, 7065 SNPs (~82.3%) were scorable and evenly distributed over the genome with a mean distance of 54.4 Kb between markers. With this array, we improved the Robusta high-density genetic map by adding 1307 SNP markers, whereas 945 SNPs were found segregating in the Arabica mapping progeny. A panel of C. canephora accessions was successfully discriminated and over 70% of the SNP markers were transferable across the three species. Furthermore, the canephora-derived subgenome of C. arabica was shown to be more closely related to C. canephora accessions from northern Uganda than to other current populations. These validated SNP markers and high-density genetic maps will be useful to molecular genetics and for innovative approaches in coffee breeding.

Rice Genome Resequencing Reveals a Major Quantitative Trait Locus for Resistance to Bakanae Disease Caused by Fusarium fujikuroi

Bakanae disease (BD), caused by the fungal pathogen Fusarium fujikuroi, has become a serious threat in rice-cultivating regions worldwide. In the present study, quantitative trait locus (QTL) mapping was performed using F2 and F3 plants derived after crossing a BD-resistant and a BD-susceptible Korean japonica rice variety, ‘Samgwang’ and ‘Junam’, respectively. Resequencing of ‘Junam’ and ‘Samgwang’ genomes revealed 151,916 DNA polymorphisms between the two varieties. After genotyping 188 F2 plants, we constructed a genetic map comprising 184 markers, including 175 kompetitive allele-specific PCR markers, eight cleaved amplified polymorphic sequence (CAPS) markers, and a derived CAPS (dCAPS) marker. The degree of BD susceptibility of each F2 plant was evaluated on the basis of the mortality rate measured with corresponding F3 progeny seedlings by in vitro screening. Consequently, qFfR9, a major QTL, was discovered at 30.1 centimorgan (cM) on chromosome 9 with a logarithm of the odds score of 60.3. For the QTL interval, 95% probability lay within a 7.24–7.56 Mbp interval. In this interval, we found that eight genes exhibited non-synonymous single nucleotide polymorphisms (SNPs) by comparing the ‘Junam’ and ‘Samgwang’ genome sequence data, and are possibly candidate genes for qFfR9; therefore, qFfR9 could be utilized as a valuable resource for breeding BD-resistant rice varieties.

An Ultra-Dense Haploid Genetic Map for Evaluating the Highly Fragmented Genome Assembly of Norway Spruce (Picea abies)

Norway spruce (Picea abies (L.) Karst.) is a conifer species of substanital economic and ecological importance. In common with most conifers, the P. abies genome is very large (∼20 Gbp) and contains a high fraction of repetitive DNA. The current P. abies genome assembly (v1.0) covers approximately 60% of the total genome size but is highly fragmented, consisting of >10 million scaffolds. The genome annotation contains 66,632 gene models that are at least partially validated (www.congenie.org), however, the fragmented nature of the assembly means that there is currently little information available on how these genes are physically distributed over the 12 P. abies chromosomes. By creating an ultra-dense genetic linkage map, we anchored and ordered scaffolds into linkage groups, which complements the fine-scale information available in assembly contigs. Our ultra-dense haploid consensus genetic map consists of 21,056 markers derived from 14,336 scaffolds that contain 17,079 gene models (25.6% of the validated gene models) that we have anchored to the 12 linkage groups. We used data from three independent component maps, as well as comparisons with previously published Picea maps to evaluate the accuracy and marker ordering of the linkage groups. We demonstrate that approximately 3.8% of the anchored scaffolds and 1.6% of the gene models covered by the consensus map have likely assembly errors as they contain genetic markers that map to different regions within or between linkage groups. We further evaluate the utility of the genetic map for the conifer research community by using an independent data set of unrelated individuals to assess genome-wide variation in genetic diversity using the genomic regions anchored to linkage groups. The results show that our map is sufficiently dense to enable detailed evolutionary analyses across the P. abies genome.

Chromosome-Wide Evolution and Sex Determination in the Three-Sexed Nematode Auanema rhodensis

Trioecy, a mating system in which males, females and hermaphrodites co-exist, is a useful system to investigate the origin and maintenance of alternative mating strategies. In the trioecious nematode Auanema rhodensis, males have one X chromosome (XO), whereas females and hermaphrodites have two (XX). The female vs. hermaphrodite sex determination mechanisms have remained elusive. In this study, RNA-seq analyses show a 20% difference between the L2 hermaphrodite and female gene expression profiles. RNAi experiments targeting the DM (doublesex/mab-3) domain transcription factor dmd-10/11 suggest that the hermaphrodite sexual fate requires the upregulation of this gene. The genetic linkage map (GLM) shows that there is chromosome-wide heterozygosity for the X chromosome in F2 hermaphrodite-derived lines originated from crosses between two parental inbred strains. These results confirm the lack of recombination of the X chromosome in hermaphrodites, as previously reported. We also describe conserved chromosome elements (Nigon elements), which have been mostly maintained throughout the evolution of Rhabditina nematodes. The seven-chromosome karyotype of A. rhodensis, instead of the typical six found in other rhabditine species, derives from fusion/rearrangements events involving three Nigon elements. The A. rhodensis X chromosome is the smallest and most polymorphic with the least proportion of conserved genes. This may reflect its atypical mode of father-to-son transmission and its lack of recombination in hermaphrodites and males. In conclusion, this study provides a framework for studying the evolution of chromosomes in rhabditine nematodes, as well as possible mechanisms for the sex determination in a three-sexed species.

Genetic Mapping of the Incompatibility Locus in Olive and Development of a Linked Sequence-Tagged Site Marker

The genetic control of self-incompatibility (SI) has been recently disclosed in olive. Inter-varietal crossing confirmed the presence of only two incompatibility groups (G1 and G2), suggesting a simple Mendelian inheritance of the trait. A double digest restriction associated DNA (ddRAD) sequencing of a biparental population segregating for incompatibility groups has been performed and high-density linkage maps were constructed in order to map the SI locus and identify gene candidates and linked markers. The progeny consisted of a full-sib family of 229 individuals derived from the cross 'Leccino' (G1) × 'Dolce Agogia' (G2) varieties, segregating 1:1 (G1:G2), in accordance with a diallelic self-incompatibility (DSI) model. A total of 16,743 single nucleotide polymorphisms was identified, 7,006 in the female parent 'Leccino' and 9,737 in the male parent 'Dolce Agogia.' Each parental map consisted of 23 linkage groups and showed an unusual large size (5,680 cM in 'Leccino' and 3,538 cM in 'Dolce Agogia'). Recombination was decreased across all linkage groups in pollen mother cells of 'Dolce Agogia,' the parent with higher heterozygosity, compared to megaspore mother cells of 'Leccino,' in a context of a species that showed exceptionally high recombination rates. A subset of 109 adult plants was assigned to either incompatibility group by a stigma test and the diallelic self-incompatibility (DSI) locus was mapped to an interval of 5.4 cM on linkage group 18. This region spanned a size of approximately 300 Kb in the olive genome assembly. We developed a sequence-tagged site marker in the DSI locus and identified five haplotypes in 57 cultivars with known incompatibility group assignment. A combination of two single-nucleotide polymorphisms (SNPs) was sufficient to predict G1 or G2 phenotypes in olive cultivars, enabling early marker-assisted selection of compatible genotypes and allowing for a rapid screening of inter-compatibility among cultivars in order to guarantee effective fertilization and increase olive production. The construction of high-density linkage maps has led to the development of the first functional marker in olive and provided positional candidate genes in the SI locus.

Recent polyploidization events in three Saccharum founding species

The complexity of polyploid Saccharum genomes hindered progress of genome research and crop improvement in sugarcane. To understand their genome structure, transcriptomes of 59 F₁ individuals derived from S. officinarumLA Purple and S. robustum Molokai 5829 (2n = 80, x = 10 for both) were sequenced, yielding 11 157 and 8998 SNPs and 83 and 105 linkage groups, respectively. Most markers in each linkage group aligned to single sorghum chromosome. However, 71 interchromosomal rearrangements were detected between sorghum and S. officinarum or S. robustum, and 24 (33.8%) of them were shared between S. officinarum and S. robustum, indicating their occurrence before the speciation event that separated these two species. More than 2000 gene pairs from S. spontaneum, S. officinarum and S. robustum were analysed to estimate their divergence time. Saccharum officinarum and S. robustum diverged about 385 thousand years ago, and the whole-genome duplication events occurred after the speciation event because of shared interchromosomal rearrangements. The ancestor of these two species diverged from S. spontaneum about 769 thousand years ago, and the reduction in basic chromosome number from 10 to 8 in S. spontaneum occurred after the speciation event but before the two rounds of whole-genome duplication. Our results proved that S. officinarum is a legitimate species in its own right and not a selection from S. robustum during the domestication process in the past 10 000 years. Our findings rejected a long-standing hypothesis and clarified the timing of speciation and whole-genome duplication events in Saccharum.

Construction of a High-Density Genetic Map and Identification of Loci Related to Hollow Stem Trait in Broccoli (Brassic oleracea L. italica)

A high-quality genetic map is important for mapping of compound traits. In this study, a genetic map was constructed based on the reference genome TO1000 after specific locus amplified fragment (SLAF) sequencing in a double-haploid segregation population of broccoli, and loci controlling hollow stem trait were identified in the genetic map. The genetic map contains 4,787 SLAF markers, with a mean marker distance of 0.22 cM and the mean sequencing depths of 91.14-fold in the maternal line, 88.97-fold in the paternal line and 17.11-fold in each DH progeny. A locus controlling the hollow stem trait, QHS.C09-2, which could explain 14.1% of the phenotypic variation, was steadily detected on the linkage group nine in the indicated data of 3 years' trials and BLUE analysis. The genetic map could lay an important foundation for mapping of compound traits, and mapping of hollow stem trait would be basis to clone the genes related to hollow stems in broccoli.

Genetic Variation for Protein Content and Yield-Related Traits in a Durum Population Derived From an Inter-Specific Cross Between Hexaploid and Tetraploid Wheat Cultivars

Wheat grain protein content (GPC) and yield components are complex quantitative traits influenced by a multi-factorial system consisting of both genetic and environmental factors. Although seed storage proteins represent less than 15% of mature kernels, they are crucial in determining end-use properties of wheat, as well as the nutritional value of derived products. Yield and GPC are negatively correlated, and this hampers breeding programs of commercially valuable wheat varieties. The goal of this work was the evaluation of genetic variability for quantity and composition of seed storage proteins, together with yield components [grain yield per spike (GYS) and thousand-kernel weight (TKW)] in a durum wheat population obtained by an inter-specific cross between a common wheat accession and the durum cv. Saragolla. Quantitative trait loci (QTL) analysis was conducted and closely associated markers identified on a genetic map composed of 4,366 SNP markers previously obtained in the same durum population genotyped with the 90K iSelect SNP assay. A total of 22 QTL were detected for traits related to durum wheat quality. Six genomic regions responsible for GPC control were mapped on chromosomes 2B, 3A, 4A, 4B, 5B, and 7B, with major QTL on chromosomes 2B, 4A, and 5B. Nine loci were detected for GYS: two on chromosome 5B and 7A and one on chromosomes 2A, 2B, 4A, 4B, 7B, with the strongest QTL on 2B. Eight QTL were identified for TKW, three of which located on chromosome 3A, two on 1B and one on 4B, 5A, and 5B. Only small overlapping was found among QTL for GYS, TKW, and GPC, and increasing alleles coming from both parents on different chromosomes. Good candidate genes were identified in the QTL confidence intervals for GYS and TKW.

A Major and Stable QTL for Bacterial Wilt Resistance on Chromosome B02 Identified Using a High-Density SNP-Based Genetic Linkage Map in Cultivated Peanut Yuanza 9102 Derived Population

Bacterial wilt (BW) is one of the important diseases limiting the production of peanut (Arachis hypogaea L.) worldwide. The sufficient precise information on the quantitative trait loci (QTL) for BW resistance is essential for facilitating gene mining and applying in molecular breeding. Cultivar Yuanza 9102 is BW resistant, bred from wide cross between cultivated peanut Baisha 1016 and a wild diploid peanut species A. chacoense with BW resistance. In this study, we aim to map the major QTLs related to BW-resistance in Yuanza 9102. A high density SNP-based genetic linkage map was constructed through double-digest restriction-site-associated DNA sequencing (ddRADseq) technique based on Yuanza 9102 derived recombinant inbred lines (RILs) population. The map contained 2,187 SNP markers distributed on 20 linkage groups (LGs) spanning 1566.10 cM, and showed good synteny with AA genome from A. duranensis and BB genome from A. ipaensis. Phenotypic frequencies of BW resistance among RIL population showed two-peak distribution in four environments. Four QTLs explaining 5.49 to 23.22% phenotypic variance were identified to be all located on chromosome B02. The major QTL, qBWB02.1 (12.17–23.33% phenotypic variation explained), was detected in three environments showing consistent and stable expression. Furthermore, there was positive additive effect among these major and minor QTLs. The major QTL region was mapped to a region covering 2.3 Mb of the pseudomolecule B02 of A. ipaensis which resides in 21 nucleotide-binding site -leucine-rich repeat (NBS-LRR) encoding genes. The result of the major stable QTL (qBWB02.1) not only offers good foundation for discovery of BW resistant gene but also provide opportunity for deployment of the QTL in marker-assisted breeding in peanut.

Construction of High-Density Genetic Maps and Detection of QTLs Associated With Huanglongbing Tolerance in Citrus

Huanglongbing (HLB), or citrus greening, is the most devastating disease in citrus worldwide. Commercial citrus varieties including sweet orange (Citrus sinensis) are highly susceptible to HLB, and trifoliate orange (Poncirus trifoliata, a close Citrus relative) is widely considered resistant or highly tolerant to HLB. In this study, an intergeneric F1 population of sweet orange and trifoliate orange was genotyped by Genotyping-by-Sequencing, and high-density SNP-based genetic maps were constructed separately for trifoliate orange and sweet orange. The two genetic maps exhibited high synteny and high coverage of the citrus genome. Progenies of the F1 population and their parents were planted in a replicated field trial, exposed to intense HLB pressure for 3 years, and then evaluated for susceptibility to HLB over 2 years. The F1 population exhibited a wide range in severity of HLB foliar symptom and canopy damage. Genome-wide QTL analysis based on the phenotypic data of foliar symptom and canopy damage in 2 years identified three clusters of repeatable QTLs in trifoliate orange linkage groups LG-t6, LG-t8 and LG-t9. Co-localization of QTLs for two traits was observed within all three regions. Additionally, one cluster of QTLs in sweet orange (linkage group LG-s7) was also detected. The majority of the identified QTLs each explained 18-30% of the phenotypic variation, indicating their major role in determining HLB responses. These results show, for the first time, a quantitative genetic nature yet the presence of major loci for the HLB tolerance in trifoliate orange. The results suggest that sweet orange also contains useful genetic factor(s) for improving HLB tolerance in commercial citrus varieties. Findings from this study should be very valuable and timely to researchers worldwide as they are hastily searching for genetic solutions to the devastating HLB crisis through breeding, genetic engineering, or genome editing.

Construction of A High-Density Genetic Map and Mapping of Fruit Traits in Watermelon (Citrullus Lanatus L.) Based on Whole-Genome Resequencing

Watermelon (Citrullus lanatus L.) is an important horticultural crop that is grown worldwide and has a high economic value. To dissect the loci associated with important horticultural traits and to analyze the genetic and genomic information of this species, a high-density genetic map was constructed based on whole-genome resequencing (WGR), a powerful high-resolution method for single-nucleotide polymorphism (SNP) marker development, genetic map construction, and gene mapping. Resequencing of both parental lines and 126 recombinant inbred lines (RIL) resulted in the detection of 178,762 single-nucleotide polymorphism (SNP) markers in the parental lines at a sequencing depth greater than four-fold. Additionally, 2132 recombination bin markers comprising 103,029 SNP markers were mapped onto 11 linkage groups (LGs). Substantially more SNP markers were mapped to the genetic map compared with other recent studies. The total length of the linkage map was 1508.94 cM, with an average distance of 0.74 cM between adjacent bin markers. Based on this genetic map, one locus for fruit bitterness, one locus for rind color, and one locus for seed coat color with high LOD scores (58.361, 18.353, 26.852) were identified on chromosome 1, chromosome 8, and chromosome 3, respectively. These prominent loci were identified in a region of 6.16 Mb, 2.07 Mb, and 0.37 Mb, respectively. On the basis of current research, the high-density map and mapping results will provide a valuable tool for identifying candidate genes, map-based gene cloning, comparative mapping, and marker-assisted selection (MAS) in watermelon breeding.

A High-Density EST-SSR-Based Genetic Map and QTL Analysis of Dwarf Trait in Cucurbita pepo L

As one of the earliest domesticated species, Cucurbita pepo (including squash and pumpkin) is rich in phenotypic polymorphism and has huge economic value. In this research, using 1660 expressed sequence tags-simple sequence repeats (EST-SSRs) and 632 genomic simple sequence repeats (gSSRs), we constructed the highest-density EST-SSR-based genetic map in Cucurbita genus, which spanned 2199.1 cM in total and harbored 623 loci distributed in 20 linkage groups. Using this map as a bridge, the two previous gSSR maps were integrated by common gSSRs and the corresponding relationships around chromosomes in three sets of genomes were also collated. Meanwhile, one large segmental inversion that existed between our map and the C. pepo genome was detected. Furthermore, three Quantitative Trait Loci (QTLs) of the dwarf trait (gibberellin-sensitive dwarf type) in C. pepo were located, and the candidate region that covered the major QTL spanned 1.39 Mb, which harbored a predicted gibberellin 2-β-oxidase gene. Considering the rich phenotypic polymorphism, the important economic value in the Cucurbita genus species and several advantages of the SSR marker were identified; thus, this high-density EST-SSR-based genetic map will be useful in Pumpkin and Squash breeding work in the future.

Construction of genetic linkage map and identification of QTLs related to agronomic traits in maize using DNA transposon-based markers

Transposable elements (TEs), are a rich source for molecular marker development as they constitute a significant fraction of the eukaryotic genome and impact the overall genome structure. Here, we utilize Mutator-based transposon display (Mu-TD), and CACTA-derived sequence-characterized amplified regions (SCAR) anchored by simple sequence repeats and single nucleotide polymorphisms to locate quantitative trait loci (QTLs) linked to agriculturally important traits on a genetic map. Specifically, we studied recombinant inbred line populations derived from a cross between dent corn and waxy corn. The resulting linkage map included 259 Mu-anchored fragments, 34 SCARs, and 614 SSR markers distributed throughout the ten maize chromosomes. Linkage analysis revealed three SNP loci associated with kernel starch synthesis genes (sh2, su1, wx1) linked to either Mu-TD loci or SSR markers, which may be useful for maize breeding programs. In addition, we used QTL analysis to determine the chromosomal location of traits related to grain yield and kernel quality. We identified 24 QTLs associated with nine traits located on nine out of ten maize chromosomes. Among these, 13 QTLs involved Mu loci and two involved SCARs. This study demonstrates the potential use of DNA transposon-based markers to construct linkage maps and identify QTLs linked to agronomic traits.

Integration of Genetic and Cytogenetic Maps and Identification of Sex Chromosome in Garden Asparagus (Asparagus officinalis L.)

A genetic linkage map of dioecious garden asparagus (Asparagus officinalis L., 2n = 2x = 20) was constructed using F₁ population, simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers. In total, 1376 SNPs and 27 SSRs were used for genetic mapping. Two resulting parental maps contained 907 and 678 markers spanning 1947 and 1814 cM, for female and male parent, respectively, over ten linkage groups representing ten haploid chromosomes of the species. With the aim to anchor the ten genetic linkage groups to individual chromosomes and develop a tool to facilitate genome analysis and gene cloning, we have optimized a protocol for flow cytometric chromosome analysis and sorting in asparagus. The analysis of DAPI-stained suspensions of intact mitotic chromosomes by flow cytometry resulted in histograms of relative fluorescence intensity (flow karyotypes) comprising eight major peaks. The analysis of chromosome morphology and localization of 5S and 45S rDNA by FISH on flow-sorted chromosomes, revealed that four chromosomes (IV, V, VI, VIII) could be discriminated and sorted. Seventy-two SSR markers were used to characterize chromosome content of individual peaks on the flow karyotype. Out of them, 27 were included in the genetic linkage map and anchored genetic linkage groups to chromosomes. The sex determining locus was located on LG5, which was associated with peak V representing a chromosome with 5S rDNA locus. The results obtained in this study will support asparagus improvement by facilitating targeted marker development and gene isolation using flow-sorted chromosomes.

NMR Metabolomics Defining Genetic Variation in Pea Seed Metabolites

Nuclear magnetic resonance (NMR) spectroscopy profiling was used to provide an unbiased assessment of changes to the metabolite composition of seeds and to define genetic variation for a range of pea seed metabolites. Mature seeds from recombinant inbred lines, derived from three mapping populations for which there is substantial genetic marker linkage information, were grown in two environments/years and analyzed by non-targeted NMR. Adaptive binning of the NMR metabolite data, followed by analysis of quantitative variation among lines for individual bins, identified the main genomic regions determining this metabolic variability and the variability for selected compounds was investigated. Analysis by t-tests identified a set of bins with highly significant associations to genetic map regions, based on probability (p) values that were appreciably lower than those determined for randomized data. The correlation between bins showing high mean absolute deviation and those showing low p-values for marker association provided an indication of the extent to which the genetics of bin variation might be explained by one or a few loci. Variation in compounds related to aromatic amino acids, branched-chain amino acids, sucrose-derived metabolites, secondary metabolites and some unidentified compounds was associated with one or more genetic loci. The combined analysis shows that there are multiple loci throughout the genome that together impact on the abundance of many compounds through a network of interactions, where individual loci may affect more than one compound and vice versa. This work therefore provides a framework for the genetic analysis of the seed metabolome, and the use of genetic marker data in the breeding and selection of seeds for specific seed quality traits and compounds that have high commercial value.

Dissecting Key Adaptation Traits in the Polyploid Perennial Medicago sativa Using GBS-SNP Mapping

Understanding key adaptation traits is crucial to developing new cultivars with broad adaptations. The main objective of this research is to understand the genetic basis of winter hardiness (WH) and fall dormancy (FD) in alfalfa and the association between the two traits. QTL analysis was conducted in a pseudo-testcross F1 population developed from two cultivars contrasting in FD (3010 with FD = 2 and CW 1010 with FD = 10). The mapping population was evaluated in three replications at two locations (Watkinsville and Blairsville, GA). FD levels showed low to moderate correlations with WH (0.22-0.57). Assessing dormancy in winter is more reliable than in the fall in southern regions with warm winters. The mapping population was genotyped using Genotyping-by-sequencing (GBS). Single dose allele SNPs (SDA) were used for constructing linkage maps. The parental map (CW 1010) consisted of 32 linkage groups spanning 2127.5 cM with 1377 markers and an average marker density of 1.5 cM/SNP. The maternal map (3010) had 32 linkage groups spanning 2788.4 cM with 1837 SDA SNPs with an average marker density of 1.5 cM/SNP. Forty-five significant (P < 0.05) QTLs for FD and 35 QTLs for WH were detected on both male and female linkage maps. More than 75% (22/28) of the dormancy QTL detected from the 3010 parent did not share genomic regions with WH QTLs and more than 70% (12/17) dormancy QTLs detected from CW 1010 parent were localized in different genomic regions than WH QTLs. These results suggest that the two traits have independent inheritance and therefore can be improved separately in breeding programs.