Accurate construction of consensus genetic maps via integer linear programming

Construction of a high-density linkage map and QTL detection for growth traits in South African abalone (Haliotis midae)

Haliotis midae is one of the most important molluscs in South African commercial aquaculture. In this study, a high-resolution integrated linkage map was constructed, and QTL identified using 2b-RADseq for genotyping SNPs in three families. The final integrated linkage map was composed by merging the individual family maps, resulting in 3290 informative SNPs mapping to 18 linkage groups, conforming to the known haploid chromosome number for H. midae. The total map spanned 1798.25 cM with an average marker interval of 0.55 cM, representing a genome coverage of 98.76%. QTL analysis, across all three families, resulted in a total of five QTL identified for growth-related traits, shell width, shell length, and total body weight. For shell width and total body weight, one QTL was identified for each trait respectively, whilst three QTL were identified for shell length. The identified QTL respectively explained between 7.20% and 11.40% of the observed phenotypic variance. All three traits were significantly correlated (r = 0.862-0.970; p < 0.01) and shared overlapping QTL. The QTL for growth traits were mapped back to the H. midae draft genome and BLAST searches revealed the identity of candidate genes, such as egf-1, megf10, megf6, tnx, sevp1, kcp, notch1, and scube2 with possible functional roles in H. midae growth. The constructed high-density linkage map and mapped QTL have given valuable insights regarding the genetic architecture of growth-related traits and will be important genetic resources for marker-assisted selection. It remains, however, important to validate causal variants through linkage disequilibrium fine mapping in future.

Meta-QTL analysis in wheat: progress, challenges and opportunities

Wheat, an important cereal crop globally, faces major challenges due to increasing global population and changing climates. The production and productivity are challenged by several biotic and abiotic stresses. There is also a pressing demand to enhance grain yield and quality/nutrition to ensure global food and nutritional security. To address these multifaceted concerns, researchers have conducted numerous meta-QTL (MQTL) studies in wheat, resulting in the identification of candidate genes that govern these complex quantitative traits. MQTL analysis has successfully unraveled the complex genetic architecture of polygenic quantitative traits in wheat. Candidate genes associated with stress adaptation have been pinpointed for abiotic and biotic traits, facilitating targeted breeding efforts to enhance stress tolerance. Furthermore, high-confidence candidate genes (CGs) and flanking markers to MQTLs will help in marker-assisted breeding programs aimed at enhancing stress tolerance, yield, quality and nutrition. Functional analysis of these CGs can enhance our understanding of intricate trait-related genetics. The discovery of orthologous MQTLs shared between wheat and other crops sheds light on common evolutionary pathways governing these traits. Breeders can leverage the most promising MQTLs and CGs associated with multiple traits to develop superior next-generation wheat cultivars with improved trait performance. This review provides a comprehensive overview of MQTL analysis in wheat, highlighting progress, challenges, validation methods and future opportunities in wheat genetics and breeding, contributing to global food security and sustainable agriculture.

A chromosome-level genome assembly enables the identification of the follicule stimulating hormone receptor as the master sex-determining gene in the flatfish Solea senegalensis

Sex determination (SD) shows huge variation among fish and a high evolutionary rate, as illustrated by the Pleuronectiformes (flatfishes). This order is characterized by its adaptation to demersal life, compact genomes and diversity of SD mechanisms. Here, we assembled the Solea senegalensis genome, a flatfish of great commercial value, into 82 contigs (614 Mb) combining long- and short-read sequencing, which were next scaffolded using a highly dense genetic map (28,838 markers, 21 linkage groups), representing 98.9% of the assembly. Further, we established the correspondence between the assembly and the 21 chromosomes by using BAC-FISH. Whole genome resequencing of six males and six females enabled the identification of 41 single nucleotide polymorphism variants in the follicle stimulating hormone receptor (fshr) consistent with an XX/XY SD system. The observed sex association was validated in a broader independent sample, providing a novel molecular sexing tool. The fshr gene displayed differential expression between male and female gonads from 86 days post-fertilization, when the gonad is still an undifferentiated primordium, concomitant with the activation of amh and cyp19a1a, testis and ovary marker genes, respectively, in males and females. The Y-linked fshr allele, which included 24 nonsynonymous variants and showed a highly divergent 3D protein structure, was overexpressed in males compared to the X-linked allele at all stages of gonadal differentiation. We hypothesize a mechanism hampering the action of the follicle stimulating hormone driving the undifferentiated gonad toward testis.

A high-density integrated map for grapevine based on three mapping populations genotyped by the Vitis18K SNP chip

We present a high-density integrated map for grapevine, allowing refinement and improved understanding of the grapevine genome, while demonstrating the applicability of the Vitis18K SNP chip for linkage mapping. The improvement of grapevine through biotechnology requires identification of the molecular bases of target traits by studying marker-trait associations. The Vitis18K SNP chip provides a useful genotyping tool for genome-wide marker analysis. Most linkage maps are based on single mapping populations, but an integrated map can increase marker density and show order conservation. Here we present an integrated map based on three mapping populations. The parents consist of the well-known wine cultivars 'Cabernet Sauvignon', 'Corvina' and 'Rhine Riesling', the lesser-known wine variety 'Deckrot', and a table grape selection, G1-7720. Three high-density population maps with an average inter-locus gap ranging from 0.74 to 0.99 cM were developed. These maps show high correlations (0.9965-0.9971) with the reference assembly, containing only 93 markers with large order discrepancies compared to expected physical positions, of which a third is consistent across multiple populations. Moreover, the genetic data aid the further refinement of the grapevine genome assembly, by anchoring 104 yet unanchored scaffolds. From these population maps, an integrated map was constructed which includes 6697 molecular markers and reduces the inter-locus gap distance to 0.60 cM, resulting in the densest integrated map for grapevine thus far. A small number of discrepancies, mainly of short distance, involve 88 markers that remain conflictual across maps. The integrated map shows similar collinearity to the reference assembly (0.9974) as the single maps. This high-density map increases our understanding of the grapevine genome and provides a useful tool for its further characterization and the dissection of complex traits.

The first high-density genetic map of common cockle (Cerastoderma edule) reveals a major QTL controlling shell color variation

Shell color shows broad variation within mollusc species and despite information on the genetic pathways involved in shell construction and color has recently increased, more studies are needed to understand its genetic architecture. The common cockle (Cerastoderma edule) is a valuable species from ecological and commercial perspectives which shows important variation in shell color across Northeast Atlantic. In this study, we constructed a high-density genetic map, as a tool for screening common cockle genome, which was applied to ascertain the genetic basis of color variation in the species. The consensus genetic map comprised 19 linkage groups (LGs) in accordance with the cockle karyotype (2n = 38) and spanned 1073 cM, including 730 markers per LG and an inter-marker distance of 0.13 cM. Five full-sib families showing segregation for several color-associated traits were used for a genome-wide association study and a major QTL on chromosome 13 associated to different color-traits was detected. Mining on this genomic region revealed several candidate genes related to shell construction and color. A genomic region previously reported associated with divergent selection in cockle distribution overlapped with this QTL suggesting its putative role on adaptation.

Current status and trends in forest genomics

Forests are not only the most predominant of the Earth's terrestrial ecosystems, but are also the core supply for essential products for human use. However, global climate change and ongoing population explosion severely threatens the health of the forest ecosystem and aggravtes the deforestation and forest degradation. Forest genomics has great potential of increasing forest productivity and adaptation to the changing climate. In the last two decades, the field of forest genomics has advanced quickly owing to the advent of multiple high-throughput sequencing technologies, single cell RNA-seq, clustered regularly interspaced short palindromic repeats (CRISPR)-mediated genome editing, and spatial transcriptomes, as well as bioinformatics analysis technologies, which have led to the generation of multidimensional, multilayered, and spatiotemporal gene expression data. These technologies, together with basic technologies routinely used in plant biotechnology, enable us to tackle many important or unique issues in forest biology, and provide a panoramic view and an integrative elucidation of molecular regulatory mechanisms underlying phenotypic changes and variations. In this review, we recapitulated the advancement and current status of 12 research branches of forest genomics, and then provided future research directions and focuses for each area. Evidently, a shift from simple biotechnology-based research to advanced and integrative genomics research, and a setup for investigation and interpretation of many spatiotemporal development and differentiation issues in forest genomics have just begun to emerge.

Linkage mapping and QTL analysis of growth traits in Rhopilema esculentum

