Genome-wide association mapping of canopy wilting in diverse soybean genotypes

Key message

Genome-wide association analysis identified 61 SNP markers for canopy wilting, which likely tagged 51 different loci. Based on the allelic effects of the significant SNPs, the slowest and fastest wilting genotypes were identified.

Abstract

Drought stress is a major global constraint for crop production, and slow canopy wilting is a promising trait for improving drought tolerance. The objective of this study was to identify genetic loci associated with canopy wilting and to confirm those loci with previously reported canopy wilting QTLs. A panel of 373 maturity group (MG) IV soybean genotypes was grown in four environments to evaluate canopy wilting. Statistical analysis of phenotype indicated wide variation for the trait, with significant effects of genotype (G), environment (E), and G × E interaction. Over 42,000 SNP markers were obtained from the Illumina Infinium SoySNP50K iSelect SNP Beadchip. After filtration for quality control, 31,260 SNPs with a minor allele frequency (MAF) ≥5% were used for association mapping using the Fixed and random model Circulating Probability Unification (FarmCPU) model. There were 61 environment-specific significant SNP-canopy wilting associations, and 21 SNPs that associated with canopy wilting in more than one environment. There were 34 significant SNPs associated with canopy wilting when averaged across environments. Together, these SNPs tagged 23 putative loci associated with canopy wilting. Six of the putative loci were located within previously reported chromosomal regions that were associated with canopy wilting through bi-parental mapping. Several significant SNPs were located within a gene or very close to genes that had a reported biological connection to transpiration or water transport. Favorable alleles from significant SNPs may be an important resource for pyramiding genes to improve drought tolerance and for identifying parental genotypes for use in breeding programs.

     

Association Mapping of Total Carotenoids in Diverse Soybean Genotypes Based on Leaf Extracts and High-Throughput Canopy Spectral Reflectance Measurements

Carotenoids are organic pigments that are produced predominantly by photosynthetic organisms and provide antioxidant activity to a wide variety of plants, animals, bacteria, and fungi. The carotenoid biosynthetic pathway is highly conserved in plants and occurs mostly in chromoplasts and chloroplasts. Leaf carotenoids play important photoprotective roles and targeted selection for leaf carotenoids may offer avenues to improve abiotic stress tolerance. A collection of 332 soybean [Glycine max (L.) Merr.] genotypes was grown in two years and total leaf carotenoid content was determined using three different methods. The first method was based on extraction and spectrophotometric determination of carotenoid content (eCaro) in leaf tissue, whereas the other two methods were derived from high-throughput canopy spectral reflectance measurements using wavelet transformed reflectance spectra (tCaro) and a spectral reflectance index (iCaro). An association mapping approach was employed using 31,253 single nucleotide polymorphisms (SNPs) to identify SNPs associated with total carotenoid content using a mixed linear model based on data from two growing seasons. A total of 28 SNPs showed a significant association with total carotenoid content in at least one of the three approaches. These 28 SNPs likely tagged 14 putative loci for carotenoid content. Six putative loci were identified using eCaro, five loci with tCaro, and nine loci with iCaro. Three of these putative loci were detected by all three carotenoid determination methods. All but four putative loci were located near a known carotenoid-related gene. These results showed that carotenoid markers can be identified in soybean using extract-based as well as by high-throughput canopy spectral reflectance-based approaches, demonstrating the utility of field-based canopy spectral reflectance phenotypes for association mapping.     

Time series canopy phenotyping enables the identification of genetic variants controlling dynamic phenotypes in soybean

