High-density genetic map construction and QTL mapping of a zigzag-shaped stem trait in tea plant (Camellia sinensis)

The highly unique zigzag-shaped stem phenotype in tea plants boasts significant ornamental value and is exceptionally rare. To investigate the genetic mechanism behind this trait, we developed BC1 artificial hybrid populations. Our genetic analysis revealed the zigzag-shaped trait as a qualitative trait. Utilizing whole-genome resequencing, we constructed a high-density genetic map from the BC1 population, incorporating 5,250 SNP markers across 15 linkage groups, covering 3,328.51 cM with an average marker interval distance of 0.68 cM. A quantitative trait locus (QTL) for the zigzag-shaped trait was identified on chromosome 4, within a 61.2 to 97.2 Mb range, accounting for a phenotypic variation explained (PVE) value of 13.62%. Within this QTL, six candidate genes were pinpointed. To better understand their roles, we analyzed gene expression in various tissues and individuals with erect and zigzag-shaped stems. The results implicated CsXTH (CSS0035625) and CsCIPK14 (CSS0044366) as potential key contributors to the zigzag-shaped stem formation. These discoveries lay a robust foundation for future functional genetic mapping and tea plant genetic enhancement.

Genome-wide QTL and eQTL mapping reveal genes associated with growth rate trait of the Pacific white shrimp (Litopenaeus vannamei)

BACKGROUND

Growth rate is a crucial economic trait for farmed animals, but the genetic regulation of this trait is largely unknown in non-model organisms such as shrimp.

RESULTS

In this study, we performed genome-wide phenotypic quantitative trait loci (QTL) and expression quantitative trait loci (eQTL) mapping analyses to identify genes affecting the growth rate of Pacific white shrimp (Litopenaeus vannamei), which is the most commercially-farmed crustacean worldwide. We used RNA-sequencing of 268 individuals in a mapping population, and subsequently validated our findings through gene silencing and shrimp growth experiments. We constructed a high-density genetic linkage map comprising 5533 markers spanning 44 linkage groups, with a total distance of 6205.75 cM and an average marker interval of 1.12 cM. Our analyses identified 11 QTLs significantly correlated with growth rate, and 117,525 eQTLs. By integrating QTL and eQTL data, we identified a gene (metalloreductase STEAP4) highly associated with shrimp growth rate. RNA interference (RNAi) analysis and growth experiments confirmed that STEAP4 was significantly correlated with growth rate in L. vannamei.

CONCLUSIONS

Our results indicate that the comprehensive analysis of QTL and eQTL can effectively identify genes involved in complex animal traits. This is important for marker-assisted selection (MAS) of animals. Our work contributes to the development of shrimp breeding and available genetic resources.

Refinement of Rf1-gene localization and development of the new molecular markers for fertility restoration in sunflower

BACKGROUND

Ability to restore male fertility is important trait for sunflower breeding. The most commonly used fertility restoration gene in the production of sunflower hybrids is Rf1. The localization of Rf1 on the linkage group 13 has been previously shown, however, its exact position, its sequence and molecular mechanism for fertility restoration remain unknown. Therefore, several markers linked to Rf1 gene, commonly used for MAS, don't always allow to identify the genotype of plants. For this reason, the search for new markers and precise localization of the Rf1 gene is an urgent task.

METHODS AND RESULTS

Based on previously identified single nucleotide polymorphisms (SNPs) at LG13, significantly associated with the ability to restore male fertility, two markers have been developed that have performed well after careful evaluation. These markers, together with other Rf1 markers, were applied for genotyping 72 diversity panel accessions and 291 individuals of F2 segregating population, obtained from crossing the cytoplasmic male sterility (CMS) AHO33 and restorer RT085HO lines. The analysis revealed no recombinants between Rf1 gene and SRF833 marker, the distance between Rf1 and SRF122 marker was 1.0 cM.

CONCLUSIONS

Data obtained made it possible to specify the localization of the Rf1 gene and reduce the list of candidate genes to the 3 closely linked PPR-genes spanning a total of 59 Kb. However, it cannot be ruled out that analysis of the candidate region in the genome of fertility restorer lines can reveal new candidate genes in this locus that are absent in the cytoplasmic male sterility maintainer reference sequence.

Genomics-assisted genetics of complex region from chickpea chromosome 4 reveals two candidate genes for Ascochyta blight resistance

