QTL mapping of yield-related traits in sesame

Improving yield is one of the most important targets of sesame breeding. Identifying quantitative trait loci (QTLs) of yield-related traits is a prerequisite for marker-assisted selection (MAS) and QTL/gene cloning. In this study, a BC1 population was developed and genotyped with the specific-locus amplified fragment (SLAF) sequencing technology, and a high-density genetic map was constructed. The map consisted of 13 linkage groups, contained 3528 SLAF markers, and covered a total of 1312.52 cM genetic distance, with an average distance of 0.37 cM between adjacent markers. Based on the map, 46 significant QTLs were identified for seven yield-related traits across three environments. These QTLs distributed on 11 linkage groups, each explaining 2.34-71.41% of the phenotypic variation. Of the QTLs, 23 were stable QTLs that were detected in more than one environment, and 20 were major QTLs that explained more than 10% of the corresponding phenotypic variation in at least one environment. Favorable alleles of 38 QTLs originated from the locally adapted variety, Yuzhi 4; the exotic germplasm line, BS, contributed favorable alleles to only 8 QTLs. The results should provide useful information for future molecular breeding and functional gene cloning.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01236-x.

AFLAP: assembly-free linkage analysis pipeline using k-mers from genome sequencing data

Our assembly-free linkage analysis pipeline (AFLAP) identifies segregating markers as k-mers in the raw reads without using a reference genome assembly for calling variants and provides genotype tables for the construction of unbiased, high-density genetic maps without a genome assembly. AFLAP is validated and contrasted to a conventional workflow using simulated data. AFLAP is applied to whole genome sequencing and genotype-by-sequencing data of F1, F2, and recombinant inbred populations of two different plant species, producing genetic maps that are concordant with genome assemblies. The AFLAP-based genetic map for Bremia lactucae enables the production of a chromosome-scale genome assembly.

Epidemiology of ataxia and hereditary spastic paraplegia in Spain: a cross-sectional study

INTRODUCTION

Ataxia and hereditary spastic paraplegia are rare neurodegenerative syndromes. We aimed to determine the prevalence of these disorders in Spain in 2019.

PATIENTS AND METHODS

We conducted a cross-sectional, multicentre, retrospective, descriptive study of patients with ataxia and hereditary spastic paraplegia in Spain between March 2018 and December 2019.

RESULTS

We gathered data from a total of 1.809 patients from 11 autonomous communities, provided by 47 neurologists or geneticists. Mean (SD) age in our sample was 53.64 (20.51) years; 920 patients were men (50.8%) and 889 were women (49.2%). The genetic defect was unidentified in 920 patients (47.6%). A total of 1371 patients (70.9%) had ataxia and 562 (29.1%) had hereditary spastic paraplegia. Prevalence rates for ataxia and hereditary spastic paraplegia were estimated at 5.48 and 2.24 cases per 100 000 population, respectively. The most frequent type of dominant ataxia in our sample was SCA3, and the most frequent recessive ataxia was Friedreich ataxia. The most frequent type of dominant hereditary spastic paraplegia in our sample was SPG4, and the most frequent recessive type was SPG7.

CONCLUSIONS

In our sample, the estimated prevalence of ataxia and hereditary spastic paraplegia was 7.73 cases per 100 000 population. This rate is similar to those reported for other countries. Genetic diagnosis was not available in 47.6% of cases. Despite these limitations, our study provides useful data for estimating the necessary healthcare resources for these patients, raising awareness of these diseases, determining the most frequent causal mutations for local screening programmes, and promoting the development of clinical trials.

Development and validation of genome-wide InDel markers with high levels of polymorphism in bitter gourd (Momordica charantia)

Background

The preferred choice for molecular marker development is identifying existing variation in populations through DNA sequencing. With the genome resources currently available for bitter gourd (Momordica charantia), it is now possible to detect genome-wide insertion-deletion (InDel) polymorphisms among bitter gourd populations, which guides the efficient development of InDel markers.

Results

Here, using bioinformatics technology, we detected 389,487 InDels from 61 Chinese bitter gourd accessions with an average density of approximately 1298 InDels/Mb. Then we developed a total of 2502 unique InDel primer pairs with a polymorphism information content (PIC) ≥0.6 distributed across the whole genome. Amplification of InDels in two bitter gourd lines ‘47–2–1-1-3’ and ‘04–17,’ indicated that the InDel markers were reliable and accurate. To highlight their utilization, the InDel markers were employed to construct a genetic map using 113 ‘47–2–1-1-3’ × ‘04–17’ F₂ individuals. This InDel genetic map of bitter gourd consisted of 164 new InDel markers distributed on 15 linkage groups with a coverage of approximately half of the genome.