Advances in plant phenotyping technologies are dramatically reducing the marginal costs of collecting multiple phenotypic measurements across several time points. Yet, most current approaches and best statistical practices implemented to link genetic and phenotypic variation in plants have been developed in an era of single-time-point data. Here, we used time-series phenotypic data collected with an unmanned aircraft system for a large panel of soybean (Glycine max (L.) Merr.) varieties to identify previously uncharacterized loci. Specifically, we focused on the dissection of canopy coverage (CC) variation from this rich data set. We also inferred the speed of canopy closure, an additional dimension of CC, from the time-series data, as it may represent an important trait for weed control. Genome-wide association studies (GWASs) identified 35 loci exhibiting dynamic associations with CC across developmental stages. The time-series data enabled the identification of 10 known flowering time and plant height quantitative trait loci (QTLs) detected in previous studies of adult plants and the identification of novel QTLs influencing CC. These novel QTLs were disproportionately likely to act earlier in development, which may explain why they were missed in previous single-time-point studies. Moreover, this time-series data set contributed to the high accuracy of the GWASs, which we evaluated by permutation tests, as evidenced by the repeated identification of loci across multiple time points. Two novel loci showed evidence of adaptive selection during domestication, with different genotypes/haplotypes favored in different geographic regions. In summary, the time-series data, with soybean CC as an example, improved the accuracy and statistical power to dissect the genetic basis of traits and offered a promising opportunity for crop breeding with quantitative growth curves.     

Genome-wide association study for flowering time,maturity dates and plant height in early maturing soybean (Glycine max) germplasm

Background

Soybean (Glycine max) is a photoperiod-sensitive and self-pollinated species. Days to flowering (DTF) and maturity (DTM), duration of flowering-to-maturity (DFTM) and plant height (PH) are crucial for soybean adaptability and yield. To dissect the genetic architecture of these agronomically important traits, a population consisting of 309 early maturity soybean germplasm accessions was genotyped with the Illumina Infinium SoySNP50K BeadChip and phenotyped in multiple environments. A genome-wide association study (GWAS) was conducted using a mixed linear model that involves both relative kinship and population structure.

Results

The linkage disequilibrium (LD) decayed slowly in soybean, and a substantial difference in LD pattern was observed between euchromatic and heterochromatic regions. A total of 27, 6, 18 and 27 loci for DTF, DTM, DFTM and PH were detected via GWAS, respectively. The Dt1 gene was identified in the locus strongly associated with both DTM and PH. Ten candidate genes homologous to Arabidopsis flowering genes were identified near the peak single nucleotide polymorphisms (SNPs) associated with DTF. Four of them encode MADS-domain containing proteins. Additionally, a pectin lyase-like gene was also identified in a major-effect locus for PH where LD decayed rapidly.

Conclusions

This study identified multiple new loci and refined chromosomal regions of known loci associated with DTF, DTM, DFTM and/or PH in soybean. It demonstrates that GWAS is powerful in dissecting complex traits and identifying candidate genes although LD decayed slowly in soybean. The loci and trait-associated SNPs identified in this study can be used for soybean genetic improvement, especially the major-effect loci associated with PH could be used to improve soybean yield potential. The candidate genes may serve as promising targets for studies of molecular mechanisms underlying the related traits in soybean.

Electronic supplementary material

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

Genome-wide association mapping of quantitative resistance to sudden death syndrome in soybean

Background

Sudden death syndrome (SDS) is a serious threat to soybean production that can be managed with host plant resistance. To dissect the genetic architecture of quantitative resistance to the disease in soybean, two independent association panels of elite soybean cultivars, consisting of 392 and 300 unique accessions, respectively, were evaluated for SDS resistance in multiple environments and years. The two association panels were genotyped with 52,041 and 5,361 single nucleotide polymorphisms (SNPs), respectively. Genome-wide association mapping was carried out using a mixed linear model that accounted for population structure and cryptic relatedness.

Result

A total of 20 loci underlying SDS resistance were identified in the two independent studies, including 7 loci localized in previously mapped QTL intervals and 13 novel loci. One strong peak of association on chromosome 18, associated with all disease assessment criteria across the two panels, spanned a physical region of 1.2 Mb around a previously cloned SDS resistance gene (GmRLK18-1) in locus Rfs2. An additional variant independently associated with SDS resistance was also found in this genomic region. Other peaks were within, or close to, sequences annotated as homologous to genes previously shown to be involved in plant disease resistance. The identified loci explained an average of 54.5% of the phenotypic variance measured by different disease assessment criteria.

Conclusions

This study identified multiple novel loci and refined the map locations of known loci related to SDS resistance. These insights into the genetic basis of SDS resistance can now be used to further enhance durable resistance to SDS in soybean. Additionally, the associations identified here provide a basis for further efforts to pinpoint causal variants and to clarify how the implicated genes affect SDS resistance in soybean.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-809) contains supplementary material, which is available to authorized users.     