Ascochyta blight (AB) disease caused by the fungus Ascochyta rabiei is a major threat to global chickpea production. Molecular breeding for improved AB resistance requires the identification of robust fine-mapped QTLs/candidate genes and associated markers. Earlier, we identified three QTLs (qABR4.1, qABR4.2, and qABR4.3) for AB resistance on chickpea chromosome 4 by employing multiple quantitative trait loci sequencing strategy on an intra-specific (FLIP84-92C x PI359075) and an inter-specific (FLIP84-92C x PI599072) crosses derived recombinant inbred lines. Here, we report the identification of AB resistance providing candidate genes under the fine mapped qABR4.2 and qABR4.3 genomic region by combining genetic mapping, haplotype block inheritance, and expression analysis. The qABR4.2 region was narrowed down from 5.94Mb to ~800kb. Among 34 predicted gene models, a secreted class III peroxidase encoding gene showed higher expression in AB resistant parent after A. rabiei conidia inoculation. Under qABR4.3, we identified a frame-shift mutation in a cyclic nucleotide-gated channel CaCNGC1 gene leading to the truncated N-terminal domain in resistant accession of chickpea. This N-terminal domain of CaCNGC1 interacts with chickpea calmodulin. Thus, our analysis has revealed narrowed genomic regions and their associated polymorphic markers, CaNIP43 and CaCNGCPD1. These co-dominant markers significantly associate with AB resistance on qABR4.2 and qABR4.3 regions. Our genetic analysis revealed that the presence of resistant alleles for two major QTLs (qABR4.1 and qABR4.2) together provide AB resistance in the field while minor QTL qABR4.3 determines the degree of resistance. The identified candidate genes and their diagnostic markers will help in biotechnological and AB resistance introgression into farmers adapted local chickpea varieties.

Identification of an important QTL for seed oil content in soybean

Seed oil content is one of the most important quantitative traits in soybean (Glycine max) breeding. Here, we constructed a high-density single nucleotide polymorphism linkage map using two genetically similar parents, Heinong 84 and Kenfeng 17, that differ dramatically in their seed oil contents, and performed quantitative trait loci (QTL) mapping of seed oil content in a recombinant inbred line (RIL) population derived from their cross. We detected five QTL related to seed oil content distributed on five chromosomes. The QTL for seed oil content explained over 10% of the phenotypic variation over two years. This QTL was mapped to an interval containing 20 candidate genes, including a previously reported gene, soybean RING Finger 1a (RNF1a) encoding an E3 ubiquitin ligase. Notably, two short sequences were inserted in the GmRNF1a coding region of KF 17 compared to that of HN 84, resulting in a longer protein variant in KF 17. Our results thus provide information for uncovering the genetic mechanisms determining seed oil content in soybean, as well as identifying an additional QTL and highlighting GmRNF1a as candidate gene for modulating seed oil content in soybean.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01384-2.

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 1933 patients from 11 autonomous communities, provided by 47 neurologists or geneticists. Mean (SD) age in our sample was 53.64 (20.51) years; 982 patients were men (50.8%) and 951 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.

Integrated consensus genetic map and genomic scaffold re-ordering of oil palm (Elaeis guineensis) genome

A high-quality reference genome is an important resource that can help decipher the genetic basis of traits in combination with linkage or association analyses. The publicly available oil palm draft genome sequence of AVROS pisifera (EG5) accounts for 1.535 Gb of the 1.8 Gb oil palm genome. However, the assemblies are fragmented, and the earlier assembly only had 43% of the sequences placed on pseudo-chromosomes. By integrating a number of SNP and SSR-based genetic maps, a consensus map (AM_EG5.1), comprising of 828.243 Mb genomic scaffolds anchored to 16 pseudo-chromosomes, was generated. This accounted for 54% of the genome assembly, which is a significant improvement to the original assembly. The total length of N50 scaffolds anchored to the pseudo-chromosomes increased by ∼18% compared to the previous assembly. A total of 139 quantitative trait loci for agronomically important quantitative traits, sourced from literature, were successfully mapped on the new pseudo-chromosomes. The improved assembly could also be used as a reference to identify potential errors in placement of specific markers in the linkage groups of the genetic maps used to assemble the consensus map. The 3422 unique markers from five genetic maps, anchored to the pseudo-chromosomes of AM_EG5.1, are an important resource that can be used preferentially to either construct new maps or fill gaps in existing genetic maps. Synteny analysis further revealed that the AM_EG5.1 had high collinearity with the date palm genome cultivar 'Barhee BC4' and shared most of its segmental duplications. This improved chromosomal-level genome is a valuable resource for genetic research in oil palm.

Construction of the first high-density SNP genetic map and identification of QTLs for the natural rubber content in Taraxacum kok-saghyz Rodin

BACKGROUND

