High‐throughput genomics in sorghum: from whole‐genome resequencing to a SNP screening array

With its small, diploid and completely sequenced genome, sorghum (Sorghum bicolor L. Moench) is highly amenable to genomics‐based breeding approaches. Here, we describe the development and testing of a robust single‐nucleotide polymorphism (SNP) array platform that enables polymorphism screening for genome‐wide and trait‐linked polymorphisms in genetically diverse S. bicolor populations. Whole‐genome sequences with 6× to 12× coverage from five genetically diverse S. bicolor genotypes, including three sweet sorghums and two grain sorghums, were aligned to the sorghum reference genome. From over 1 million high‐quality SNPs, we selected 2124 Infinium Type II SNPs that were informative in all six source genomes, gave an optimal Assay Design Tool (ADT) score, had allele frequencies of 50% in the six genotypes and were evenly spaced throughout the S. bicolor genome. Furthermore, by phenotype‐based pool sequencing, we selected an additional 876 SNPs with a phenotypic association to early‐stage chilling tolerance, a key trait for European sorghum breeding. The 3000 attempted bead types were used to populate half of a dual‐species Illumina iSelect SNP array. The array was tested using 564 Sorghum spp. genotypes, including offspring from four unrelated recombinant inbred line (RIL) and F₂ populations and a genetic diversity collection. A high call rate of over 80% enabled validation of 2620 robust and polymorphic sorghum SNPs, underlining the efficiency of the array development scheme for whole‐genome SNP selection and screening, with diverse applications including genetic mapping, genome‐wide association studies and genomic selection.     

Genome‐wide association study for ham weight loss at first salting in Italian Large White pigs: towards the genetic dissection of a key trait for dry‐cured ham production

Protected designation of origin dry‐cured hams are the most important productions of the Italian heavy pig industry. Hams capable of minimal seasoning losses produce better quality dry‐cured hams. Ham weight loss during the first 7 days in brine (first salting) is highly correlated with the total loss of weight up to the end of seasoning, and it has quite high heritability (0.30–0.61). For these reasons, ham weight loss at first salting has been included as a meat quality trait in the Italian heavy pig selection program. In this work, we carried out a genome‐wide association study for this parameter in the Italian Large White pig breed by genotyping 1365 animals with the Illumina BeadChip PorcineSNP60 chip. A total of 44 single nucleotide polymorphisms (SNPs) had a P_{nominal value} below 5.0E‐04, five of which were below 5.0E‐05 and one of them (ALGA0057985 on chromosome 10) was associated with this trait at a P_Bonferroni threshold of 0.10. These SNPs identified a total of at least 29 putative QTLs that were located on most porcine autosomal chromosomes. This study provides genomic information that could be useful in dissecting this complex trait by identifying potential candidate genes whose function could contribute to understanding the biological mechanisms affecting meat quality for seasoning aptitude.     

The OsmiR396c‐OsGRF4‐OsGIF1 regulatory module determines grain size and yield in rice

Grain weight is the most important component of rice yield and is mainly determined by grain size, which is generally controlled by quantitative trait loci (QTLs). Although numerous QTLs that regulate grain weight have been identified, the genetic network that controls grain size remains unclear. Herein, we report the cloning and functional analysis of a dominant QTL, grain length and width 2 (GLW2), which positively regulates grain weight by simultaneously increasing grain length and width. The GLW2 locus encodes OsGRF4 (growth‐regulating factor 4) and is regulated by the microRNA miR396c in vivo. The mutation in OsGRF4 perturbs the OsmiR396 target regulation of OsGRF4, generating a larger grain size and enhanced grain yield. We also demonstrate that OsGIF1 (GRF‐interacting factors 1) directly interacts with OsGRF4, and increasing its expression improves grain size. Our results suggest that the miR396c‐OsGRF4‐OsGIF1 regulatory module plays an important role in grain size determination and holds implications for rice yield improvement.     