Genome-wide association study (GWAS) of salt tolerance in worldwide soybean germplasm lines

Salt is a severe abiotic stress causing soybean yield loss in saline soils and irrigated fields. Marker-assisted selection (MAS) is a powerful genomic tool for improving the efficiency of breeding salt-tolerant soybean varieties. The objectives of this study were to uncover novel single-nucleotide polymorphisms (SNPs) and quantitative trait loci (QTLs) associated with salt tolerance and to confirm the previously identified genomic regions and SNPs for salt tolerance. A total of 283 diverse soybean plant introductions (PIs) were screened for salt tolerance in the greenhouse based on leaf chloride concentrations and leaf chlorophyll concentrations after 12–18 days of 120-mM NaCl treatment. A total of 33,009 SNPs across 283 genotypes from the Illumina Infinium SoySNP50K BeadChip database were employed in the association analysis with leaf chloride concentrations and leaf chlorophyll concentrations. Genome-wide association mapping showed that 45 SNPs representing nine genomic regions on chromosomes (Chr.) 2, 3, 7, 8, 10, 13, 14, 16, and 20 were significantly associated with both leaf chloride concentrations and leaf chlorophyll concentrations in 2014, 2015, and combined years. A total of 31 SNPs on Chr. 3 were mapped at or near the previously reported major salt tolerance QTL. The significant SNP on Chr. 2 was also in proximity to the previously reported SNP for salt tolerance. The other significant SNPs represent seven putative novel QTLs for salt tolerance. The significant SNP markers on Chr. 2, 3, 14, 16, and 20, which were identified in both general linear model and mixed linear model, were highly recommended for MAS in breeding salt-tolerant soybean varieties.     

Identification of loci governing eight agronomic traits using a GBS‐GWAS approach and validation by QTL mapping in soya bean

Soya bean is a major source of edible oil and protein for human consumption as well as animal feed. Understanding the genetic basis of different traits in soya bean will provide important insights for improving breeding strategies for this crop. A genome‐wide association study (GWAS) was conducted to accelerate molecular breeding for the improvement of agronomic traits in soya bean. A genotyping‐by‐sequencing (GBS) approach was used to provide dense genome‐wide marker coverage (>47 000 SNPs) for a panel of 304 short‐season soya bean lines. A subset of 139 lines, representative of the diversity among these, was characterized phenotypically for eight traits under six environments (3 sites × 2 years). Marker coverage proved sufficient to ensure highly significant associations between the genes known to control simple traits (flower, hilum and pubescence colour) and flanking SNPs. Between one and eight genomic loci associated with more complex traits (maturity, plant height, seed weight, seed oil and protein) were also identified. Importantly, most of these GWAS loci were located within genomic regions identified by previously reported quantitative trait locus (QTL) for these traits. In some cases, the reported QTLs were also successfully validated by additional QTL mapping in a biparental population. This study demonstrates that integrating GBS and GWAS can be used as a powerful complementary approach to classical biparental mapping for dissecting complex traits in soya bean.     

Evaluation of Clustering and Genotype Distribution for Replication in Genome Wide Association Studies: The Age-Related Eye Disease Study

Genome-wide association studies (GWASs) assess correlation between traits and DNA sequence variation using large numbers of genetic variants such as single nucleotide polymorphisms (SNPs) distributed across the genome. A GWAS produces many trait-SNP associations with low p-values, but few are replicated in subsequent studies. We sought to determine if characteristics of the genomic loci associated with a trait could be used to identify initial associations with a higher chance of replication in a second cohort. Data from the age-related eye disease study (AREDS) of 100,000 SNPs on 395 subjects with and 198 without age-related macular degeneration (AMD) were employed. Loci highly associated with AMD were characterized based on the distribution of genotypes, level of significance, and clustering of adjacent SNPs also associated with AMD suggesting linkage disequilibrium or multiple effects. Forty nine loci were highly associated with AMD, including 3 loci (CFH, C2/BF, LOC387715/HTRA1) already known to contain important genetic risks for AMD. One additional locus (C3) reported during the course of this study was identified and replicated in an additional study group. Tag-SNPs and haplotypes for each locus were evaluated for association with AMD in additional cohorts to account for population differences between discovery and replication subjects, but no additional clearly significant associations were identified. Relying on a significant genotype tests using a log-additive model would have excluded 57% of the non-replicated and none of the replicated loci, while use of other SNP features and clustering might have missed true associations.     