R. esculentum is a popular seafood in Asian countries and an economic marine fishery resource in China. However, the genetic linkage map and growth-related molecular markers are still lacking, hindering marker assisted selection (MAS) for genetic improvement of R. esculentum. Therefore, we firstly used 2b-restriction site-associated DNA (2b-RAD) method to sequence 152 R. esculentum specimens and obtained 9100 single nucleotide polymorphism (SNP) markers. A 1456.34 cM linkage map was constructed using 2508 SNP markers with an average interval of 0.58 cM. Then, six quantitative trait loci (QTLs) for umbrella diameter and body weight were detected by QTL analysis based on the new linkage map. The six QTLs are located on four linkage groups (LGs), LG4, LG13, LG14 and LG15, explaining 9.4% to 13.4% of the phenotypic variation. Finally, 27 candidate genes in QTLs regions of LG 14 and 15 were found associated with growth and one gene named RE13670 (sushi, von Willebrand factor type A, EGF and pentraxin domain-containing protein 1-like) may play an important role in controlling the growth of R. esculentum. This study provides valuable information for investigating the growth mechanism and MAS breeding in R. esculentum.

High-density genetic linkage mapping reveals low stability of QTLs across environments for economic traits in Eucalyptus

Introduction

Eucalyptus urophylla, E. tereticornis and their hybrids are the most important commercial forest tree species in South China where they are grown for pulpwood and solid wood production. Construction of a fine-scale genetic linkage map and detecting quantitative trait loci (QTL) for economically important traits linked to these end-uses will facilitate identification of the main candidate genes and elucidate the regulatory mechanisms.

Method

A high-density consensus map (a total of 2754 SNPs with 1359.18 cM) was constructed using genotyping by sequencing (GBS) on clonal progenies of E. urophylla × tereticornis hybrids. QTL mapping of growth and wood property traits were conducted in three common garden experiments, resulting in a total of 108 QTLs. A total of 1052 candidate genes were screened by the efficient combination of QTL mapping and transcriptome analysis.

Results

Only ten QTLs were found to be stable across two environments, and only one (qSG10Stable mapped on chromosome 10, and associated with lignin syringyl-to-guaiacyl ratio) was stable across all three environments. Compared to other QTLs, qSG10Stable explained a very high level of phenotypic variation (18.4-23.6%), perhaps suggesting that QTLs with strong effects may be more stably inherited across multiple environments. Screened candidate genes were associated with some transcription factor families, such as TALE, which play an important role in the secondary growth of plant cell walls and the regulation of wood formation.

Discussion

While QTLs such as qSG10Stable, found to be stable across three sites, appear to be comparatively uncommon, their identification is likely to be a key to practical QTL-based breeding. Further research involving clonally-replicated populations, deployed across multiple target planting sites, will be required to further elucidate QTL-by-environment interactions.

Construction of Consensus Genetic Map With Applications in Gene Mapping of Wheat (Triticum aestivum L.) Using 90K SNP Array

Wheat is one of the most important cereal crops worldwide. A consensus map combines genetic information from multiple populations, providing an effective alternative to improve the genome coverage and marker density. In this study, we constructed a consensus map from three populations of recombinant inbred lines (RILs) of wheat using a 90K single nucleotide polymorphism (SNP) array. Phenotypic data on plant height (PH), spike length (SL), and thousand-kernel weight (TKW) was collected in six, four, and four environments in the three populations, and then used for quantitative trait locus (QTL) mapping. The mapping results obtained using the constructed consensus map were compared with previous results obtained using individual maps and previous studies on other populations. A simulation experiment was also conducted to assess the performance of QTL mapping with the consensus map. The constructed consensus map from the three populations spanned 4558.55 cM in length, with 25,667 SNPs, having high collinearity with physical map and individual maps. Based on the consensus map, 21, 27, and 19 stable QTLs were identified for PH, SL, and TKW, much more than those detected with individual maps. Four PH QTLs and six SL QTLs were likely to be novel. A putative gene called TraesCS4D02G076400 encoding gibberellin-regulated protein was identified to be the candidate gene for one major PH QTL located on 4DS, which may enrich genetic resources in wheat semi-dwarfing breeding. The simulation results indicated that the length of the confidence interval and standard errors of the QTLs detected using the consensus map were much smaller than those detected using individual maps. The consensus map constructed in this study provides the underlying genetic information for systematic mapping, comparison, and clustering of QTL, and gene discovery in wheat genetic study. The QTLs detected in this study had stable effects across environments and can be used to improve the wide adaptation of wheat cultivars through marker-assisted breeding.

Genome survey and high-resolution genetic map provide valuable genetic resources for Fenneropenaeus chinensis

Fenneropenaeus chinensis is one of the most important aquaculture species in China. Research on its genomic and genetic structure not only helps us comprehend the genetic basis of complex economic traits, but also offers theoretical guidance in selective breeding. In the present study, a genome survey sequencing was performed to generate a rough reference genome utilized for groping preliminary genome characteristics and facilitate linkage and quantitative trait locus (QTL) mapping. Linkage mapping was conducted using a reduced-representation sequencing method 2b-RAD. In total, 36,762 SNPs were genotyped from 273 progenies in a mapping family, and a high-resolution linkage map was constructed. The consensus map contained 12,884 markers and spanned 5257.81 cM with an average marker interval of 0.41 cM, which was the first high-resolution genetic map in F. chinensis to our knowledge. QTL mapping and association analysis were carried out in 29 characters including body size, sex and disease resistance. 87 significant QTLs were detected in several traits and they were also evaluated by association analysis. Results of this study provide us valuable suggestions in genetic improvement and breeding of new varieties and also lay a basic foundation for further application of cloning of economic genes in selective breeding program and marker-assisted selection.

Construction of a High-Density Genetic Map of Acca sellowiana (Berg.) Burret, an Outcrossing Species, Based on Two Connected Mapping Populations

Acca sellowiana, known as feijoa or pineapple guava, is a diploid, (2n = 2x = 22) outcrossing fruit tree species native to Uruguay and Brazil. The species stands out for its highly aromatic fruits, with nutraceutical and therapeutic value. Despite its promising agronomical value, genetic studies on this species are limited. Linkage genetic maps are valuable tools for genetic and genomic studies, and constitute essential tools in breeding programs to support the development of molecular breeding strategies. A high-density composite genetic linkage map of A. sellowiana was constructed using two genetically connected populations: H5 (TCO × BR, N = 160) and H6 (TCO × DP, N = 184). Genotyping by sequencing (GBS) approach was successfully applied for developing single nucleotide polymorphism (SNP) markers. A total of 4,921 SNP markers were identified using the reference genome of the closely related species Eucalyptus grandis, whereas other 4,656 SNPs were discovered using a de novo pipeline. The individual H5 and H6 maps comprised 1,236 and 1,302 markers distributed over the expected 11 linkage groups, respectively. These two maps spanned a map length of 1,593 and 1,572 cM, with an average inter-marker distance of 1.29 and 1.21 cM, respectively. A large proportion of markers were common to both maps and showed a high degree of collinearity. The composite map consisted of 1,897 SNPs markers with a total map length of 1,314 cM and an average inter-marker distance of 0.69. A novel approach for the construction of composite maps where the meiosis information of individuals of two connected populations is captured in a single estimator is described. A high-density, accurate composite map based on a consensus ordering of markers provides a valuable contribution for future genetic research and breeding efforts in A. sellowiana. A novel mapping approach based on an estimation of multipopulation recombination fraction described here may be applied in the construction of dense composite genetic maps for any other outcrossing diploid species.

Genetic architecture of agronomic and quality traits in a nested association mapping population of spring wheat

Germplasm collections are rich sources of genetic variation to improve crops for many valuable traits. Nested association mapping (NAM) populations can overcome the limitations of genome-wide association studies (GWAS) in germplasm collections by reducing the effect of population structure. We exploited the genetic diversity of the USDA-ARS wheat (Triticum aestivum L.) core collection by developing the Spring Wheat Multiparent Introgression Population (SWMIP). To develop this population, twenty-five core parents were crossed and backcrossed to the Minnesota spring wheat cultivar RB07. The NAM population and 26 founder parents were genotyped using genotyping-by-sequencing and phenotyped for heading date, height, test weight, and grain protein content. After quality control, 20,312 markers with physical map positions were generated for 2,038 recombinant inbred lines (RILs). The number of RILs in each family varied between 58 and 96. Three GWAS models were utilized for quantitative trait loci (QTL) detection and accounted for known family stratification, genetic kinship, and both covariates. GWAS was performed on the whole population and also by bootstrap sampling of an equal number of RILs from each family. Greater power of QTL detection was achieved by treating families equally through bootstrapping. In total 16, 15, 12, and 13 marker-trait associations (MTAs) were identified for heading date, height, test weight, and grain protein content, respectively. Some of these MTAs were coincident with major genes known to control the traits, but others were novel and contributed by the wheat core parents. The SWMIP will be a valuable source of genetic variation for spring wheat breeding.