Conclusions

This is the first report on the development of genome-wide InDel markers for bitter gourd. The validation of the amplification and genetic map construction suggests that these unique InDel markers may enhance the efficiency of genetic studies and marker-assisted selection for bitter gourd.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07499-0.

QTLs and candidate genes analyses for fruit size under domestication and differentiation in melon (Cucumis melo L.) based on high resolution maps

BACKGROUND

Melon is a very important horticultural crop produced worldwide with high phenotypic diversity. Fruit size is among the most important domestication and differentiation traits in melon. The molecular mechanisms of fruit size in melon are largely unknown.

RESULTS

Two high-density genetic maps were constructed by whole-genome resequencing with two F₂ segregating populations (WAP and MAP) derived from two crosses (cultivated agrestis × wild agrestis and cultivated melo × cultivated agrestis). We obtained 1,871,671 and 1,976,589 high quality SNPs that show differences between parents in WAP and MAP. A total of 5138 and 5839 recombination events generated 954 bins in WAP and 1027 bins in MAP with the average size of 321.3 Kb and 301.4 Kb respectively. All bins were mapped onto 12 linkage groups in WAP and MAP. The total lengths of two linkage maps were 904.4 cM (WAP) and 874.5 cM (MAP), covering 86.6% and 87.4% of the melon genome. Two loci for fruit size were identified on chromosome 11 in WAP and chromosome 5 in MAP, respectively. An auxin response factor and a YABBY transcription factor were inferred to be the candidate genes for both loci.

CONCLUSION

The high-resolution genetic maps and QTLs analyses for fruit size described here will provide a better understanding the genetic basis of domestication and differentiation, and provide a valuable tool for map-based cloning and molecular marker assisted breeding.

Mapping QTLs for 1000-grain weight and genes controlling hull type using SNP marker in Tartary buckwheat (Fagopyrum tataricum)

BACKGROUND

Tartary buckwheat (Fagopyrum tataricum), an important pseudocereal crop, has high economic value due to its nutritional and medicinal properties. However, dehulling of Tartary buckwheat is difficult owing to its thick and tough hull, which has greatly limited the development of the Tartary buckwheat processing industry. The construction of high-resolution genetic maps serves as a basis for identifying quantitative trait loci (QTLs) and qualitative trait genes for agronomic traits. In this study, a recombinant inbred lines (XJ-RILs) population derived from a cross between the easily dehulled Rice-Tartary type and Tartary buckwheat type was genotyped using restriction site-associated DNA (RAD) sequencing to construct a high-density SNP genetic map. Furthermore, QTLs for 1000-grain weight (TGW) and genes controlling hull type were mapped in multiple environments.

RESULTS

In total, 4151 bin markers comprising 122,185 SNPs were used to construct the genetic linkage map. The map consisted of 8 linkage groups and covered 1444.15 cM, with an average distance of 0.35 cM between adjacent bin markers. Nine QTLs for TGW were detected and distributed on four loci on chromosome 1 and 4. A major locus detected in all three trials was mapped in 38.2-39.8 cM region on chromosome 1, with an LOD score of 18.1-37.0, and explained for 23.6-47.5% of the phenotypic variation. The genes controlling hull type were mapped to chromosome 1 between marker Block330 and Block331, which was closely followed by the major locus for TGW. The expression levels of the seven candidate genes controlling hull type present in the region between Block330 and Block336 was low during grain development, and no significant difference was observed between the parental lines. Six non-synonymous coding SNPs were found between the two parents in the region.

CONCLUSIONS

We constructed a high-density SNP genetic map for the first time in Tartary buckwheat. The mapped major loci controlling TGW and hull type will be valuable for gene cloning and revealing the mechanism underlying grain development and easy dehulling, and marker-assisted selection in Tartary buckwheat.

QTL mapping of agronomic traits in wheat using the UK Avalon × Cadenza reference mapping population grown in Kazakhstan

Background

The success of wheat production is largely dependent on local breeding projects that focus on the development of high-yielding cultivars with the use of novel molecular tools. One strategy for improving wheat productivity involves the deployment of diverse germplasms with a high potential yield. An important factor for achieving success involves the dissection of quantitative trait loci (QTLs) for complex agronomic traits, such as grain yield components, in targeted environments for wheat growth.

Methods

In this study, we tested the United Kingdom (UK) spring set of the doubled haploid (DH) reference population derived from the cross between two British cultivars, Avalon (winter wheat) and Cadenza (spring wheat), in the Northern, Central, and Southern regions (Karabalyk, Karaganda, Kyzylorda) of Kazakhstan over three years (2013–2015). The DH population has previously been genotyped by UK scientists using 3647 polymorphic DNA markers. The list of tested traits includes the heading time, seed maturation time, plant height, spike length, productive tillering, number of kernels per spike, number of kernels per meter, thousand kernel weight, and yield per square meter. Windows QTL Cartographer was applied for QTL mapping using the composite interval mapping method.