Association between Expression Quantitative Trait Loci and Metabolic Traits in Two Korean Populations

Most genome-wide association studies consider genes that are located closest to single nucleotide polymorphisms (SNPs) that are highly significant for those studies. However, the significance of the associations between SNPs and candidate genes has not been fully determined. An alternative approach that used SNPs in expression quantitative trait loci (eQTL) was reported previously for Crohn’s disease; it was shown that eQTL-based preselection for follow-up studies was a useful approach for identifying risk loci from the results of moderately sized GWAS. In this study, we propose an approach that uses eQTL SNPs to support the functional relationships between an SNP and a candidate gene in a genome-wide association study. The genome-wide SNP genotypes and 10 biochemical measures (fasting glucose levels, BUN, serum albumin levels, AST, ALT, gamma GTP, total cholesterol, HDL cholesterol, triglycerides, and LDL cholesterol) were obtained from the Korean Association Resource (KARE) consortium. The eQTL SNPs were isolated from the SNP dataset based on the RegulomeDB eQTL-SNP data from the ENCODE projects and two recent eQTL reports. A total of 25,658 eQTL SNPs were tested for their association with the 10 metabolic traits in 2 Korean populations (Ansung and Ansan). The proportion of phenotypic variance explained by eQTL and non-eQTL SNPs showed that eQTL SNPs were more likely to be associated with the metabolic traits genetically compared with non-eQTL SNPs. Finally, via a meta-analysis of the two Korean populations, we identified 14 eQTL SNPs that were significantly associated with metabolic traits. These results suggest that our approach can be expanded to other genome-wide association studies.     

Genomic analysis of genetic markers associated with inherited cardiomyopathy (round heart disease) in the turkey (Meleagris gallopavo)

In turkeys, spontaneous cardiomyopathy or round heart (RH) disease is characterised by dilated ventricles and cardiac muscle hypertrophy. Although the aetiology of RH is still unknown, the disease can have a significant economic impact on turkey producers. In an initial attempt to identify genomic regions associated with RH, we utilised the chicken genome sequence to target short DNA sequences (sequence-characterised amplified regions, SCARs) identified in previous studies that had significant differences in frequency distribution between RH+ and RH- turkeys. SCARs were comparatively aligned with the chicken whole-genome sequence to identify flanking regions for primer design. Primers from 32 alignments were tested and target sequences were successfully amplified for 30 loci (94%). Comparative re-sequencing identified putative SNPs in 20 of the 30 loci (67%). Genetically informative SNPs at 16 loci were genotyped in the UMN/NTBF turkey mapping population. As a result of this study, 34 markers were placed on the turkey/chicken comparative map and 15 markers were added to the turkey genetic linkage map. The position of these markers relative to cardiac-related genes is presented. In addition, analysis of genotypes at 109 microsatellite loci presumed to flank the SCAR sequences in the turkey genome identified four significant associations with RH.     

Identification of single nucleotide polymorphisms and haplotypes associated with yield and yield components in soybean (<Emphasis Type="Italic">Glycine max</Emphasis>) landraces across multiple environments

Genome-wide association analysis is a powerful approach to identify the causal genetic polymorphisms underlying complex traits. In this study, we evaluated a population of 191 soybean landraces in five environments to detect molecular markers associated with soybean yield and its components using 1,536 single-nucleotide polymorphisms (SNPs) and 209 haplotypes. The analysis revealed that abundant phenotypic and genetic diversity existed in the studied population. This soybean population could be divided into two subpopulations and no or weak relatedness was detected between pair-wise landraces. The level of intra-chromosomal linkage disequilibrium was about 500 kb. Genome-wide association analysis based on the unified mixed model identified 19 SNPs and 5 haplotypes associated with soybean yield and yield components in three or more environments. Nine markers were found co-associated with two or more traits. Many markers were located in or close to previously reported quantitative trait loci mapped by linkage analysis. The SNPs and haplotypes identified in this study will help to further understand the genetic basis of soybean yield and its components, and may facilitate future high-yield breeding by marker-assisted selection in soybean.     

