Genetic dissection of wheat panicle traits using linkage analysis and a genome-wide association study

Key message

Coincident regions on chromosome 4B for GW, on 5A for SD and TSS, and on 3A for SL and GNS were detected through an integration of a linkage analysis and a genome-wide association study (GWAS). In addition, six stable QTL clusters on chromosomes 2D, 3A, 4B, 5A and 6A were identified with high PVE% on a composite map.

Abstract

The panicle traits of wheat, such as grain number per spike and 1000-grain weight, are closely correlated with grain yield. Superior and effective alleles at loci related to panicles developments play a crucial role in the progress of molecular improvement in wheat yield breeding. Here, we revealed several notable allelic variations of seven panicle-related traits through an integration of genome-wide association mapping and a linkage analysis. The linkage analysis was performed using a recombinant inbred line (RIL) population (173 lines of F_8:9) with a high-density genetic map constructed with 90K SNP arrays, Diversity Arrays Technology (DArT) and simple sequence repeat (SSR) markers in five environments. Thirty-five additive quantitative trait loci (QTL) were discovered, including eleven stable QTLs on chromosomes 1A, 2D, 4B, 5B, 6B, and 6D. The marker interval between EX_C101685 and RAC875_C27536 on chromosome 4B exhibited pleiotropic effects for GW, SL, GNS, FSN, SSN, and TSS, with the phenotypic variation explained (PVE) ranging from 5.40 to 37.70%. In addition, an association analysis was conducted using a diverse panel of 205 elite wheat lines with a composite map (24,355 SNPs) based on the Illumina Infinium assay in four environments. A total of 73 significant marker-trait associations (MTAs) were detected for panicle traits, which were distributed across all wheat chromosomes except for 4D, 5D, and 6D. Consensus regions between RAC875_C27536_611 and Tdurum_contig4974_355 on chromosome 4B for GW in multiple environments, between QTSS5A.7-43 and BS00021805_51 on 5A for SD and TSS, and between QSD3A.2-164 and RAC875_c17479_359 on 3A for SL and GNS in multiple environments were detected through linkage analysis and a genome-wide association study (GWAS). In addition, six stable QTL clusters on chromosomes 2D, 3A, 4B, 5A, and 6A were identified with high PVE% on a composite map. This study provides potentially valuable information on the dissection of yield-component traits and valuable genetic alleles for molecular-design breeding or functional gene exploration.

     

QTL mapping of potato chip color and tuber traits within an autotetraploid family

Cultivated potato (Solanum tuberosum L.) is a highly heterozygous autotetraploid crop species, and this creates challenges for traditional line development and molecular breeding. Recent availability of a single-nucleotide polymorphism (SNP) array with 8303 features and software packages for linkage and association mapping in autotetraploid species present new opportunities for the identification of genomic regions that contribute to high-value traits in cultivated potato. A biparental tetraploid potato population was evaluated across three field seasons and storage trials in order to identify quantitative trait loci (QTL) for multiple tuber traits including fried chip color after 5.5–7.2 °C storage. Tetra-allelic dosage information was used to construct a genetic linkage map that covered 1041 cM and contained 2095 SNP markers with a median marker interval of 0.4 cM. A total of 41 QTL were identified for flower color, tuber yield, tuber number per plant, tuber weight, tuber size, and chip color after various storage regimes. Moderate effect QTL for chip color at 3 months were identified that co-localized with candidate genes vacuolar invertase (VInv), invertase inhibitor (INH2), and apoplastic invertase (Inv _ap -b). A separate QTL for chip color after 6 months of storage was identified in the short arm of chromosome 2, and this locus may contribute to variation in senescent sweetening resistance. QTL for tuber weight, length, and width co-localized with a known QTL for plant maturity on chromosome 5. Genome-wide association mapping using a polyploid model detected the tuber size QTL and identified a number of candidate SNPs, but was unable to detect markers significantly associated with chip color.     

Molecular dissection of QTL governing grain size traits employing association and linkage mapping in Basmati rice