Results

In total, 83 out of 232 QTLs were identified as stable QTLs from at least two environments. A literature survey suggests that 40 QTLs had previously been reported elsewhere, indicating that this study identified 43 QTLs that are presumably novel marker-trait associations (MTA) for these environments. Hence, the phenotyping of the DH population in new environments led to the discovery of novel MTAs. The identified SNP markers associated with agronomic traits in the DH population could be successfully used in local Kazakh breeding projects for the improvement of wheat productivity.

Identification of a quantitative trait loci (QTL) associated with ammonia tolerance in the Pacific white shrimp (Litopenaeus vannamei)

Background

Ammonia is one of the most common toxicological environment factors affecting shrimp health. Although ammonia tolerance in shrimp is closely related to successful industrial production, few genetic studies of this trait are available.

Results

In this study, we constructed a high-density genetic map of the Pacific white shrimp (Litopenaeus vannamei) using specific length amplified fragment sequencing (SLAF-seq). The constructed genetic map contained 17,338 polymorphic markers spanning 44 linkage groups, with a total distance of 6360.12 centimorgans (cM) and an average distance of 0.37 cM. Using this genetic map, we identified a quantitative trait locus (QTL) that explained 7.41–8.46% of the phenotypic variance in L. vannamei survival time under acute ammonia stress. We then sequenced the transcriptomes of the most ammonia-tolerant and the most ammonia-sensitive individuals from each of four genetically distinct L. vannamei families. We found that 7546 genes were differentially expressed between the ammonia-tolerant and ammonia-sensitive individuals. Using QTL analysis and the transcriptomes, we identified one candidate gene (annotated as an ATP synthase g subunit) associated with ammonia tolerance.

Conclusions

In this study, we constructed a high-density genetic map of L. vannamei and identified a QTL for ammonia tolerance. By combining QTL and transcriptome analyses, we identified a candidate gene associated with ammonia tolerance. Our work provides the basis for future genetic studies focused on molecular marker-assisted selective breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07254-x.

Genetic mapping of male sterility and pollen fertility QTLs in triticale with sterilizing Triticum timopheevii cytoplasm

Cytoplasmic male sterility (CMS) phenomenon is widely exploited in commercial hybrid seed production in economically important crop species, including rye, wheat, maize, rice, sorghum, cotton, sugar beets, and many vegetables. Although some commercial successes, little is known about QTLs responsible for the trait in case of triticale with sterilizing Triticum timopheevii (Tt) cytoplasm. Recombinant inbred line (RIL) F6 mapping population encompassing 182 individuals derived from the cross of individual plants representing the HT352 line and cv Borwo was employed for genetic map construction using SNP markers and identification of QTLs conferring pollen sterility in triticale with CMS Tt. The phenotypes of the F1 lines resulting from crossing of the HT352 (Tt) with HT352 (maintainer) × Borwo were determined by assessing the number of the F2 seeds per spike. A genetic map with 21 linkage groups encompasses 29,737 markers and spanned over the distance of 2549 cM. Composite (CIM) and multiple (MIM) interval mappings delivered comparable results. Single QTLs mapped to the 1A, 1B, 2A, 2R, 3B, 3R, 4B, and 5B chromosomes, whereas the 5R and 6B chromosomes shared 3 and 2 QTLs, respectively. The QTLs with the highest LOD score mapped to the 5R, 3R, 1B, and 4B chromosomes; however, the QRft-5R.3 has the highest explained variance of the trait.

Genetic map and heritability of Aspergillus flavus

The carcinogenic aflatoxins are a human health concern as well as an economic burden to corn, peanut and other crops grown within the United States and globally. Aflatoxins are produced by fungi species in Aspergillus section Flavi, primarily Aspergillus flavus. Though previously thought of as only asexual, A. flavus has recently been found to undergo sexual reproduction both in laboratory crosses and in the field. To elucidate the consequences of genetic exchange through a single generation of the sexual cycle within A. flavus, we constructed genetic maps based on three mapping populations, each composed of the parental strains and approximately 70 F1 progeny. Genome-wide data using double digest Restriction Associated DNA sequencing identified 496, 811, and 576 significant polymorphisms differentiating parents across eight linkage groups; these polymorphisms served as markers. Average spacing between marker loci was 3.1, 2.1, and 3.5 map units and overall map length was 1504.4, 1669.2, and 2001.3 cM. Recombination was non-randomly distributed across chromosomes with an average rate of recombination of about 46.81 cM per Mbp. We showed inheritance of mitochondrial loci from the sclerotial (female) parent in crosses, whereas nuclear loci showed a 1:1 segregation ratio from both parents. The linkage map will be useful in QTL analyses to identify traits that increase sexual fertility in A. flavus and modulate aflatoxin production, both of which have significant implications for sustainable reduction of aflatoxin contamination using biological control agents.