Confirmation and fine‐mapping of clinical mastitis and somatic cell score QTL in Nordic Holstein cattle

A genome‐wide association study of 2098 progeny‐tested Nordic Holstein bulls genotyped for 36 387 SNPs on 29 autosomes was conducted to confirm and fine‐map quantitative trait loci (QTL) for mastitis traits identified earlier using linkage analysis with sparse microsatellite markers in the same population. We used linear mixed model analysis where a polygenic genetic effect was fitted as a random effect and single SNPs were successively included as fixed effects in the model. We detected 143 SNP‐by‐trait significant associations (P < 0.0001) on 20 chromosomes affecting mastitis‐related traits. Among them, 21 SNP‐by‐trait combinations exceeded the genome‐wide significant threshold. For 12 chromosomes, both the present association study and the previous linkage study detected QTL, and of these, six were in the same chromosomal locations. Strong associations of SNPs with mastitis traits were observed on bovine autosomes 6, 13, 14 and 20. Possible candidate genes for these QTL were identified. Identification of SNPs in linkage disequilibrium with QTL will enable marker‐based selection for mastitis resistance. The candidate genes identified should be further studied to detect candidate polymorphisms underlying these QTL.     

Large‐scale GWAS in sorghum reveals common genetic control of grain size among cereals

Grain size is a key yield component of cereal crops and a major quality attribute. It is determined by a genotype’s genetic potential and its capacity to fill the grains. This study aims to dissect the genetic architecture of grain size in sorghum. An integrated genome‐wide association study (GWAS) was conducted using a diversity panel (n = 837) and a BC‐NAM population (n = 1421). To isolate genetic effects associated with genetic potential of grain size, rather than the genotype’s capacity to fill the grains, a treatment of removing half of the panicle was imposed during flowering. Extensive and highly heritable variation in grain size was observed in both populations in 5 field trials, and 81 grain size QTL were identified in subsequent GWAS. These QTL were enriched for orthologues of known grain size genes in rice and maize, and had significant overlap with SNPs associated with grain size in rice and maize, supporting common genetic control of this trait among cereals. Grain size genes with opposite effect on grain number were less likely to overlap with the grain size QTL from this study, indicating the treatment facilitated identification of genetic regions related to the genetic potential of grain size. These results enhance understanding of the genetic architecture of grain size in cereal, and pave the way for exploration of underlying molecular mechanisms and manipulation of this trait in breeding practices.     

A whole‐genome SNP array (RICE6K) for genomic breeding in rice

The advances in genotyping technology provide an opportunity to use genomic tools in crop breeding. As compared to field selections performed in conventional breeding programmes, genomics‐based genotype screen can potentially reduce number of breeding cycles and more precisely integrate target genes for particular traits into an ideal genetic background. We developed a whole‐genome single nucleotide polymorphism (SNP) array, RICE6K, based on Infinium technology, using representative SNPs selected from more than four million SNPs identified from resequencing data of more than 500 rice landraces. RICE6K contains 5102 SNP and insertion–deletion (InDel) markers, about 4500 of which were of high quality in the tested rice lines producing highly repeatable results. Forty‐five functional markers that are located inside 28 characterized genes of important traits can be detected using RICE6K. The SNP markers are evenly distributed on the 12 chromosomes of rice with the average density of 12 SNPs per 1 Mb and can provide information for polymorphisms between indica and japonica subspecies as well as varieties within indica and japonica groups. Application tests of RICE6K showed that the array is suitable for rice germplasm fingerprinting, genotyping bulked segregating pools, seed authenticity check and genetic background selection. These results suggest that RICE6K provides an efficient and reliable genotyping tool for rice genomic breeding.     

Effects of hydroxy‐delta‐5‐steroid dehydrogenase, 3 beta‐ and steroid delta‐isomerase 1 polymorphisms on fat androstenone level and gene expression in Duroc pigs

