Negative dominance and dominance-by-dominance epistatic effects reduce grain-yield heterosis in wide crosses in wheat.

The genetics underlying heterosis, the difference in performance of crosses compared with midparents, is hypothesized to vary with relatedness between parents. We established a unique germplasm comprising three hybrid wheat sets differing in the degree of divergence between parents and devised a genetic distance measure giving weight to heterotic loci. Heterosis increased steadily with heterotic genetic distance for all 1903 hybrids. Midparent heterosis, however, was significantly lower in the hybrids including crosses between elite and exotic lines than in crosses among elite lines. The analysis of the genetic architecture of heterosis revealed this to be caused by a higher portion of negative dominance and dominance-by-dominance epistatic effects. Collectively, these results expand our understanding of heterosis in crops, an important pillar toward global food security.

Whole-Genome Association Mapping and Genomic Prediction for Iron Concentration in Wheat Grains.

Malnutrition of iron (Fe) affects two billion people worldwide. Therefore, enhancing grain Fe concentration (GFeC) in wheat (Triticum aestivum L.) is an important goal for breeding. Here we study the genetic factors underlying GFeC trait by genome-wide association studies (GWAS) and the prediction abilities using genomic prediction (GP) in a panel of 369 European elite wheat varieties which was genotyped with 15,523 mapped single-nucleotide polymorphism markers (SNP) and a subpanel of 183 genotypes with 44,233 SNP markers. The resulting means of GFeC from three field experiments ranged from 24.42 to 52.42 µg·g-1 with a broad-sense heritability (H²) equaling 0.59 over the years. GWAS revealed 41 and 137 significant SNPs in the whole and subpanel, respectively, including significant marker-trait associations (MTAs) for best linear unbiased estimates (BLUEs) of GFeC on chromosomes 2A, 3B and 5A. Putative candidate genes such as NAC transcription factors and transmembrane proteins were present on chromosome 2A (763,689,738⁻765,710,113 bp). The GP for a GFeC trait ranged from low to moderate values. The current study reported GWAS of GFeC for the first time in hexaploid wheat varieties. These findings confirm the utility of GWAS and GP to explore the genetic architecture of GFeC for breeding programs aiming at the improvement of wheat grain quality.

Identifying Candidate Genes for Enhancing Grain Zn Concentration in Wheat.

Wheat (Triticum aestivum L.) is one of the major staple food crops worldwide. Despite efforts in improving wheat quality, micronutrient levels are still below the optimal range for human nutrition. In particular, zinc (Zn) deficiency is a widespread problem in human nutrition in countries relying mainly on a cereal diet; hence improving Zn accumulation in grains is an imperative need. This study was designed to understand the genetic architecture of Zn grain concentrations in wheat grains. We performed a genome-wide association study (GWAS) for grain Zn concentrations in 369 European wheat genotypes, using field data from 3 years. The complete wheat panel was genotyped by high-density arrays of single nucleotide polymorphic (SNP) markers (90k iSELECT Infinium and 35k Affymetrix arrays) resulting in 15,523 polymorphic markers. Additionally, a subpanel of 183 genotypes was analyzed with a novel 135k Affymetrix marker array resulting in 28,710 polymorphic SNPs for high-resolution mapping of the potential genomic regions. The mean grain Zn concentration of the genotypes ranged from 25.05-52.67 µg g-1 dry weight across years with a moderate heritability value. Notably, 40 marker-trait associations (MTAs) were detected in the complete panel of varieties on chromosomes 2A, 3A, 3B, 4A, 4D, 5A, 5B, 5D, 6D, 7A, 7B, and 7D. The number of MTAs in the subpanel was increased to 161 MTAs whereas the most significant and consistent associations were located on chromosomes 3B (723,504,241-723,611,488 bp) and 5A (462,763,758-466,582,184 bp) having major effects. These genomic regions include newly identified putative candidate genes, which are related to Zn uptake and transport or represent bZIP and mitogen-activated protein kinase genes. These findings provide the basis for understanding the genetic background of Zn concentration in wheat grains that in turn may help breeders to select high Zn-containing genotypes to improve human health and grain quality.

Genome-wide association mapping of resistance to eyespot disease (Pseudocercosporella herpotrichoides) in European winter wheat (Triticum aestivum L.) and fine-mapping of Pch1.