Improved reconstruction and comparative analysis of chromosome 12 to rectify Mis-assemblies in Gossypium arboreum

BACKGROUND

Genome sequencing technologies have been improved at an exponential pace but precise chromosome-scale genome assembly still remains a great challenge. The draft genome of cultivated G. arboreum was sequenced and assembled with shotgun sequencing approach, however, it contains several misassemblies. To address this issue, we generated an improved reassembly of G. arboreum chromosome 12 using genetic mapping and reference-assisted approaches and evaluated this reconstruction by comparing with homologous chromosomes of G. raimondii and G. hirsutum.

RESULTS

In this study, we generated a high quality assembly of the 94.64 Mb length of G. arboreum chromosome 12 (A_A12) which comprised of 144 scaffolds and contained 3361 protein coding genes. Evaluation of results using syntenic and collinear analysis of reconstructed G. arboreum chromosome A_A12 with its homologous chromosomes of G. raimondii (D_D08) and G. hirsutum (AD_A12 and AD_D12) confirmed the significant improved quality of current reassembly as compared to previous one. We found major misassemblies in previously assembled chromosome 12 (A_Ca9) of G. arboreum particularly in anchoring and orienting of scaffolds into a pseudo-chromosome. Further, homologous chromosomes 12 of G. raimondii (D_D08) and G. arboreum (A_A12) contained almost equal number of transcription factor (TF) related genes, and showed good collinear relationship with each other. As well, a higher rate of gene loss was found in corresponding homologous chromosomes of tetraploid (AD_A12 and AD_D12) than diploid (A_A12 and D_D08) cotton, signifying that gene loss is likely a continuing process in chromosomal evolution of tetraploid cotton.

CONCLUSION

This study offers a more accurate strategy to correct misassemblies in sequenced draft genomes of cotton which will provide further insights towards its genome organization.

High resolution mapping of restoration of fertility (Rf) by combining large population and high density genetic map in pigeonpea [Cajanus cajan (L.) Millsp]

Background

Restoration of fertility (Rf) is an important trait for pigeonpea hybrid breeding. Few coarse quantitative trait locus (QTL) studies conducted in the past identified QTLs with large confidence intervals on the genetic map and could not provide any information on possible genes responsible for Rf in pigeonpea. Therefore, a larger population comprising of 369 F₂s derived from ICPA 2039 × ICPL 87119 was genotyped with high density Axiom Cajanus SNP Array with 56 K single nucleotide polymorphism (SNPs) for high resolution mapping of Rf.

Results

A genetic map with 4867 markers was developed and a total of four QTLs for Rf were identified. While one major effect QTL (qRf8.1) was co-localized with the QTL identified in two previous studies and its size was refined from 1.2 Mb to 0.41 Mb. Further analysis of qRf8.1 QTL with genome sequence provided 20 genes including two genes namely flowering locus protein T and 2-oxoglutarate/Fe (II)-dependent dioxygenases (2-ODDs) superfamily protein with known function in the restoration of fertility.

Conclusion

The qRf8.1 QTL and the potential candidate genes present in this QTL will be valuable for genomics-assisted breeding and identification of causal genes/nucleotides for the restoration of fertility in the hybrid breeding program of pigeonpea.

Identification of QTLs associated with curd architecture in cauliflower

BACKGROUND

Curd architecture is one of the most important characters determining the curd morphology of cauliflower. However, the genetic mechanism dissection of this complex trait at molecular level is lacking. Genes/QTLs responsible for the morphological differences between present-day loose-curd and compact-curd cauliflower haven't been well revealed.

RESULTS

Herein, by using a common compact-curd parent and two loose-curd parents, we developed two double haploid (DH) populations including 122 and 79 lines, respectively. For each population, we decomposed the curd architecture concept into four parameters (basal diameter, stalk length, stalk angle and curd solidity), and collected corresponding phenotypic data for each parameter across two environments. The Kosambi function and composite interval mapping algorithm were conducted to construct the linkage map and analyze the QTLs associated with curd architecture parameters. A total of 20 QTLs were detected with the minimum likelihood of odd (LOD) values ranging from 2.61 to 8.38 and the percentage of the phenotypic variance explained by each QTL (PVE) varying between 7.69 and 25.10%. Of these, two QTLs controlling stalk length (qSL.C6-1, qSL.C6-2) and two QTLs controlling curd solidity (qCS.C6-1 and qCS.C6-2) were steadily expressed in both environments. Further, qSL.C6-1, qSL.C6-2, qCS.C6-1 and qCS.C6-4 fell into the same chromosomal region of the reference genome, indicating that these loci are involved in pleiotropic effects or are tightly linked.

