Development of F1 hybrid population and the high-density linkage map for European aspen (<Emphasis Type="Italic">Populus tremula</Emphasis> L.) using RADseq technology

Background

Restriction-site associated DNA sequencing (RADseq) technology was recently employed to identify a large number of single nucleotide polymorphisms (SNP) for linkage mapping of a North American and Eastern Asian Populus species. However, there is also the need for high-density genetic linkage maps for the European aspen (P. tremula) as a tool for further mapping of quantitative trait loci (QTLs) and marker-assisted selection of the Populus species native to Europe.

Results

We established a hybrid F1 population from the cross of two aspen parental genotypes diverged in their phenological and morphological traits. We performed RADseq of 122 F1 progenies and two parents yielding 15,732 high-quality SNPs that were successfully identified using the reference genome of P. trichocarpa. 2055 SNPs were employed for the construction of maternal and paternal linkage maps. The maternal linkage map was assembled with 1000 SNPs, containing 19 linkage groups and spanning 3054.9 cM of the genome, with an average distance of 3.05 cM between adjacent markers. The paternal map consisted of 1055 SNPs and the same number of linkage groups with a total length of 3090.56 cM and average interval distance of 2.93 cM. The linkage maps were employed for QTL mapping of one-year-old seedlings height variation. The most significant QTL (LOD = 5.73) was localized to LG5 (96.94 cM) of the male linkage map, explaining 18% of the phenotypic variation.

Conclusions

The set of 15,732 SNPs polymorphic in aspen and high-density genetic linkage maps constructed for the P. tremula intra-specific cross will provide a valuable source for QTL mapping and identification of candidate genes facilitating marker-assisted selection in European aspen.

     

Quantitative trait locus mapping under irrigated and drought treatments based on a novel genetic linkage map in mungbean (<Emphasis Type="Italic">Vigna radiata</Emphasis> L.)

Key message

A novel genetic linkage map was constructed using SSR markers and stable QTLs were identified for six drought tolerance related-traits using single-environment analysis under irrigation and drought treatments.

Abstract

Mungbean (Vigna radiata L.) is one of the most important leguminous food crops. However, mungbean production is seriously constrained by drought. Isolation of drought-responsive genetic elements and marker-assisted selection breeding will benefit from the detection of quantitative trait locus (QTLs) for traits related to drought tolerance. In this study, we developed a full-coverage genetic linkage map based on simple sequence repeat (SSR) markers using a recombinant inbred line (RIL) population derived from an intra-specific cross between two drought-resistant varieties. This novel map was anchored with 313 markers. The total map length was 1010.18 cM across 11 linkage groups, covering the entire genome of mungbean with a saturation of one marker every 3.23 cM. We subsequently detected 58 QTLs for plant height (PH), maximum leaf area (MLA), biomass (BM), relative water content, days to first flowering, and seed yield (Yield) and 5 for the drought tolerance index of 3 traits in irrigated and drought environments at 2 locations. Thirty-eight of these QTLs were consistently detected two or more times at similar linkage positions. Notably, qPH5A and qMLA2A were consistently identified in marker intervals from GMES5773 to MUS128 in LG05 and from Mchr11-34 to the HAAS_VR_1812 region in LG02 in four environments, contributing 6.40–20.06% and 6.97–7.94% of the observed phenotypic variation, respectively. None of these QTLs shared loci with previously identified drought-related loci from mungbean. The results of these analyses might facilitate the isolation of drought-related genes and help to clarify the mechanism of drought tolerance in mungbean.

     

Fine-mapping of QTLs for individual and total isoflavone content in soybean (<Emphasis Type="Italic">Glycine max L.</Emphasis>) using a high-density genetic map

Key message

Fifteen stable QTLs were identified using a high-density soybean genetic map across multiple environments. One major QTL, qIF5-1, contributing to total isoflavone content explained phenotypic variance 49.38, 43.27, 46.59, 45.15 and 52.50%, respectively.

Abstract