KEY MESSAGE: Genotypes with recombination events in the Triticum ventricosum introgression on chromosome 7D allowed to fine-map resistance gene Pch1, the main source of eyespot resistance in European winter wheat cultivars. Eyespot (also called Strawbreaker) is a common and serious fungal disease of winter wheat caused by the necrotrophic fungi Oculimacula yallundae and Oculimacula acuformis (former name Pseudocercosporella herpotrichoides). A genome-wide association study (GWAS) for eyespot was performed with 732 microsatellite markers (SSR) and 7761 mapped SNP markers derived from the 90 K iSELECT wheat array using a panel of 168 European winter wheat varieties as well as three spring wheat varieties and phenotypic evaluation of eyespot in field tests in three environments. Best linear unbiased estimations (BLUEs) were calculated across all trials and ranged from 1.20 (most resistant) to 5.73 (most susceptible) with an average value of 4.24 and a heritability of H 2 = 0.91. A total of 108 SSR and 235 SNP marker-trait associations (MTAs) were identified by considering associations with a -log10 (P value) ≥3.0. Significant MTAs for eyespot-score BLUEs were found on chromosomes 1D, 2A, 2D, 3D, 5A, 5D, 6A, 7A and 7D for the SSR markers and chromosomes 1B, 2A, 2B, 2D, 3B and 7D for the SNP markers. For 18 varieties (10.5%), a highly resistant phenotype was detected that was linked to the presence of the resistance gene Pch1 on chromosome 7D. The identification of genotypes with recombination events in the introgressed genomic segment from Triticum ventricosum harboring the Pch1 resistance gene on chromosome 7DL allowed the fine-mapping of this gene using additional SNP markers and a potential candidate gene Traes_7DL_973A33763 coding for a CC-NBS-LRR class protein was identified.

Diversity analysis and genomic prediction of Sclerotinia resistance in sunflower using a new 25 K SNP genotyping array.

KEY MESSAGE: We have developed a SNP array for sunflower containing more than 25 K markers, representing single loci mostly in or near transcribed regions of the genome. The array was successfully applied to genotype a diversity panel of lines, hybrids, and mapping populations and represented well the genetic diversity of cultivated sunflower. Results of PCoA and population substructure analysis underlined the complexity of the genetic composition of current elite breeding material. The performance of this genotyping platform for genome-based prediction of phenotypes and detection of QTL with improved resolution could be demonstrated based on the re-evaluation of a population segregating for resistance to Sclerotinia midstalk rot. Given our results, the newly developed 25 K SNP array is expected to be of great utility for the most important applications in genome-based sunflower breeding and research. ABSTRACT: Genotyping with a large number of molecular markers is a prerequisite to conduct genome-based genetic analyses with high precision. Here, we report the design and performance of a 25 K SNP genotyping array for sunflower (Helianthus annuus L.). SNPs were discovered based on variant calling in de novo assembled, UniGene-based contigs of sunflower derived from whole genome sequencing and amplicon sequences originating from four and 48 inbred lines, respectively. After inclusion of publically available transcriptome-derived SNPs, in silico design of the Illumina(®) Infinium iSelect HD BeadChip yielded successful assays for 22,299 predominantly haplotype-specific SNPs. The array was validated in a sunflower diversity panel including inbred lines, open-pollinated varieties, introgression lines, landraces, recombinant inbred lines, and F2 populations. Validation provided 20,502 high-quality bi-allelic SNPs with stable cluster performance whereby each SNP marker represents a single locus mostly in or near transcribed regions of the sunflower genome. Analyses of population structure and quantitative resistance to Sclerotinia midstalk rot demonstrate that this array represents a significant improvement over currently available genomic tools for genetic diversity analyses, genome-wide marker-trait association studies, and genetic mapping in sunflower.

Analysis of main effect QTL for thousand grain weight in European winter wheat (Triticum aestivum L.) by genome-wide association mapping.