CONCLUSION

The current study identified a series of QTLs associated with curd architecture parameters, which might contribute essentially to the formation of present-day loose-curd cauliflower that is widely cultivated in China. These results may pave the way for intensive deciphering the molecular mechanisms of curd development and for marker-assisted selection of curd morphology in cauliflower breeding.

Ethylene-responsive factor 4 is associated with the desirable rind hardness trait conferring cracking resistance in fresh fruits of watermelon

Fruit rind plays a pivotal role in alleviating water loss and disease and particularly in cracking resistance as well as the transportability, storability and shelf-life quality of the fruit. High susceptibility to cracking due to low rind hardness is largely responsible for severe annual yield losses of fresh fruits such as watermelon in the field and during the postharvest process. However, the candidate gene controlling the rind hardness phenotype remains unclear to date. Herein, we report, for the first time, an ethylene-responsive transcription factor 4 (ClERF4) associated with variation in rind hardness via a combinatory genetic map with bulk segregant analysis (BSA). Strikingly, our fine-mapping approach revealed an InDel of 11 bp and a neighbouring SNP in the ClERF4 gene on chromosome 10, conferring cracking resistance in F₂ populations with variable rind hardness. Furthermore, the concomitant kompetitive/competitive allele-specific PCR (KASP) genotyping data sets of 104 germplasm accessions strongly supported candidate ClERF4 as a causative gene associated with fruit rind hardness variability. In conclusion, our results provide new insight into the underlying mechanism controlling rind hardness, a desirable trait in fresh fruit. Moreover, the findings will further enable the molecular improvement of fruit cracking resistance in watermelon via precisely targeting the causative gene relevant to rind hardness, ClERF4.

Chromosomal features revealed by comparison of genetic maps of Glycine max and Glycine soja

Recombination is a crucial component of evolution and breeding. New combinations of variation on chromosomes are shaped by recombination. Recombination is also involved in chromosomal rearrangements. However, recombination rates vary tremendously among chromosome segments. Genome-wide genetic maps are one of the best tools to study variation of recombination. Here, we describe high density genetic maps of Glycine max and Glycine soja constructed from four segregating populations. The maps were used to identify chromosomal rearrangements and find the highly predictable pattern of cross-overs on the broad scale in soybean. Markers on these genetic maps were used to evaluate assembly quality of the current soybean reference genome sequence. We find a strong inversion candidate larger than 3 Mb based on patterns of cross-overs. We also identify quantitative trait loci (QTL) that control number of cross-overs. This study provides fundamental insights relevant to practical strategy for breeding programs and for pan-genome researches.

Comparative analysis of chromosome 2A molecular organization in diploid and hexaploid wheat

Diploid A genome wheat species harbor immense genetic variability which has been targeted and proven useful in wheat improvement. Development and deployment of sequence-based markers has opened avenues for comparative analysis, gene transfer and marker assisted selection (MAS) using high throughput cost effective genotyping techniques. Chromosome 2A of wheat is known to harbor several economically important genes. The present study aimed at identification of genic sequences corresponding to full length cDNAs and mining of SSRs and ISBPs from 2A draft sequence assembly of hexaploid wheat cv. Chinese Spring for marker development. In total, 1029 primer pairs including 478 gene derived, 501 SSRs and 50 ISBPs were amplified in diploid A genome species Triticum monococcum and T. boeoticum identifying 221 polymorphic loci. Out of these, 119 markers were mapped onto a pre-existing chromosome 2A genetic map consisting of 42 mapped markers. The enriched genetic map constituted 161 mapped markers with final map length of 549.6 cM. Further, 2A genetic map of T. monococcum was anchored to the physical map of 2A of cv. Chinese Spring which revealed several rearrangements between the two species. The present study generated a highly saturated genetic map of 2A and physical anchoring of genetically mapped markers revealed a complex genetic architecture of chromosome 2A that needs to be investigated further.

Construction of a high-density genetic map and QTL mapping of leaf traits and plant growth in an interspecific F1 population of Catalpa bungei × Catalpa duclouxii Dode

BACKGROUND

Catalpa bungei is an important tree species used for timber in China and widely cultivated for economic and ornamental purposes. A high-density linkage map of C. bungei would be an efficient tool not only for identifying key quantitative trait loci (QTLs) that affect important traits, such as plant growth and leaf traits, but also for other genetic studies.