Taraxacum kok-saghyz Rodin (TKS) is a promising commercial alternative natural rubber (NR) yielding plant. Cultivating TKS with a high NR content is an important breeding target, and developing molecular markers related to NR content can effectively accelerate the breeding process of TKS.

RESULTS

To construct a high-density SNP genetic map and uncover genomic regions related to the NR content in TKS, an F₁ mapping population of TKS was constructed by crossing two parents (l66 and X51) with significant differences in NR contents. The NR content of the F₁ plants ranged from 0.30 to 15.14% and was distributed normally with a coefficient of variation of 47.61%, indicating quantitative trait inheritance. Then, employing whole-genome resequencing (WGR), a TKS genetic linkage map of 12,680 bin markers comprising 322,439 SNPs was generated. Based on the genetic map and NR content of the F₁ population, six quantitative trait loci (QTLs) for NR content with LOD > 4.0 were identified on LG01/Chr01 and LG06/Chr06. Of them, the 2.17 Mb genomic region between qHRC-C6-1 and qHRC-C6-2 on ChrA06, with 65.62% PVE in total, was the major QTL region. In addition, the six QTLs have significant additive genetic effects on NR content and could be used to develop markers for marker-assisted selection (MAS) in TKS with a high NR content.

CONCLUSION

This work constructed the first high-density TKS genetic map and identified the QTLs and genomic regions controlling the NR content, which provides useful information for fine mapping, map-based cloning, and MAS in TKS.

A technical guide to TRITEX, a computational pipeline for chromosome-scale sequence assembly of plant genomes

BACKGROUND

As complete and accurate genome sequences are becoming easier to obtain, more researchers wish to get one or more of them to support their research endeavors. Reliable and well-documented sequence assembly workflows find use in reference or pangenome projects.

RESULTS

We describe modifications to the TRITEX genome assembly workflow motivated by the rise of fast and easy long-read contig assembly of inbred plant genomes and the routine deployment of the toolchains in pangenome projects. New features include the use as surrogates of or complements to dense genetic maps and the introduction of user-editable tables to make the curation of contig placements easier and more intuitive.

CONCLUSION

Even maximally contiguous sequence assemblies of the telomere-to-telomere sort, and to a yet greater extent, the fragmented kind require validation, correction, and comparison to reference standards. As pangenomics is burgeoning, these tasks are bound to become more widespread and TRITEX is one tool to get them done. This technical guide is supported by a step-by-step computational tutorial accessible under https://tritexassembly.bitbucket.io/ . The TRITEX source code is hosted under this URL: https://bitbucket.org/tritexassembly .

Genetic analysis of a DROOPING LEAF mutant allele dl-6 associated with a twisted and folded leaf base caused by a deficiency in midrib development in rice

The failure of midrib formation in rice leaf blades results in the drooping leaf (dl) phenotype. A normal DROOPING LEAF (DL) gene is necessary for leaf homeotic transformation, which affects midrib and pistil development. Genetic analysis was performed on a new drooping leaf (dl) mutant named dl-6 in rice. The dl-6 allelic mutant exhibited drooping leaves that were severely folded and twisted at the base but had normal flower structure. The dl-6 allele is a nuclear recessive trait that fits a 3:1 Mendelian segregation ratio. The dl-6 mutant leaves displayed an abnormal main vein (midrib-less) with undeveloped aerenchyma and vascular bundles, resulting in severe leaf drooping. The lack of a midrib in dl-6 caused weak mechanical support, which resulted in folding at the collar junction of the leaf base and downward bending. Through genetic mapping, the dl-6 allele was identified at approximately 28.2 cM on rice chromosome 3. The allele was caused by mutations within the DL (LOC_Os03g11600.1) gene, with specific amino acid substitutions and additions in the encoded protein of the YABBY transcription factor. The dl-6 mutant is a recessive allele encoding a dysfunctional YABBY transcription factor that regulates leaf midrib development and aerenchymatous clear cell structures, leading to a drooping leaf phenotype in rice.

Construction of a highly saturated genetic map and identification of quantitative trait loci for leaf traits in jujube