Soybeans (Glycine max L.) are a major source of dietary isoflavones. To identify novel quantitative trait loci (QTL) underlying isoflavone content, and to improve the accuracy of marker-assisted breeding in soybean, a valuable mapping population comprised of 196 F_7:8–10 recombinant inbred lines (RILs, Huachun 2 × Wayao) was utilized to evaluate individual and total isoflavone content in plants grown in four different environments in Guangdong. A high-density genetic linkage map containing 3469 recombination bin markers based on 0.2 × restriction site-associated DNA tag sequencing (RAD-seq) technology was used to finely map QTLs for both individual and total isoflavone contents. Correlation analyses showed that total isoflavone content, and that of five individual isoflavone, was significantly correlated across the four environments. Based on the high-density genetic linkage map, a total of 15 stable quantitative trait loci (QTLs) associated with isoflavone content across multiple environments were mapped onto chromosomes 02, 05, 07, 09, 10, 11, 13, 16, 17, and 19. Further, one of them, qIF5-1, localized to chromosomes 05 (38,434,171–39,045,620 bp) contributed to almost all isoflavone components across all environments, and explained 6.37–59.95% of the phenotypic variance, especially explained 49.38, 43.27, 46.59, 45.15 and 52.50% for total isoflavone. The results obtained in the present study will pave the way for a better understanding of the genetics of isoflavone accumulation and reveals the scope available for improvement of isoflavone content through marker-assisted selection.

     

Identification of candidate genes of QTLs for seed weight in <Emphasis Type="Italic">Brassica napus</Emphasis> through comparative mapping among <Emphasis Type="Italic">Arabidopsis</Emphasis> and <Emphasis Type="Italic">Brassica</Emphasis>species

Background

Map-based cloning of quantitative trait loci (QTLs) in polyploidy crop species remains a challenge due to the complexity of their genome structures. QTLs for seed weight in B. napus have been identified, but information on candidate genes for identified QTLs of this important trait is still rare.

Results

In this study, a whole genome genetic linkage map for B. napus was constructed using simple sequence repeat (SSR) markers that covered a genetic distance of 2,126.4 cM with an average distance of 5.36 cM between markers. A procedure was developed to establish colinearity of SSR loci on B. napus with its two progenitor diploid species B. rapa and B. oleracea through extensive bioinformatics analysis. With the aid of B. rapa and B. oleracea genome sequences, the 421 homologous colinear loci deduced from the SSR loci of B. napus were shown to correspond to 398 homologous loci in Arabidopsis thaliana. Through comparative mapping of Arabidopsis and the three Brassica species, 227 homologous genes for seed size/weight were mapped on the B. napus genetic map, establishing the genetic bases for the important agronomic trait in this amphidiploid species. Furthermore, 12 candidate genes underlying 8 QTLs for seed weight were identified, and a gene-specific marker for BnAP2 was developed through molecular cloning using the seed weight/size gene distribution map in B. napus.

Conclusions

Our study showed that it is feasible to identify candidate genes of QTLs using a SSR-based B. napus genetic map through comparative mapping among Arabidopsis and B. napus and its two progenitor species B. rapa and B. oleracea. Identification of candidate genes for seed weight in amphidiploid B. napus will accelerate the process of isolating the mapped QTLs for this important trait, and this approach may be useful for QTL identification of other traits of agronomic significance.

     

Comparison of single-nucleotide polymorphisms and microsatellites in detecting quantitative trait loci for alcoholism: The Collaborative Study on the Genetics of Alcoholism

Background

The feasibility of effectively analyzing high-density single nucleotide polymorphism (SNP) maps in whole genome scans of complex traits is not known. The purpose of this study was to compare variance components linkage results using different density marker maps in data from the Collaborative Study on the Genetics of Alcoholism (COGA). Marker maps having an average spacing of 10 cM (microsatellite), 0.78 cM (SNP1), and 0.31 cM (SNP2) were used to identify quantitative trait loci (QTLs) affecting maximum number of alcoholic drinks consumed in a 24-hour period (lnmaxalc).

Results