RESULTS

Restriction site-associated DNA sequencing (RAD-seq) was used to identify molecular markers and construct a genetic map. Approximately 280.77 Gb of clean data were obtained after sequencing, and in total, 25,614,295 single nucleotide polymorphisms (SNPs) and 2,871,647 insertions-deletions (InDels) were initially identified in the genomes of 200 individuals of a C. bungei (7080) × Catalpa duclouxii (16-PJ-3) F₁ population and their parents. Finally, 9072 SNP and 521 InDel markers that satisfied the requirements for constructing a genetic map were obtained. The integrated genetic map contained 9593 pleomorphic markers in 20 linkage groups and spanned 3151.63 cM, with an average distance between adjacent markers of 0.32 cM. Twenty QTLs for seven leaf traits and 13 QTLs for plant height at five successive time points were identified using our genetic map by inclusive composite interval mapping (ICIM). Q16-60 was identified as a QTL for five leaf traits, and three significant QTLs (Q9-1, Q18-66 and Q18-73) associated with plant growth were detected at least twice. Genome annotation suggested that a cyclin gene participates in leaf trait development, while the growth of C. bungei may be influenced by CDC48C and genes associated with phytohormone synthesis.

CONCLUSIONS

This is the first genetic map constructed in C. bungei and will be a useful tool for further genetic study, molecular marker-assisted breeding and genome assembly.

A Molecular Cytogenetic Map of Scallop (Patinopecten yessoensis)

To consolidate the genetic, physical, and cytogenetic maps of scallop (Patinopecten yessoensis), we constructed a molecular cytogenetic map by localizing 84 fosmid clones that contain different SNP markers from 19 linkage groups (LGs) using fluorescence in situ hybridization (FISH). Among these 84 SNP-anchored clones, 56 clones produced specific and stable signals on one pair of chromosomes. Dual-color FISH assigned 19 LGs to their corresponding chromosomes with 38 SNP-anchored clones as probes. Among these 19 LGs, 17 LGs were assigned to their corresponding one pair of chromosomes, while two clones containing SNPs from LG10 and LG19 were located on two different pairs of chromosomes separately. The orientation of 7 LGs was corrected according to the chromosome location of SNPs within the same LG. In addition, a probe panel of SNP-anchored clones was developed to identify each chromosome of P. yessoensis. The molecular cytogenetic map will facilitate molecular breeding in scallop and enable comparative studies on chromosome evolution of bivalve mollusk.

Mapping and Validating QTL for Fatty Acid Compositions and Growth Traits in Asian Seabass

Asian seabass is an important food fish species. While improving growth, increasing the nutritional value is important, omega-3 fatty acids are indispensable to human health. Identifying and validating DNA markers associated with traits is the first step towards marker-assisted selection (MAS). We quantified 13 different fatty acids and three growth traits in 213 F₂ Asian seabass from a family at the age 270 days post hatch, and screened QTL for these traits. The content of total fatty acids in 100 g flesh was 2.57 ± 0.80 g, while the proportions of docosahexaenoic acid (DHA) and eicosapentaenoic acid (EPA) were 16.96 ± 2.20% and 5.42 ± 0.90%, respectively. A linkage map with 2424 SNPs was constructed and used for QTL mapping. For fatty acid compositions, 14 significant QTL were identified on three linkage groups (LG5, LG11 and LG14), with phenotypic variance explained (PVE) from 12.8 to 24.6%. Thirty-nine suggestive QTL were detected on 16 LGs. Two significant QTL for EPA were identified on LG5 and LG14, with PVE of 15.2% and 15.1%, respectively. No significant QTL was identified for DHA. For growth traits, six significant and 13 suggestive QTL were identified on two and seven LGs, respectively. Only a few significant QTL for fatty acids overlapped with previously mapped QTL for these traits, suggesting that most QTL detected in a family are family-specific and could only be used in MAS in the family per se. To facilitate population-wide molecular breeding, more powerful methods (e.g. GWAS) should be used to identify SNPs for genomic selection.

A genomic resource derived from the integration of genome sequences, expressed transcripts and genetic markers in ramie

Background

The redundancy of genomic resources, including transcript and molecular markers, and their uncertain position in the genome have dramatically hindered the study of traits in ramie, an important natural fiber crop.

Results

We obtained a high-quality transcriptome consisting of 30,591 non-redundant transcripts using single-molecule long-read sequencing and proposed it as a universal ramie transcriptome. Additionally, 55,882 single nucleotide polymorphisms (SNPs) were identified and a high-density genetic map was developed. Based on this genetic map, 181.7 Mb ramie genome sequences were assembled into 14 chromosomes. For the convenient use of these resources, 29,286 (~ 95.7%) of the transcripts and all 55,882 SNPs, along with 1827 previously reported sequence repeat markers (SSRs), were mapped into the ramie genome, and 22,343 (~ 73.0%) transcripts, 50,154 (~ 89.7%) SNPs, and 1466 (~ 80.3%) SSRs were assigned to a specific location in the corresponding chromosome.