Chinese jujube (Ziziphus jujuba Mill.), a member of the genus Ziziphus, which comes under the family Rhamnaceae, is the most important species in terms of its economic, ecological, and social benefits. To dissect the loci associated with important phenotypical traits and analyze their genetic and genomic information in jujube, a whole-genome resequencing (WGR) based highly saturated genetic map was constructed using an F1 hybrid population of 140 progeny individuals derived from the cross of 'JMS2' × 'Jiaocheng 5'. The average sequencing depth of the parents was 14.09× and that of the progeny was 2.62×, and the average comparison efficiency between the sample and the reference genome was 97.09%. Three sets of genetic maps were constructed for a female parent, a male parent, and integrated. A total of 8,684 markers, including 8,158 SNP and 526 InDel markers, were evenly distributed across all 12 linkage groups (LGs) in the integrated map, spanning 1,713.22 cM with an average marker interval of 0.2 cM. In terms of marker number and density, this is the most saturated genetic map of jujube to date, nearly doubling that of the best ones previously reported. Based on this genetic map and phenotype data from 2019 to 2021, 31 leaf trait QTLs were identified in the linkage groups (LG1, 15; LG3, 1; LG5, 8; LG7, 4; LG8, 1, and LG11, 2), including 17 major QTLs. There were 4, 8, 14, and 5 QTLs that contributed to leaf length, leaf width, leaf shape index, and leaf area, respectively. Six QTLs clusters were detected on LG1 (8.05 cM-9.52 cM; 13.12 cM-13.99 cM; 123.84 cM-126.09 cM), LG5 (50.58 cM-50.86 cM; 80.10 cM-81.76 cM) and LG11 (35.98 cM-48.62 cM). Eight candidate genes were identified within the QTLs cluster regions. Annotation information showed that 4 genes (LOC107418196, LOC107418241, LOC107417968, and LOC112492570) in these QTLs are related to cell division and cell wall integrity. This research will provide a valuable tool for further QTL analysis, candidate gene identification, map-based gene cloning, comparative mapping, and marker-assisted selection (MAS) in jujube.

The recombination landscape of the Khoe-San likely represents the upper limits of recombination divergence in humans

Background

Recombination maps are  important resources for epidemiological and evolutionary analyses; however, there are currently no recombination maps representing any African population outside of those with West African ancestry. We infer the demographic history for the Nama, an indigenous Khoe-San population of southern Africa, and derive a novel, population-specific recombination map from the whole genome sequencing of 54 Nama individuals. We hypothesise that there are no publicly available recombination maps representative of the Nama, considering the deep population divergence and subsequent isolation of the Khoe-San from other African groups.

Results

We show that the recombination landscape of the Nama does not cluster with any continental groups with publicly available representative recombination maps. Finally, we use selection scans as an example of how fine-scale differences between the Nama recombination map and the combined Phase II HapMap recombination map can impact the outcome of selection scans.

Conclusions

Fine-scale differences in recombination can meaningfully alter the results of a selection scan. The recombination map we infer likely represents an upper bound on the extent of divergence we expect to see for a recombination map in humans and would be of interest to any researcher that wants to test the sensitivity of population genetic or GWAS analysis to recombination map input.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02744-5.

Accurate recombination estimation from pooled genotyping and sequencing: a case study on barley

Sexual reproduction involves meiotic recombination and the creation of crossing over between homologous chromosomes, which leads to new allele combinations. We present a new approach that uses the allele frequency differences and the physical distance of neighboring polymorphisms to estimate the recombination rate from pool genotyping or sequencing. This allows a considerable cost reduction compared to conventional mapping based on genotyping or sequencing data of single individuals. We evaluated the approach based on computer simulations at various genotyping depths and population sizes as well as applied it to experimental data of 45 barley populations, comprising 4182 RIL. High correlations between the recombination rates from this new pool genetic mapping approach and conventional mapping in simulated and experimental barley populations were observed. The proposed method therefore provides a reliable genetic map position and recombination rate estimation in defined genomic windows.

Quantitative trait loci mapping of stem sugar content and stem diameter in sorghum recombinant inbred lines using genotyping-by-sequencing

BACKGROUND

Sweet sorghum is an important crop with sugary stem that can mainly be used for syrup, fodder and bio-fuel. Many sugar content QTLs have been discovered from different sources through breeding worldwide. Most of these QTLs are detected using exotic germplasm as a mapping population. This study aimed to detect and map QTLs for stem sugar content and stem diameter targeting Ethiopian recombinant inbred lines of sorghum using genotyping-by-sequencing.

METHODS AND RESULT

Genotyping-by-sequencing and phenotyping using 139 recombinant inbred lines of sorghum as mapping populations were conducted. A total of 1082 polymorphic and high quality SNP markers that are evenly distributed across the ten linkage groups of sorghum were selected to detect and map the trait of interest. A genetic linkage map using 1082 SNP markers was constructed and several QTLs associated with stem sugar content and stem diameter were identified. Phenotypic variation explained by qBrix4-1 and qBrix2-1 ranged from 6.33 to 14%, respectively. Over two seasons, four QTLs for stem sugar content (qBrix1-1, qBrix2-1, qBrix4-1 and qBrix4-2) and three QTLs for stem diameter (qSD1-1, qSD8-1 and qSD9-1) were detected.