Genome-Wide Association of Stem Water Soluble Carbohydrates in Bread Wheat

Water soluble carbohydrates (WSC) in stems play an important role in buffering grain yield in wheat against biotic and abiotic stresses; however, knowledge of genes controlling WSC is very limited. We conducted a genome-wide association study (GWAS) using a high-density 90K SNP array to better understand the genetic basis underlying WSC, and to explore marker-based breeding approaches. WSC was evaluated in an association panel comprising 166 Chinese bread wheat cultivars planted in four environments. Fifty two marker-trait associations (MTAs) distributed across 23 loci were identified for phenotypic best linear unbiased estimates (BLUEs), and 11 MTAs were identified in two or more environments. Liner regression showed a clear dependence of WSC BLUE scores on numbers of favorable (increasing WSC content) and unfavorable alleles (decreasing WSC), indicating that genotypes with higher numbers of favorable or lower numbers of unfavorable alleles had higher WSC content. In silico analysis of flanking sequences of trait-associated SNPs revealed eight candidate genes related to WSC content grouped into two categories based on the type of encoding proteins, namely, defense response proteins and proteins triggered by environmental stresses. The identified SNPs and candidate genes related to WSC provide opportunities for breeding higher WSC wheat cultivars.     

A polymorphism in the HLA-DPB1 gene is associated with susceptibility to multiple sclerosis

We conducted an association study across the human leukocyte antigen (HLA) complex to identify loci associated with multiple sclerosis (MS). Comparing 1927 SNPs in 1618 MS cases and 3413 controls of European ancestry, we identified seven SNPs that were independently associated with MS conditional on the others (each P ≤ 4 x 10(-6)). All associations were significant in an independent replication cohort of 2212 cases and 2251 controls (P ≤ 0.001) and were highly significant in the combined dataset (P ≤ 6 x 10(-8)). The associated SNPs included proxies for HLA-DRB1*15:01 and HLA-DRB1*03:01, and SNPs in moderate linkage disequilibrium (LD) with HLA-A*02:01, HLA-DRB1*04:01 and HLA-DRB1*13:03. We also found a strong association with rs9277535 in the class II gene HLA-DPB1 (discovery set P = 9 x 10(-9), replication set P = 7 x 10(-4), combined P = 2 x 10(-10)). HLA-DPB1 is located centromeric of the more commonly typed class II genes HLA-DRB1, -DQA1 and -DQB1. It is separated from these genes by a recombination hotspot, and the association is not affected by conditioning on genotypes at DRB1, DQA1 and DQB1. Hence rs9277535 represents an independent MS-susceptibility locus of genome-wide significance. It is correlated with the HLA-DPB1*03:01 allele, which has been implicated previously in MS in smaller studies. Further genotyping in large datasets is required to confirm and resolve this association.     

Identification of SNPs tightly linked to the QTL for pod shattering in soybean

The pod shattering or dehiscence is essential for the propagation of pod-bearing plant species in the wild, but it causes significant yield losses during harvest of domesticated crop plants. Identifying novel molecular makers, which are linked to seed-shattering genes, is needed to employ the molecular marker-assisted selection for efficiently developing shattering-resistant soybean varieties. In this study, a genetic linkage map was constructed using 115 recombinant inbred lines (RILs) developed from crosses between the pod shattering susceptible variety, Keunol, and resistant variety, Sinpaldal. A 180 K Axiom® SoyaSNPs data and pod shattering data from two environments in 2001 and 2015 were used to identify quantitative trait loci (QTL) for pod shattering. A major QTL was identified between two flanking single nucleotide polymorphism (SNP) markers, AX-90320801 and AX-90306327 on chromosome 16 with 1.3 cM interval, 857 kb of physical range. In sequence, genotype distribution analysis was conducted using extreme phenotype RILs. This could narrow down the QTL down to 153 kb on the physical map and was designated as qPDH1-KS with 6 annotated gene models. All exons within qPDH1-KS were sequenced and the 6 polymorphic SNPs affecting the amino acid sequence were identified. To develop universally available molecular markers, 38 Korean soybean cultivars were investigated by the association study using the 6 identified SNPs. Only two SNPs were strongly associated with the pod shattering. These two identified SNPs will help to identify the pod shattering responsible gene and to develop pod shattering-resistant soybean plants using marker-assisted selection.     

