Genome‐wide association mapping for phenotypic plasticity in rice

Phenotypic plasticity of plants in response to environmental changes is important for adapting to changing climate. Less attention has been paid to exploring the advantages of phenotypic plasticity in resource‐rich environments to enhance the productivity of agricultural crops. Here, we examined genetic variation for phenotypic plasticity in indica rice (Oryza sativa L.) across two diverse panels: (1) a Phenomics of Rice Adaptation and Yield (PRAY) population comprising 301 accessions; and (2) a Multi‐parent Advanced Generation Inter‐Cross (MAGIC) indica population comprising 151 accessions. Altered planting density was used as a proxy for elevated atmospheric CO₂ response. Low planting density significantly increased panicle weight per plant compared with normal density, and the magnitude of the increase ranged from 1.10 to 2.78 times among accessions for the PRAY population and from 1.05 to 2.45 times for the MAGIC population. Genome‐wide‐association studies validate three E nvironmental R esponsiveness (ER) candidate alleles (qER1–3) that were associated with relative response of panicle weight to low density. Two of these alleles were tested in 13 genotypes to clarify their biomass responses during vegetative growth under elevated CO₂ in Japan. Our study provides evidence for polymorphisms that control rice phenotypic plasticity in environments that are rich in resources such as light and CO₂.     

Genome‐wide association analyses reveal the genetic basis of combining ability in rice

Combining ability is a measure for selecting elite parents and predicting hybrid performance in plant breeding. However, the genetic basis of combining ability remains unclear and a global view of combining ability from diverse mating designs is lacking. We developed a North Carolina II (NCII) population of 96 Oryza sativa and four male sterile lines to identify parents of greatest value for hybrid rice production. Statistical analyses indicated that general combining ability (GCA) and specific combining ability (SCA) contributed variously to different agronomic traits. In a genome‐wide association study (GWAS) of agronomic traits, GCA and SCA, we identified 34 significant associations (P < 2.39 × 10^?7). The superior alleles of GCA loci (Ghd8, GS3 and qSSR4) accumulated in parental lines with high GCA and explained 30.03% of GCA variance in grain yield, indicating that molecular breeding of high GCA parental lines is feasible. The distinct distributions of these QTLs contributed to the differentiation of parental GCA in subpopulations. GWAS of SCA identified 12 more loci that showed dominance on corresponding agronomic traits. We conclude that the accumulation of superior GCA and SCA alleles is an important contributor to heterosis and QTLs that greatly contributed to combining ability in our study would accelerate the identification of elite inbred lines and breeding of super hybrids.     

Genome‐wide association mapping of heritable temperament variation in the Tennessee Walking Horse

Temperament is a key criterion in the selection of horses for both leisure and competitive riding to ensure optimal performance and safety. The Tennessee Walking Horse (TWH) is described as a calm, docile breed and is often used as a trail, show and pleasure horse. However, among horse owners and caretakers, there are anecdotes supporting familial and disciplinal typical behaviors and personalities. To investigate the contribution of genetics to temperament, we collected a behavior questionnaire, brief training history and identifying information for 276 TWH, as well as blood or hair samples for DNA. Factor analysis was conducted on the 20‐item questionnaire for the set of 216 horses that met inclusion thresholds. Factor analysis identified four temperament factors in TWH: ‘anxious’, ‘tractable’, ‘agonistic’ and ‘gregarious’. These four factors account for 64% of the total trait variance. DNA from 113 TWHs were selected and genotyped using the Equine SNP70 bead chip for three separate genome‐wide association studies (GWAs) using the factor 1‐anxious, factor 2‐tractable and factor 3‐agonistic scores as the phenotype. Quantitative association analysis identified significant candidate loci for each factor that warrant further investigation.     

Genome‐wide association study identifies an NLR gene that confers partial resistance to Magnaporthe oryzae in rice

Because of the frequent breakdown of major resistance (R) genes, identification of new partial R genes against rice blast disease is an important goal of rice breeding. In this study, we used a core collection of the Rice Diversity Panel II (C‐RDP‐II), which contains 584 rice accessions and are genotyped with 700 000 single‐nucleotide polymorphism (SNP) markers. The C‐RDP‐II accessions were inoculated with three blast strains collected from different rice‐growing regions in China. Genome‐wide association study identified 27 loci associated with rice blast resistance (LABRs). Among them, 22 LABRs were not associated with any known blast R genes or QTLs. Interestingly, a nucleotide‐binding site leucine‐rich repeat (NLR) gene cluster exists in the LABR12 region on chromosome 4. One of the NLR genes is highly conserved in multiple partially resistant rice cultivars, and its expression is significantly up‐regulated at the early stages of rice blast infection. Knockout of this gene via CRISPR‐Cas9 in transgenic plants partially reduced blast resistance to four blast strains. The identification of this new non‐strain specific partial R gene, tentatively named rice blast Partial Resistance gene 1 (PiPR1), provides genetic material that will be useful for understanding the partial resistance mechanism and for breeding durably resistant cultivars against blast disease of rice.     