CONCLUSION

QTLs that significantly associated with stem sugar content and stem diameter have been detected and mapped. This will help sorghum breeding program to develop superior sweet sorghum varieties through the use of appropriate crop improvement approaches like marker assisted breeding. This ultimately contributes to the current development plan to considerably improve food, feed and bio-fuel supply in developing countries like Ethiopia.

Genome resequencing-based high-density genetic map and QTL detection for yield and fiber quality traits in diploid Asiatic cotton (Gossypium arboreum)

Cotton is the most important fiber crop in the world. Asiatic cotton (Gossypium arboreum, genome A₂) is a diploid cotton species producing spinnable fibers and important germplasm for cotton breeding and a significant model for fiber biology. However, the genetic map of Asiatic cotton has been lagging behind tetraploid cottons, as well as other stable crops. This study aimed to construct a high-density SNP genetic map and to map QTLs for important yield and fiber quality traits. Using a recombinant inbred line (RIL) population and genome resequencing technology, we constructed a high-density genetic map that covered 1980.17 cM with an average distance of 0.61 cM between adjacent markers. QTL analysis revealed a total of 297 QTLs for 13 yield and fiber quality traits in three environments, explaining 5.0-37.4% of the phenotypic variance, among which 75 were stably detected in two or three environments. Besides, 47 QTL clusters, comprising 131 QTLs for representative traits, were identified. Our works laid solid foundation for fine mapping and cloning of QTL for yield and fiber quality traits in Asiatic cotton.

Construction of a high-density genetic map based on specific-locus amplified fragment sequencing and identification of loci controlling anthocyanin pigmentation in Yunnan red radish

Radish (Raphanus sativus L.) belongs to the family Brassicaceae. The Yunnan red radish variety contains fairly relatively large amounts of anthocyanins, making them important raw materials for producing edible red pigment. However, the genetic mechanism underlying this pigmentation has not been fully characterized. Herein, the radish inbred line YAAS-WR1 (white root-skin and white root-flesh) was crossed with the inbred line YAAS-RR1 (red root-skin and red root-flesh) to produce F1, F2, BC1P1, and BC1P2 populations. Genetic analyses revealed that the pigmented/non-pigmented (PiN) and purple/red (PR) traits were controlled by two genetic loci. The F2 population and the specific-locus amplified fragment sequencing (SLAF-seq) technique were used to construct a high-density genetic map (1230.16 cM), which contained 4032 markers distributed in nine linkage groups, with a mean distance between markers of 0.31 cM. Additionally, two QTL (QAC1 and QAC2) considerably affecting radish pigmentation were detected. A bioinformatics analysis of the QAC1 region identified 58 predicted protein-coding genes. Of these genes, RsF3'H, which is related to anthocyanin biosynthesis, was revealed as a likely candidate gene responsible for the PR trait. The results were further verified by analyzing gene structure and expression. Regarding QAC2, RsMYB1.3 was determined to be a likely candidate gene important for the PiN trait, with a 4-bp insertion in the first exon that introduced a premature termination codon in the YAAS-WR1 sequence. Assays demonstrated that RsMYB1.3 interacted with RsTT8 and activates RsTT8 and RsUFGT expression. These findings may help clarify the complex regulatory mechanism underlying radish anthocyanin synthesis. Furthermore, this study's results may be relevant for the molecular breeding of radish to improve the anthocyanin content and appearance of the taproots.

Construction of a high-density genetic map using specific-length amplified fragment markers and identification of QTLs for branching angle in poplar

Branching angle is a critical factor that determines the morphological establishment and is a typical quantitative trait controlled by multiple genes. In this study, we used SLAF-seq to construct a high-density genetic map, to investigate the genetic architecture of branching angle in poplar (Populus leucopyramidalis). A total of 240,672 SLAF tags were obtained, including 103,691 polymorphic SLAF tags. After filtering, 53,407 polymorphic markers were sorted into eight segregation types, and 11,162 of them were used to construct the genetic map. 8447 were on the female parent map, 8532 were on the male parent map, and 11,162 were on the integrated map. The marker coverage was 4820.84 and 5044.80 cM for the female and male maps, and 3142.61 cM for the integrated map. The average intervals between two adjacent mapped markers were 0.55, 0.59, and 0.28 cM for the three maps, respectively. Two quantitative trait loci (QTLs) were detected. Seven markers that exceeded the threshold in these two regions were considered as being associated with branching angle and the phenotypic variance explained by each of these marker was 10.64-11.66%. After functional annotation, we identified 15 candidate genes and analyzed the expression of candidate genes in narrow and wide crown progenies by qRT-PCR. These results show that the combination of QTL and SLAF-seq will contribute to future breeding plans in poplar breeding, especially in narrow crown poplar breeding.