A high level of androstenone in porcine adipose tissue is a major factor contributing to boar taint. Porcine hydroxy‐delta‐5‐steroid dehydrogenase, 3 beta‐ and steroid delta‐isomerase 1 (3β‐HSD, also known as HSD3B1) plays a key role in the hepatic metabolism that catalyzes androstenone to β‐androstenol. Therefore, 3β‐HSD is a candidate gene for boar taint. This study aimed to investigate functional 3β‐HSD polymorphisms in Duroc pigs. We found eight single nucleotide polymorphisms (SNPs) in the full‐length porcine 3β‐HSD. Four of the SNPs had restriction enzyme sites, and we genotyped them in 147 uncastrated male Duroc pigs using a polymerase chain reaction–restriction fragment length polymorphism method. Pigs with the GG genotype at the g.165262G>A locus (SNP5) had significantly lower androstenone levels than did those with other genotypes (P = 0.030). SNP5 also was associated with differences in 3β‐HSD mRNA levels: pigs with the GG genotype had higher levels than those with other genotypes (P = 0.019). The SNP5 polymorphism could affect the hepatic catabolism of androstenone and consequently impact androstenone accumulation in the adipose tissue. Therefore, SNP5 in the 3β‐HSD of Duroc pigs could be a useful selective marker for decreasing boar taint.     

LIN28A gene polymorphisms modify neuroblastoma susceptibility: A four‐centre case‐control study

Neuroblastoma ranks the most common seen solid tumour in childhood. Overexpression of LIN28A gene has been linked to the development of multiple human malignancies, but the relationship between LIN28A single nucleotide polymorphisms (SNPs) and neuroblastoma susceptibility is still under debate. Herein, we evaluated the correlation of four potentially functional LIN28A SNPs (rs3811464 G>A, rs3811463 T>C, rs34787247 G>A, and rs11247957 G>A) and neuroblastoma susceptibility in 505 neuroblastoma patients and 1070 controls from four independent hospitals in China. The correlation strengths were determined by using odds ratios (ORs) and corresponding 95% confidence intervals (CIs). Among these SNPs, rs34787247 G>A exhibited a significant association with increased susceptibility in neuroblastoma (GA vs GG: adjusted OR = 1.30, 95% CI = 1.03‐1.64; AA vs GG: adjusted OR = 2.51, 95% CI = 1.36‐4.64, AA/GA vs GG: adjusted OR = 1.42, 95% CI = 1.12‐1.80, AA vs GG/GA: adjusted OR = 2.39, 95% CI = 1.29‐4.42). Furthermore, the combined analysis of risk genotypes revealed that subjects carrying three risk genotypes (adjusted OR = 1.64, 95% CI = 1.02‐2.63) are more inclined to develop neuroblastoma than those without risk genotype, and so do carriers of 1‐4 risk genotypes (adjusted OR = 1.26, 95% CI = 1.01‐1.56). Stratification analysis further revealed risk effect of rs3811464 G>A, rs34787247 G>A and 1‐4 risk genotypes in some subgroups. Haplotype analysis of these four SNPs yields two haplotypes significantly correlated with increased neuroblastoma susceptibility. Overall, our finding indicated that LIN28A SNPs, especially rs34787247 G>A, may increase neuroblastoma risk.     

High‐density genetic map of Miscanthus sinensis reveals inheritance of zebra stripe