Heritability of lnmaxalc was estimated to be 15%. Multipoint variance components linkage analysis revealed similar linkage patterns among the three marker panels, with the SNP maps consistently yielding higher LOD scores. Robust LOD scores > 1.0 were observed on chromosomes 1 and 13 for all three marker maps. Additional LODs > 1.0 were observed on chromosome 4 with both SNP maps and on chromosomes 18 and 21 with the SNP2 map. Peak LOD scores for lnmaxalc were observed on chromosome 1, although none reached genome-wide statistical significance. Quantile-quantile plots revealed that the multipoint distribution of SNP results appeared to fit the asymptotic null distribution better than the twopoint results.

Conclusion

In conclusion, variance-components linkage analysis using high-density SNP maps provided higher LOD scores compared with the standard microsatellite map, similar to studies using nonparametric linkage methods. Widespread application of SNP maps will depend on further improvements in the computational methods implemented in current software packages.

     

Quantitative Trait Locus (QTL) meta-analysis and comparative genomics for candidate gene prediction in perennial ryegrass (<Emphasis Type="Italic">Lolium perenne</Emphasis> L.)

Background

In crop species, QTL analysis is commonly used for identification of factors contributing to variation of agronomically important traits. As an important pasture species, a large number of QTLs have been reported for perennial ryegrass based on analysis of biparental mapping populations. Further characterisation of those QTLs is, however, essential for utilisation in varietal improvement programs.

Results

A bibliographic survey of perennial ryegrass trait-dissection studies identified a total of 560 QTLs from previously published papers, of which 189, 270 and 101 were classified as morphology-, physiology- and resistance/tolerance-related loci, respectively. The collected dataset permitted a subsequent meta-QTL study and implementation of a cross-species candidate gene identification approach. A meta-QTL analysis based on use of the BioMercator software was performed to identify two consensus regions for pathogen resistance traits. Genes that are candidates for causal polymorphism underpinning perennial ryegrass QTLs were identified through in silico comparative mapping using rice databases, and 7 genes were assigned to the p150/112 reference map. Markers linked to the Lp DGL1, Lp Ph1 and Lp PIPK1 genes were located close to plant size, leaf extension time and heading date-related QTLs, respectively, suggesting that these genes may be functionally associated with important agronomic traits in perennial ryegrass.

Conclusions

Functional markers are valuable for QTL meta-analysis and comparative genomics. Enrichment of such genetic markers may permit further detailed characterisation of QTLs. The outcomes of QTL meta-analysis and comparative genomics studies may be useful for accelerated development of novel perennial ryegrass cultivars with desirable traits.

     

QTL mapping of flag leaf-related traits in wheat (<Emphasis Type="Italic">Triticum aestivum</Emphasis> L.)

Key message

QTL controlling flag leaf length, flag leaf width, flag leaf area and flag leaf angle were mapped in wheat.

Abstract

This study aimed to advance our understanding of the genetic mechanisms underlying morphological traits of the flag leaves of wheat (Triticum aestivum L.). A recombinant inbred line (RIL) population derived from ND3331 and the Tibetan semi-wild wheat Zang1817 was used to identify quantitative trait loci (QTLs) controlling flag leaf length (FLL), flag leaf width (FLW), flag leaf area (FLA), and flag leaf angle (FLANG). Using an available simple sequence repeat genetic linkage map, 23 putative QTLs for FLL, FLW, FLA, and FLANG were detected on chromosomes 1B, 2B, 3A, 3D, 4B, 5A, 6B, 7B, and 7D. Individual QTL explained 4.3–68.52% of the phenotypic variance in different environments. Four QTLs for FLL, two for FLW, four for FLA, and five for FLANG were detected in at least two environments. Positive alleles of 17 QTLs for flag leaf-related traits originated from ND3331 and 6 originated from Zang1817. QTLs with pleiotropic effects or multiple linked QTL were also identified on chromosomes 1B, 4B, and 5A; these are potential target regions for fine-mapping and marker-assisted selection in wheat breeding programs.

     

Nucleotide binding domain 1 pharmacophore modeling for visualization and analysis of P-glycoprotein–flavonoid molecular interactions