Genome‐wide association study identifies variation of glucosidase being linked to natural variation of the maximal quantum yield of photosystem II

The maximum quantum yield of photosystem II (as reflected by variable to maximum chlorophyll a fluorescence, F_v/F_m) is regarded as one of the most important photosynthetic parameters. The genetic basis underlying natural variation in F_v/F_m, which shows low level of variations in plants under non‐stress conditions, is not easy to be exploited using the conventional gene cloning approaches. Thus, in order to answer this question, we have followed another strategy: we used genome‐wide association study (GWAS) and transgenic analysis in a rice mini‐core collection. We report here that four single‐nucleotide polymorphisms, located in the promoter region of β‐glucosidase 5 (BGlu‐5), are associated with observed variation in F_v/F_m. Indeed, our transgenic analysis showed a good correlation between BGlu‐5 and F_v/F_m. Thus, our work demonstrates the feasibility of using GWAS to study natural variation in F_v/F_m, suggesting that cis‐element polymorphism, affecting the BGlu‐5 expression level, may, indirectly, contribute to F_v/F_m variation in rice through the gibberellin signaling pathway. Further research is needed to understand the mechanism of our novel observation.     

Genome‐wide association mapping of natural variation in odour‐guided behaviour in Drosophila

A defining goal in the field of behavioural genetics is to identify the key genes or genetic networks that shape behaviour. A corollary to this goal is the goal of identifying genetic variants that are responsible for variation in the behaviour. These goals are achieved by measuring behavioural responses to controlled stimuli, in the present case the responses of Drosophila melanogaster to olfactory stimuli. We used a high‐throughput behavioural assay system to test a panel of 157 Drosophila inbred lines derived from a natural population for both temporal and spatial dynamics of odour‐guided behaviour. We observed significant variation in response to the odourant 2,3‐butanedione, a volatile compound present in fermenting fruit. The recent whole genome sequencing of these inbred lines allowed us to then perform genome‐wide association analyses in order to identify genetic polymorphisms underlying variation in responses. These analyses revealed numerous single nucleotide polymorphisms associated with variation in responses. Among the candidate genes identified were both novel and previously identified olfaction‐related genes. Further, gene network analyses suggest that genes influencing variation in odour‐guided behaviour are enriched for functions involving neural processing and that these genes form a pleiotropic interaction network. We examined several of these candidate genes that were highly connected in the protein‐ and genetic interaction networks using RNA interference. Our results showed that subtle changes influencing nervous system function can result in marked differences in behaviour .     

Genome‐wide association study of body morphological traits in Sudanese goats

Long‐term selection of goats for a certain production system and/or different environmental conditions will be reflected in the body morphology of the animals under selection. To investigate the variation contributing to different morphological traits and to identify genomic regions that are associated with body morphological traits in Sudanese goats, we genotyped 96 females belonging to four Sudanese goat breeds with the SNP52 BeadChip. After quality control of the data, the genome‐wide association study was performed using 95 goats and 24 027 informative single nucleotide polymorphisms (SNPs). Bicoastal diameter was significantly associated (LOD = 6.32) with snp10185‐scaffold1365‐620922 on chromosome 2. The minor allele has an additive effect, increasing the bicoastal diameter by 2.6 cm. A second significant association was found between body length and snp56482‐scaffold89‐467312 on chromosome 3 (LOD = 5.65). The minor allele is associated with increased body length. Additionally, five regions were suggestive for cannon bone, head width, rump length and withers height (LOD > 5). Only one gene (CNTNAP5) is located within the 1‐Mb region surrounding the significant SNP for bicoastal diameter on chromosome 2. The body length QTL on chromosome 3 harbors 49 genes. Further research is required to validate the observed associations and to prioritize candidate genes.     

Genome‐wide association study uncovers four QTL predisposing to supernumerary teats in cattle