MG2C: a user-friendly online tool for drawing genetic maps

Genetic map is a linear arrangement of the relative positions of sites in the chromosome or genome based on the recombination frequency between genetic markers. It is the important basis for genetic analysis. Several kinds of software have been designed for genetic mapping, but all these tools require users to write or edit code, making it time-costing and difficult for researchers without programming skills to handle with. Here, MG2C, a new online tool was designed, based on PERL and SVG languages.Users can get a standard genetic map, only by providing the location of genes (or quantitative trait loci) and the length of the chromosome, without writing additional code. The operation interface of MG2C contains three sections: data input, data output and parameters. There are 33 attribute parameters in MG2C, which are further divided into 8 modules. Values of the parameters can be changed according to the users' requirements. The information submitted by users will be transformed into the genetic map in SVG file, which can be further modified by other image processing tools.MG2C is a user-friendly and time-saving online tool for drawing genetic maps, especially for those without programming skills. The tool has been running smoothly since 2015, and updated to version 2.1. It significantly lowers the technical barriers for the users, and provides great convenience for the researchers.

A linkage disequilibrium-based approach to position unmapped SNPs in crop species

BACKGROUND

High-density SNP arrays are now available for a wide range of crop species. Despite the development of many tools for generating genetic maps, the genome position of many SNPs from these arrays is unknown. Here we propose a linkage disequilibrium (LD)-based algorithm to allocate unassigned SNPs to chromosome regions from sparse genetic maps. This algorithm was tested on sugarcane, wheat, and barley data sets. We calculated the algorithm's efficiency by masking SNPs with known locations, then assigning their position to the map with the algorithm, and finally comparing the assigned and true positions.

RESULTS

In the 20-fold cross-validation, the mean proportion of masked mapped SNPs that were placed by the algorithm to a chromosome was 89.53, 94.25, and 97.23% for sugarcane, wheat, and barley, respectively. Of the markers that were placed in the genome, 98.73, 96.45 and 98.53% of the SNPs were positioned on the correct chromosome. The mean correlations between known and new estimated SNP positions were 0.97, 0.98, and 0.97 for sugarcane, wheat, and barley. The LD-based algorithm was used to assign 5920 out of 21,251 unpositioned markers to the current Q208 sugarcane genetic map, representing the highest density genetic map for this species to date.

CONCLUSIONS

Our LD-based approach can be used to accurately assign unpositioned SNPs to existing genetic maps, improving genome-wide association studies and genomic prediction in crop species with fragmented and incomplete genome assemblies. This approach will facilitate genomic-assisted breeding for many orphan crops that lack genetic and genomic resources.

Shape, firmness and fruit quality QTLs shared in two non-related strawberry populations

The cultivated strawberry (Fragaria x ananassa) is an octoploid species (2n = 8x = 56), appreciated widely for its fruit. There have been very few studies on fruit quality traits, which are known to be mostly polygenic and environmentally dependent. To identify higher genetic variability, two non-related populations were genotyped: an F1 population cross between 'FC50' and 'FD54' and an F2 population cross between 'Camarosa' and 'Dover', hybridizing both with IStraw35k and IStraw90k SNP arrays, respectively. The F1 genetic map was constructed with 14595 SNPs and the F2 map with 7977 SNPs. High collinearity was observed when comparing one genetic map with the other and on comparing both with the octoploid genome. To assess fruit variability, both populations were phenotyped for shape, firmness, taste and other fruit traits over the 2016-2019 period. With QTL analyses, 33 stable QTLs were mapped in the 'FC50xFD54' population, and three hotspot regions were found for shape traits in LG3A, LG4D and LG6D. In the '21AF' population, only eight stable QTLs were detected. Despite that, two major and stable QTLs were mapped in the same interval of confidence for both populations. A shared fruit shape ratio QTL which explained around 25 % of trait variance was mapped in LG3A, and a shared firmness QTL explaining 26.9 % of trait variance in LG7C. For the first time, two QTLs were discovered in LG3A and LG4A for a phenotype neck without achenes. When analysing two different mapping populations, in addition to finding specific QTL regions for the studied traits, a narrowing down of the interval of confidence for the shared QTLs is achieved. As a result of this study, a new set of SNPs for fruit firmness and shape is now available for use in MAS in strawberry breeding programs.