Grain size is one of the key traits that determines the quality of Basmati rice from the consumers’ as well as the traders’ point of view. Though many genes governing grain size have been identified in indica and japonica, little work has been done in Basmati rice. The present study aims at dissection of a QTL region governing grain size traits in Basmati employing association and linkage mapping approaches. Association mapping revealed that three markers, i.e., RM 6024 (grain breadth), RM1237 and RM18582 (grain length-breadth ratio), which cover 889 kb in the targeted QTL region have been significantly associated with grain size traits. Using linkage mapping, the targeted QTL region has been further delimited to a physical distance of 268 kb that comprises 24 annotated genes. The gene expression analysis of parents, revealed 19 genes differentially expressing within the QTL. Of them, 15 genes showed high expression in Basmati370, while four were expressed in Jaya, and whereas five genes did not show any differential expression between parents. Among differentially expressed genes, a highly expressed gene in Basmati370, Os05g0374200 (Cytokinin dehydrogenase 1 precursor) seems to be involved in accumulation of cytokinins, thus affecting the grain size. Therefore, our findings demonstrated that by complimenting association and linkage mapping, it is likely to dissect a QTL governing grain size traits in Basmati rice and also the QTL could be a potential target for marker-assisted breeding and map-based cloning studies.     

Genetic dissection of sorghum grain quality traits using diverse and segregating populations

Key message

Coordinated association and linkage mapping identified 25 grain quality QTLs in multiple environments, and fine mapping of the Wx locus supports the use of high-density genetic markers in linkage mapping.

Abstract

There is a wide range of end-use products made from cereal grains, and these products often demand different grain characteristics. Fortunately, cereal crop species including sorghum [Sorghum bicolor (L.) Moench] contain high phenotypic variation for traits influencing grain quality. Identifying genetic variants underlying this phenotypic variation allows plant breeders to develop genotypes with grain attributes optimized for their intended usage. Multiple sorghum mapping populations were rigorously phenotyped across two environments (SC Coastal Plain and Central TX) in 2 years for five major grain quality traits: amylose, starch, crude protein, crude fat, and gross energy. Coordinated association and linkage mapping revealed several robust QTLs that make prime targets to improve grain quality for food, feed, and fuel products. Although the amylose QTL interval spanned many megabases, the marker with greatest significance was located just 12 kb from waxy (Wx), the primary gene regulating amylose production in cereal grains. This suggests higher resolution mapping in recombinant inbred line (RIL) populations can be obtained when genotyped at a high marker density. The major QTL for crude fat content, identified in both a RIL population and grain sorghum diversity panel, encompassed the DGAT1 locus, a critical gene involved in maize lipid biosynthesis. Another QTL on chromosome 1 was consistently mapped in both RIL populations for multiple grain quality traits including starch, crude protein, and gross energy. Collectively, these genetic regions offer excellent opportunities to manipulate grain composition and set up future studies for gene validation.

     

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.

     

Characterization of aphid resistance loci in black raspberry (<Emphasis Type="Italic">Rubus occidentalis</Emphasis> L.)

Viruses vectored by the aphid Amphorophora agathonica cause decline in black raspberry plant health resulting in a shortened life and poor fruit quality of the infected plantings. New aphid resistant cultivars could increase the longevity of plantings providing growers and processors with consistent fruit production. Recent exploration of the native range of black raspberry identified three sources of aphid resistance: Ag₄ from Ontario (ON), Canada, Ag₅ from Maine (ME), and a third source from Michigan (MI) with no formal designation. The objectives of this study were to assess segregation of these three sources of aphid resistance in populations with single and combined sources and develop markers that can identify each source of resistance. A genetic linkage map constructed for ORUS 4305 placed the ON aphid resistance locus on Rubus linkage group (RLG) 6. Segregation ratios in populations with single and combined sources, and linkage mapping in two populations (ORUS 4304 and ORUS 4812) segregating for the Ag₅ and MI sources, respectively, indicated that these three sources of resistance are each conferred by single dominant genes/alleles that are linked on RLG6. Confirmation of marker association in 16 validation populations identified four markers that could be used to predict resistance; however, none could distinguish between the ON and MI sources. These four markers may be useful for screening populations to enrich the field-planted progeny for aphid resistance. Fine mapping of the resistance loci is needed to develop functional markers at each of the resistance loci to enable pyramiding and durable aphid resistance.     

Revealing of allelic polymorphism in the populations of parthenogenetic lizards <Emphasis Type="Italic">Darevskia dahli</Emphasis> (Lacertidae) using locus-specific PCR