Association Mapping and Development of Marker-Assisted Selection Tools for the Resistance to White Pine Blister Rust in the Alberta Limber Pine Populations

Since its introduction to North America in the early 1900s, white pine blister rust (WPBR) caused by the fungal pathogen Cronartium ribicola has resulted in substantial economic losses and ecological damage to native North American five-needle pine species. The high susceptibility and mortality of these species, including limber pine (Pinus flexilis), creates an urgent need for the development and deployment of resistant germplasm to support recovery of impacted populations. Extensive screening for genetic resistance to WPBR has been underway for decades in some species but has only started recently in limber pine using seed families collected from wild parental trees in the USA and Canada. This study was conducted to characterize Alberta limber pine seed families for WPBR resistance and to develop reliable molecular tools for marker-assisted selection (MAS). Open-pollinated seed families were evaluated for host reaction following controlled infection using C. ribicola basidiospores. Phenotypic segregation for presence/absence of stem symptoms was observed in four seed families. The segregation ratios of these families were consistent with expression of major gene resistance (MGR) controlled by a dominant R locus. Based on linkage disequilibrium (LD)-based association mapping used to detect single nucleotide polymorphism (SNP) markers associated with MGR against C. ribicola, MGR in these seed families appears to be controlled by Cr4 or other R genes in very close proximity to Cr4. These associated SNPs were located in genes involved in multiple molecular mechanisms potentially underlying limber pine MGR to C. ribicola, including NBS-LRR genes for recognition of C. ribicola effectors, signaling components, and a large set of defense-responsive genes with potential functions in plant effector-triggered immunity (ETI). Interactions of associated loci were identified for MGR selection in trees with complex genetic backgrounds. SNPs with tight Cr4-linkage were further converted to TaqMan assays to confirm their effectiveness as MAS tools. This work demonstrates the successful translation and deployment of molecular genetic knowledge into specific MAS tools that can be easily applied in a selection or breeding program to efficiently screen MGR against WPBR in Alberta limber pine populations.

Construction of a high-density genetic map and mapping of growth related QTLs in the grass carp (Ctenopharyngodon idellus)

BACKGROUND

Grass carp (Ctenopharyngodon idellus) are important species in Asian aquaculture. A draft genome for grass carp has already been published in 2015. However, there is still a requirement for a suitable genetic linkage map to arrange scaffolds on chromosomal frameworks. QTL analysis is a powerful tool to detect key locations for quantitative traits, especially in aquaculture. There no growth related QTLs of grass carp have been published yet. Even the growth trait is one of the focuses in grass carp culture.

RESULTS

In this study, a pair of distantly related parent grass carps and their 100 six-month-old full-sib offspring were used to construct a high-density genetic map with 6429 single nucleotide polymorphisms (SNPs) by 2b-RAD technology. The total length of the consensus map is 5553.43 cM with the average marker interval of 1.92 cM. The map has a good collinearity with both the grass carp draft genome and the zebrafish genome, and it assembled 89.91% of the draft genome to a chromosomal level. Additionally, according to the growth-related traits of progenies, 30 quantitative trait loci (QTLs), including 7 for body weight, 9 for body length, 5 for body height and 9 for total length, were identified in 16 locations on 5 linkage groups. The phenotypic variance explained for these QTLs varies from 13.4 to 21.6%. Finally, 17 genes located in these regions were considered to be growth-related because they either had functional mutations predicted from the resequencing data of the parents.

CONCLUSION

A high density genetic linkage map of grass carp was built and it assembled the draft genome to a chromosomal level. Thirty growth related QTLs were detected. After the cross analysis of Parents resequencing data, 17 candidate genes were obtained for further researches.

A novel high-density grapevine (Vitis vinifera L.) integrated linkage map using GBS in a half-diallel population

A half-diallel population involving five elite grapevine cultivars was generated and genotyped by GBS, and highly-informative segregation data was used to construct a high-density genetic map for Vitis vinifera L. Grapevine is one of the most relevant fruit crops in the world. Deeper genetic knowledge could assist modern grapevine breeding programs to develop new wine grape varieties able to face climate change effects. To assist in the rapid identification of markers for crop yield components, grape quality traits and adaptation potential, we generated a large Vitis vinifera L. population (N = 624) by crossing five red wine cultivars in a half-diallel scheme, which was subsequently sequenced by an efficient GBS procedure. A high number of fully informative genetic variants was detected using a novel mapping approach capable of reconstructing local haplotypes from adjacent biallelic SNPs, which were subsequently used to construct the densest consensus genetic map available for the cultivated grapevine to date. This 1378.3-cM map integrates 10 bi-parental consensus maps and orders 4437 markers in 3353 unique positions on 19 chromosomes. Markers are well distributed all along the grapevine reference genome, covering up to 98.8% of its genomic sequence. Additionally, a good agreement was observed between genetic and physical orders, adding confidence in the quality of this map. Collectively, our results pave the way for future genetic studies (such as fine QTL mapping) aimed to understand the complex relationship between genotypic and phenotypic variation in the cultivated grapevine. In addition, the method used (which efficiently delivers a high number of fully informative markers) could be of interest to other outbred organisms, notably perennial fruit crops.

Discovering consensus genomic regions in wheat for root-related traits by QTL meta-analysis

Root system architecture is crucial for wheat adaptation to drought stress, but phenotyping for root traits in breeding programmes is difficult and time-consuming owing to the belowground characteristics of the system. Identifying quantitative trait loci (QTLs) and linked molecular markers and using marker-assisted selection is an efficient way to increase selection efficiency and boost genetic gains in breeding programmes. Hundreds of QTLs have been identified for different root traits in the last few years. In the current study, consensus QTL regions were identified through QTL meta-analysis. First, a consensus map comprising 7352 markers was constructed. For the meta-analysis, 754 QTLs were retrieved from the literature and 634 of them were projected onto the consensus map. Meta-analysis grouped 557 QTLs in 94 consensus QTL regions, or meta-QTLs (MQTLs), and 18 QTLs remained as singletons. The recently published genome sequence of wheat was used to search for gene models within the MQTL peaks. As a result, gene models for 68 of the 94 Root_MQTLs were found, 35 of them related to root architecture and/or drought stress response. This work will facilitate QTL cloning and pyramiding to develop new cultivars with specific root architecture for coping with environmental constraints.

A physical and genetic map of Cannabis sativa identifies extensive rearrangements at the THC/CBD acid synthase loci

Cannabis sativa is widely cultivated for medicinal, food, industrial, and recreational use, but much remains unknown regarding its genetics, including the molecular determinants of cannabinoid content. Here, we describe a combined physical and genetic map derived from a cross between the drug-type strain Purple Kush and the hemp variety "Finola." The map reveals that cannabinoid biosynthesis genes are generally unlinked but that aromatic prenyltransferase (AP), which produces the substrate for THCA and CBDA synthases (THCAS and CBDAS), is tightly linked to a known marker for total cannabinoid content. We further identify the gene encoding CBCA synthase (CBCAS) and characterize its catalytic activity, providing insight into how cannabinoid diversity arises in cannabis. THCAS and CBDAS (which determine the drug vs. hemp chemotype) are contained within large (>250 kb) retrotransposon-rich regions that are highly nonhomologous between drug- and hemp-type alleles and are furthermore embedded within ∼40 Mb of minimally recombining repetitive DNA. The chromosome structures are similar to those in grains such as wheat, with recombination focused in gene-rich, repeat-depleted regions near chromosome ends. The physical and genetic map should facilitate further dissection of genetic and molecular mechanisms in this commercially and medically important plant.

Efficient algorithms for Longest Common Subsequence of two bucket orders to speed up pairwise genetic map comparison