Supernumerary teats (hyperthelia, SNTs) are a common abnormality of the bovine udder with a medium to high heritability and a postulated oligogenic or polygenic inheritance pattern. SNTs not only negatively affect machine milking ability but also act as a reservoir for bacteria. A genome‐wide association study was carried out to identify genes involved in the development of SNTs in the dual‐purpose Fleckvieh breed. A total of 2467 progeny‐tested bulls were genotyped at 43 698 single nucleotide polymorphisms, and daughter yield deviations (DYDs) for ‘udder clearness' (UC) were used as high‐heritability phenotypes. Massive structuring of the study population was accounted for by principal components analysis‐based and mixed model‐based approaches. Four loci on BTA5, BTA6, BTA11 and BTA17 were significantly associated with the UC DYD. Three associated regions contain genes of the highly conserved Wnt signalling pathway. The four QTL together account for 10.7% of the variance of the UC DYD, whereas the major fraction of the DYD variance is attributable to chromosomes with no identified QTL. Our results support both an oligogenic and a polygenic inheritance pattern of SNTs in cattle. The identified candidate genes permit insights into the genetic architecture of teat malformations in cattle and provide clues to unravel the molecular mechanisms of mammary gland alterations in cattle and other species.     

Genome‐wide association study for pigmentation traits in Chinese Holstein population

With the Illumina BovineSNP50K BeadChip, we performed a genome‐wide association study (GWAS) for two pigmentation traits in a Chinese Holstein population: proportion of black (PB) and teat colour (TC). A case–control design was used. Cases were the cows with PB <0.30 (n = 129) and TC <2 points (n = 140); controls were those with PB >0.90 (n = 58) and TC >4 points (n = 281). The RM test of roadtrips (version 1.2) was applied to detect SNPs for the two traits with 42 883 and 42 741 SNPs respectively. A total of nine and 12 genome‐wide significant (P < 0.05) SNPs associated with PB and TC respectively were identified. Of these, two SNPs for PB were located within the KIT and IGFBP7 genes, and the other four SNPs were 23~212 kb away from the PDGFRA gene on BTA6; nine SNPs associated with TC were located within or 21~78.8 kb away from known genes on chromosomes 4, 11, 22, 23 and 24. By combing through our GWAS results and the biological functions of the genes, we suggest that the KIT, IGFBP7, PDGFRA, MITF, ING3 and WNT16 genes are promising candidates for PB and TC in Holstein cattle, providing a basis for further investigation on the genetic mechanism of pigmentation formation.     

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.     

Genome‐wide association study for semen traits of the bulls in Chinese Holstein

A genome‐wide association study (GWAS) was performed to identify markers and candidate genes for five semen traits in the Holstein bull population in China. The analyzed dataset consisted of records from 692 bulls from eight bull stations; each bull was genotyped using the Illumina BovineSNP50 BeadChip. Association tests between each trait and the 41 188 informative high‐quality SNPs were achieved with gapit software. In total, 19 suggestive significant SNPs, partly located within the reported QTL regions or within or close to the reported candidate genes, associated with five semen traits were detected. By combining our GWAS results with the biological functions of these genes, eight novel promising candidate genes, including ETNK1, PDE3A, PDGFRB, CSF1R, WT1, DSCAML1, SOD1 and RUNX2, were identified that potentially relate to semen traits. Our findings may provide a basis for further research on the genetic mechanism of semen traits and marker‐assisted selection of such traits in Holstein bulls.     

Genome‐wide association with delayed puberty in swine

An improvement in the proportion of gilts entering the herd that farrow a litter would increase overall herd performance and profitability. A significant proportion (10–30%) of gilts that enter the herd never farrow a litter; reproductive reasons account for approximately a third of gilt removals, with anestrous and failure to conceive the most common reasons for culling. Tools to select gilts for reproductive longevity through genomics or alternative phenotypes would be of great benefit to the producer. Ninety‐one gilts that failed to display behavioral estrus by 240 days (cases) and 127 pubertal littermates (controls) were genotyped with the Illumina Porcine SNP60 Beadchip. One hundred and seventy‐four SNPs with the most significant associations were genotyped in an additional 86 cases and 103 controls. Twelve of these associations were significant in the final analysis. The most significant (P < 1.5 × 10^?14) region associated with failure to attain puberty was on chromosome 4 surrounding the NHLH2 gene. Delayed pubertal development and age at first estrus have been associated with NHLH2 in mice. Because attainment of puberty is a complex trait, identifying genes that affect pubertal age would greatly contribute to our knowledge of reproductive development as well as overall fertility.     