Background

P-glycoprotein (P-gp) is a 170-kDa membrane protein. It provides a barrier function and help to excrete toxins from the body as a transporter. Some bioflavonoids have been shown to block P-gp activity.

Objective

To evaluate the important amino acid residues within nucleotide binding domain 1 (NBD1) of P-gp that play a key role in molecular interactions with flavonoids using structure-based pharmacophore model.

Methods

In the molecular docking with NBD1 models, a putative binding site of flavonoids was proposed and compared with the site for ATP. The binding modes for ligands were achieved using LigandScout to generate the P-gp–flavonoid pharmacophore models.

Results

The binding pocket for flavonoids was investigated and found these inhibitors compete with the ATP for binding site in NBD1 including the NBD1 amino acid residues identified by the in silico techniques to be involved in the hydrogen bonding and van der Waals (hydrophobic) interactions with flavonoids.

Conclusion

These flavonoids occupy with the same binding site of ATP in NBD1 proffering that they may act as an ATP competitive inhibitor.

     

Multi-environment QTL mapping reveals genetic architecture of fruit cracking in a tomato RIL <Emphasis Type="Italic">Solanum lycopersicum</Emphasis> × <Emphasis Type="Italic">S. pimpinellifolium</Emphasis> population

Key message

QTL and codominant genetic markers for fruit cracking have been identified in a tomato genetic map derived from a RIL population, providing molecular tools for marker-assisted breeding of this trait.

Abstract

In tomato, as well as in other fleshy fruits, one of the main disorders that widely limit quality and production is fruit cracking or splitting of the epidermis that is observed on the fruit skin and flesh at any stage of fruit growth and maturation. To elucidate the genetic basis of fruit cracking, a quantitative trait loci (QTL) analysis was conducted in a recombinant inbred line (RIL) population derived from a cross between tomato (Solanum lycopersicum) and the wild-relative species S. pimpinellifolium. The RIL population was evaluated for fruit cracking during three consecutive growing seasons. Construction of a high-density linkage map based on codominant markers, covering more than 1000 cM of the whole genome, led to the identification of both main and epistatic QTL controlling fruit cracking on the basis of a single-environment as well as multiple-environment analysis. This information will enhance molecular breeding for novel cracking resistant varieties and simultaneously assist the identification of genes underlying these QTL, helping to reveal the genetic basis of fruit cracking in tomato.

     

Construction of an ultra-high density consensus genetic map,and enhancement of the physical map from genome sequencing in <Emphasis Type="Italic">Lupinus angustifolius</Emphasis>

Key message

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.

Abstract

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.

     

Genome-wide linkage analysis of age at onset of alcohol dependence: a comparison between microsatellites and single-nucleotide polymorphisms

Background

Using the dataset provided for Genetic Analysis Workshop 14 by the Collaborative Study on the Genetics of Alcoholism, we performed genome-wide linkage analysis of age at onset of alcoholism to compare the utility of microsatellites and single-nucleotide polymorphisms (SNPs) in genetic linkage study.

Methods

A multipoint nonparametric variance component linkage analysis method was applied to the survival distribution function obtained from semiparametric proportional hazards model of the age at onset phenotype of alcoholism. Three separate linkage analyses were carried out using 315 microsatellites, 2,467 and 9,467 SNPs, spanning the 22 autosomal chromosomes.

Results

Heritability of age at onset was estimated to be approximately 12% (p < 0.001). We observed weak correlation, both in trend and strength, of genome-wide linkage signals between microsatellites and SNPs. Results from SNPs revealed more and stronger linkage signals across the genome compared with those from microsatellites. The only suggestive evidence of linkage from microsatellites was on chromosome 1 (LOD of 1.43). Differences in map densities between the two sets of SNPs used in this study did not appear to confer an advantage in terms of strength of linkage signals.

Conclusion

Our study provided support for better performance of dense SNP maps compared with the sparse mirosatellite maps currently available for linkage analysis of quantitative traits. This better performance could be attributable to precise definition and high map resolutions achievable with dense SNP maps, thus resulting in increased power to detect possible loci affecting given trait or disease.

     