Genetic maps order genetic markers along chromosomes. They are, for instance, extensively used in marker-assisted selection to accelerate breeding programs. Even for the same species, people often have to deal with several alternative maps obtained using different ordering methods or different datasets, e.g. resulting from different segregating populations. Having efficient tools to identify the consistency and discrepancy of alternative maps is thus essential to facilitate genetic map comparisons. We propose to encode genetic maps by bucket order, a kind of order, which takes into account the blurred parts of the marker order while being an efficient data structure to achieve low complexity algorithms. The main result of this paper is an O(n log(n)) procedure to identify the largest agreements between two bucket orders of n elements, their Longest Common Subsequence (LCS), providing an efficient solution to highlight discrepancies between two genetic maps. The LCS of two maps, being the largest set of their collinear markers, is used as a building block to compute pairwise map congruence, to visually emphasize maker collinearity and in some scaffolding methods relying on genetic maps to improve genome assembly. As the LCS computation is a key subroutine of all these genetic map related tools, replacing the current LCS subroutine of those methods by ours -to do the exact same work but faster- could significantly speed up those methods without changing their accuracy. To ease such transition we provide all required algorithmic details in this self contained paper as well as an R package implementing them, named LCSLCIS, which is freely available at: https://github.com/holtzy/LCSLCIS.

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.

High-density multi-population consensus genetic linkage map for peach

Highly saturated genetic linkage maps are extremely helpful to breeders and are an essential prerequisite for many biological applications such as the identification of marker-trait associations, mapping quantitative trait loci (QTL), candidate gene identification, development of molecular markers for marker-assisted selection (MAS) and comparative genetic studies. Several high-density genetic maps, constructed using the 9K SNP peach array, are available for peach. However, each of these maps is based on a single mapping population and has limited use for QTL discovery and comparative studies. A consensus genetic linkage map developed from multiple populations provides not only a higher marker density and a greater genome coverage when compared to the individual maps, but also serves as a valuable tool for estimating genetic positions of unmapped markers. In this study, a previously developed linkage map from the cross between two peach cultivars 'Zin Dai' and 'Crimson Lady' (ZC2) was improved by genotyping additional progenies. In addition, a peach consensus map was developed based on the combination of the improved ZC2 genetic linkage map with three existing high-density genetic maps of peach and a reference map of Prunus. A total of 1,476 SNPs representing 351 unique marker positions were mapped across eight linkage groups on the ZC2 genetic map. The ZC2 linkage map spans 483.3 cM with an average distance between markers of 1.38 cM/marker. The MergeMap and LPmerge tools were used for the construction of a consensus map based on markers shared across five genetic linkage maps. The consensus linkage map contains a total of 3,092 molecular markers, consisting of 2,975 SNPs, 116 SSRs and 1 morphological marker associated with slow ripening in peach (SR). The consensus map provides valuable information on marker order and genetic position for QTL identification in peach and other genetic studies within Prunus and Rosaceae.

A physical and genetic map of Cannabis sativa identifies extensive rearrangements at the THC/CBD acid synthase loci

Cannabis sativa is widely cultivated for medicinal, food, industrial, and recreational use, but much remains unknown regarding its genetics, including the molecular determinants of cannabinoid content. Here, we describe a combined physical and genetic map derived from a cross between the drug-type strain Purple Kush and the hemp variety “Finola.” The map reveals that cannabinoid biosynthesis genes are generally unlinked but that aromatic prenyltransferase (AP), which produces the substrate for THCA and CBDA synthases (THCAS and CBDAS), is tightly linked to a known marker for total cannabinoid content. We further identify the gene encoding CBCA synthase (CBCAS) and characterize its catalytic activity, providing insight into how cannabinoid diversity arises in cannabis. THCAS and CBDAS (which determine the drug vs. hemp chemotype) are contained within large (>250 kb) retrotransposon-rich regions that are highly nonhomologous between drug- and hemp-type alleles and are furthermore embedded within ∼40 Mb of minimally recombining repetitive DNA. The chromosome structures are similar to those in grains such as wheat, with recombination focused in gene-rich, repeat-depleted regions near chromosome ends. The physical and genetic map should facilitate further dissection of genetic and molecular mechanisms in this commercially and medically important plant.

Highly dense linkage maps from 31 full-sibling families of turbot (Scophthalmus maximus) provide insights into recombination patterns and chromosome rearrangements throughout a newly refined genome assembly

Highly dense linkage maps enable positioning thousands of landmarks useful for anchoring the whole genome and for analysing genome properties. Turbot is the most important cultured flatfish worldwide and breeding programs in the fifth generation of selection are targeted to improve growth rate, obtain disease resistant broodstock and understand sex determination to control sex ratio. Using a Restriction-site Associated DNA approach, we genotyped 18,214 single nucleotide polymorphism in 1,268 turbot individuals from 31 full-sibling families. Individual linkage maps were combined to obtain a male, female and species consensus maps. The turbot consensus map contained 11,845 markers distributed across 22 linkage groups representing a total normalised length of 3,753.9 cM. The turbot genome was anchored to this map, and scaffolds representing 96% of the assembly were ordered and oriented to obtain the expected 22 megascaffolds according to its karyotype. Recombination rate was lower in males, especially around centromeres, and pairwise comparison of 44 individual maps suggested chromosome polymorphism at specific genomic regions. Genome comparison across flatfish provided new evidence on karyotype reorganisations occurring across the evolution of this fish group.

Development and characterization of the first dsRNA-resistant insect population from western corn rootworm, Diabrotica virgifera virgifera LeConte

The use of dsRNA to control insect pests via the RNA interference (RNAi) pathway is being explored by researchers globally. However, with every new class of insect control compounds, the evolution of insect resistance needs to be considered, and understanding resistance mechanisms is essential in designing durable technologies and effective resistance management strategies. To gain insight into insect resistance to dsRNA, a field screen with subsequent laboratory selection was used to establish a population of DvSnf7 dsRNA-resistant western corn rootworm, Diabrotica virgifera virgifera, a major maize insect pest. WCR resistant to ingested DvSnf7 dsRNA had impaired luminal uptake and resistance was not DvSnf7 dsRNA-specific, as indicated by cross resistance to all other dsRNAs tested. No resistance to the Bacillus thuringiensis Cry3Bb1 protein was observed. DvSnf7 dsRNA resistance was inherited recessively, located on a single locus, and autosomal. Together these findings will provide insights for dsRNA deployment for insect pest control.

Populations of doubled haploids for genetic mapping in hexaploid winter triticale

To create a framework for genetic dissection of hexaploid triticale, six populations of doubled haploid (DH) lines were developed from pairwise hybrids of high-yielding winter triticale cultivars. The six populations comprise between 97 and 231 genotyped DH lines each, totaling 957 DH lines. A consensus genetic map spans 4593.9 cM is composed of 1576 unique DArT markers. The maps reveal several structural rearrangements in triticale genomes. In preliminary tests of the populations and maps, markers specific to wheat segments of the engineered rye chromosome 1R (RM1B) were identified. Example QTL mapping of days to heading in cv. Krakowiak revealed loci on chromosomes 2BL and 2R responsible for extended vernalization requirement, and candidate genes were identified. The material is available to all parties interested in triticale genetics.

Genetic mapping and comparative genomics to inform restoration enhancement and culture of southern flounder, Paralichthys lethostigma

BACKGROUND

Southern flounder, Paralichthys lethostigma, historically support a substantial fishery along the Atlantic and Gulf coasts of the southern United States. Low year-class strengths over the past few years in the western Gulf of Mexico have raised concern that spawning stocks may be overfished. Current management of the resource includes releasing hatchery-raised juveniles to restock bays and estuaries; additionally, there is a growing interest in the potential for commercial aquaculture of the species. Currently, genomic resources for southern flounder do not exist. Here, we used two hatchery-reared families and double-digest, restriction-site-associated DNA (ddRAD) sequencing to create a reduced-representation genomic library consisting of several thousand single nucleotide polymorphisms (SNPs) located throughout the genome.

RESULTS

The relative position of each SNP-containing locus was determined to create a high-density genetic map spanning the 24 linkage groups of the southern flounder genome. The consensus map was used to identify regions of shared synteny between southern flounder and seven other fish species for which genome assemblies are available. Finally, syntenic blocks were used to localize genes identified from transcripts in European flounder as potentially being involved in ecotoxicological and osmoregulatory responses, as well as QTLs associated with growth and disease resistance in Japanese flounder, on the southern flounder linkage map.

CONCLUSIONS

The information provided by the linkage map will enrich restoration efforts by providing a foundation for interpreting spatial genetic variation within the species, ultimately furthering an understanding of the adaptive potential and resilience of southern flounder to future changes in local environmental conditions. Further, the map will facilitate the use of genetic markers to enhance restoration and commercial aquaculture.

Construction of an ultra-high density consensus genetic map, and enhancement of the physical map from genome sequencing in Lupinus angustifolius