Genetic mapping of adult-plant resistance genes to powdery mildew in triticale

Triticale is a cereal of high economic importance; however, along with the increase in the area of this cereal, it is more often infected by the fungal pathogen Blumeria graminis, which causes powdery mildew. The rapid development of molecular biology techniques, in particular methods based on molecular markers may be an important tool used in modern plant breeding. Development of genetic maps, location of the QTLs defining the region of the genome associated with resistance and selection of markers linked to particular trait can be used to select resistant genotypes as well as to pyramidize several resistance genes in one variety. In this paper, we present a new, high-density genetic map of triticale doubled haploids (DH) population “Grenado” × “Zorro” composed of DArT, silicoDArT, and SNP markers. Composite interval mapping method was used to detect eight QTL regions associated with the area under disease progress curve (AUDPC) and 15 regions with the average value of powdery mildew infection (avPM) based on observation conducted in 3-year period in three different locations across the Poland. Two regions on rye chromosome 4R, and single loci on 5R and 6R were reported for the first time as regions associated with powdery mildew resistance. Among all QTLs, 14 candidate genes were identified coded cyclin-dependent kinase, serine/threonine-protein kinase-like protein as well as AMEIOTIC 1 homolog DYAD-like protein, DETOXIFICATION 16-like protein, and putative disease resistance protein RGA3. Three of identified candidate genes were found among newly described QTL regions associated with powdery mildew resistance in triticale.

Quantitative trait loci and candidate genes associated with freezing tolerance of winter triticale (× Triticosecale Wittmack)

Freezing tolerance of triticale is a major trait contributing to its winter hardiness. The identification of genomic regions — quantitative trait loci (QTL) and molecular markers associated with freezing tolerance in winter hexaploid triticale — was the aim of this study. For that purpose, a new genetic linkage map was developed for the population of 92 doubled haploid lines derived from ‘Hewo’ × ‘Magnat’ F₁ hybrid. Those lines, together with parents were subjected to freezing tolerance test three times during two winter seasons. Plants were grown and cold-hardened under natural fall/winter conditions and then subjected to freezing in controlled conditions. Freezing tolerance was assessed as the plants recovery (REC), the electrolyte leakage (EL) from leaves and chlorophyll fluorescence parameters (JIP) after freezing. Three consistent QTL for several fluorescence parameters, electrolyte leakage, and the percentage of the survived plants were identified with composite interval mapping (CIM) and single marker analysis (SMA). The first locus Qfr.hm-7A.1 explained 9% of variation of both electrolyte leakage and plants recovery after freezing. Two QTL explaining up to 12% of variation in plants recovery and shared by selected chlorophyll fluorescence parameters were found on 4R and 5R chromosomes. Finally, main locus Qchl.hm-5A.1 was detected for chlorophyll fluorescence parameters that explained up to 19.6% of phenotypic variation. The co-located QTL on chromosomes 7A.1, 4R and 5R, clearly indicated physiological and genetic relationship of the plant survival after freezing with the ability to maintain optimal photochemical activity of the photosystem II and preservation of the cell membranes integrity. The genes located in silico within the identified QTL include those encoding BTR1-like protein, transmembrane helix proteins like potassium channel, and phosphoric ester hydrolase involved in response to osmotic stress as well as proteins involved in the regulation of the gene expression, chloroplast RNA processing, and pyrimidine salvage pathway. Additionally, our results confirm that the JIP test is a valuable tool to evaluate freezing tolerance of triticale under unstable winter environments.

A consensus and saturated genetic map provides insight into genome anchoring, synteny of Solanaceae and leaf- and fruit-related QTLs in wolfberry (Lycium Linn.)

BACKGROUND

Lycium Linn. (Solanaceae) is a genus of economically important plants producing fruits and leaves with high nutritional value and medicinal benefits. However, genetic analysis of this plant and molecular breeding for quality improvement are limited by the lack of sufficient molecular markers.

RESULTS