Miscanthus is a perennial C₄ grass that has recently become an important bioenergy crop. The efficiency of breeding improved Miscanthus biomass cultivars could be greatly increased by marker‐assisted selection. Thus, a high‐density genetic map is critical to Miscanthus improvement. In this study, a mapping population of 261 F₁ progeny was developed from a cross between two diploid M. sinensis cultivars, ‘Strictus’ and ‘Kaskade’. High‐density genetic maps for the two parents were produced with 3044 newly developed single nucleotide polymorphisms (SNPs) obtained from restriction site‐associated DNA sequencing, and 138 previously mapped GoldenGate SNPs. The female parent (‘Strictus’) map spanned 1599 cM, with 1989 SNPs on 19 linkage groups, and an average intermarker spacing of 0.8 cM. The length of the male parent (‘Kaskade’) map was 1612 cM, with 1821 SNPs, and an average intermarker spacing of 0.9 cM. The utility of the map was confirmed by locating quantitative trait loci (QTL) for the zebra‐striped trait, which was segregating in this population. Three QTL for zebra‐striped presence/absence (zb1, zb2 on LG 7, and zb3 on LG 10) and three for zebra‐striped intensity (zbi1, zbi2, zbi3 on LGs 7, 10, 3) were identified. Each allele that caused striping was recessive. Incomplete penetrance was observed for each zb QTL, but penetrance was greatest when two or more zb QTL were homozygous for the causative alleles. Similarly, the intensity of striping was greatest when two or more zbi QTL were homozygous for alleles that conferred the trait. Comparative mapping indicated putative correspondence between zb3 and/or zbi2 on LG 10 to previously sequenced genes conferring zebra stripe in maize and rice. These results demonstrate that the new map is useful for identifying marker–trait associations. The mapped markers will become a valuable community resource, facilitating comparisons among studies and the breeding of Miscanthus.     

Genetic dissection of pre‐anthesis sub‐phase durations during the reproductive spike development of wheat

Flowering time is an important factor affecting grain yield in wheat. In this study, we divided reproductive spike development into eight sub‐phases. These sub‐phases have the potential to be delicately manipulated to increase grain yield. We measured 36 traits with regard to sub‐phase durations, determined three grain yield‐related traits in eight field environments and mapped 15 696 single nucleotide polymorphism (SNP, based on 90k Infinium chip and 35k Affymetrix chip) markers in 210 wheat genotypes. Phenotypic and genetic associations between grain yield traits and sub‐phase durations showed significant consistency (Mantel test; r = 0.5377, P < 0.001). The shared quantitative trait loci (QTLs) revealed by the genome‐wide association study suggested a close association between grain yield and sub‐phase duration, which may be attributed to effects on spikelet initiation/spikelet number (double ridge to terminal spikelet stage, DR‐TS) and assimilate accumulation (green anther to anthesis stage, GA‐AN). Moreover, we observed that the photoperiod‐sensitivity allele at the Ppd‐D1 locus on chromosome 2D markedly extended all sub‐phase durations, which may contribute to its positive effects on grain yield traits. The dwarfing allele at the Rht‐D1 (chromosome 4D) locus altered the sub‐phase duration and displayed positive effects on grain yield traits. Data for 30 selected genotypes (from among the original 210 genotypes) in the field displayed a close association with that from the greenhouse. Most importantly, this study demonstrated specific connections to grain yield in narrower time windows (i.e. the eight sub‐phases), rather than the entire stem elongation phase as a whole.     

Single nucleotide polymorphisms in the insulin‐like growth factor 1 (IGF1) gene are associated with growth‐related traits in farmed Atlantic salmon