An ultra-high density genetic map containing 34,574 sequence-defined markers was developed in Lupinus angustifolius. Markers closely linked to nine genes of agronomic traits were identified. A physical map was improved to cover 560.5 Mb genome sequence. Lupin (Lupinus angustifolius L.) is a recently domesticated legume grain crop. In this study, we applied the restriction-site associated DNA sequencing (RADseq) method to genotype an F₉ recombinant inbred line population derived from a wild type × domesticated cultivar (W × D) cross. A high density linkage map was developed based on the W × D population. By integrating sequence-defined DNA markers reported in previous mapping studies, we established an ultra-high density consensus genetic map, which contains 34,574 markers consisting of 3508 loci covering 2399 cM on 20 linkage groups. The largest gap in the entire consensus map was 4.73 cM. The high density W × D map and the consensus map were used to develop an improved physical map, which covered 560.5 Mb of genome sequence data. The ultra-high density consensus linkage map, the improved physical map and the markers linked to genes of breeding interest reported in this study provide a common tool for genome sequence assembly, structural genomics, comparative genomics, functional genomics, QTL mapping, and molecular plant breeding in lupin.

Increased genomic prediction accuracy in wheat breeding using a large Australian panel

KEY MESSAGE

Genomic prediction accuracy within a large panel was found to be substantially higher than that previously observed in smaller populations, and also higher than QTL-based prediction. In recent years, genomic selection for wheat breeding has been widely studied, but this has typically been restricted to population sizes under 1000 individuals. To assess its efficacy in germplasm representative of commercial breeding programmes, we used a panel of 10,375 Australian wheat breeding lines to investigate the accuracy of genomic prediction for grain yield, physical grain quality and other physiological traits. To achieve this, the complete panel was phenotyped in a dedicated field trial and genotyped using a custom AxiomTM Affymetrix SNP array. A high-quality consensus map was also constructed, allowing the linkage disequilibrium present in the germplasm to be investigated. Using the complete SNP array, genomic prediction accuracies were found to be substantially higher than those previously observed in smaller populations and also more accurate compared to prediction approaches using a finite number of selected quantitative trait loci. Multi-trait genetic correlations were also assessed at an additive and residual genetic level, identifying a negative genetic correlation between grain yield and protein as well as a positive genetic correlation between grain size and test weight.

Integrated high-density consensus genetic map of Pyrus and anchoring of the 'Bartlett' v1.0 (Pyrus communis) genome

Genetic maps are essential tools for pear genetics and genomics research. In this study, we first constructed an integrated simple sequence repeat (SSR) and single nucleotide polymorphism (SNP)-based consensus genetic map for pear based on common SSR markers between nine published maps. A total of 5,085 markers, including 1,232 SSRs and 3,853 SNPs, were localized on a consensus map spanning 3,266.0 cM in total, with an average marker interval of 0.64 cM, which represents the highest density consensus map of pear to date. Using three sets of high-density SNP-based genetic maps with European pear genetic backgrounds, we anchored a total of 291.5 Mb of the ‘Bartlett’ v1.0 (Pyrus communis L.) genome scaffolds into 17 pseudo-chromosomes. This accounted for 50.5% of the genome assembly, which was a great improvement on the 29.7% achieved originally. Intra-genome and inter-genome synteny analyses of the new ‘Bartlett’ v1.1 genome assembly with the Asian pear ‘Dangshansuli’ (Pyrus bretschneideri Rehd.) and apple (Malus × domestica Borkh.) genomes uncovered four new segmental duplication regions. The integrated high-density SSR and SNP-based consensus genetic map provided new insights into the genetic structure patterns of pear and assisted in the genome assembly of ‘Bartlett’ through further exploration of different pear genetic maps.

Linkage Map of Lissotriton Newts Provides Insight into the Genetic Basis of Reproductive Isolation

Linkage maps are widely used to investigate structure, function, and evolution of genomes. In speciation research, maps facilitate the study of the genetic architecture of reproductive isolation by allowing identification of genomic regions underlying reduced fitness of hybrids. Here we present a linkage map for European newts of the Lissotriton vulgaris species complex, constructed using two families of F2 L. montandoni × L. vulgaris hybrids. The map consists of 1146 protein-coding genes on 12 linkage groups, equal to the haploid chromosome number, with a total length of 1484 cM (1.29 cM per marker). It is notably shorter than two other maps available for salamanders, but the differences in map length are consistent with cytogenetic estimates of the number of chiasmata per chromosomal arm. Thus, large salamander genomes do not necessarily translate into long linkage maps, as previously suggested. Consequently, salamanders are an excellent model to study evolutionary consequences of recombination rate variation in taxa with large genomes and a similar number of chromosomes. A complex pattern of transmission ratio distortion (TRD) was detected: TRD occurred mostly in one family, in one breeding season, and was clustered in two genomic segments. This is consistent with environment-dependent mortality of individuals carrying L. montandoni alleles in these two segments and suggests a role of TRD blocks in reproductive isolation. The reported linkage map will empower studies on the genomic architecture of divergence and interactions between the genomes of hybridizing newts.

Scallop genome provides insights into evolution of bilaterian karyotype and development

Reconstructing the genomes of bilaterian ancestors is central to our understanding of animal evolution, where knowledge from ancient and/or slow-evolving bilaterian lineages is critical. Here we report a high-quality, chromosome-anchored reference genome for the scallop Patinopecten yessoensis, a bivalve mollusc that has a slow-evolving genome with many ancestral features. Chromosome-based macrosynteny analysis reveals a striking correspondence between the 19 scallop chromosomes and the 17 presumed ancestral bilaterian linkage groups at a level of conservation previously unseen, suggesting that the scallop may have a karyotype close to that of the bilaterian ancestor. Scallop Hox gene expression follows a new mode of subcluster temporal co-linearity that is possibly ancestral and may provide great potential in supporting diverse bilaterian body plans. Transcriptome analysis of scallop mantle eyes finds unexpected diversity in phototransduction cascades and a potentially ancient Pax2/5/8-dependent pathway for noncephalic eyes. The outstanding preservation of ancestral karyotype and developmental control makes the scallop genome a valuable resource for understanding early bilaterian evolution and biology.

SSR-enriched genetic linkage maps of bermudagrass (Cynodon dactylon × transvaalensis), and their comparison with allied plant genomes

We report SSR-enriched genetic maps of bermudagrass that: (1) reveal partial residual polysomic inheritance in the tetraploid species, and (2) provide insights into the evolution of chloridoid genomes. This study describes genetic linkage maps of two bermudagrass species, Cynodon dactylon (T89) and Cynodon transvaalensis (T574), that integrate heterologous microsatellite markers from sugarcane into frameworks built with single-dose restriction fragments (SDRFs). A maximum likelihood approach was used to construct two separate parental maps from a population of 110 F₁ progeny of a cross between the two parents. The T89 map is based on 291 loci on 34 cosegregating groups (CGs), with an average marker spacing of 12.5 cM. The T574 map is based on 125 loci on 14 CGs, with an average marker spacing of 10.7 cM. Six T89 and one T574 CG(s) deviated from disomic inheritance. Furthermore, marker segregation data and linkage phase analysis revealed partial residual polysomic inheritance in T89, suggesting that common bermudagrass is undergoing diploidization following whole genome duplication (WGD). Twenty-six T89 CGs were coalesced into 9 homo(eo)logous linkage groups (LGs), while 12 T574 CGs were assembled into 9 LGs, both putatively representing the basic chromosome complement (x = 9) of the species. Eight T89 and two T574 CGs remain unassigned. The marker composition of bermudagrass ancestral chromosomes was inferred by aligning T89 and T574 homologs, and used in comparisons to sorghum and rice genome sequences based on 108 and 91 significant blast hits, respectively. Two nested chromosome fusions (NCFs) shared by two other chloridoids (i.e., zoysiagrass and finger millet) and at least three independent translocation events were evident during chromosome number reduction from 14 in the polyploid common ancestor of Poaceae to 9 in Cynodon.

Linkage Mapping and Comparative Genomics of Red Drum (Sciaenops ocellatus) Using Next-Generation Sequencing

Developments in next-generation sequencing allow genotyping of thousands of genetic markers across hundreds of individuals in a cost-effective manner. Because of this, it is now possible to rapidly produce dense genetic linkage maps for nonmodel species. Here, we report a dense genetic linkage map for red drum, a marine fish species of considerable economic importance in the southeastern United States and elsewhere. We used a prior microsatellite-based linkage map as a framework and incorporated 1794 haplotyped contigs derived from high-throughput, reduced representation DNA sequencing to produce a linkage map containing 1794 haplotyped restriction-site associated DNA (RAD) contigs, 437 anonymous microsatellites, and 44 expressed sequence-tag-linked microsatellites (EST-SSRs). A total of 274 candidate genes, identified from transcripts from a preliminary hydrocarbon exposure study, were localized to specific chromosomes, using a shared synteny approach. The linkage map will be a useful resource for red drum commercial and restoration aquaculture, and for better understanding and managing populations of red drum in the wild.