In this study, two parental strains, 'Ningqi No. 1' (Lycium barbarum L.) and 'Yunnan Gouqi' (Lycium yunnanense Kuang et A.M. Lu), and 200 F₁ hybrid individuals were resequenced for genetic analysis. In total, 8,507 well-selected SNPs were developed, and a high-density genetic map (NY map) was constructed with a total genetic distance of 2,122.24 cM. A consensus genetic map was established by integrating the NY map and a previously published genetic map (NC map) containing 15,240 SNPs, with a total genetic distance of 3,058.19 cM and an average map distance of 0.21 cM. The 12 pseudochromosomes of the Lycium reference genome were anchored using this consensus genetic map, with an anchoring rate of 64.3%. Moreover, weak collinearities between the consensus map and the pepper, potato, and tomato genomes were observed. Twenty-five stable QTLs were identified for leaf- and fruit-related phenotypes, including fruit weight, fruit longitude, leaf length, the fruit index, and the leaf index; these stable QTLs were mapped to four different linkage groups, with LOD scores ranging from 2.51 to 19.37 and amounts of phenotypic variance explained from 6.2% to 51.9%. Finally, 82 out of 188 predicted genes underlying stable QTLs for fruit-related traits were differentially expressed according to RNA-seq analysis.

CONCLUSIONS

A chromosome-level assembly can provide a foundation for further functional genomics research for wolfberry. The genomic regions of these stably expressed QTLs could be used as targets for further fine mapping and development of molecular markers for marker-assisted selection (MAS). The present study provided valuable information on saturated SNP markers and reliable QTLs for map-based cloning of functional genes related to yield and morphological traits in Lycium spp.

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.

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.

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.

High-density SNP linkage map construction and QTL mapping for flavonoid-related traits in a tea plant (Camellia sinensis) using 2b-RAD sequencing

BACKGROUND

Flavonoids are important components that confer upon tea plants a unique flavour and health functions. However, the traditional breeding method for selecting a cultivar with a high or unique flavonoid content is time consuming and labour intensive. High-density genetic map construction associated with quantitative trait locus (QTL) mapping provides an effective way to facilitate trait improvement in plant breeding. In this study, an F1 population (LJ43×BHZ) was genotyped using 2b-restriction site-associated DNA (2b-RAD) sequencing to obtain massive single nucleotide polymorphism (SNP) markers to construct a high-density genetic map for a tea plant. Furthermore, QTLs related to flavonoids were identified using our new genetic map.

RESULTS

A total of 13,446 polymorphic SNP markers were developed using 2b-RAD sequencing, and 4,463 of these markers were available for constructing the genetic linkage map. A 1,678.52-cM high-density map at an average interval of 0.40 cM with 4,217 markers, including 427 frameset simple sequence repeats (SSRs) and 3,800 novel SNPs, mapped into 15 linkage groups was successfully constructed. After QTL analysis, a total of 27 QTLs related to flavonoids or caffeine content (CAF) were mapped to 8 different linkage groups, LG01, LG03, LG06, LG08, LG10, LG11, LG12, and LG13, with an LOD from 3.14 to 39.54, constituting 7.5% to 42.8% of the phenotypic variation.

CONCLUSIONS

To our knowledge, the highest density genetic map ever reported was constructed since the largest mapping population of tea plants was adopted in present study. Moreover, novel QTLs related to flavonoids and CAF were identified based on the new high-density genetic map. In addition, two markers were located in candidate genes that may be involved in flavonoid metabolism. The present study provides valuable information for gene discovery, marker-assisted selection breeding and map-based cloning for functional genes that are related to flavonoid content in tea plants.

Construction of a highly saturated Genetic Map for Vitis by Next-generation Restriction Site-associated DNA Sequencing

BACKGROUND

High-saturate molecular linkage maps are an important tool in studies on plant molecular biology and assisted breeding. Development of a large set of single nucleotide polymorphisms (SNPs) via next-generation sequencing (NGS)-based methods, restriction-site associated DNA sequencing (RAD-seq), and the generation of a highly saturated genetic map help improve fine mapping of quantitative trait loci (QTL).

RESULTS

We generated a highly saturated genetic map to identify significant traits in two elite grape cultivars and 176 F₁ plants. In total, 1,426,967 high-quality restriction site-associated DNA tags were detected; 51,365, 23,683, and 70,061 markers were assessed in 19 linkage groups (LGs) for the maternal, paternal, and integrated maps, respectively. Our map was highly saturated in terms of marker density and average "Gap ≤ 5 cM" percentage.

CONCLUSIONS

In this study, RAD-seq of 176 F₁ plants and their parents yielded 8,481,484 SNPs and 1,646,131 InDel markers, of which 65,229 and 4832, respectively, were used to construct a highly saturated genetic map for grapevine. This map is expected to facilitate genetic studies on grapevine, including an evaluation of grapevine and deciphering the genetic basis of economically and agronomically important traits. Our findings provide basic essential genetic data the grapevine genetic research community, which will lead to improvements in grapevine breeding.