Understanding the genetic basis of variation in traits related to growth and fillet quality in Atlantic salmon is of importance to the aquaculture industry. Several growth‐related QTL have been identified via the application of genetic markers. The IGF1 gene is considered a highly conserved and crucial growth‐regulating gene in salmonid species. However, the association between polymorphisms in the IGF1 gene and growth‐related traits in Atlantic salmon is unknown. Therefore, in this study, regions of the Atlantic salmon IGF1 gene were sequenced, aligned and compared across individuals. Three SNPs were identified in the putative promoter (SNP1, g.5763G>T; GenBank no. AGKD01012745 ), intron 1 (SNP2, g.7292C>T; GenBank no. AGKD01012745 ) and intron 3 (SNP3, g.4671A>C; GenBank no. AGKD01133398 ) regions respectively. These SNPs were genotyped in a population of 4800 commercial Atlantic salmon with data on several weight and fillet traits measured at harvest (at approximately 3 years of age). In a mixed model, association analysis of individual SNPs, SNP1 and SNP3 were both significantly associated with several weight traits (P < 0.05). The estimated additive effect on overall harvest weight was approximately 35 and 110 g for SNPs 1 and 3 respectively. A haplotype analysis confirmed the association between genetic variation in the IGF1 gene with overall body weight (P < 0.05) and fillet component traits (P < 0.05). Our findings suggest the identified nucleotide polymorphisms of the IGF1 gene may either affect farmed Atlantic salmon growth directly or be in population‐wide linkage disequilibrium with causal variation, highlighting their possible utility as candidates for marker‐assisted selection in the aquaculture industry.     

Genome‐wide association study for jumping performances in French sport horses

A genome‐wide association study was performed to identify single nucleotide polymorphisms (SNPs) associated with jumping performances of warmbloods in France. The 999 horses included in the study for jumping performances were sport horses [mostly Selle Français (68%), Anglo‐Arabians (13%) and horses from the other European studbooks]. Horses were genotyped using the Illumina EquineSNP50 BeadChip. Of the 54 602 SNPs available on this chip, 44 424 were retained after quality testing. Phenotypes were obtained by deregressing official breeding values for jumping competitions to use all available information, that is, the performances of each horse as well as those of its relatives. Two models were used to test the effects of the genotypes on deregressed phenotypes: a single‐marker mixed model and a haplotype‐based mixed model (significant: P < 1E‐05; suggestive: P < 1E‐04). Both models included a polygenic effect to take into account familial structures. For jumping performances, one suggestive quantitative trait locus (QTL) located on chromosome 1 (BIEC2_31196 and BIEC2_31198) was detected with both models. This QTL explains 0.7% of the phenotypic variance. RYR2, a gene encoding a major calcium channel in cardiac muscle in humans and mice, is located 0.55 Mb from this potential QTL.     

Genome‐wide association study to identify genomic regions affecting prolificacy in Lori‐Bakhtiari sheep

Several causative mutations in candidate genes affecting prolificacy have been detected in various sheep breeds. A genome‐wide association study was performed on estimated breeding values for litter size in Lori‐Bakhtiari sheep. Prolific ewes with twinning records and others with only singleton records were genotyped using the medium‐density Illumina Ovine SNP50 array. Four single nucleotide polymorphisms (SNPs) associated with litter size were identified on chromosomes 3, 6 and 22. The region on sheep chromosome 3 between 75 739 167 and 75 745 152 bp included two significant SNPs (s52383.1 and OAR3_80038014_X.1) in high linkage disequilibrium with each other. The region that surrounds these SNPs contains a novel putative candidate gene: luteinizing hormone/choriogonadotropin receptor (LHCGR), known to be involved in ovarian steroidogenesis and organism‐specific biosystem pathways in sheep. Known prolificacy genes BMPR1B, BMP15 and GDF9 were not associated with litter size in Lori‐Bakhtiari sheep, suggesting that other biological mechanisms could be responsible for the trait's variation in this breed.     

Analysis of the genetic architecture of maize kernel size traits by combined linkage and association mapping