Understanding genetic control of root system architecture in soybean: Insights into the genetic basis of lateral root number

Developing crops with better root systems is a promising strategy to ensure productivity in both optimum and stress environments. Root system architectural traits in 397 soybean accessions were characterized and a high‐density single nucleotide polymorphisms (SNPs)‐based genome‐wide association study was performed to identify the underlying genes associated with root structure. SNPs associated with root architectural traits specific to landraces and elite germplasm pools were detected. Four loci were detected in landraces for lateral root number (LRN) and distribution of root thickness in diameter Class I with a major locus on chromosome 16. This major loci was detected in the coding region of unknown protein, and subsequent analyses demonstrated that root traits are affected with mutated haplotypes of the gene. In elite germplasm pool, 3 significant SNPs in alanine‐glyoxalate aminotransferase, Leucine‐Rich Repeat receptor/No apical meristem, and unknown functional genes were found to govern multiple traits including root surface area and volume. However, no major loci were detected for LRN in elite germplasm. Nucleotide diversity analysis found evidence of selective sweeps around the landraces LRN gene. Soybean accessions with minor and mutated allelic variants of LRN gene were found to perform better in both water‐limited and optimal field conditions.     

Variants in Notch signalling pathway genes,PSEN1 and MAML2, predict overall survival in Chinese patients with epithelial ovarian cancer

To identify genetic variants in Notch signalling pathway genes that may predict survival of Han Chinese patients with epithelial ovarian cancer (EOC), we analysed a total of 1273 single nucleotide polymorphisms (SNPs) within 75 Notch genes in 480 patients from a published EOC genomewide association study (GWAS). We found that PSEN1 rs165934 and MAML2 rs76032516 were associated with overall survival (OS) of patients by multivariate Cox proportional hazards regression analysis. Specifically, the PSEN1 rs165934 AA genotype was associated with a poorer survival (adjusted hazards ratio [adjHR] = 1.41, 95% CI = 1.07‐1.84, and P = .014), compared with the CC + CA genotype, while MAML2 rs76032516 AA + AC genotypes were associated with a poorer survival (adjHR = 1.58, 95% CI = 1.16‐2.14, P = .004), compared with the CC genotype. The combined analysis of these two SNPs revealed that the death risk increased as the number of unfavourable genotypes increased in a dose‐dependent manner (P_trend < .001). Additionally, the expression quantitative trait loci analysis revealed that the SNP rs165932 in the rs165934 LD block (r² = .946) was associated with expression levels of PSEN1, which might be responsible for the observed association with SNP rs165934. The associations of PSEN1 rs165934 and MAML2 rs76032516 of the Notch signalling pathway genes with OS in Chinese EOC patients are novel findings, which need to be validated in other large and independent studies.     

Genome-wide association study for soybean cyst nematode resistance in Chinese elite soybean cultivars

Soybean cyst nematode (SCN) (Heterodera glycines Ichinohe) is a highly recalcitrant endoparasite of soybean roots, causing more yield loss than any other pest. To identify quantitative trait loci (QTL) controlling resistance to SCN (HG type 2.5.7, race 1), a genome-wide association study (GWAS) was performed. The association panel, consisting of 120 Chinese soybean cultivars, was genotyped with 7189 single nucleotide polymorphism (SNPs). A total of 6204 SNPs with minor allele frequency >0.05 were used to estimate linkage disequilibrium (LD) and population structure. The mean level of LD measured by r ² declined very rapidly to half its maximum value (0.51) at 220 kb. The overall population structure was approximately coincident with geographic origin. The GWAS results identified 13 SNPs in 7 different genomic regions significantly associated with SCN resistance. Of these, three SNPs were localized in previously mapped QTL intervals, including rhg1 and Rhg4. The GWAS results also detected 10 SNPs in 5 different genomic regions associated with SCN resistance. The identified loci explained an average of 95.5% of the phenotypic variance. The proportion of phenotypic variance was due to additive genetic variance of the validated SNPs. The present study identified multiple new loci and refined chromosomal regions of known loci associated with SCN resistance. The loci and trait-associated SNPs identified in this study can be used for developing soybean cultivars with durable resistance against SCN.     