The potential of probabilistic graphical models in linkage map construction

Key message

Probabilistic graphical models show great potential for robust and reliable construction of linkage maps. We show how to use probabilistic graphical models to construct high-quality linkage maps in the face of data perturbations caused by genotyping errors and reciprocal translocations.

Abstract

It has been shown that linkage map construction can be hampered by the presence of genotyping errors and chromosomal rearrangements such as inversions and translocations. Here, we report a novel method for linkage map construction using probabilistic graphical models. The method is proven, both theoretically and practically, to be effective in filtering out markers that contain genotyping errors. In particular, it carries out marker filtering and ordering simultaneously, and is therefore superior to the standard post hoc filtering using nearest-neighbour stress. Furthermore, we demonstrate empirically that the proposed method offers a promising solution to linkage map construction in the case of a reciprocal translocation.

     

Construction of the First High-Density Genetic Linkage Map and Analysis of Quantitative Trait Loci for Growth-Related Traits in <Emphasis Type="Italic">Sinonovacula constricta</Emphasis>

The razor clam (Sinonovacula constricta) is an important aquaculture species, for which a high-density genetic linkage map would play an important role in marker-assisted selection (MAS). In this study, we constructed a high-density genetic map and detected quantitative trait loci (QTLs) for Sinonovacula constricta with an F₁ cross population by using the specific locus amplified fragment sequencing (SLAF-seq) method. A total of 315,553 SLAF markers out of 467.71 Mreads were developed. The final linkage map was composed of 7516 SLAFs (156.60-fold in the parents and 20.80-fold in each F₁ population on average). The total distance of the linkage map was 2383.85 cM, covering 19 linkage groups with an average inter-marker distance of 0.32 cM. The proportion of gaps less than 5.0 cM was on average 96.90%. A total of 16 suggestive QTLs for five growth-related traits (five QTLs for shell height, six QTLs for shell length, three QTLs for shell width, one QTL for total body weight, and one QTL for soft body weight) were identified. These QTLs were distributed on five linkage groups, and the regions showed overlapping on LG9 and LG13. In conclusion, the high-density genetic map and QTLs for S. constricta provide a valuable genetic resource and a basis for MAS.     

Bayesian QTL mapping using genome-wide SSR markers and segregating population derived from a cross of two commercial F<Subscript>1</Subscript> hybrids of tomato

Key message

Using newly developed euchromatin-derived genomic SSR markers and a flexible Bayesian mapping method, 13 significant agricultural QTLs were identified in a segregating population derived from a four-way cross of tomato.

Abstract

So far, many QTL mapping studies in tomato have been performed for progeny obtained from crosses between two genetically distant parents, e.g., domesticated tomatoes and wild relatives. However, QTL information of quantitative traits related to yield (e.g., flower or fruit number, and total or average weight of fruits) in such intercross populations would be of limited use for breeding commercial tomato cultivars because individuals in the populations have specific genetic backgrounds underlying extremely different phenotypes between the parents such as large fruit in domesticated tomatoes and small fruit in wild relatives, which may not be reflective of the genetic variation in tomato breeding populations. In this study, we constructed F₂ population derived from a cross between two commercial F₁ cultivars in tomato to extract QTL information practical for tomato breeding. This cross corresponded to a four-way cross, because the four parental lines of the two F₁ cultivars were considered to be the founders. We developed 2510 new expressed sequence tag (EST)-based (euchromatin-derived) genomic SSR markers and selected 262 markers from these new SSR markers and publicly available SSR markers to construct a linkage map. QTL analysis for ten agricultural traits of tomato was performed based on the phenotypes and marker genotypes of F₂ plants using a flexible Bayesian method. As results, 13 QTL regions were detected for six traits by the Bayesian method developed in this study.

     

Epistatic determinism of durum wheat resistance to the wheat spindle streak mosaic virus

Key message

The resistance of durum wheat to the Wheat spindle streak mosaic virus (WSSMV) is controlled by two main QTLs on chromosomes 7A and 7B, with a huge epistatic effect.