Kernel size‐related traits are the most direct traits correlating with grain yield. The genetic basis of three kernel traits of maize, kernel length (KL), kernel width (KW) and kernel thickness (KT), was investigated in an association panel and a biparental population. A total of 21 single nucleotide polymorphisms (SNPs) were detected to be most significantly (P < 2.25 × 10^?6) associated with these three traits in the association panel under four environments. Furthermore, 50 quantitative trait loci (QTL) controlling these traits were detected in seven environments in the intermated B73 × Mo17 (IBM) Syn10 doubled haploid (DH) population, of which eight were repetitively identified in at least three environments. Combining the two mapping populations revealed that 56 SNPs (P < 1 × 10^?3) fell within 18 of the QTL confidence intervals. According to the top significant SNPs, stable‐effect SNPs and the co‐localized SNPs by association analysis and linkage mapping, a total of 73 candidate genes were identified, regulating seed development. Additionally, seven miRNAs were found to situate within the linkage disequilibrium (LD) regions of the co‐localized SNPs, of which zma‐miR164e was demonstrated to cleave the mRNAs of Arabidopsis CUC1, CUC2 and NAC6 in vitro. Overexpression of zma‐miR164e resulted in the down‐regulation of these genes above and the failure of seed formation in Arabidopsis pods, with the increased branch number. These findings provide insights into the mechanism of seed development and the improvement of molecular marker‐assisted selection (MAS) for high‐yield breeding in maize.     

Unlocking the secondary gene‐pool of barley with next‐generation sequencing

Crop wild relatives (CWR) provide an important source of allelic diversity for any given crop plant species for counteracting the erosion of genetic diversity caused by domestication and elite breeding bottlenecks. Hordeum bulbosum L. is representing the secondary gene pool of the genus Hordeum. It has been used as a source of genetic introgressions for improving elite barley germplasm (Hordeum vulgare L.). However, genetic introgressions from H. bulbosum have yet not been broadly applied, due to a lack of suitable molecular tools for locating, characterizing, and decreasing by recombination and marker‐assisted backcrossing the size of introgressed segments. We applied next‐generation sequencing (NGS) based strategies for unlocking genetic diversity of three diploid introgression lines of cultivated barley containing chromosomal segments of its close relative H. bulbosum. Firstly, exome capture‐based (re)‐sequencing revealed large numbers of single nucleotide polymorphisms (SNPs) enabling the precise allocation of H. bulbosum introgressions. This SNP resource was further exploited by designing a custom multiplex SNP genotyping assay. Secondly, two‐enzyme‐based genotyping‐by‐sequencing (GBS) was employed to allocate the introgressed H. bulbosum segments and to genotype a mapping population. Both methods provided fast and reliable detection and mapping of the introgressed segments and enabled the identification of recombinant plants. Thus, the utilization of H. bulbosum as a resource of natural genetic diversity in barley crop improvement will be greatly facilitated by these tools in the future.     

Gene editing of the wheat homologs of TONNEAU1‐recruiting motif encoding gene affects grain shape and weight in wheat

Grain size and weight are important components of a suite of yield‐related traits in crops. Here, we showed that the CRISPR‐Cas9 gene editing of TaGW7, a homolog of rice OsGW7 encoding a TONNEAU1‐recruiting motif (TRM) protein, affects grain shape and weight in allohexaploid wheat. By editing the TaGW7 homoeologs in the B and D genomes, we showed that mutations in either of the two or both genomes increased the grain width and weight but reduced the grain length. The effect sizes of mutations in the TaGW7 gene homoeologs coincided with the relative levels of their expression in the B and D genomes. The effects of gene editing on grain morphology and weight traits were dosage dependent with the double‐copy mutant showing larger effect than the respective single copy mutants. The TaGW7‐centered gene co‐expression network indicated that this gene is involved in the pathways regulating cell division and organ growth, also confirmed by the cellular co‐localization of TaGW7 with α‐ and β‐tubulin proteins, the building blocks of microtubule arrays. The analyses of exome capture data in tetraploid domesticated and wild emmer, and hexaploid wheat revealed the loss of diversity around TaGW7‐associated with domestication selection, suggesting that TaGW7 is likely to play an important role in the evolution of yield component traits in wheat. Our study showed how integrating CRISPR‐Cas9 system with cross‐species comparison can help to uncover the function of a gene fixed in wheat for allelic variants targeted by domestication selection and select targets for engineering new gene variants for crop improvement.