Locus-specific PCR was used to study the genetic polymorphism in three populations of parthenogenetic lizard species Darevskia dahli. The analysis was carried at the two (GATA) _n -containing loci (Du215 and Du281) using the sample of 26 individuals. A total of eight Du215 and three Du281 allelic variants were detected. It was demonstrated that all the lizards examined were heterozygous at these loci. In 12 animals, unusual Du215 allelic variant was revealed, the origin of which was thought to be associated with different types of genomic rearrangements, or segmental duplication. The populations studied were substantially different relative to the levels of allelic polymorphism, which could be explained by different habitation conditions, leading to accumulation of mutations in noncoding genome regions.     

Analysis of the genetic architecture of maize ear and grain morphological traits by combined linkage and association mapping

Key message

Using combined linkage and association mapping, 26 stable QTL and six stable SNPs were detected across multiple environments for eight ear and grain morphological traits in maize. One QTL, PKS2, might play an important role in maize yield improvement.

Abstract

In the present study, one bi-parental population and an association panel were used to identify quantitative trait loci (QTL) for eight ear and grain morphological traits. A total of 108 QTL related to these traits were detected across four environments using an ultra-high density bin map constructed using recombinant inbred lines (RILs) derived from a cross between Ye478 and Qi319, and 26 QTL were identified in more than two environments. Furthermore, 64 single nucleotide polymorphisms (SNPs) were found to be significantly associated with the eight ear and grain morphological traits (?log₁₀(P)?>?4) in an association panel of 240 maize inbred lines. Combining the two mapping populations, a total of 17 pleiotropic QTL/SNPs (pQTL/SNPs) were associated with various traits across multiple environments. PKS2, a stable locus influencing kernel shape identified on chromosome 2 in a genome-wide association study (GWAS), was within the QTL confidence interval defined by the RILs. The candidate region harbored a short 13-Kb LD block encompassing four SNPs (SYN11386, PHM14783.16, SYN11392, and SYN11378). In the association panel, 13 lines derived from the hybrid PI78599 possessed the same allele as Qi319 at the PHM14783.16 (GG) locus, with an average value of 0.21 for KS, significantly lower than that of the 34 lines derived from Ye478 that carried a different allele (0.25, P?<?0.05). Therefore, further fine mapping of PKS2 will provide valuable information for understanding the genetic components of grain yield and improving molecular marker-assisted selection (MAS) in maize.

     

Genetic analysis of genetic basis of a physiological disorder “straighthead” in rice (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.)

Straighthead is a physiological disorder in rice that causes yield losses and is a serious threat to rice production worldwide. Identification of QTL conferring resistance will help develop resistant cultivars for straighthead control. We conducted linkage mapping to identify QTL involved with straighthead. The study was based on a F₂ population developed from a cross between ‘Zhe733(resistant)/R312(susceptible)’. Using phenotypic data of F₂ plants and their F_2:3 families, two major QTL, qSTH-2 and qSTH-8, were identified using bulked segregant analysis, explaining 11.1 and 28.1 % of the phenotypic variation on chromosome 2 and 8, respectively. The qSTH-2 for straighthead resistance was identified by linkage mapping. qSTH-2 was situated near a QTL “AsS” responsible for arsenic accumulation. Straighthead is frequently observed on land where As has accumulated. The result suggests a kind of internal connection between qSTH-2 and AsS. Additionally, the QTL qSTH-8 was located close to HD5 related with heading date. The close location may be associated with the observation of early heading among straighthead resistant varieties. These findings should be useful for further genetic study of straighthead.     

A morphophysiological analysis of the effects of drought and shade on <Emphasis Type="Italic">Catalpa bungei</Emphasis> plantlets