Genome resources for climate-resilient cowpea, an essential crop for food security

Cowpea (Vigna unguiculata L. Walp.) is a legume crop that is resilient to hot and drought-prone climates, and a primary source of protein in sub-Saharan Africa and other parts of the developing world. However, genome resources for cowpea have lagged behind most other major crops. Here we describe foundational genome resources and their application to the analysis of germplasm currently in use in West African breeding programs. Resources developed from the African cultivar IT97K-499-35 include a whole-genome shotgun (WGS) assembly, a bacterial artificial chromosome (BAC) physical map, and assembled sequences from 4355 BACs. These resources and WGS sequences of an additional 36 diverse cowpea accessions supported the development of a genotyping assay for 51 128 SNPs, which was then applied to five bi-parental RIL populations to produce a consensus genetic map containing 37 372 SNPs. This genetic map enabled the anchoring of 100 Mb of WGS and 420 Mb of BAC sequences, an exploration of genetic diversity along each linkage group, and clarification of macrosynteny between cowpea and common bean. The SNP assay enabled a diversity analysis of materials from West African breeding programs. Two major subpopulations exist within those materials, one of which has significant parentage from South and East Africa and more diversity. There are genomic regions of high differentiation between subpopulations, one of which coincides with a cluster of nodulin genes. The new resources and knowledge help to define goals and accelerate the breeding of improved varieties to address food security issues related to limited-input small-holder farming and climate stress.

Wheat Landrace Genome Diversity

Understanding the genomic complexity of bread wheat is important for unraveling domestication processes, environmental adaptation, and for future of...

Understanding the genomic complexity of bread wheat (Triticum aestivum L.) is a cornerstone in the quest to unravel the processes of domestication and the following adaptation of domesticated wheat to a wide variety of environments across the globe. Additionally, it is of importance for future improvement of the crop, particularly in the light of climate change. Focusing on the adaptation after domestication, a nested association mapping (NAM) panel of 60 segregating biparental populations was developed, mainly involving landrace accessions from the core set of the Watkins hexaploid wheat collection optimized for genetic diversity. A modern spring elite variety, “Paragon,” was used as common reference parent. Genetic maps were constructed following identical rules to make them comparable. In total, 1611 linkage groups were identified, based on recombination from an estimated 126,300 crossover events over the whole NAM panel. A consensus map, named landrace consensus map (LRC), was constructed and contained 2498 genetic loci. These newly developed genetics tools were used to investigate the rules underlying genome fluidity or rigidity, e.g., by comparing marker distances and marker orders. In general, marker order was highly correlated, which provides support for strong synteny between bread wheat accessions. However, many exceptional cases of incongruent linkage groups and increased marker distances were also found. Segregation distortion was detected for many markers, sometimes as hot spots present in different populations. Furthermore, evidence for translocations in at least 36 of the maps was found. These translocations fell, in general, into many different translocation classes, but a few translocation classes were found in several accessions, the most frequent one being the well-known T5B:7B translocation. Loci involved in recombination rate, which is an interesting trait for plant breeding, were identified by QTL analyses using the crossover counts as a trait. In total, 114 significant QTL were detected, nearly half of them with increasing effect from the nonreference parents.

QTL mapping identifies a major locus for resistance in wheat to Sunn pest (Eurygaster integriceps) feeding at the vegetative growth stage

This research provides the first report of a major locus controlling wheat resistance to Sunn pest. It developed and validated SNP markers that will be useful for marker-assisted selection. Sunn pest (Eurygaster integriceps Puton) is the most destructive insect pest of bread wheat and durum wheat in West and Central Asia and East Europe. Breeding for resistance at the vegetative stage of growth is vital in reducing the damage caused by overwintered adult populations that feed on shoot and leaves of seedlings, and in reducing the next generation of pest populations (nymphs and adults), which can cause damage to grain quality by feeding on spikes. In the present study, two doubled haploid (DH) populations involving resistant landraces from Afghanistan were genotyped with the 90k SNP iSelect assay and candidate gene-based KASP markers. The DH lines and parents were phenotyped for resistance to Sunn pest feeding, using artificial infestation cages at Terbol station, in Lebanon, over three years. Quantitative trait locus (QTL) analysis identified a single major locus on chromosome 4BS in the two populations, with the resistance allele derived from the landrace accessions, IG139431 and IG139883. The QTL explained a maximum of 42 % of the phenotypic variation in the Cham6 × IG139431 and 56 % in the Cham6 × IG139883 populations. SNP markers closest to the QTL showed high similarity to rice genes that putatively encode proteins for defense response to herbivory and wounding. The markers were validated in a large, unrelated population of parental wheat genotypes. All wheat lines carrying the 'C-G' haplotype at the identified SNPs were resistant, suggesting that selection based on a haplotype of favourable alleles would be effective in predicting resistance status of unknown genotypes.

Construction of high-resolution recombination maps in Asian seabass

Background

A high-density genetic map is essential for de novo genome assembly, fine mapping QTL for important complex traits, comparative genomic studies and understanding the mechanisms of genome evolution. Although a number of genomic resources are available in Asian seabass (Lates calcarifer), a high-density linkage map is still lacking. To facilitate QTL mapping for marker-assisted selection and genome assembly, and to understand the genome-wide recombination rates, we constructed high density linkage maps using three families and genotyping by sequencing.

Results

A high-density consensus linkage map consisting of 8, 274 markers was constructed based on sex-averaged genetic maps. The genetic maps were then aligned and integrated with the current genome assembly of Asian seabass. More than 90% of the genome contig sequences were anchored onto the consensus genetic map. Evidence of assembly errors in the current genome assembly was identified. A fragment of up to 2.5 Mb belonging to LG14 was assembled into Chr15. The length of family-specific sex-averaged maps ranged from 1348.96 to 1624.65 cM. Female maps were slightly longer than male maps using common markers. Female-to-male ratios were highly variable both across chromosomes within each family and throughout three families for each chromosome. However, the distribution patterns of recombination along chromosomes were similar between sexes across the whole genome. The overall recombination rates were significantly correlated with genome-wide GC content and the correlations were revealed to be stronger in females than in males.

Conclusions

These high-density genetic maps provide not only essential tools for facilitating de novo genome assembly and comparative genomic studies in teleosts, but also critical resources for fine mapping QTL and genome-wide association mapping for economically important traits in Asian seabass.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3462-z) contains supplementary material, which is available to authorized users.

SNP-Based Linkage Mapping for Validation of QTLs for Resistance to Ascochyta Blight in Lentil

Lentil (Lens culinaris Medik.) is a self-pollinating, diploid, annual, cool-season, food legume crop that is cultivated throughout the world. Ascochyta blight (AB), caused by Ascochyta lentis Vassilievsky, is an economically important and widespread disease of lentil. Development of cultivars with high levels of durable resistance provides an environmentally acceptable and economically feasible method for AB control. A detailed understanding of the genetic basis of AB resistance is hence highly desirable, in order to obtain insight into the number and influence of resistance genes. Genetic linkage maps based on single nucleotide polymorphisms (SNP) and simple sequence repeat (SSR) markers have been developed from three recombinant inbred line (RIL) populations. The IH × NF map contained 460 loci across 1461.6 cM, while the IH × DIG map contained 329 loci across 1302.5 cM and the third map, NF × DIG contained 330 loci across 1914.1 cM. Data from these maps were combined with a map from a previously published study through use of bridging markers to generate a consensus linkage map containing 689 loci distributed across seven linkage groups (LGs), with a cumulative length of 2429.61 cM at an average density of one marker per 3.5 cM. Trait dissection of AB resistance was performed for the RIL populations, identifying totals of two and three quantitative trait loci (QTLs) explaining 52 and 69% of phenotypic variation for resistance to infection in the IH × DIG and IH × NF populations, respectively. Presence of common markers in the vicinity of the AB_IH1- and AB_IH2.1/AB_IH2.2-containing regions on both maps supports the inference that a common genomic region is responsible for conferring resistance and is associated with the resistant parent, Indianhead. The third QTL was derived from Northfield. Evaluation of markers associated with AB resistance across a diverse lentil germplasm panel revealed that the identity of alleles associated with AB_IH1 predicted the phenotypic responses with high levels of accuracy (~86%), and therefore have the potential to be widely adopted in lentil breeding programs. The availability of RIL-based maps, a consensus map, and validated markers linked to AB resistance provide important resources for lentil improvement.