Grain weight, an essential yield component, is under strong genetic control and at the same time markedly influenced by the environment. Genetic analysis of the thousand grain weight (TGW) by genome-wide association study (GWAS) was performed with a panel of 358 European winter wheat (Triticum aestivum L.) varieties and 14 spring wheat varieties using phenotypic data of field tests in eight environments. Wide phenotypic variations were indicated for the TGW with BLUEs (best linear unbiased estimations) values ranging from 35.9 to 58.2 g with a mean value of 45.4 g and a heritability of H(2) = 0.89. A total of 12 candidate genes for plant height, photoperiodism and grain weight were genotyped on all varieties. Only three candidates, the photoperiodism gene Ppd-D1, dwarfing gene Rht-B1and the TaGW-6A gene were significant explaining up to 14.4, 2.3, and 3.4% of phenotypic variation, respectively. For a comprehensive genome-wide analysis of TGW-QTL genotyping data from 732 microsatellite markers and a set of 7769 mapped SNP-markers genotyped with the 90k iSELECT array were analyzed. In total, 342 significant (-log10 (P-value) ≥ 3.0) marker trait associations (MTAs) were detected for SSR-markers and 1195 MTAs (-log10(P-value) ≥ 3.0) for SNP-markers in all single environments plus the BLUEs. After Bonferroni correction, 28 MTAs remained significant for SSR-markers (-log10 (P-value) ≥ 4.82) and 58 MTAs for SNP-markers (-log10 (P-value) ≥ 5.89). Apart from chromosomes 4B and 6B for SSR-markers and chromosomes 4D and 5D for SNP-markers, MTAs were detected on all chromosomes. The highest number of significant SNP-markers was found on chromosomes 3B and 1B, while for the SSRs most markers were significant on chromosomes 6D and 3D. Overall, TGW was determined by many markers with small effects. Only three SNP-markers had R(2) values above 6%.

Large SNP arrays for genotyping in crop plants.

Genotyping with large numbers of molecular markers is now an indispensable tool within plant genetics and breeding. Especially through the identification of large numbers of single nucleotide polymorphism (SNP) markers using the novel high-throughput sequencing technologies, it is now possible to reliably identify many thousands of SNPs at many different loci in a given plant genome. For a number of important crop plants, SNP markers are now being used to design genotyping arrays containing thousands of markers spread over the entire genome and to analyse large numbers of samples. In this article, we discuss aspects that should be considered during the design of such large genotyping arrays and the analysis of individuals. The fact that crop plants are also often autopolyploid or allopolyploid is given due consideration. Furthermore, we outline some potential applications of large genotyping arrays including high-density genetic mapping, characterization (fingerprinting) of genetic material and breeding-related aspects such as association studies and genomic selection.

A large maize (Zea mays L.) SNP genotyping array: development and germplasm genotyping, and genetic mapping to compare with the B73 reference genome.

SNP genotyping arrays have been useful for many applications that require a large number of molecular markers such as high-density genetic mapping, genome-wide association studies (GWAS), and genomic selection. We report the establishment of a large maize SNP array and its use for diversity analysis and high density linkage mapping. The markers, taken from more than 800,000 SNPs, were selected to be preferentially located in genes and evenly distributed across the genome. The array was tested with a set of maize germplasm including North American and European inbred lines, parent/F1 combinations, and distantly related teosinte material. A total of 49,585 markers, including 33,417 within 17,520 different genes and 16,168 outside genes, were of good quality for genotyping, with an average failure rate of 4% and rates up to 8% in specific germplasm. To demonstrate this array's use in genetic mapping and for the independent validation of the B73 sequence assembly, two intermated maize recombinant inbred line populations - IBM (B73×Mo17) and LHRF (F2×F252) - were genotyped to establish two high density linkage maps with 20,913 and 14,524 markers respectively. 172 mapped markers were absent in the current B73 assembly and their placement can be used for future improvements of the B73 reference sequence. Colinearity of the genetic and physical maps was mostly conserved with some exceptions that suggest errors in the B73 assembly. Five major regions containing non-colinearities were identified on chromosomes 2, 3, 6, 7 and 9, and are supported by both independent genetic maps. Four additional non-colinear regions were found on the LHRF map only; they may be due to a lower density of IBM markers in those regions or to true structural rearrangements between lines. Given the array's high quality, it will be a valuable resource for maize genetics and many aspects of maize breeding.

Analysis of DNA polymorphisms in sugar beet (Beta vulgaris L.) and development of an SNP-based map of expressed genes.