The interactive effects of shade and drought on the morphological and physiological traits of Catalpa bungei plantlets were assessed. Seedling growth, biomass, biomass allocation, leaf morphology, chlorophyll (Chl) content and gas-exchange parameters were measured in plants raised for 3 months under three light levels [80% (HI), 50% (MI), 30% (LI)] and two water levels [moisture (M) and drought (D)]. The results showed that shade greatly decreased growth, biomass, leaf area (LA) and Chl a/b; increased specific leaf area (SLA) and Chl content; and reduced photosynthetic rate (P _n). Drought reduced the growth, biomass, LA, SLA, Chl a/b, P _n, stomatal conductance (G _s), transpiration rate (T _r) and intercellular carbon dioxide concentration (C _i) and increased the Chl content. Stomatal closure was an early physiological response to water stress. Light, water and their interaction significantly affected plant traits and their bivariate relationships. The phenotypic plasticity index of light (0.47) was much higher than that of water (0.21), indicating that light was the main driver of the variations observed. Under drought stress, growth, biomass, leaf and stem biomass allocation significantly decreased in the HI and MI environments, whereas no significant difference was observed in growth or biomass parameters under the LI condition. Furthermore, no significant difference was observed in P _n, G _s, or T _r under the LI condition under water stress. Our results showed that shade did not alter the negative effects caused by drought stress in MI but did alleviate the negative effects of the LI condition. In summary, the effect of drought on C. bungei plantlets depends on the irradiance conditions.     

Genetic mapping of local adaptation along the altitudinal gradient in <Emphasis Type="Italic">Abies sachalinensis</Emphasis>

Understanding the genetic bases of local adaptation in dominant conifer species is critical in predicting the impacts of rapid climate change on forest ecosystems. However, the genetic basis of adaptation is not yet fully understood due to the huge and complex genomes of conifers and the unavailability to date of suitable crossing material. In this study, we constructed a linkage map for Abies sachalinensis (2n = 24) and investigated quantitative trait loci (QTLs) associated with local adaptation along an altitudinal gradient. A segregating population of 239 seedlings was produced from a cross between two F₁ hybrids (high-altitude × low-altitude genotypes). QTL mapping of phenological and growth traits was performed using a pseudo-testcross strategy with linkage maps based on 1251 single-nucleotide polymorphism (SNP) and three simple sequence repeat (SSR) markers. Two maps consisting of 12 linkage groups with an average marker interval of ca. 3 cM were constructed for each parent. The total lengths of the maps were 1861 and 1949 cM. A permutation test identified four significant QTLs and 11 additional suggestive QTLs, with high logarithm of odds (LOD) scores (> 3.0). This is the first highly saturated linkage map produced for Abies taxa. Our results suggest that spring bud phenology is controlled by several QTLs with moderate effects. The use of the mapping population created by crossing two hybrids (high × low altitude genotypes) and numerous SNP markers enabled us to investigate the genetic basis of adaptive traits in conifer species.     

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.

     

Molecular mapping of the genomic region conferring resistance to soybean stem canker in Hutcheson soybean

Genetic resistance to soybean stem canker, caused by the fungus Diaporthe phaseolorum var. meridionalis (Dpm), is controlled by five major, dominant, nonallelic genes Rdm1 to Rdm5. A genomic region containing the Rdm4 and Rdm5 genes was first described in Hutcheson soybean, where they were found to confer specific resistance to Argentinean physiological races of Dpm. Here, we report the genetic mapping of Rdm4 and Rdm5 loci using two pheno- and genotypically characterized F_2:3 populations derived from Hutcheson cultivar. The mapping populations were screened with amplified fragment length polymorphism (AFLP) markers using bulk segregant analysis, and with simple sequence repeat (SSR) markers. Linkage analysis indicated that the Rdm4 and Rdm5 resistance loci were located in a genomic region collinear with the molecular linkage group (MLG) A2 (chromosome 8) of the soybean genetic map. The linkage group contains two SSR markers, Sat_162 and Satt233, flanking the Rdm4 and Rdm5 loci. These SSR will be useful to increase the efficiency of selection in breeding programs aimed to incorporate Rdm4 and Rdm5 genes into soybean elite germplasm.     

Association of <Emphasis Type="Italic">ACYP2</Emphasis> and <Emphasis Type="Italic">TSPYL6</Emphasis> Genetic Polymorphisms with Risk of Ischemic Stroke in Han Chinese Population