Combined linkage mapping and association analysis reveals genetic control of maize kernel moisture content

The free moisture in crop kernels after being naturally dried is referred to as kernel moisture content (KMC). Maize KMC reflects grain quality and influences transportation and storage of seeds. We used an IBM Syn10 DH maize population consisting of 249 lines and an association panel comprising 310 maize inbred lines to identify the genetic loci affecting maize KMC in three environments. Using the IBM population detected 13 QTL on seven chromosomes, which were clustered into nine common QTL. Genome-wide association analysis (GWAS) identified 16 significant SNPs across the 3 environments, which were linked to 158 genes across the three environments. Combined QTL mapping and GWAS found two SNPs that were located in two of the mapped QTL, respectively. Twenty-three genes were linked with the loci co-localized in both populations. Of these 181 genes, five have previously been reported to be associated with KMC or to regulate seed development. These associations were verified by candidate gene association analysis. Two superior alleles and one favorable haplotype for Zm00001d007774 and Zm00001d047868 were found to influence KMC. These findings provide insights into molecular mechanisms underlying maize KMC and contribute to the use of marker-assisted selection for breeding low-KMC maize.     

单羧酸转运蛋白基因SNP 与肝细胞肝癌预后的关联研究

目的:探讨单羧酸转运蛋白基因(monocarboxylate transporter,MCT)单核苷酸多态性(single nucleotide polymorphism,SNPs)与肝细胞肝癌(hepatocellular carcinoma,HCC)根治术患者预后的关系。方法:运用Sequenom i PLEX分型技术对830例原发性HCC患者MCT家族(MCT1、MCT2和MCT4)基因上的8个功能性SNP位点进行基因分型,并分析这些SNP与HCC患者预后的相关性。结果:MCT1基因rs1049434位点和MCT2基因rs995343位点基因型与HCC患者总体生存期及无复发生存期均显著相关(P0.05)。携带MCT1 AT基因型或TT基因型的患者死亡及复发风险均显著低于携带AA基因型的患者(HR=0.72;P=0.042或HR=0.64;P=0.002);携带MCT2 CT基因型或TT基因型的患者死亡及复发风险均显著高于携带CC基因型的患者(HR=1.64;P=0.018或HR=1.52;P=0.026)。而且,MCT1基因rs1049434位点和MCT2基因rs995343位点对HCC预后存在显著的累积效应,携带2个危险基因型的患者死亡及复发风险分别是没有危险基因型患者的2.16倍和2.54倍。此外,携带2个危险基因型的HCC患者在术后行TACE辅助治疗后死亡及复发风险均显著降低(P0.05)。结论:MCT1和MCT2基因上的功能性SNP位点有可能作为HCC根治术后预后评估和TACE辅助治疗反应预测的独立标志物。     

Coherent synthesis of genomic associations with phenotypes and home environments

Local adaptation is often studied via (i) multiple common garden experiments comparing performance of genotypes in different environments and (ii) sequencing genotypes from multiple locations and characterizing geographic patterns in allele frequency. Both approaches aim to characterize the same pattern (local adaptation), yet the complementary information from each has not yet been coherently integrated. Here, we develop a genome‐wide association model of genotype interactions with continuous environmental gradients (G × E), that is reaction norms. We present an approach to impute relative fitness, allowing us to coherently synthesize evidence from common garden and genome–environment associations. Our approach identifies loci exhibiting environmental clines where alleles are associated with higher fitness in home environments. Simulations show our approach can increase power to detect loci causing local adaptation. In a case study on Arabidopsis thaliana, most identified SNPs exhibited home allele advantage and fitness trade‐offs along climate gradients, suggesting selective gradients can maintain allelic clines. SNPs exhibiting G × E associations with fitness were enriched in genic regions, putative partial selective sweeps and associations with an adaptive phenotype (flowering time plasticity). We discuss extensions for situations where only adaptive phenotypes other than fitness are available. Many types of data may point towards the loci underlying G × E and local adaptation; coherent models of diverse data provide a principled basis for synthesis.