Genetic mapping of Pinus flexilis major gene (Cr4) for resistance to white pine blister rust using transcriptome-based SNP genotyping

Background

Linkage of DNA markers with phenotypic traits provides essential information to dissect clustered genes with potential phenotypic contributions in a target genome region. Pinus flexilis E. James (limber pine) is a keystone five-needle pine species in mountain-top ecosystems of North America. White pine blister rust (WPBR), caused by a non-native fungal pathogen Cronartium ribicola (J.C. Fisch.), has resulted in mortality in this conifer species and is still spreading through the distribution. The objective of this research was to develop P. flexilis transcriptome-wide single nucleotide polymorphism (SNP) markers using RNA-seq analysis for genetic mapping of the major gene (Cr4) that confers complete resistance to C. ribicola.

Results

Needle tissues of one resistant and two susceptible seedling families were subjected to RNA-seq analysis. In silico SNP markers were uncovered by mapping the RNA-seq reads back to the de novo assembled transcriptomes. A total of 110,573 in silico SNPs and 2,870 indels were identified with an average of 3.7 SNPs per Kb. These SNPs were distributed in 17,041 unigenes. Of these polymorphic P. flexilis unigenes, 6,584 were highly conserved as compared to the genome sequence of P. taeda L (loblolly pine). High-throughput genotyping arrays were designed and were used to search for Cr4-linked genic SNPs in megagametophyte populations of four maternal trees by haploid-segregation analysis. A total of 32 SNP markers in 25 genes were localized on the Cr4 linkage group (LG). Syntenic relationships of this Cr4-LG map with the model conifer species P. taeda anchored Cr4 on Pinus consensus LG8, indicating that R genes against C. ribicola have evolved independently in different five-needle pines. Functional genes close to Cr4 were annotated and their potential roles in Cr4-mediated resistance were further discussed.

Conclusions

We demonstrated a very effective, low-cost strategy for developing a SNP genetic map of a phenotypic trait of interest. SNP discovery through transcriptome comparison was integrated with high-throughput genotyping of a small set of in silico SNPs. This strategy may be applied to mapping any trait in non-model plant species that have complex genomes. Whole transcriptome sequencing provides a powerful tool for SNP discovery in conifers and other species with complex genomes, for which sequencing and annotation of complex genomes is still challenging. The genic SNP map for the consensus Cr4-LG may help future molecular breeding efforts by enabling both Cr4 positional characterization and selection of this gene against WPBR.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3079-2) contains supplementary material, which is available to authorized users.

Comparison of Spinach Sex Chromosomes with Sugar Beet Autosomes Reveals Extensive Synteny and Low Recombination at the Male-Determining Locus

Spinach (Spinacia oleracea, 2n = 12) and sugar beet (Beta vulgaris, 2n = 18) are important crop members of the family Chenopodiaceae ss Sugar beet has a basic chromosome number of 9 and a cosexual breeding system, as do most members of the Chenopodiaceae ss. family. By contrast, spinach has a basic chromosome number of 6 and, although certain cultivars and genotypes produce monoecious plants, is considered to be a dioecious species. The loci determining male and monoecious sexual expression were mapped to different loci on the spinach sex chromosomes. In this study, a linkage map with 46 mapped protein-coding sequences was constructed for the spinach sex chromosomes. Comparison of the linkage map with a reference genome sequence of sugar beet revealed that the spinach sex chromosomes exhibited extensive synteny with sugar beet chromosomes 4 and 9. Tightly linked protein-coding genes linked to the male-determining locus in spinach corresponded to genes located in or around the putative pericentromeric and centromeric regions of sugar beet chromosomes 4 and 9, supporting the observation that recombination rates were low in the vicinity of the male-determining locus. The locus for monoecism was confined to a chromosomal segment corresponding to a region of approximately 1.7Mb on sugar beet chromosome 9, which may facilitate future positional cloning of the locus.

Quantitative Trait Loci Associated with Phenological Development, Low-Temperature Tolerance, Grain Quality, and Agronomic Characters in Wheat (Triticum aestivum L.)

Plants must respond to environmental cues and schedule their development in order to react to periods of abiotic stress and commit fully to growth and reproduction under favorable conditions. This study was initiated to identify SNP markers for characters expressed from the seedling stage to plant maturity in spring and winter wheat (Triticum aestivum L.) genotypes adapted to western Canada. Three doubled haploid populations with the winter cultivar ‘Norstar’ as a common parent were developed and genotyped with a 90K Illumina iSelect SNP assay and a 2,998.9 cM consensus map with 17,541 markers constructed. High heritability’s reflected large differences among the parents and relatively low genotype by environment interactions for all characters considered. Significant QTL were detected for the 15 traits examined. However, different QTL for days to heading in controlled environments and the field provided a strong reminder that growth and development are being orchestrated by environmental cues and caution should be exercised when extrapolating conclusions from different experiments. A QTL on chromosome 6A for minimum final leaf number, which determines the rate of phenological development in the seedling stage, was closely linked to QTL for low-temperature tolerance, grain quality, and agronomic characters expressed up to the time of maturity. This suggests phenological development plays a critical role in programming subsequent outcomes for many traits. Transgressive segregation was observed for the lines in each population and QTL with additive effects were identified suggesting that genes for desirable traits could be stacked using Marker Assisted Selection. QTL were identified for characters that could be transferred between the largely isolated western Canadian spring and winter wheat gene pools demonstrating the opportunities offered by Marker Assisted Selection to act as bridges in the identification and transfer of useful genes among related genetic islands while minimizing the drag created by less desirable genes.

Determining the order of resistance genes against Stagonospora nodorum blotch, Fusarium head blight and stem rust on wheat chromosome arm 3BS

BACKGROUND

Stagonospora nodorum blotch (SNB), Fusarium head blight (FHB) and stem rust (SR), caused by the fungi Parastagonospora (synonym Stagonospora) nodorum, Fusarium graminearum and Puccinia graminis, respectively, significantly reduce yield and quality of wheat. Three resistance factors, QSng.sfr-3BS, Fhb1 and Sr2, conferring resistance, respectively, to SNB, FHB and SR, each from a unique donor line, were mapped previously to the short arm of wheat chromosome 3B. Based on published reports, our hypothesis was that Sr2 is the most distal, Fhb1 the most proximal and QSng.sfr-3BS is in between Sr2 and Fhb1 on wheat chromosome arm 3BS.

RESULTS

To test this hypothesis, 1600 F2 plants from crosses between parental lines Arina, Alsen and Ocoroni86, conferring resistance genes QSng.sfr-3BS, Fhb1 and Sr2, respectively, were genotyped and phenotyped for SNB along with the parental lines. Five closely linked single-nucleotide polymorphism (SNP) markers were used to make the genetic map and determine the gene order.

CONCLUSIONS

The results indicate that QSng.sfr-3BS is located between the other two resistance genes on chromosome 3BS. Knowing the positional order of these resistance genes will aid in developing a wheat line with all three genes in coupling, which has the potential to provide broad-spectrum resistance preventing grain yield and quality losses.

Genome-Wide QTL Mapping for Wheat Processing Quality Parameters in a Gaocheng 8901/Zhoumai 16 Recombinant Inbred Line Population

Dough rheological and starch pasting properties play an important role in determining processing quality in bread wheat (Triticum aestivum L.). In the present study, a recombinant inbred line (RIL) population derived from a Gaocheng 8901/Zhoumai 16 cross grown in three environments was used to identify quantitative trait loci (QTLs) for dough rheological and starch pasting properties evaluated by Mixograph, Rapid Visco-Analyzer (RVA), and Mixolab parameters using the wheat 90 and 660 K single nucleotide polymorphism (SNP) chip assays. A high-density linkage map constructed with 46,961 polymorphic SNP markers from the wheat 90 and 660 K SNP assays spanned a total length of 4121 cM, with an average chromosome length of 196.2 cM and marker density of 0.09 cM/marker; 6596 new SNP markers were anchored to the bread wheat linkage map, with 1046 and 5550 markers from the 90 and 660 K SNP assays, respectively. Composite interval mapping identified 119 additive QTLs on 20 chromosomes except 4D; among them, 15 accounted for more than 10% of the phenotypic variation across two or three environments. Twelve QTLs for Mixograph parameters, 17 for RVA parameters and 55 for Mixolab parameters were new. Eleven QTL clusters were identified. The closely linked SNP markers can be used in marker-assisted wheat breeding in combination with the Kompetitive Allele Specific PCR (KASP) technique for improvement of processing quality in bread wheat.