Conclusion

This is the first study to characterize the ramie transcriptome by long-read sequencing, and the substantial number of transcripts of significant length obtained will accelerate our understanding of ramie growth and development. This integration of genome sequences, expressed transcripts, and genetic markers will provide an extremely useful resource for genetic, molecular, and breeding studies of ramie.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5878-8) contains supplementary material, which is available to authorized users.

Constructing High-Density Genetic Maps and Developing Sexing Markers in Northern Snakehead (Channa argus)

High-density genetic maps are essential for mapping QTL, improving genome assembly, comparative genomics, and studying sex chromosome evolution. The northern snakehead (Channa argus) is an economically important foodfish species with significant sexual dimorphism, where the males grow much faster and bigger than the females. However, to date, the sex determination pattern is still not clear, limiting identification of sex chromosomes, even sex determination genes and development of monosex populations that are valuable for both sex evolution of vertebrates and aquaculture practices. Here, a sex-averaged map and two sex-specific genetic maps were constructed with 2974, 2323, and 2338 SNPs, respectively. Little difference was observed in the pattern of sex-specific recombination between female- and male-specific genetic maps. Genome scan identified a major locus for sex determination at LG16. Females and males are, respectively, homogametic and heterogametic, suggesting an XY sex determination system for this species. By resequencing genomes, InDels in the sex-associated QTL region were discovered and used for developing sex-specific PCR assays for fast sexing of snakehead. These high-density genetic maps provide useful resources for future genomic studies in snakehead and its related species. The PCR assays for sexing are of importance in developing all male populations for aquaculture.

Construction of a high-density genetic map and identification of loci controlling purple sepal trait of flower head in Brassica oleracea L. italica

BACKGROUND

Some broccoli (Brassica oleracea L. italic) accessions have purple sepals and cold weather would deepen the purple color, while the sepals of other broccoli lines are always green even in cold winter. The related locus or gene is still unknown. In this study, a high-density genetic map was constructed based on specific locus amplified fragment (SLAF) sequencing in a doubled-haploid segregation population with 127 individuals. And mapping of the purple sepal trait in flower heads based on phenotypic data collected during three seasons was performed.

RESULTS

A genetic map was constructed, which contained 6694 SLAF markers with an average sequencing depth of 81.37-fold in the maternal line, 84-fold in the paternal line, and 15.76-fold in each individual population studied. In all of the annual data recorded, three quantitative trait loci (QTLs) were identified that were all distributed within the linkage group (LG) 1. Among them, a major locus, qPH.C01-2, located at 36.393 cM LG1, was consistently detected in all analysis. Besides this locus, another two minor loci, qPH.C01-4 and qPH.C01-5, were identified near qPH.C01-2, based on the phenotypic data from spring of 2018.

CONCLUSION

The purple sepal trait could be controlled by a major single locus and two minor loci. The genetic map and location of the purple sepal trait of flower heads provide an important foundation for mapping other compound traits and the identification of the genes related to purple sepal trait in broccoli.

A high-density genetic map developed by specific-locus amplified fragment (SLAF) sequencing and identification of a locus controlling anthocyanin pigmentation in stalk of Zicaitai (Brassica rapa L. ssp. chinensis var. purpurea)

BACKGROUND

Caixin and Zicaitai (Brassica rapa) belong to Southern and Central China respectively. Zicaitai contains high amount of anthocyanin in leaf and stalk resulting to the purple color. Stalk is the major edible part and stalk color is an economically important trait for the two vegetables. The aim of this study is to construct a high density genetic map using the specific length amplified fragment sequencing (SLAF-seq) technique to explore genetic basis for anthocyanin pigmentation traits via quantitative trait loci (QTL) mapping.

RESULTS

We constructed a high generation linkage map with a mapping panel of F2 populations derived from 150 individuals of parental lines "Xianghongtai 01" and "Yinong 50D" with purple and green stalk respectively. The map was constructed containing 4253 loci, representing 10,940 single nucleotide polymorphism (SNP) markers spanning 1030.04 centiMorgans (cM) over 10 linkage groups (LGs), with an average distance between markers of 0.27 cM. Quantitative trait loci (QTL) analysis revealed that a major locus on chromosome 7 and 4 minor QTLs explaining 2.69-61.21% of phenotypic variation (PVE) were strongly responsible for variation in stalk color trait. Bioinformatics analysis of the major locus identified 62 protein-coding genes. Among the major locus, there were no biosynthetic genes related to anthocyanin. However, there were several transcription factors like helix-loop-helix (bHLH) bHLH, MYB in the locus. Seven predicted candidate genes were selected for the transcription level analysis. Only bHLH49 transcription factor, was significantly higher expressed in both stalks and young leaves of Xianghongtai01 than Yinong50D. An insertion and deletion (InDel) marker developed from deletion/insertion in the promoter region of bHLH49 showed significant correlation with the stalk color trait in the F2 population.