Abstract

Wheat spindle streak mosaic virus (WSSMV) is a major disease of durum wheat in Europe and North America. Breeding WSSMV-resistant cultivars is currently the only way to control the virus since no treatment is available. This paper reports studies of the inheritance of WSSMV resistance using two related durum wheat populations obtained by crossing two elite cultivars with a WSSMV-resistant emmer cultivar. In 2012 and 2015, 354 recombinant inbred lines (RIL) were phenotyped using visual notations, ELISA and qPCR and genotyped using locus targeted capture and sequencing. This allowed us to build a consensus genetic map of 8568 markers and identify three chromosomal regions involved in WSSMV resistance. Two major regions (located on chromosomes 7A and 7B) jointly explain, on the basis of epistatic interactions, up to 43% of the phenotypic variation. Flanking sequences of our genetic markers are provided to facilitate future marker-assisted selection of WSSMV-resistant cultivars.

     

An ultra-dense integrated linkage map for hexaploid chrysanthemum enables multi-allelic QTL analysis

Key message

We constructed the first integrated genetic linkage map in a polysomic hexaploid. This enabled us to estimate inheritance of parental haplotypes in the offspring and detect multi-allelic QTL.

Abstract

Construction and use of linkage maps are challenging in hexaploids with polysomic inheritance. Full map integration requires calculations of recombination frequency between markers with complex segregation types. In addition, detection of QTL in hexaploids requires information on all six alleles at one locus for each individual. We describe a method that we used to construct a fully integrated linkage map for chrysanthemum (Chrysanthemum × morifolium, 2n = 6x = 54). A bi-parental F1 population of 406 individuals was genotyped with an 183,000 SNP genotyping array. The resulting linkage map consisted of 30,312 segregating SNP markers of all possible marker dosage types, representing nine chromosomal linkage groups and 107 out of 108 expected homologues. Synteny with lettuce (Lactuca sativa) showed local colinearity. Overall, it was high enough to number the chrysanthemum chromosomal linkage groups according to those in lettuce. We used the integrated and phased linkage map to reconstruct inheritance of parental haplotypes in the F1 population. Estimated probabilities for the parental haplotypes were used for multi-allelic QTL analyses on four traits with different underlying genetic architectures. This resulted in the identification of major QTL that were affected by multiple alleles having a differential effect on the phenotype. The presented linkage map sets a standard for future genetic mapping analyses in chrysanthemum and closely related species. Moreover, the described methods are a major step forward for linkage mapping and QTL analysis in hexaploids.

     

Targeted metabolomic profiling of low and high grade serous epithelial ovarian cancer tissues: a pilot study

Introduction

Epithelial ovarian cancer (EOC) remains the leading cause of death from gynecologic malignancies and has an alarming global fatality rate. Besides the differences in underlying pathogenesis, distinguishing between high grade (HG) and low grade (LG) EOC is imperative for the prediction of disease progression and responsiveness to chemotherapy.

Objectives

The aim of this study was to investigate, the tissue metabolome associated with HG and LG serous epithelial ovarian cancer.

Methods

A combination of one dimensional proton nuclear magnetic resonance (1D H NMR) spectroscopy and targeted mass spectrometry (MS) was employed to profile the tissue metabolome of HG, LG serous EOCs, and controls.

Results

Using partial least squares-discriminant analysis, we observed significant separation between all groups (p?<?0.05) following cross validation. We identified which metabolites were significantly perturbed in each EOC grade as compared with controls and report the biochemical pathways which were perturbed due to the disease. Among these metabolic pathways, ascorbate and aldarate metabolism was identified, for the first time, as being significantly altered in both LG and HG serous cancers. Further, we have identified potential biomarkers of EOC and generated predictive algorithms with AUC (CI)?=?0.940 and 0.929 for HG and LG, respectively.

Conclusion

These previously unreported biochemical changes provide a framework for future metabolomic studies for the development of EOC biomarkers. Finally, pharmacologic targeting of the key metabolic pathways identified herein could lead to novel and effective treatments of EOC.