The development of ischemic stroke is associated with advanced age. Telomere length, as a marker of biological aging, has been reported to influence the risk of several age-related diseases, including ischemic stroke. Recent studies have identified the genetic variant within ACYP2 and TSPYL6 associated with shorter telomere length. The objective of this study is to investigate the putative association of ischemic stroke with common polymorphisms in ACYP2 and TSPYL6 genes in a Chinese Han population. We found that the risk alleles of six single nucleotide polymorphisms (SNPs), including rs11125529, rs12615793, rs843711, rs11896604, and rs843706 within both ACYP2 and TSPYL6, and rs17045754 in ACYP2 gene, were related with increased risk of ischemic stroke according to both allelic and genotype association analyses. The significant correlations between ACYP2 and TSPYL6 SNPs and ischemic stroke risk were also observed in dominant, recessive, and additive models, respectively. Two blocks in high linkage disequilibrium were identified in this study, and two haplotypes were associated with higher ischemic stroke susceptibility. In conclusion, the genetic polymorphisms of ACYP2 and TSPYL6 are associated with increased risk of developing ischemic stroke. Further studies with larger sample sizes are required to validate our findings.     

Association of SNPs in dopamine and serotonin pathway genes and their interacting genes with temperament traits in Charolais cows

Cattle temperament is a complex trait, and molecular studies aimed at defining this trait are scarce. We used an interaction networks approach to identify new genes (interacting genes) and to estimate their effects and those of 19 dopamine- and serotonin-related genes on the temperament traits of Charolais cattle. The genes proopiomelanocortin (POMC), neuropeptide Y (NPY), solute carrier family 18, member 2 (SLC18A2) and FBJ murine osteosarcoma viral oncogene homologue (FOSFBJ) were identified as new candidates. Their potential to be associated with temperament was estimated according to their reported biological activities, which included interactions with neural activity, receptor function, targeting or synthesis of neurotransmitters and association with behaviour. Pen score (PS) and exit velocity (EV) measures were determined from 412 Charolais cows to calculate their temperament score (TS). Based on the TS, calm (n?=?55; TS, 1.09?±?0.33) and temperamental (n?=?58; TS, 2.27?±?0.639) cows were selected and genotyped using a 248 single-nucleotide variation (SNV) panel. Of the 248 variations in the panel, only 151 were confirmed to be polymorphic (single-nucleotide polymorphisms; SNPs) in the tested population. Single-marker association analyses between genotypes and temperament measures (EV, PS and/or TS) indicated significant associations of six SNPs from four candidate genes. The markers rs109576799 and rs43696138, located in the DRD3 and HTR2A genes, respectively, were significantly associated with both EV and TS traits. Four markers, rs110365063 and rs137756569 from the POMC gene and rs110365063 and rs135155082 located in SLC18A2 and DRD2, respectively, were associated with PS. The variant rs110365063 located in bovine SLC18A2 causes a change in the amino acid sequence from Ala to Thr. Further studies are needed to confirm the association of genetic profile with cattle temperament; however, our study represents important progress in understanding the regulation of cattle temperament by different genes with divergent functions.     

Genetic diversity near the <Emphasis Type="Italic">DGAT1-2</Emphasis> gene for high oleic acid content and kernel trait variation in a maize breeding collection

Previous studies in maize have identified three quantitative trait loci (QTLs) coding for high oleic acid content (HOAC) (QTLs oleic6-1, oleic6-2, and oleic6-3) at bins 6.04–6.05, proximal to the DGAT1-2 gene. The aims of this work were (i) to discover new markers for linkage disequilibrium (LD) and haplotype analysis distal to DGAT1-2, (ii) to develop a new DGAT1-2 PCR probe to detect the allele determining HOAC (F469 insertion), (iii) to conduct cluster analysis for kernel traits, and (iv) to assess genetic diversity, LD, and association analysis for kernel traits with the DGAT1-2 PCR probe and 13 markers previously mapped near DGAT1-2 in the Argentine temperate maize collection of 111 inbred lines. The results showed high haplotype diversity distal to DGAT1-2 and relatedness between the inbred line LP199 (with HOAC) and the Non-Stiff Stalk line W22 (reference genome). The frequency of F469 was low (20%). F469 was clustered with flint-grain type characteristics, whereas HOAC was associated with F469 across linear models. Genetic diversity at bins 6.04–6.05 was high (0.62), whereas LD extent was low (r²?≤?0.45). This low extent of LD indicates a high level of recombination and no LD between DGAT1-2 and markers flanking QTLs oleic6-1 to 3. Nevertheless, the significant LD between markers flanking those QTLs and the cosegregation of F469 with nc009 (markers flanking QTLs oleic6-1 and 2) during inbred line conversion suggest that these QTLs might contribute to HOAC in the breeding collection. However, further studies are needed to precise mapping at bins 6.04–6.05 for breeding purposes.