CONCLUSION

Using the constructed high-qualified linkage map, this study successfully identified QTLs for stalk color trait. The identified valuable markers and candidate genes for anthocyanin accumulation in stalk will provide useful information for molecular regulation of anthocyanin biosynthesis. Overall our findings will lay a foundation for functional gene cloning, marker-assisted selection (MAS) and molecular breeding of important economic traits in B. rapa.

Construction of genetic linkage map and identification of QTLs related to agronomic traits in DH population of maize (Zea mays L.) using SSR markers

BACKGROUND

In this study, we used phenotypic and genetic analysis to investigate Double haploid (DH) lines derived from normal corn parents (HF1 and 11S6169). DH technology offers an array of advantages in maize genetics and breeding as follows: first, it significantly shortens the breeding cycle by development of completely homozygous lines in two or three generations; and second, it simplifies logistics, including requiring less time, labor, and financial resources for developing new DH lines compared with the conventional RIL population development process.

OBJECTIVES

In our study, we constructed a maize genetic linkage map using SSR markers and a DH population derived from a cross of normal corn (HF1) and normal corn (11S6169).

METHODS

The DH population used in this study was developed by the following methods: we crossed normal corn (HF1) and normal corn (11S6169), which are parent lines of a normal corn cultivar, in 2014; and the next year, the F₁ hybrids were crossed with a tropicalized haploid inducer line (TAIL), which is homozygous for the dominant marker gene R1-nj (Nanda and Chase in Crop Sci 6:213-215, 1966), and we harvested seeds of the haploid lines.

RESULTS

A total of 200 SSR markers were assigned to 10 linkage groups that spanned 1145.4 cM with an average genetic distance between markers of 5.7 cM. 68 SSR markers showed Mendelian segregation ratios in the DH population at a 5% significance threshold. A total of 15 quantitative trait loci (QTLs) for plant height (PH), ear height (EH), ear height ratio (ER), leaf length (LL), ear length (EL), set ear length (SEL), set ear ratio (SER), ear width (EW), 100 kernel weight (100 KW), and cob color (CC) were found in the 121 lines in the DH population.

CONCLUSION

The results of this study may help to improve the detection and characterization of agronomic traits and provide great opportunities for maize breeders and researchers using a DH population in maize breeding programs.

A multi-parent recombinant inbred line population of C. elegans allows identification of novel QTLs for complex life history traits

BACKGROUND

The nematode Caenorhabditis elegans has been extensively used to explore the relationships between complex traits, genotypes, and environments. Complex traits can vary across different genotypes of a species, and the genetic regulators of trait variation can be mapped on the genome using quantitative trait locus (QTL) analysis of recombinant inbred lines (RILs) derived from genetically and phenotypically divergent parents. Most RILs have been derived from crossing two parents from globally distant locations. However, the genetic diversity between local C. elegans populations can be as diverse as between global populations and could thus provide means of identifying genetic variation associated with complex traits relevant on a broader scale.

RESULTS

To investigate the effect of local genetic variation on heritable traits, we developed a new RIL population derived from 4 parental wild isolates collected from 2 closely located sites in France: Orsay and Santeuil. We crossed these 4 genetically diverse parental isolates to generate a population of 200 multi-parental RILs and used RNA-seq to obtain sequence polymorphisms identifying almost 9000 SNPs variable between the 4 genotypes with an average spacing of 11 kb, doubling the mapping resolution relative to currently available RIL panels for many loci. The SNPs were used to construct a genetic map to facilitate QTL analysis. We measured life history traits such as lifespan, stress resistance, developmental speed, and population growth in different environments, and found substantial variation for most traits. We detected multiple QTLs for most traits, including novel QTLs not found in previous QTL analysis, including those for lifespan and pathogen responses. This shows that recombining genetic variation across C. elegans populations that are in geographical close proximity provides ample variation for QTL mapping.

CONCLUSION

Taken together, we show that using more parents than the classical two parental genotypes to construct a RIL population facilitates the detection of QTLs and that the use of wild isolates facilitates the detection of QTLs. The use of multi-parent RIL populations can further enhance our understanding of local adaptation and life history trade-offs.