Genome‐wide association study identified genes in the response to Salmonella pullorum infection in chickens

Pullorum is a bacterial disease that threatens the modern poultry industry. Over the years, research on this topic has focused mainly on its epidemiology, whereas the hosts’ genetic basis of infection is still vague. In order to identify chickens’ genes associated with pullorum, we sequenced 300 New Pudong chicken by double digest genotyping‐by‐sequencing. We obtained 1 527 953 SNPs for a genome‐wide association analysis, which identified 43 genome‐wide significant markers. Most of the significant SNPs were in the interval of 57.7–59.0 Mb on chromosome 5. The gene set enrichment analysis suggests a potential manner for bacterial infection and remaining inside the host. This work provides basic data for the purification, prevention and treatment of pullorum disease.     

Genome‐wide association study for body weight in cattle populations from Siberia

Body weight is a complex trait in cattle associated with commonly used commercial breeding measurements related to growth. Although many quantitative trait loci (QTL) for body weight have been identified in cattle so far, searching for genetic determinants in different breeds or environments is promising. Therefore, we carried out a genome‐wide association study (GWAS) in two cattle populations from the Russian Federation (Siberian region) using the GGP HD150K array containing 139 376 single nucleotide polymorphism (SNP) markers. Association tests for 107 550 SNPs left after filtering revealed five statistically significant SNPs on BTA5, considering a false discovery rate of less than 0.05. The chromosomal region containing these five SNPs contains the CCND2 gene, which was previously associated with average daily weight gain and body mass index in US beef cattle populations and in humans respectively. Our study is the first GWAS for body weight in beef cattle populations from the Russian Federation. The results provided here suggest that, despite the existence of breed‐ and species‐specific QTL, the genetic architecture of body weight could be evolutionarily conserved in mammals.     

Genome‐wide association reveals the locus responsible for four‐horned ruminant

Phenotypic variability in horn characteristics, such as their size, number and shape, offers the opportunity to elucidate the molecular basis of horn development. The objective of this study was to map the genetic determinant controlling the production of four horns in two breeds, Jacob sheep and Navajo‐Churro, and examine whether an eyelid abnormality occurring in the same populations is related. Genome‐wide association mapping was performed using 125 animals from the two breeds that contain two‐ and four‐horned individuals. A case–control design analysis of 570 712 SNPs genotyped with the ovine HD SNP Beadchip revealed a strong association signal on sheep chromosome 2. The 10 most strongly associated SNPs were all located in a region spanning Mb positions 131.9–132.6, indicating the genetic architecture underpinning the production of four horns is likely to involve a single gene. The closest genes to the most strongly associated marker (OAR2_132568092) were MTX2 and the HOXD cluster, located approximately 93 Kb and 251 Kb upstream respectively. The occurrence of an eyelid malformation across both breeds was restricted to polled animals and those carrying more than two horns. This suggests the eyelid abnormality may be associated with departures from the normal developmental production of two‐horned animals and that the two conditions are developmentally linked. This study demonstrated the presence of separate loci responsible for the polled and four‐horned phenotypes in sheep and advanced our understanding of the complexity that underpins horn morphology in ruminants.     

Genome‐wide association studies for hematological traits in swine

Improving immune capacity may increase the profitability of animal production if it enables animals to better cope with infections. Hematological traits play pivotal roles in animal immune capacity and disease resistance. Thus far, few studies have been conducted using a high‐density swine SNP chip panel to unravel the genetic mechanism of the immune capability in domestic animals. In this study, using mixed model‐based single‐locus regression analyses, we carried out genome‐wide association studies, using the Porcine SNP60 BeadChip, for immune responses in piglets for 18 hematological traits (seven leukocyte traits, seven erythrocyte traits, and four platelet traits) after being immunized with classical swine fever vaccine. After adjusting for multiple testing based on permutations, 10, 24, and 77 chromosome‐wise significant SNPs were identified for the leukocyte traits, erythrocyte traits, and platelet traits respectively, of which 10 reached genome‐wise significance level. Among the 53 SNPs for mean platelet volume, 29 are located in a linkage disequilibrium block between 32.77 and 40.59 Mb on SSC6. Four genes of interest are located within the block, providing genetic evidence that this genomic segment may be considered a candidate region relevant to the platelet traits. Other candidate genes of interest for red blood cell, hemoglobin, and red blood cell volume distribution width also have been found near the significant SNPs. Our genome‐wide association study provides a list of significant SNPs and candidate genes that offer valuable information for future dissection of molecular mechanisms regulating hematological traits.