A high-density, multi-parental SNP genetic map on apple validates a new mapping approach for outcrossing species

Quantitative trait loci (QTL) mapping approaches rely on the correct ordering of molecular markers along the chromosomes, which can be obtained from genetic linkage maps or a reference genome sequence. For apple (Malus domestica Borkh), the genome sequence v1 and v2 could not meet this need; therefore, a novel approach was devised to develop a dense genetic linkage map, providing the most reliable marker-loci order for the highest possible number of markers. The approach was based on four strategies: (i) the use of multiple full-sib families, (ii) the reduction of missing information through the use of HaploBlocks and alternative calling procedures for single-nucleotide polymorphism (SNP) markers, (iii) the construction of a single backcross-type data set including all families, and (iv) a two-step map generation procedure based on the sequential inclusion of markers. The map comprises 15 417 SNP markers, clustered in 3 K HaploBlock markers spanning 1 267 cM, with an average distance between adjacent markers of 0.37 cM and a maximum distance of 3.29 cM. Moreover, chromosome 5 was oriented according to its homoeologous chromosome 10. This map was useful to improve the apple genome sequence, design the Axiom Apple 480 K SNP array and perform multifamily-based QTL studies. Its collinearity with the genome sequences v1 and v3 are reported. To our knowledge, this is the shortest published SNP map in apple, while including the largest number of markers, families and individuals. This result validates our methodology, proving its value for the construction of integrated linkage maps for any outbreeding species.

A SNP-based consensus genetic map for synteny-based trait targeting in faba bean (Vicia faba L.)

Faba bean (Vicia faba L.) is a globally important nitrogen-fixing legume, which is widely grown in a diverse range of environments. In this work, we mine and validate a set of 845 SNPs from the aligned transcriptomes of two contrasting inbred lines. Each V. faba SNP is assigned by BLAST analysis to a single Medicago orthologue. This set of syntenically anchored polymorphisms were then validated as individual KASP assays, classified according to their informativeness and performance on a panel of 37 inbred lines, and the best performing 757 markers used to genotype six mapping populations. The six resulting linkage maps were merged into a single consensus map on which 687 SNPs were placed on six linkage groups, each presumed to correspond to one of the six V. faba chromosomes. This sequence-based consensus map was used to explore synteny with the most closely related crop species, lentil and the most closely related fully sequenced genome, Medicago. Large tracts of uninterrupted colinearity were found between faba bean and Medicago, making it relatively straightforward to predict gene content and order in mapped genetic interval. As a demonstration of this, we mapped a flower colour gene to a 2-cM interval of Vf chromosome 2 which was highly colinear with Mt3. The obvious candidate gene from 78 gene models in the collinear Medicago chromosome segment was the previously characterized MtWD40-1 gene controlling anthocyanin production in Medicago and resequencing of the Vf orthologue showed a putative causative deletion of the entire 5' end of the gene.

High-density interspecific genetic linkage mapping provides insights into genomic incompatibility between channel catfish and blue catfish

Catfish is the leading aquaculture species in the United States. The interspecific hybrid catfish produced by mating female channel catfish with male blue catfish outperform both of their parent species in a number of traits. However, mass production of the hybrids has been difficult because of reproductive isolation. Investigations of genome structure and organization of the hybrids provide insights into the genetic basis for maintenance of species divergence in the face of gene flow, thereby helping develop strategies for introgression and efficient production of the hybrids for aquaculture. In this study, we constructed a high-density genetic linkage map using the hybrid catfish system with the catfish 250K SNP array. A total of 26,238 SNPs were mapped to 29 linkage groups, with 12,776 unique marker positions. The linkage map spans approximately 3240 cM with an average intermarker distance of 0.25 cM. A fraction of markers (986 of 12,776) exhibited significant deviation from the expected Mendelian ratio of segregation, and they were clustered in major genomic blocks across 15 LGs, most notably LG9 and LG15. The distorted markers exhibited significant bias for maternal alleles among the backcross progenies, suggesting strong selection against the blue catfish alleles. The clustering of distorted markers within genomic blocks should lend insights into speciation as marked by incompatibilities between the two species. Such findings should also have profound implications for understanding the genomic evolution of closely related species as well as the introgression of hybrid production programs in aquaculture.

Second-Generation Linkage Maps for the Pacific Oyster Crassostrea gigas Reveal Errors in Assembly of Genome Scaffolds

The Pacific oyster Crassostrea gigas, a widely cultivated marine bivalve mollusc, is becoming a genetically and genomically enabled model for highly fecund marine metazoans with complex life-histories. A genome sequence is available for the Pacific oyster, as are first-generation, low-density, linkage and gene-centromere maps mostly constructed from microsatellite DNA makers. Here, higher density, second-generation, linkage maps are constructed from more than 1100 coding (exonic) single-nucleotide polymorphisms (SNPs), as well as 66 previously mapped microsatellite DNA markers, all typed in five families of Pacific oysters (nearly 172,000 genotypes). The map comprises 10 linkage groups, as expected, has an average total length of 588 cM, an average marker-spacing of 1.0 cM, and covers 86% of a genome estimated to be 616 cM. All but seven of the mapped SNPs map to 618 genome scaffolds; 260 scaffolds contain two or more mapped SNPs, but for 100 of these scaffolds (38.5%), the contained SNPs map to different linkage groups, suggesting widespread errors in scaffold assemblies. The 100 misassembled scaffolds are significantly longer than those that map to a single linkage group. On the genetic maps, marker orders and intermarker distances vary across families and mapping methods, owing to an abundance of markers segregating from only one parent, to widespread distortions of segregation ratios caused by early mortality, as previously observed for oysters, and to genotyping errors. Maps made from framework markers provide stronger support for marker orders and reasonable map lengths and are used to produce a consensus high-density linkage map containing 656 markers.

A Consensus Genetic Map for Pinus taeda and Pinus elliottii and Extent of Linkage Disequilibrium in Two Genotype-Phenotype Discovery Populations of Pinus taeda

A consensus genetic map for Pinus taeda (loblolly pine) and Pinus elliottii (slash pine) was constructed by merging three previously published P. taeda maps with a map from a pseudo-backcross between P. elliottii and P. taeda. The consensus map positioned 3856 markers via genotyping of 1251 individuals from four pedigrees. It is the densest linkage map for a conifer to date. Average marker spacing was 0.6 cM and total map length was 2305 cM. Functional predictions of mapped genes were improved by aligning expressed sequence tags used for marker discovery to full-length P. taeda transcripts. Alignments to the P. taeda genome mapped 3305 scaffold sequences onto 12 linkage groups. The consensus genetic map was used to compare the genome-wide linkage disequilibrium in a population of distantly related P. taeda individuals (ADEPT2) used for association genetic studies and a multiple-family pedigree used for genomic selection (CCLONES). The prevalence and extent of LD was greater in CCLONES as compared to ADEPT2; however, extended LD with LGs or between LGs was rare in both populations. The average squared correlations, r², between SNP alleles less than 1 cM apart were less than 0.05 in both populations and r² did not decay substantially with genetic distance. The consensus map and analysis of linkage disequilibrium establish a foundation for comparative association mapping and genomic selection in P. taeda and P. elliottii.

Genetic architecture of carbon isotope composition and growth in Eucalyptus across multiple environments

In the context of climate change, the water-use efficiency (WUE) of highly productive tree varieties, such as eucalypts, has become a major issue for breeding programmes. This study set out to dissect the genetic architecture of carbon isotope composition (δ(13) C), a proxy of WUE, across several environments. A family of Eucalyptus urophylla × E. grandis was planted in three trials and phenotyped for δ(13) C and growth traits. High-resolution genetic maps enabled us to target genomic regions underlying δ(13) C quantitative trait loci (QTLs) on the E. grandis genome. Of the 15 QTLs identified for δ(13) C, nine were stable across the environments and three displayed significant QTL-by-environment interaction, suggesting medium to high genetic determinism for this trait. Only one colocalization was found between growth and δ(13) C. Gene ontology (GO) term enrichment analysis suggested candidate genes related to foliar δ(13) C, including two involved in the regulation of stomatal movements. This study provides the first report of the genetic architecture of δ(13) C and its relation to growth in Eucalyptus. The low correlations found between the two traits at phenotypic and genetic levels suggest the possibility of improving the WUE of Eucalyptus varieties without having an impact on breeding for growth.