A panel of 13 sugar beet lines and one genotype each of the Beta vulgaris cultivars red beet and Swiss chard, and B. vulgaris ssp. maritima were used to identify polymorphisms in alignments of genomic DNA sequences derived from 315 EST- and 43 non-coding RFLP-derived loci. In sugar beet lines, loci of expressed genes showed an average SNP frequency of 1/72 bp, 1 in 58 bp in non-coding sequences, increasing to 1/47 bp upon the addition of the remaining genotypes. Within analysed DNA fragments, alleles at different SNP positions displayed linkage disequilibrium indicative of haplotype structures. On average 2.7 haplotypes were found in sugar beet lines, and haplotype conservation in expressed genes appeared to exceed 500 bp in length. Seven different genotyping techniques including SNP detection by MALDI-TOF mass spectrometry, pyrosequencing and fluorescence scanning of labelled nucleotides were employed to perform 712 segregation analyses for 538 markers in three F(2) populations. Functions were predicted for 492 mapped sequences. Genetic maps comprised 305 loci covering 599.8 cM in population K1, 241 loci distributed over 636.6 cM in population D2, and 166 loci over 507.1 cM in population K2, respectively. Based on 156 markers common to more than one population an integrated map was constructed with 524 loci covering 664.3 cM. For 377 loci the genome positions of the most similar sequences from A. thaliana were identified, but little evidence for previously presented ancestral genome structures was found.

Serial analysis of gene expression identifies metallothionein-II as major neuroprotective gene in mouse focal cerebral ischemia.

We applied serial analysis of gene expression (SAGE) to study differentially expressed genes in mouse brain 14 hr after the induction of focal cerebral ischemia. Analysis of >60,000 transcripts revealed 83 upregulated and 94 downregulated transcripts (more than or equal to eightfold). Reproducibility was demonstrated by performing SAGE in duplicate on the same starting material. Metallothionein-II (MT-II) was the most significantly upregulated transcript in the ischemic hemisphere. MT-I and MT-II are assumed to be induced by metals, glucocorticoids, and inflammatory signals in a coordinated manner, yet their function remains elusive. Upregulation of both MT-I and MT-II was confirmed by Northern blotting. MT-I and MT-II mRNA expression increased as early as 2 hr after 2 hr of transient ischemia, with a maximum after 16 hr. Western blotting and immunohistochemistry revealed MT-I/-II upregulation in the ischemic hemisphere, whereas double labeling demonstrated the colocalization of MT with markers for astrocytes as well as for monocytes/macrophages. MT-I- and MT-II-deficient mice developed approximately threefold larger infarcts than wild-type mice and a significantly worse neurological outcome. For the first time we make available a comprehensive data set on brain ischemic gene expression and underscore the important protective role of metallothioneins in ischemic damage of the brain. Our results demonstrate the usefulness of SAGE to screen functionally relevant genes and the power of knock-out models in linking function to expression data generated by high throughput techniques.

Association of familial and sporadic rheumatoid arthritis with a single corticotropin-releasing hormone genomic region (8q12.3) haplotype.

OBJECTIVE: Rheumatoid arthritis (RA) is a common disabling autoimmune disease with a complex genetic component. We have previously described linkage of a region of chromosome 8q12.3 with RA and association of the microsatellite marker CRHRA1 with RA in 295 affected sibling-pair families. In the current study we aimed to physically link the RA-associated marker with the corticotropin-releasing hormone (CRH) candidate gene, and to examine the genomic region for additional short tandem repeat (STR) genetic markers in order to clarify the association with RA. METHODS: We examined the association of 2 STR markers with disease in the original 295 multicase families and in a cohort of 131 simplex families to refine our understanding of this genetic region in disease susceptibility in sporadic and familial RA. Genomic library screening and sequencing were used to generate physical sequences in the CRH genomic region. Bioinformatic analysis of the sequence flanking the CRH structural gene was used to screen for additional STRs and other genetic features. Genotyping was carried out using a standard fluorescence approach. Estimations of haplotype frequencies were performed to assess linkage disequilibrium. The transmission disequilibrium test was performed using TRANSMIT. RESULTS: Physical cloning and sequencing analyses identified the genomic region linking the CRHRA1 marker and the CRH structural locus. Moreover, we identified a further STR, CRHRA2, which was in strong linkage disequilibrium with CRHRA1 (P = 4.0 x 10(-14)). A haplotype, CRHRA1*10;CRHRA2*14, was preferentially carried by unaffected parents at a frequency of 8.6% compared with the expected frequency of 3.1%. This haplotype was overtransmitted in the multiply affected families (P = 0.0077) and, similarly, in the simplex families (P = 0.024). Combined analysis of both family cohorts confirmed significant evidence for linkage (P = 4.9 x 10(-4)) and association (P = 5.5 x 10(-3)) for this haplotype with RA. CONCLUSION: In demonstrating significant linkage disequilibrium between these 2 markers, we have refined the disease-associated region to a single haplotype and confirmed the significance of this region in our understanding of the genetics of RA.