Genomewide rapid association using mixed model and regression: a fast and simple method for genomewide pedigree-based quantitative trait loci association analysis

Role of inflammation in depressive and anxiety disorders, affect, and cognition: genetic and non-genetic findings in the lifelines cohort study

Inflammation is associated with a range of neuropsychiatric symptoms, but the issue of causality remains unclear. We used complementary non-genetic, genetic risk score (GRS), and Mendelian randomization (MR) analyses to examine whether inflammatory markers are associated with affect, depressive and anxiety disorders, and cognition. We tested in ≈55,098 (59% female) individuals from the Dutch Lifelines cohort the concurrent/prospective associations of C-reactive protein (CRP) with: depressive and anxiety disorders; positive/negative affect; and attention, psychomotor speed, episodic memory, and executive functioning at baseline and a follow-up assessment occurring 3.91 years later (SD = 1.21). Additionally, we examined the association between inflammatory GRSs (CRP, interleukin-6 [IL-6], IL-6 receptor [IL-6R and soluble IL-6R (sIL-6R)], glycoprotein acetyls [GlycA]) on these same outcomes (Nmin = 35,300; Nmax = 57,946), followed by MR analysis examining evidence of causality of CRP on outcomes (Nmin=22,154; Nmax = 23,268). In non-genetic analyses, higher CRP was associated with depressive disorder, lower positive/higher negative affect, and worse executive function, attention, and psychomotor speed after adjusting for potential confounders. In genetic analyses, CRPGRS was associated with any anxiety disorder (β = 0.002, p = 0.037) whereas GlycAGRS was associated with major depressive disorder (β = 0.001, p = 0.036). Both CRPGRS (β = 0.006, p = 0.035) and GlycAGRS (β = 0.006, p = 0.049) were associated with greater negative affect. Inflammatory GRSs were not associated with cognition, except sIL-6RGRS which was associated with poorer memory (β = -0.009, p = 0.018). There was a non-significant CRP-anxiety association using MR (β = 0.12; p = 0.054). Genetic and non-genetic analyses provide consistent evidence for an association between CRP and negative affect. These results suggest that inflammation may impact a broad range of trans-diagnostic affective symptoms.

Genetic determinism of cortisol levels in pig

In facing the challenge of sustainability, animal breeding provides the option to improve animal robustness. In the search for new selection criteria related to robustness, the hypothalamic-pituitary-adrenocortical (HPA) axis is studied as a major neuroendocrine system involved in metabolic regulations and adaptive responses. Indeed, HPA axis activity is strongly influenced by genetic factors acting at several levels of the axis. The adrenocorticotropic hormone (ACTH) stimulation test has long been used to analyze interindividual and genetic differences in HPA axis activity in several species, including pigs. To uncover the genetic determinism of HPA activity and its influence on functional traits and robustness, a divergent selection experiment was carried out for three generations in a Large White pig population based on plasma cortisol levels measured one hour after injection of ACTH. In the present study the response to selection was very strong (confirming our previous studies), with a heritability value of cortisol level after ACTH injections reaching 0.64 (±0.03). The difference between the two divergent lines was around five genetic standard deviations after three selection steps. A genome-wide association study pointed out the importance of the glucocorticoid receptor gene (NR3C1) in this response. The measurement of plasma corticosteroid-binding globulin (CBG) binding capacity excluded any significant role of CBG in this selection process. The phenotypic effect of selection on body weight and growth rate was modest and/or inconsistent across generations. The HPA axis, a major neuroendocrine system involved in adaptation processes is highly heritable and responsive to genetic selection. The present experiment confirms the importance of glucocorticoid receptor polymorphism in genetic variation of HPA axis activity-in addition to the previously demonstrated role of CBG gene polymorphism. Further studies will explore the effect of this divergent selection on production and robustness.

Role of Inflammation in Depressive and Anxiety Disorders, Affect, and Cognition: Genetic and Non-Genetic Findings in the Lifelines Cohort Study

Inflammation is associated with a range of neuropsychiatric symptoms; however, the nature of the causal relationship is unclear. We used complementary non-genetic, genetic risk score (GRS), and Mendelian randomization (MR) analyses to examine whether inflammatory markers are associated with affect, depressive and anxiety disorders, and cognition. We tested in ≈ 55,098 (59% female) individuals from the Dutch Lifelines cohort the concurrent/prospective associations of C-reactive protein (CRP) with: depressive and anxiety disorders; positive/negative affect; and attention, psychomotor speed, episodic memory, and executive functioning. Additionally, we examined the association between inflammatory GRSs (CRP, interleukin-6 [IL-6], IL-6 receptor [IL-6R and soluble IL-6R (sIL-6R)], glycoprotein acetyls [GlycA]) on these same outcomes (Nmax=57,946), followed by MR analysis examining evidence of causality of CRP on outcomes (Nmax=23,268). In non-genetic analyses, higher CRP was associated with a depressive disorder, lower positive/higher negative affect, and worse executive function, attention, and psychomotor speed after adjusting for potential confounders. In genetic analyses, CRPgrs was associated with any anxiety disorder (β = 0.002, p = 0.037) whereas GlycAGRS was associated with major depressive disorder (β = 0.001, p = 0.036). Both CRPgrs (β = 0.006, p = 0.035) and GlycAGRS (β = 0.006, p = 0.049) were associated with greater negative affect. Inflammatory GRSs were not associated with cognition, except slL-6RGRS which was associated with poorer memory (β=-0.009, p = 0.018). There was weak evidence for a CRP-anxiety association using MR (β = 0.12; p = 0.054). Genetic and non-genetic analyses provide consistent evidence for an association between CRP and negative affect. These results suggest that dysregulated immune physiology may impact a broad range of trans-diagnostic affective symptoms.

Phenotypic correlations of carpal gland diverticular number with production traits and its genome-wide association analysis in multiple pig populations

Pig carpal glands play crucial roles in territorial recognition, reproductive behavior, and information exchange; however, their effects on production traits and underlying genetic mechanisms remain unclear. In this study, 1028 pigs from six populations were counted for the carpal gland diverticular numbers (CGDNs) on the left (CGDNL) and right (CGDNR) legs, and their carcass and meat quality traits were assessed. The CGDNs were significantly different among the populations, and Licha Black pigs had a lower CGDN than the Bama Xiang breed. It was also significantly different between sexes, with males having more diverticula than females (p ≤ 0.0391). Moreover, the number was asymmetric, with CGDNR being significantly higher than CGDNL. Notably, CGDNs was significantly correlated with each other in phenotype and genetics and with 24-h pH, 24-h meat color score, 24-h marbling score, fat content, moisture content, sodium salt content, and saturated fatty acid content in phenotype. Furthermore, genome-wide association analyses identified seven SNPs in association with CGDNs at a 5% genome-wide significance level, all of which were located in a 1.78-Mb (35.347-37.129 Mb) region on chromosome 1. CNC10010837 and CNC10010840 were the top SNPs: both had an additive effect of 0.789 ± 0.120 on CGDNR with p = 8.31E-10. These findings provide important insights into the functions and underlying genetic mechanisms of swine carpal glands.

Role of Inflammation in Depressive and Anxiety Disorders, Affect, and Cognition: Genetic and Non-Genetic Findings in the Lifelines Cohort Study

Background

Low-grade systemic inflammation is implicated in the pathogenesis of various neuropsychiatric conditions affecting mood and cognition. While much of the evidence concerns depression, large-scale population studies of anxiety, affect, and cognitive function are scarce. Importantly, causality remains unclear. We used complementary non-genetic, genetic risk score (GRS), and Mendelian randomization (MR) analyses to examine whether inflammatory markers are associated with affect, depressive and anxiety disorders, and cognitive performance in the Lifelines Cohort; and whether associations are likely to be causal.

Methods

Using data from up to 55,098 (59% female) individuals from the Dutch Lifelines cohort, we tested the cross-sectional and longitudinal associations of C-reactive protein (CRP) with (i) depressive and anxiety disorders; (ii) positive and negative affect scores, and (iii) five cognitive measures assessing attention, psychomotor speed, episodic memory, and executive functioning (figural fluency and working memory). Additionally, we examined the association between inflammatory marker GRSs (CRP, interleukin-6 [IL-6], IL-6 receptor [IL-6R and soluble IL-6R (sIL-6R)], glycoprotein acetyls [GlycA]) on these same outcomes (Nmax=57,946), followed by MR analysis examining evidence of causality of CRP on outcomes (Nmax=23,268). In genetic analyses, all GRSs and outcomes were z-transformed.

Results

In non-genetic analyses, higher CRP was associated with diagnosis of any depressive disorder, lower positive and higher negative affect scores, and worse performance on tests of figural fluency, attention, and psychomotor speed after adjusting for potential confounders, although the magnitude of these associations was small. In genetic analyses, CRPGRS was associated with any anxiety disorder (β=0.002, p=0.037, N=57,047) whereas GlycAGRS was associated with major depressive disorder (β=0.001, p=0.036; N=57,047). Both CRPGRS (β=0.006, p=0.035, N=57,946) and GlycAGRS (β=0.006, p=0.049; N=57,946) were associated with higher negative affect score. Inflammatory marker GRSs were not associated with cognitive performance, except sIL-6RGRS which was associated with poorer memory performance (β=-0.009, p=0.018, N=36,783). Further examination of the CRP-anxiety association using MR provided some weak evidence of causality (β=0.12; p=0.054).

Conclusions

Genetic and non-genetic analyses provide consistent evidence for an association between CRP and negative affect. Genetic analyses suggest that IL-6 signaling could be relevant for memory, and that the association between CRP and anxiety disorders could be causal. These results suggest that dysregulated immune physiology may impact a broad range of trans-diagnostic affective symptoms. However, given the small effect sizes and multiple tests conducted, future studies are required to investigate whether effects are moderated by sub-groups and whether these findings replicate in other cohorts.

Genome-Wide Association Studies for Body Conformation Traits in Korean Holstein Population

The objective of this study was to identify quantitative trait loci (QTL) and nearby candidate genes that influence body conformation traits. Phenotypic data for 24 body conformation traits were collected from a population of 2329 Korean Holstein cattle, and all animals were genotyped using the 50 K Illumina bovine SNP chip. A total of 24 genome-wide significant SNPs associated with 24 body conformation traits were identified by genome-wide association analysis. The selection of the most promising candidate genes was based on gene ontology (GO) terms and the previously identified functions that influence various body conformation traits as determined in our study. These genes include KCNA1, RYBP, PTH1R, TMIE, and GNAI3 for body traits; ANGPT1 for rump traits; MALRD1, INHBA, and HOXA13 for feet and leg traits; and CDK1, RHOBTB1, and SLC17A1 for udder traits, respectively. These findings contribute to our understanding of the genetic basis of body conformation traits in this population and pave the way for future breeding strategies aimed at enhancing desirable traits in dairy cattle.

Using residual regressions to quantify and map signal leakage in genomic prediction

BACKGROUND

Most genomic prediction applications in animal breeding use genotypes with tens of thousands of single nucleotide polymorphisms (SNPs). However, modern sequencing technologies and imputation algorithms can generate ultra-high-density genotypes (including millions of SNPs) at an affordable cost. Empirical studies have not produced clear evidence that using ultra-high-density genotypes can significantly improve prediction accuracy. However, (whole-genome) prediction accuracy is not very informative about the ability of a model to capture the genetic signals from specific genomic regions. To address this problem, we propose a simple methodology that detects chromosome regions for which a specific model (e.g., single-step genomic best linear unbiased prediction (ssGBLUP)) may fail to fully capture the genetic signal present in such segments-a phenomenon that we refer to as signal leakage. We propose to detect regions with evidence of signal leakage by testing the association of residuals from a pedigree or a genomic model with SNP genotypes. We discuss how this approach can be used to map regions with signals that are poorly captured by a model and to identify strategies to fix those problems (e.g., using a different prior or increasing marker density). Finally, we explored the proposed approach to scan for signal leakage of different models (pedigree-based, ssGBLUP, and various Bayesian models) applied to growth-related phenotypes (average daily gain and backfat thickness) in pigs.

RESULTS

We report widespread evidence of signal leakage for pedigree-based models. Including a percentage of animals with SNP data in ssGBLUP reduced the extent of signal leakage. However, local peaks of missed signals remained in some regions, even when all animals were genotyped. Using variable selection priors solves leakage points that are caused by excessive shrinkage of marker effects. Nevertheless, these models still miss signals in some regions due to low linkage disequilibrium between the SNPs on the array used and causal variants. Thus, we discuss how such problems could be addressed by adding sequence SNPs from those regions to the prediction model.

CONCLUSIONS

Residual single-marker regression analysis is a simple approach that can be used to detect regional genomic signals that are poorly captured by a model and to indicate ways to fix such problems.

Genetic association models are robust to common population kinship estimation biases

Common genetic association models for structured populations, including principal component analysis (PCA) and linear mixed-effects models (LMMs), model the correlation structure between individuals using population kinship matrices, also known as genetic relatedness matrices. However, the most common kinship estimators can have severe biases that were only recently determined. Here we characterize the effect of these kinship biases on genetic association. We employ a large simulated admixed family and genotypes from the 1000 Genomes Project, both with simulated traits, to evaluate key kinship estimators. Remarkably, we find practically invariant association statistics for kinship matrices of different bias types (matching all other features). We then prove using statistical theory and linear algebra that LMM association tests are invariant to these kinship biases, and PCA approximately so. Our proof shows that the intercept and relatedness effect coefficients compensate for the kinship bias, an argument that extends to generalized linear models. As a corollary, association testing is also invariant to changing the reference ancestral population of the kinship matrix. Lastly, we observed that all kinship estimators, except for popkin ratio-of-means, can give improper non-positive semidefinite matrices, which can be problematic although some LMMs handle them surprisingly well, and condition numbers can be used to choose kinship estimators. Overall, we find that existing association studies are robust to kinship estimation bias, and our calculations may help improve association methods by taking advantage of this unexpected robustness, as well as help determine the effects of kinship bias in related problems.

From SNP to pathway-based GWAS meta-analysis: do current meta-analysis approaches resolve power and replication in genetic association studies?

Genome-wide association studies (GWAS) have benefited greatly from enhanced high-throughput technology in recent decades. GWAS meta-analysis has become increasingly popular to highlight the genetic architecture of complex traits, informing about the replicability and variability of effect estimations across human ancestries. A wealth of GWAS meta-analysis methodologies have been developed depending on the input data and the outcome information of interest. We present a survey of current approaches from SNP to pathway-based meta-analysis by acknowledging the range of resources and methodologies in the field, and we provide a comprehensive review of different categories of Genome-Wide Meta-analysis methods employed. These methods highlight different levels at which GWAS meta-analysis may be done, including Single Nucleotide Polymorphisms, Genes and Pathways, for which we describe their framework outline. We also discuss the strengths and pitfalls of each approach and make suggestions regarding each of them.

Advances in sequencing technologies for amyotrophic lateral sclerosis research

Amyotrophic lateral sclerosis (ALS) is caused by upper and lower motor neuron loss and has a fairly rapid disease progression, leading to fatality in an average of 2-5 years after symptom onset. Numerous genes have been implicated in this disease; however, many cases remain unexplained. Several technologies are being used to identify regions of interest and investigate candidate genes. Initial approaches to detect ALS genes include, among others, linkage analysis, Sanger sequencing, and genome-wide association studies. More recently, next-generation sequencing methods, such as whole-exome and whole-genome sequencing, have been introduced. While those methods have been particularly useful in discovering new ALS-linked genes, methodological advances are becoming increasingly important, especially given the complex genetics of ALS. Novel sequencing technologies, like long-read sequencing, are beginning to be used to uncover the contribution of repeat expansions and other types of structural variation, which may help explain missing heritability in ALS. In this review, we discuss how popular and/or upcoming methods are being used to discover ALS genes, highlighting emerging long-read sequencing platforms and their role in aiding our understanding of this challenging disease.

Two decades of association mapping: Insights on disease resistance in major crops

Climate change across the globe has an impact on the occurrence, prevalence, and severity of plant diseases. About 30% of yield losses in major crops are due to plant diseases; emerging diseases are likely to worsen the sustainable production in the coming years. Plant diseases have led to increased hunger and mass migration of human populations in the past, thus a serious threat to global food security. Equipping the modern varieties/hybrids with enhanced genetic resistance is the most economic, sustainable and environmentally friendly solution. Plant geneticists have done tremendous work in identifying stable resistance in primary genepools and many times other than primary genepools to breed resistant varieties in different major crops. Over the last two decades, the availability of crop and pathogen genomes due to advances in next generation sequencing technologies improved our understanding of trait genetics using different approaches. Genome-wide association studies have been effectively used to identify candidate genes and map loci associated with different diseases in crop plants. In this review, we highlight successful examples for the discovery of resistance genes to many important diseases. In addition, major developments in association studies, statistical models and bioinformatic tools that improve the power, resolution and the efficiency of identifying marker-trait associations. Overall this review provides comprehensive insights into the two decades of advances in GWAS studies and discusses the challenges and opportunities this research area provides for breeding resistant varieties.

Dissecting Complex Traits Using Omics Data: A Review on the Linear Mixed Models and Their Application in GWAS

Genome-wide association study (GWAS) is the most popular approach to dissecting complex traits in plants, humans, and animals. Numerous methods and tools have been proposed to discover the causal variants for GWAS data analysis. Among them, linear mixed models (LMMs) are widely used statistical methods for regulating confounding factors, including population structure, resulting in increased computational proficiency and statistical power in GWAS studies. Recently more attention has been paid to pleiotropy, multi-trait, gene-gene interaction, gene-environment interaction, and multi-locus methods with the growing availability of large-scale GWAS data and relevant phenotype samples. In this review, we have demonstrated all possible LMMs-based methods available in the literature for GWAS. We briefly discuss the different LMM methods, software packages, and available open-source applications in GWAS. Then, we include the advantages and weaknesses of the LMMs in GWAS. Finally, we discuss the future perspective and conclusion. The present review paper would be helpful to the researchers for selecting appropriate LMM models and methods quickly for GWAS data analysis and would benefit the scientific society.

Gene-Based Genome-Wide Association Study Identified Genes for Agronomic Traits in Maize

A gene integrates the effects of all SNPs in its sequence span, which benefits the genome-wide association study. To explore gene-level variations affecting economic traits in maize, we extended the SNP-based GWAS analysis software Single-RunKing developed by our team to gene-based GWAS, which used the FaST-LMM algorithm to convert the linear mixed model into simple linear model association analysis. An F-test statistic was formulated to test and identify candidate genes. We compared the statistical efficiency of using 80% principal components (EPC), the first principal component (FPC), and all SNP markers (ALLSNP) as independent variables, which predecessors commonly used to integrate SNPs and represent genes. With a Huazhong Agricultural University (HAU) genomic dataset of 2.65M SNPs from 540 maize plants, 34,774 genes were annotated across the whole genome. Genome-wide association studies with 20 agronomic traits were performed using the software developed here. Another maize dataset from the Ames panel (AP) was also analyzed. The EPC method fits the model well and has good statistical efficiency. It not only overcomes the false negative problem when using all SNP markers for analysis (ALLSNP) but also solves the false positive problem of its corresponding simple linear model method EPCLM. Compared with FPC, the EPC method has higher statistical efficiency. A total of 132 quantitative trait genes (QTG) were identified for the 20 traits from HAU maize dataset and one trait of AP maize.

Identification of QTLs for wheat heading time across multiple-environments

KEY MESSAGE

The genetic response to changing climatic factors selects consistent across the tested environments and location-specific thermo-sensitive and photoperiod susceptible alleles in lower and higher altitudes, respectively, for starting flowering in winter wheat. Wheat breeders select heading date to match the most favorable conditions for their target environments and this is favored by the extensive genetic variation for this trait that has the potential to be further explored. In this study, we used a germplasm with broad geographic distribution and tested it in multi-location field trials across Germany over three years. The genotypic response to the variation in the climatic parameters depending on location and year uncovered the effect of photoperiod and spring temperatures in accelerating heading date in higher and lower latitudes, respectively. Spring temperature dominates other factors in inducing heading, whereas the higher amount of solar radiation delays it. A genome-wide scan of marker-trait associations with heading date detected two QTL: an adapted allele at locus TaHd102 on chromosome 5A that has a consistent effect on HD in German cultivars in multiple environments and a non-adapted allele at locus TaHd044 on chromosome 3A that accelerates flowering by 5.6 days. TaHd102 and TaHd044 explain 13.8% and 33% of the genetic variance, respectively. The interplay of the climatic variables led to the detection of environment specific association responding to temperature in lower latitudes and photoperiod in higher ones. Another locus TaHd098 on chromosome 5A showed epistatic interactions with 15 known regulators of flowering time when non-adapted cultivars from outside Germany were included in the analysis.

Genetic and genomic diversity in the sorghum gene bank collection of Uganda

BACKGROUND

The Plant Genetic Resources Centre at the Uganda National Gene Bank houses has over 3000 genetically diverse landraces and wild relatives of Sorghum bicolor accessions. This genetic diversity resource is untapped, under-utilized, and has not been systematically incorporated into sorghum breeding programs. In this study, we characterized the germplasm collection using whole-genome SNP markers (DArTseq). Discriminant analysis of principal components (DAPC) was implemented to study the racial ancestry of the accessions in comparison to a global sorghum diversity set and characterize the sub-groups present in the Ugandan (UG) germplasm.

RESULTS

Population structure and phylogenetic analysis revealed the presence of five subgroups among the Ugandan accessions. The samples from the highlands of the southwestern region were genetically distinct as compared to the rest of the population. This subset was predominated by the caudatum race and unique in comparison to the other sub-populations. In this study, we detected QTL for juvenile cold tolerance by genome-wide association studies (GWAS) resulting in the identification of 4 markers associated (-log10p > 3) to survival under cold stress under both field and climate chamber conditions, located on 3 chromosomes (02, 06, 09). To our best knowledge, the QTL on Sb09 with the strongest association was discovered for the first time.

CONCLUSION

This study demonstrates how genebank genomics can potentially facilitate effective and efficient usage of valuable, untapped germplasm collections for agronomic trait evaluation and subsequent allele mining. In face of adverse climate change, identification of genomic regions potentially involved in the adaptation of Ugandan sorghum accessions to cooler climatic conditions would be of interest for the expansion of sorghum production into temperate latitudes.

Detection of Genomic Imprinting for Carcass Traits in Cattle Using Imputed High-Density Genotype Data

Genomic imprinting is an epigenetic phenomenon defined as the silencing of an allele, at least partially, at a given locus based on the sex of the transmitting parent. The objective of the present study was to detect the presence of SNP-phenotype imprinting associations for carcass weight (CW), carcass conformation (CC) and carcass fat (CF) in cattle. The data used comprised carcass data, along with imputed, high-density genotype data on 618,837 single nucleotide polymorphisms (SNPs) from 23,687 cattle; all animal genotypes were phased with respect to parent of origin. Based on the phased genotypes and a series of single-locus linear models, 24, 339, and 316 SNPs demonstrated imprinting associations with CW, CC, and CF, respectively. Regardless of the trait in question, no known imprinted gene was located within 0.5 Mb of the SNPs demonstrating imprinting associations in the present study. Since all imprinting associations detected herein were at novel loci, further investigation of these regions may be warranted. Nonetheless, knowledge of these associations might be useful for improving the accuracy of genomic evaluations for these traits, as well as mate allocations systems to exploit the effects of genomic imprinting.

Deciphering the genetic basis of male fertility in Italian Brown Swiss dairy cattle

Improving reproductive performance remains a major goal in dairy cattle worldwide. Service sire has been recognized as an important factor affecting herd fertility. The main objective of this study was to reveal the genetic basis of male fertility in Italian Brown Swiss dairy cattle. Dataset included 1102 Italian Brown Swiss bulls with sire conception rate records genotyped with 454k single nucleotide polymorphisms. The analysis included whole-genome scans and gene-set analyses to identify genomic regions, individual genes and genetic mechanisms affecting Brown Swiss bull fertility. One genomic region on BTA1 showed significant additive effects. This region harbors gene RABL3 which is implicated cell proliferation and motility. Two genomic regions, located on BTA6 and BTA26, showed marked non-additive effects. These regions harbor genes, such as WDR19 and ADGRA1, that are directly involved in male fertility, including sperm motility, acrosome reaction, and embryonic development. The gene-set analysis revealed functional terms related to cell adhesion, cellular signaling, cellular transport, immune system, and embryonic development. Remarkably, a gene-set analysis also including Holstein and Jersey data, revealed significant processes that are common to the three dairy breeds, including cell migration, cell-cell interaction, GTPase activity, and the immune function. Overall, this comprehensive study contributes to a better understanding of the genetic basis of male fertility in cattle. In addition, our findings may guide the development of novel genomic strategies for improving service sire fertility in Brown Swiss cattle.

Modeling allelic diversity of multiparent mapping populations affects detection of quantitative trait loci

The search for quantitative trait loci that explain complex traits such as yield and drought tolerance has been ongoing in all crops. Methods such as biparental quantitative trait loci mapping and genome-wide association studies each have their own advantages and limitations. Multiparent advanced generation intercross populations contain more recombination events and genetic diversity than biparental mapping populations and are better able to estimate effect sizes of rare alleles than association mapping populations. Here, we discuss the results of using a multiparent advanced generation intercross population of doubled haploid maize lines created from 16 diverse founders to perform quantitative trait loci mapping. We compare 3 models that assume bi-allelic, founder, and ancestral haplotype allelic states for quantitative trait loci. The 3 methods have differing power to detect quantitative trait loci for a variety of agronomic traits. Although the founder approach finds the most quantitative trait loci, all methods are able to find unique quantitative trait loci, suggesting that each model has advantages for traits with different genetic architectures. A closer look at a well-characterized flowering time quantitative trait loci, qDTA8, which contains vgt1, highlights the strengths and weaknesses of each method and suggests a potential epistatic interaction. Overall, our results reinforce the importance of considering different approaches to analyzing genotypic datasets, and shows the limitations of binary SNP data for identifying multiallelic quantitative trait loci.

Genomic Analysis of Visceral Fat Accumulation in Holstein Cows

Visceral fat is related to important metabolic processes, including insulin sensitivity and lipid mobilization. The goal of this study was to identify individual genes, pathways, and molecular processes implicated in visceral fat deposition in dairy cows. Data from 172 genotyped Holstein cows classified at slaughterhouse as having low (n = 77; omental fold <5 mm in thickness and minimum fat deposition in omentum) or high (n = 95; omental fold ≥20 mm in thickness and marked fat deposition in omentum) omental fat were analyzed. The identification of regions with significant additive and non-additive genetic effects was performed using a two-step mixed model-based approach. Genomic scans were followed by gene-set analyses in order to reveal the genetic mechanisms controlling abdominal obesity. The association mapping revealed four regions located on BTA19, BTA20 and BTA24 with significant additive effects. These regions harbor genes, such as SMAD7, ANKRD55, and the HOXB family, that are implicated in lipolysis and insulin tolerance. Three regions located on BTA1, BTA13, and BTA24 showed marked non-additive effects. These regions harbor genes MRAP, MIS18A, PRNP and TSHZ1, that are directly implicated in adipocyte differentiation, lipid metabolism, and insulin sensitivity. The gene-set analysis revealed functional terms related to cell arrangement, cell metabolism, cell proliferation, cell signaling, immune response, lipid metabolism, and membrane permeability, among other functions. We further evaluated the genetic link between visceral fat and two metabolic disorders, ketosis, and displaced abomasum. For this, we analyzed 28k records of incidence of metabolic disorders from 14k cows across lactations using a single-step genomic BLUP approach. Notably, the region on BTA20 significantly associated with visceral fat deposition was also associated with the incidence of displaced abomasum. Overall, our findings suggest that visceral fat deposition in dairy cows is controlled by both additive and non-additive effects. We detected at least one region with marked pleiotropic effects affecting both visceral fat accumulation and displaced abomasum.

Identification of Candidate Forage Yield Genes in Sorghum (Sorghum bicolor L.) Using Integrated Genome-Wide Association Studies and RNA-Seq

Genetic dissection of forage yield traits is critical to the development of sorghum as a forage crop. In the present study, association mapping was performed with 85,585 SNP markers on four forage yield traits, namely plant height (PH), tiller number (TN), stem diameter (SD), and fresh weight per plant (FW) among 245 sorghum accessions evaluated in four environments. A total of 338 SNPs or quantitative trait nucleotides (QTNs) were associated with the four traits, and 21 of these QTNs were detected in at least two environments, including four QTNs for PH, ten for TN, six for SD, and one for FW. To identify candidate genes, dynamic transcriptome expression profiling was performed at four stages of sorghum development. One hundred and six differentially expressed genes (DEGs) that were enriched in hormone signal transduction pathways were found in all stages. Weighted gene correlation network analysis for PH and SD indicated that eight modules were significantly correlated with PH and that three modules were significantly correlated with SD. The blue module had the highest positive correlation with PH and SD, and the turquoise module had the highest negative correlation with PH and SD. Eight candidate genes were identified through the integration of genome-wide association studies (GWAS) and RNA sequencing. Sobic.004G143900, an indole-3-glycerol phosphate synthase gene that is involved in indoleacetic acid biosynthesis, was down-regulated as sorghum plants grew in height and was identified in the blue module, and Sobic.003G375100, an SD candidate gene, encoded a DNA repair RAD52-like protein 1 that plays a critical role in DNA repair-linked cell cycle progression. These findings demonstrate that the integrative analysis of omics data is a promising approach to identify candidate genes for complex traits.

Pear genetics: Recent advances, new prospects, and a roadmap for the future

Pear, belonging to the genus Pyrus, is one of the most economically important temperate fruit crops. Pyrus is an important genus of the Rosaceae family, subfamily Maloideae, and has at least 22 different species with over 5000 accessions maintained or identified worldwide. With the release of draft whole-genome sequences for Pyrus, opportunities for pursuing studies on the evolution, domestication, and molecular breeding of pear, as well as for conducting comparative genomics analyses within the Rosaceae family, have been greatly expanded. In this review, we highlight key advances in pear genetics, genomics, and breeding driven by the availability of whole-genome sequences, including whole-genome resequencing efforts, pear domestication, and evolution. We cover updates on new resources for undertaking gene identification and molecular breeding, as well as for pursuing functional validation of genes associated with desirable economic traits. We also explore future directions for "pear-omics".

Dissecting the Genetics of Early Vigour to Design Drought-Adapted Wheat

Due to the climate change and an increased frequency of drought, it is of enormous importance to identify and to develop traits that result in adaptation and in improvement of crop yield stability in drought-prone regions with low rainfall. Early vigour, defined as the rapid development of leaf area in early developmental stages, is reported to contribute to stronger plant vitality, which, in turn, can enhance resilience to erratic drought periods. Furthermore, early vigour improves weed competitiveness and nutrient uptake. Here, two sets of a multi-reference nested association mapping (MR-NAM) population of bread wheat (Triticum aestivum ssp. aestivum L.) were used to investigate early vigour in a rain-fed field environment for 3 years, and additionally assessed under controlled conditions in a greenhouse experiment. The normalised difference vegetation index (NDVI) calculated from red/infrared light reflectance was used to quantify early vigour in the field, revealing a correlation (p < 0.05; r = 0.39) between the spectral measurement and the length of the second leaf. Under controlled environmental conditions, the measured projected leaf area, using a green-pixel counter, was also correlated to the leaf area of the second leaf (p < 0.05; r = 0.38), as well as to the recorded biomass (p < 0.01; r = 0.71). Subsequently, genetic determination of early vigour was tested by conducting a genome-wide association study (GWAS) for the proxy traits, revealing 42 markers associated with vegetation index and two markers associated with projected leaf area. There are several quantitative trait loci that are collocated with loci for plant developmental traits including plant height on chromosome 2D (log₁₀ (P) = 3.19; PVE = 0.035), coleoptile length on chromosome 1B (-log₁₀ (P) = 3.24; PVE = 0.112), as well as stay-green and vernalisation on chromosome 5A (-log₁₀ (P) = 3.14; PVE = 0.115).

Genetic Bases of Complex Traits: From Quantitative Trait Loci to Prediction

Conceived as a general introduction to the book, this chapter is a reminder of the core concepts of genetic mapping and molecular marker-based prediction. It provides an overview of the principles and the evolution of methods for mapping the variation of complex traits, and methods for QTL-based prediction of human disease risk and animal and plant breeding value. The principles of linkage-based and linkage disequilibrium-based QTL mapping methods are described in the context of the simplest, single-marker, methods. Methodological evolutions are analysed in relation with their ability to account for the complexity of the genotype-phenotype relations. Main characteristics of the genetic architecture of complex traits, drawn from QTL mapping works using large populations of unrelated individuals, are presented. Methods combining marker-QTL association data into polygenic risk score that captures part of an individual's susceptibility to complex diseases are reviewed. Principles of best linear mixed model-based prediction of breeding value in animal- and plant-breeding programs using phenotypic and pedigree data, are summarized and methods for moving from BLUP to marker-QTL BLUP are presented. Factors influencing the additional genetic progress achieved by using molecular data and rules for their optimization are discussed.

Dissection of Quantitative Blackleg Resistance Reveals Novel Variants of Resistance Gene Rlm9 in Elite Brassica napus

Blackleg is one of the major fungal diseases in oilseed rape/canola worldwide. Most commercial cultivars carry R gene-mediated qualitative resistances that confer a high level of race-specific protection against Leptosphaeria maculans, the causal fungus of blackleg disease. However, monogenic resistances of this kind can potentially be rapidly overcome by mutations in the pathogen's avirulence genes. To counteract pathogen adaptation in this evolutionary arms race, there is a tremendous demand for quantitative background resistance to enhance durability and efficacy of blackleg resistance in oilseed rape. In this study, we characterized genomic regions contributing to quantitative L. maculans resistance by genome-wide association studies in a multiparental mapping population derived from six parental elite varieties exhibiting quantitative resistance, which were all crossed to one common susceptible parental elite variety. Resistance was screened using a fungal isolate with no corresponding avirulence (AvrLm) to major R genes present in the parents of the mapping population. Genome-wide association studies revealed eight significantly associated quantitative trait loci (QTL) on chromosomes A07 and A09, with small effects explaining 3-6% of the phenotypic variance. Unexpectedly, the qualitative blackleg resistance gene Rlm9 was found to be located within a resistance-associated haploblock on chromosome A07. Furthermore, long-range sequence data spanning this haploblock revealed high levels of single-nucleotide and structural variants within the Rlm9 coding sequence among the parents of the mapping population. The results suggest that novel variants of Rlm9 could play a previously unknown role in expression of quantitative disease resistance in oilseed rape.

REHE: Fast variance components estimation for linear mixed models

Linear mixed models are widely used in ecological and biological applications, especially in genetic studies. Reliable estimation of variance components is crucial for using linear mixed models. However, standard methods, such as the restricted maximum likelihood (REML), are computationally inefficient in large samples and may be unstable with small samples. Other commonly used methods, such as the Haseman-Elston (HE) regression, may yield negative estimates of variances. Utilizing regularized estimation strategies, we propose the restricted Haseman-Elston (REHE) regression and REHE with resampling (reREHE) estimators, along with an inference framework for REHE, as fast and robust alternatives that provide nonnegative estimates with comparable accuracy to REML. The merits of REHE are illustrated using real data and benchmark simulation studies.

An Efficient Score Test Integrated with Empirical Bayes for Genome-Wide Association Studies

Many methods used in multi-locus genome-wide association studies (GWAS) have been developed to improve statistical power. However, most existing multi-locus methods are not quicker than single-locus methods. To address this concern, we proposed a fast score test integrated with Empirical Bayes (ScoreEB) for multi-locus GWAS. Firstly, a score test was conducted for each single nucleotide polymorphism (SNP) under a linear mixed model (LMM) framework, taking into account the genetic relatedness and population structure. Then, all of the potentially associated SNPs were selected with a less stringent criterion. Finally, Empirical Bayes in a multi-locus model was performed for all of the selected SNPs to identify the true quantitative trait nucleotide (QTN). Our new method ScoreEB adopts the similar strategy of multi-locus random-SNP-effect mixed linear model (mrMLM) and fast multi-locus random-SNP-effect EMMA (FASTmrEMMA), and the only difference is that we use the score test to select all the potentially associated markers. Monte Carlo simulation studies demonstrate that ScoreEB significantly improved the computational efficiency compared with the popular methods mrMLM, FASTmrEMMA, iterative modified-sure independence screening EM-Bayesian lasso (ISIS EM-BLASSO), hybrid of restricted and penalized maximum likelihood (HRePML) and genome-wide efficient mixed model association (GEMMA). In addition, ScoreEB remained accurate in QTN effect estimation and effectively controlled false positive rate. Subsequently, ScoreEB was applied to re-analyze quantitative traits in plants and animals. The results show that ScoreEB not only can detect previously reported genes, but also can mine new genes.

Optimizing genomic control in mixed model associations with binary diseases

Complex computation and approximate solution hinder the application of generalized linear mixed models (GLMM) into genome-wide association studies. We extended GRAMMAR to handle binary diseases by considering genomic breeding values (GBVs) estimated in advance as a known predictor in genomic logit regression, and then reduced polygenic effects by regulating downward genomic heritability to control false negative errors produced in the association tests. Using simulations and case analyses, we showed in optimizing GRAMMAR, polygenic effects and genomic controls could be evaluated using the fewer sampling markers, which extremely simplified GLMM-based association analysis in large-scale data. Further, joint association analysis for quantitative trait nucleotide (QTN) candidates chosen by multiple testing offered significant improved statistical power to detect QTNs over existing methods.

Genomic and phenotypic insights from an atlas of genetic effects on DNA methylation

Characterizing genetic influences on DNA methylation (DNAm) provides an opportunity to understand mechanisms underpinning gene regulation and disease. In the present study, we describe results of DNAm quantitative trait locus (mQTL) analyses on 32,851 participants, identifying genetic variants associated with DNAm at 420,509 DNAm sites in blood. We present a database of >270,000 independent mQTLs, of which 8.5% comprise long-range (trans) associations. Identified mQTL associations explain 15-17% of the additive genetic variance of DNAm. We show that the genetic architecture of DNAm levels is highly polygenic. Using shared genetic control between distal DNAm sites, we constructed networks, identifying 405 discrete genomic communities enriched for genomic annotations and complex traits. Shared genetic variants are associated with both DNAm levels and complex diseases, but only in a minority of cases do these associations reflect causal relationships from DNAm to trait or vice versa, indicating a more complex genotype-phenotype map than previously anticipated.

Genome-wide hierarchical mixed model association analysis

In genome-wide mixed model association analysis, we stratified the genomic mixed model into two hierarchies to estimate genomic breeding values (GBVs) using the genomic best linear unbiased prediction and statistically infer the association of GBVs with each SNP using the generalized least square. The hierarchical mixed model (Hi-LMM) can correct confounders effectively with polygenic effects as residuals for association tests, preventing potential false-negative errors produced with genome-wide rapid association using mixed model and regression or an efficient mixed-model association expedited (EMMAX). Meanwhile, the Hi-LMM performs the same statistical power as the exact mixed model association and the same computing efficiency as EMMAX. When the GBVs have been estimated precisely, the Hi-LMM can detect more quantitative trait nucleotides (QTNs) than existing methods. Especially under the Hi-LMM framework, joint association analysis can be made straightforward to improve the statistical power of detecting QTNs.

An investigation of the effects of BMPR1B, BMP15, and GDF9 genes on litter size in Ramlıç and Dağlıç sheep

This study was carried out to determine the presence of polymorphisms in genes affecting litter size. The SNPs in bone morphogenetic protein receptor type 1B (BMPR1B), bone morphogenetic protein 15 (BMP15), and growth differentiation factor 9 (GDF9) genes were detected in 60 uniparous and 60 multiparous ewes from Ramlıç and Dağlıç breeds. The ewes are maintained in nine public herds at the breeding station of the Afyonkarahisar Sheep and Goats Breeders' Association and lambed in two consecutive breeding seasons. PCR and DNA sequencing analyses were conducted, and 36, 4, and 11 SNPs in Ramlıç and 40, 3, and 11 SNPs in Dağlıç were detected in BMPR1B, BMP15, and GDF9 genes, respectively. A total of 16 SNPs in Ramlıç and 10 SNPs in Dağlıç breeds for three genes were found to be significant (P<0.05). The resulting analyses showed that four SNPs (g.49496G>A, c.1658A>C, c.2037C>T, c.2053C>T) of the BMPR1B gene and one deletion mutation (c.28_30delCTT) in the BMP15 gene of the Ramlıç breed as well as five SNPs (c.1487C>A, c.2492C>T, c.2523G>A, c.2880A>G, and c.2763G>A) of the BMPR1B gene of the Dağlıç breed have significant positive regression coefficients in the desired direction of the rare allele. The observed mutations have potential to be used as genetic markers in the selection of prolific animals for both breeds.

A novel deletion in FLOWERING LOCUS T modulates flowering time in winter oilseed rape

A novel structural variant was discovered in the FLOWERING LOCUS T orthologue BnaFT.A02 by long-read sequencing. Nested association mapping in an elite winter oilseed rape population revealed that this 288 bp deletion associates with early flowering, putatively by modification of binding-sites for important flowering regulation genes. Perfect timing of flowering is crucial for optimal pollination and high seed yield. Extensive previous studies of flowering behavior in Brassica napus (canola, rapeseed) identified mutations in key flowering regulators which differentiate winter, semi-winter and spring ecotypes. However, because these are generally fixed in locally adapted genotypes, they have only limited relevance for fine adjustment of flowering time in elite cultivar gene pools. In crosses between ecotypes, the ecotype-specific major-effect mutations mask minor-effect loci of interest for breeding. Here, we investigated flowering time in a multiparental mapping population derived from seven elite winter oilseed rape cultivars which are fixed for major-effect mutations separating winter-type rapeseed from other ecotypes. Association mapping revealed eight genomic regions on chromosomes A02, C02 and C03 associating with fine modulation of flowering time. Long-read genomic resequencing of the seven parental lines identified seven structural variants coinciding with candidate genes for flowering time within chromosome regions associated with flowering time. Segregation patterns for these variants in the elite multiparental population and a diversity set of winter types using locus-specific assays revealed significant associations with flowering time for three deletions on chromosome A02. One of these was a previously undescribed 288 bp deletion within the second intron of FLOWERING LOCUS T on chromosome A02, emphasizing the advantage of long-read sequencing for detection of structural variants in this size range. Detailed analysis revealed the impact of this specific deletion on flowering-time modulation under extreme environments and varying day lengths in elite, winter-type oilseed rape.

Genome-Wide Association Study Identifies New Candidate Markers for Somatic Cells Score in a Local Dairy Sheep

In the Mediterranean basin countries, the dairy sheep production is usually based on local breeds, which are very well-adapted to their production systems and environments and can indeed guarantee income, employment, and economic viability in areas where production alternatives are scarce or non-existent. Mastitis is still one of the greatest problems affecting commercial milk production. However, genetic evaluation of mastitis is particularly difficult because of its low heritability and the categorical nature of the trait. The aim of this study was to identify genomic regions putatively associated with somatic cells count (SCC) in the local economically important Valle del Belice sheep breed using of deregressed breeding values (DEBV) as response variables. All the samples were genotyped using the Illumina OvineSNP50K BeadChip. Genome-wide association analysis was carried out based on regression of DEBV. A total of eight markers were found to be significantly associated with log-transformed SCC. Several candidate genes associated with SCC were identified related to immunity system and udder conformation. The results can help improving the competitiveness of the local Valle del Belìce breed. Further studies considering a higher sample size or independent population will be needed to confirm our results.

Status and prospects of genome-wide association studies in plants

Genome-wide association studies (GWAS) have developed into a powerful and ubiquitous tool for the investigation of complex traits. In large part, this was fueled by advances in genomic technology, enabling us to examine genome-wide genetic variants across diverse genetic materials. The development of the mixed model framework for GWAS dramatically reduced the number of false positives compared with naïve methods. Building on this foundation, many methods have since been developed to increase computational speed or improve statistical power in GWAS. These methods have allowed the detection of genomic variants associated with either traditional agronomic phenotypes or biochemical and molecular phenotypes. In turn, these associations enable applications in gene cloning and in accelerated crop breeding through marker assisted selection or genetic engineering. Current topics of investigation include rare-variant analysis, synthetic associations, optimizing the choice of GWAS model, and utilizing GWAS results to advance knowledge of biological processes. Ongoing research in these areas will facilitate further advances in GWAS methods and their applications.

Genome-wide association studies provide insights into the genetic determination of fruit traits of pear

Pear is a major fruit tree crop distributed worldwide, yet its breeding is a very time-consuming process. To facilitate molecular breeding and gene identification, here we have performed genome-wide association studies (GWAS) on eleven fruit traits. We identify 37 loci associated with eight fruit quality traits and five loci associated with three fruit phenological traits. Scans for selective sweeps indicate that traits including fruit stone cell content, organic acid and sugar contents might have been under continuous selection during breeding improvement. One candidate gene, PbrSTONE, identified in GWAS, has been functionally verified to be involved in the regulation of stone cell formation, one of the most important fruit quality traits in pear. Our study provides insights into the complex fruit related biology and identifies genes controlling important traits in pear through GWAS, which extends the genetic resources and basis for facilitating molecular breeding in perennial trees.

Deshrinking ridge regression for genome-wide association studies

MOTIVATION

Genome-wide association studies (GWAS) are still the primary steps toward gene discovery. The urgency is more obvious in the big data era when GWAS are conducted simultaneously for thousand traits, e.g. transcriptomic and metabolomic traits. Efficient mixed model association (EMMA) and genome-wide efficient mixed model association (GEMMA) are the widely used methods for GWAS. An algorithm with high computational efficiency is badly needed. It is interesting to note that the test statistics of the ordinary ridge regression (ORR) have the same patterns across the genome as those obtained from the EMMA method. However, ORR has never been used for GWAS due to its severe shrinkage on the estimated effects and the test statistics.

RESULTS

We introduce a degree of freedom for each marker effect obtained from ORR and use it to deshrink both the estimated effect and the standard error so that the Wald test of ORR is brought back to the same level as that of EMMA. The new method is called deshrinking ridge regression (DRR). By evaluating the methods under three different model sizes (small, medium and large), we demonstrate that DRR is more generalized for all model sizes than EMMA, which only works for medium and large models. Furthermore, DRR detect all markers in a simultaneous manner instead of scanning one marker at a time. As a result, the computational time complexity of DRR is much simpler than EMMA and about m (number of genetic variants) times simpler than that of GEMMA when the sample size is way smaller than the number of markers.

CONTACT

shizhong.xu@ucr.edu.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Estimating FST and kinship for arbitrary population structures

F_ST and kinship are key parameters often estimated in modern population genetics studies in order to quantitatively characterize structure and relatedness. Kinship matrices have also become a fundamental quantity used in genome-wide association studies and heritability estimation. The most frequently-used estimators of F_ST and kinship are method-of-moments estimators whose accuracies depend strongly on the existence of simple underlying forms of structure, such as the independent subpopulations model of non-overlapping, independently evolving subpopulations. However, modern data sets have revealed that these simple models of structure likely do not hold in many populations, including humans. In this work, we analyze the behavior of these estimators in the presence of arbitrarily-complex population structures, which results in an improved estimation framework specifically designed for arbitrary population structures. After generalizing the definition of F_ST to arbitrary population structures and establishing a framework for assessing bias and consistency of genome-wide estimators, we calculate the accuracy of existing F_ST and kinship estimators under arbitrary population structures, characterizing biases and estimation challenges unobserved under their originally-assumed models of structure. We then present our new approach, which consistently estimates kinship and F_ST when the minimum kinship value in the dataset is estimated consistently. We illustrate our results using simulated genotypes from an admixture model, constructing a one-dimensional geographic scenario that departs nontrivially from the independent subpopulations model. Our simulations reveal the potential for severe biases in estimates of existing approaches that are overcome by our new framework. This work may significantly improve future analyses that rely on accurate kinship and F_ST estimates.

Kinship coefficients and F_ST, which measure relatedness and population structure, respectively, are important quantities needed to accurately perform various analyses on genetic data, including genome-wide association studies and heritability estimation. However, existing estimators require restrictive assumptions of independence that are not met by real human and other datasets. In this work we find that existing estimators can be severely biased under reasonable scenarios, first by theoretically determining their properties, and then using an admixture simulation to illustrate our findings. In particular, we find that existing F_ST estimators are downwardly biased, and that existing kinship matrix estimators have related biases that are on average downward and of similar magnitude but vary for every pair of individuals. These insights led us to a new estimation framework for kinship and F_ST that is practically unbiased for any population structure, as demonstrated by theory and simulations. Our new approaches—available as open-source R packages—are easy to use and are more widely applicable than existing approaches, and they are likely to improve downstream analyses that require accurate kinship and F_ST estimates.

Genetic consistency between gait analysis by accelerometry and evaluation scores at breeding shows for the selection of jumping competition horses

The aim was to assess the efficiency of gaits characteristics in improving jumping performance of sport horses and confront accelerometers and judge scores for this purpose. A sample of 1,477 young jumping horses were measured using accelerometers for walk, trot, and canter. Of these, 702 were genotyped with 541,175 SNPs after quality control. Dataset of 26,914 horses scored by judges in breeding shows for gaits and dataset of 142,682 horses that performed in jumping competitions were used. Analysis of accelerometric data defined three principal components from 64% to 89% of variability explained for each gait. Animal mixed models were used to estimate genetic parameters with the inclusion to up 308,105 ancestors for the relationship matrix. Fixed effects for the accelerometric variables included velocity, gender, age, and event. A GWAS was performed on residuals with the fixed effect of each SNP. The GWAS did not reveal other QTLs for gait traits than the one related to the height at withers. The accelerometric principal components were highly heritable for the one linked to stride frequency and dorsoventral displacement at trot (0.53) and canter (0.41) and moderately for the one linked to longitudinal activities (0.33 for trot, 0.19 for canter). Low heritabilities were found for the walk traits. The genetic correlations of the accelerometric principal components with the jumping competition were essentially nil, except for a negative correlation with longitudinal activity at canter (-0.19). The genetic correlation between the judges’ scores and the jumping competition reached 0.45 for canter (0.31 for trot and 0.17 for walk). But these correlations turned negative when the scores were corrected for the known parental breeding value for competition at the time of the judging. In conclusion, gait traits were not helpful to select for jumping performances. Different gaits may be suitable for a good jumping horse.

Mixed logistic regression in genome-wide association studies

Background

Mixed linear models (MLM) have been widely used to account for population structure in case-control genome-wide association studies, the status being analyzed as a quantitative phenotype. Chen et al. proved in 2016 that this method is inappropriate in some situations and proposed GMMAT, a score test for the mixed logistic regression (MLR). However, this test does not produces an estimation of the variants’ effects. We propose two computationally efficient methods to estimate the variants’ effects. Their properties and those of other methods (MLM, logistic regression) are evaluated using both simulated and real genomic data from a recent GWAS in two geographically close population in West Africa.

Results

We show that, when the disease prevalence differs between population strata, MLM is inappropriate to analyze binary traits. MLR performs the best in all circumstances. The variants’ effects are well evaluated by our methods, with a moderate bias when the effect sizes are large. Additionally, we propose a stratified QQ-plot, enhancing the diagnosis of p values inflation or deflation when population strata are not clearly identified in the sample.

Conclusion

The two proposed methods are implemented in the R package milorGWAS available on the CRAN. Both methods scale up to at least 10,000 individuals. The same computational strategies could be applied to other models (e.g. mixed Cox model for survival analysis).

Whole Genome Sequence Data Provides Novel Insights Into the Genetic Architecture of Meat Quality Traits in Beef

Tenderness is a major quality attribute for fresh beef steaks in the United States, and meat quality traits in general are suitable candidates for genomic research. The objectives of the present analysis were to (1) perform genome-wide association (GWA) analysis for marbling, Warner-Bratzler shear force (WBSF), tenderness, and connective tissue using whole-genome data in an Angus population, (2) identify enriched pathways in each GWA analysis; (3) construct a protein-protein interaction network using the associated genes and (4) perform a μ-calpain proteolysis assessment for associated structural proteins. An Angus-sired population of 2,285 individuals was assessed. Animals were transported to a commercial packing plant and harvested at an average age of 457 ± 46 days. After 48 h postmortem, marbling was recorded by graders' visual appraisal. Two 2.54-cm steaks were sampled from each muscle for recording of WBSF, and tenderness, and connective tissue by a sensory panel. The relevance of additive effects on marbling, WBSF, tenderness, and connective tissue was evaluated on a genome-wide scale using a two-step mixed model-based approach in single-trait analysis. A tissue-restricted gene enrichment was performed for each GWA where all polymorphisms with an association p-value lower than 1 × 10-3 were included. The genes identified as associated were included in a protein-protein interaction network and a candidate structural protein assessment of proteolysis analyses. A total of 1,867, 3,181, 3,926, and 3,678 polymorphisms were significantly associated with marbling, WBSF, tenderness, and connective tissue, respectively. The associate region on BTA29 (36,432,655-44,313,046 bp) harbors 13 highly significant markers for meat quality traits. Enrichment for the GO term GO:0005634 (Nucleus), which includes transcription factors, was evident. The final protein-protein network included 431 interations between 349 genes. The 42 most important genes based on significance that encode structural proteins were included in a proteolysis analysis, and 81% of these proteins were potential μ-Calpain substrates. Overall, this comprehensive study unraveled genetic variants, genes and mechanisms of action responsible for the variation in meat quality traits. Our findings can provide opportunities for improving meat quality in beef cattle via marker-assisted selection.

Identifying blood pressure loci whose effects are modulated by multiple lifestyle exposures

Although multiple lifestyle exposures simultaneously impact blood pressure (BP) and cardiovascular health, most analysis so far has considered each single lifestyle exposure (e.g., smoking) at a time. Here, we exploit gene-multiple lifestyle exposure interactions to find novel BP loci. For each of 6,254 Framingham Heart Study participants, we computed lifestyle risk score (LRS) value by aggregating the risk of four lifestyle exposures (smoking, alcohol, education, and physical activity) on BP. Using the LRS, we performed genome-wide gene-environment interaction analysis in systolic and diastolic BP using the joint 2 degree of freedom (DF) and 1 DF interaction tests. We identified one genome-wide significant (p < 5 × 10-8 ) and 11 suggestive (p < 1 × 10-6 ) loci. Gene-environment analysis using single lifestyle exposures identified only one of the 12 loci. Nine of the 12 BP loci detected were novel. Loci detected by the LRS were located within or nearby genes with biologically plausible roles in the pathophysiology of hypertension, including KALRN, VIPR2, SNX1, and DAPK2. Our results suggest that simultaneous consideration of multiple lifestyle exposures in gene-environment interaction analysis can identify additional loci missed by single lifestyle approaches.

Physiological and morphological responses of different spring barley genotypes to water deficit and associated QTLs

Water deficit is one of the major limitations to food production worldwide and most climate change scenarios predict an aggravation of the situation. To face the expected increase in drought stress in the coming years, breeders are working to elucidate the genetic control of barley growth and productivity traits under water deficit. Barley is known as a relatively drought tolerant crop and genetic variability was observed for drought tolerance traits. The objectives of the present study were the quantification of morphological and physiological responses in a collection of 209 spring barley genotypes to drought stress, and the genetic analysis by genome-wide association study to find quantitative trait loci (QTL) and the allele contributions for each of the investigated traits. In six pot experiments, 209 spring barley genotypes were grown under a well-watered and water-limited regime. Stress phases were initiated individually for each genotype at the beginning of tillering and spiking for the vegetative- and the generative stage experiments, respectively, and terminated when the transpiration rates of stress treatments reached 10% of the well-watered control. After the stress phase, a total of 42 productivity related traits such as the dry matter of plant organs, tiller number, leaf length, leaf area, amount of water soluble carbohydrates in the stems, proline content in leaves and osmotic adjustment of corresponding well-watered and stressed plants were analysed, and QTL analyses were performed to find marker-trait associations. Significant water deficit effects were observed for almost all traits and significant genotype x treatment interactions (GxT) were observed for 37 phenotypic traits. Genome-wide association studies (GWAS) revealed 77 significant loci associated with 16 phenotypic traits during the vegetative stage experiment and a total of 85 significant loci associated with 13 phenotypic traits during the generative stage experiment for traits such as leaf area, number of green leaves, grain yield, harvest index and stem length. For traits with significant GxT interactions, genotypic differences for relative values were analysed using one way ANOVA. More than 110 loci for GxT interaction were found for 17 phenotypic traits explaining in many cases more than 50% of the genetic variance.

Gene-Based Association Testing of Dichotomous Traits With Generalized Functional Linear Mixed Models Using Extended Pedigrees: Applications to Age-Related Macular Degeneration

Genetics plays a role in age-related macular degeneration (AMD), a common cause of blindness in the elderly. There is a need for powerful methods for carrying out region-based association tests between a dichotomous trait like AMD and genetic variants on family data. Here, we apply our new generalized functional linear mixed models (GFLMM) developed to test for gene-based association in a set of AMD families. Using common and rare variants, we observe significant association with two known AMD genes: CFH and ARMS2. Using rare variants, we find suggestive signals in four genes: ASAH1, CLEC6A, TMEM63C, and SGSM1. Intriguingly, ASAH1 is down-regulated in AMD aqueous humor, and ASAH1 deficiency leads to retinal inflammation and increased vulnerability to oxidative stress. These findings were made possible by our GFLMM which model the effect of a major gene as a fixed mean, the polygenic contributions as a random variation, and the correlation of pedigree members by kinship coefficients. Simulations indicate that the GFLMM likelihood ratio tests (LRTs) accurately control the Type I error rates. The LRTs have similar or higher power than existing retrospective kernel and burden statistics. Our GFLMM-based statistics provide a new tool for conducting family-based genetic studies of complex diseases. Supplementary materials for this article, including a standardized description of the materials available for reproducing the work, are available as an online supplement.

Multi-trait multi-locus SEM model discriminates SNPs of different effects

Background

There is a plethora of methods for genome-wide association studies. However, only a few of them may be classified as multi-trait and multi-locus, i.e. consider the influence of multiple genetic variants to several correlated phenotypes.

Results

We propose a multi-trait multi-locus model which employs structural equation modeling (SEM) to describe complex associations between SNPs and traits - multi-trait multi-locus SEM (mtmlSEM). The structure of our model makes it possible to discriminate pleiotropic and single-trait SNPs of direct and indirect effect. We also propose an automatic procedure to construct the model using factor analysis and the maximum likelihood method. For estimating a large number of parameters in the model, we performed Bayesian inference and implemented Gibbs sampling. An important feature of the model is that it correctly copes with non-normally distributed variables, such as some traits and variants.

Conclusions

We applied the model to Vavilov’s collection of 404 chickpea (Cicer arietinum L.) accessions with 20-fold cross-validation. We analyzed 16 phenotypic traits which we organized into five groups and found around 230 SNPs associated with traits, 60 of which were of pleiotropic effect. The model demonstrated high accuracy in predicting trait values.

Genomic background and genetic relationships between boar taint and fertility traits in German Landrace and Large White

BACKGROUND

Due to ethical reasons, surgical castration of young male piglets in their first week of life without anesthesia will be banned in Germany from 2021. Breeding against boar taint is already implemented in sire breeds of breeding organizations but in recent years a low demand made this trait economically less important. The objective of this study was to estimate heritabilities and genetic relationships between boar taint compounds androstenone and skatole and maternal/paternal reproduction traits in 4'924 Landrace (LR) and 4'299 Large White (LW) animals from nucleus populations. Additionally, genome wide association analysis (GWAS) was performed per trait and breed to detect SNP marker with possible pleiotropic effects that are associated with boar taint and fertility.

RESULTS

Estimated heritabilities (h²) were 0.48 (±0.08) for LR (0.39 ± 0.07 for LW) for androstenone and 0.52 (±0.08) for LR (0.32 ± 0.07 for LW) for skatole. Heritabilities for reproduction did not differ between breeds except age at first insemination (LR: h² = 0.27 (±0.05), LW: h² = 0.34 (±0.05)). Estimates of genetic correlation (r_g) between boar taint and fertility were different in LR and LW breeds. In LR an unfavorable r_g of 0.31 (±0.15) was observed between androstenone and number of piglets born alive, whereas this r_g in LW (- 0.15 (±0.16)) had an opposite sign. A similar breed-specific difference is observed between skatole and sperm count. Within LR, the r_g of 0.08 (±0.13) indicates no relationship between the traits, whereas the r_g of - 0.37 (±0.14) in LW points to an unfavorable relationship. In LR GWAS identified QTL regions on SSC5 (21.1-22.3 Mb) for androstenone and on SSC6 (5.5-7.5 Mb) and SSC14 (141.1-141.6 Mb) for skatole. For LW, one marker was found on SSC17 at 48.1 Mb for androstenone and one QTL on SSC14 between 140.5 Mb and 141.6 Mb for skatole.

CONCLUSION

Knowledge about such genetic correlations could help to balance conventional breeding programs with boar taint in maternal breeds. QTL regions with unfavorable pleiotropic effects on boar taint and fertility could have deleterious consequences in genomic selection programs. Constraining the weighting of these QTL in the genomic selection formulae may be a useful strategy to avoid physiological imbalances.

Reduction of slaughter age of Hanwoo steers by early genotyping based on meat yield index

Objective: This study was conducted to determine early hereditary endowment to establish a short-term feeding program.Methods: Hanwoo steers (n = 140) were equally distributed into four groups (35/group) based on genetic meat yield index (MYI) viz. the greatest, great, low, and the lowest at Jukam Hanwoo farm, Goheung. All animals were fed in group pens (5 animals/pen) with similar feed depending on the growth stage. Rice straw was provided ad libitum, whereas concentrate was fed at 5.71 kg during the growing period (6 to 13 mo) and 9.4 kg during the fattening period (13 to 28 mo). Body weight (BW) was measured at two-month intervals, whereas carcass weight was determined at slaughtering at about 31 months of age. The Affymetrix Bovine Axiom Array 640K single nucleotide polymorphism (SNP) chip was used to determine the meat quantity-related gene in the blood.Results: After 6 months, the highest (p<0.05) BW was observed in the greatest MYI group (190.77 kg) and the lowest (p<0.05) in the lowest MYI group (173.51 kg). The great MYI group also showed significantly (p<0.05) higher BW than the lowest MYI group. After 16 and 24 months, the greatest MYI group had the highest BW gain (p<0.05) and were therefore slaughtered the earliest. Carcass weight was significantly (p<0.05) higher in the greatest and the great MYI groups followed by the low and the lowest MYI groups. Back-fat thickness in the greatest MYI group was highly correlated to carcass weight and marbling score. The SNP array analysis identified the carcass-weight related gene BTB-01280026 with an additive effect. The steers with the allele increasing carcass weight had heavier slaughter weight of about 12 kg.Conclusion: Genetic MYI is a potential tool for calf selection, which will reduce the slaughter age while simultaneously increasing carcass weight, back-fat thickness, and marbling score.

Efficient Estimation and Applications of Cross-Validated Genetic Predictions to Polygenic Risk Scores and Linear Mixed Models

Large-scale cohorts with combined genetic and phenotypic data, coupled with methodological advances, have produced increasingly accurate genetic predictors of complex human phenotypes called polygenic risk scores (PRSs). In addition to the potential translational impacts of identifying at-risk individuals, PRS are being utilized for a growing list of scientific applications, including causal inference, identifying pleiotropy and genetic correlation, and powerful gene-based and mixed-model association tests. Existing PRS approaches rely on external large-scale genetic cohorts that have also measured the phenotype of interest. They further require matching on ancestry and genotyping platform or imputation quality. In this work, we present a novel reference-free method to produce a PRS that does not rely on an external cohort. We show that naive implementations of reference-free PRS either result in substantial overfitting or prohibitive increases in computational time. We show that our algorithm avoids both of these issues and can produce informative in-sample PRSs over a single cohort without overfitting. We then demonstrate several novel applications of reference-free PRSs, including detection of pleiotropy across 246 metabolic traits and efficient mixed-model association testing.

Gene presence-absence variation associates with quantitative Verticillium longisporum disease resistance in Brassica napus

Although copy number variation (CNV) and presence-absence variation (PAV) have been discovered in selected gene families in most crop species, the global prevalence of these polymorphisms in most complex genomes is still unclear and their influence on quantitatively inherited agronomic traits is still largely unknown. Here we analyze the association of gene PAV with resistance of oilseed rape (Brassica napus) against the important fungal pathogen Verticillium longisporum, as an example for a complex, quantitative disease resistance in the strongly rearranged genome of a recent allopolyploid crop species. Using Single Nucleotide absence Polymorphism (SNaP) markers to efficiently trace PAV in breeding populations, we significantly increased the resolution of loci influencing V. longisporum resistance in biparental and multi-parental mapping populations. Gene PAV, assayed by resequencing mapping parents, was observed in 23-51% of the genes within confidence intervals of quantitative trait loci (QTL) for V. longisporum resistance, and high-priority candidate genes identified within QTL were all affected by PAV. The results demonstrate the prominent role of gene PAV in determining agronomic traits, suggesting that this important class of polymorphism should be exploited more systematically in future plant breeding.

Gene mapping and genomic prediction of bull fertility using sex chromosome markers

Service sire has been recognized as an important factor affecting dairy herd fertility. Our group has reported promising results on gene mapping and genomic prediction of dairy bull fertility using autosomal SNP markers. Little is known, however, about the genetic contribution of sex chromosomes, which are enriched in genes related to sexual development and reproduction. As such, the main goal of this study was to investigate the effect of SNP markers on X and Y chromosomes (BTAX and BTAY, respectively) on sire conception rate (SCR) in US Holstein bulls. The analysis included a total of 5,014 bulls with SCR records and genotypes for roughly 291k SNP located on the autosomes, 1.5k SNP located on the pseudoautosomal region (PAR), 13.7k BTAX-specific SNP, and 24 BTAY-specific SNP. We first performed genomic scans of the sex chromosomes, and then we evaluated the genomic prediction of SCR including BTAX SNP markers in the predictive models. Two markers located on PAR and 3 markers located on the X-specific region showed significant associations with sire fertility. Interestingly, these regions harbor genes, such as FAM9B, TBL1X, and PIH1D3, that are directly implicated in testosterone concentration, spermatogenesis, and sperm motility. On the other hand, BTAY showed very low genetic variability, and none of the segregating markers were associated with SCR. Notably, model predictive ability was largely improved by including BTAX markers. Indeed, the combination of autosomal with BTAX SNP delivered predictive correlations around 0.343, representing an increase in accuracy of about 7.5% compared with the standard whole autosomal genome approach. Overall, this study provides evidence of the importance of both PAR and X-specific regions in male fertility in dairy cattle. These findings may help to improve conception rates in dairy herds through accurate genome-guided decisions on bull fertility.

A fast-linear mixed model for genome-wide haplotype association analysis: application to agronomic traits in maize

Background

Haplotypes combine the effects of several single nucleotide polymorphisms (SNPs) with high linkage disequilibrium, which benefit the genome-wide association analysis (GWAS). In the haplotype association analysis, both haplotype alleles and blocks are tested. Haplotype alleles can be inferred with the same statistics as SNPs in the linear mixed model, while blocks require the formulation of unified statistics to fit different genetic units, such as SNPs, haplotypes, and copy number variations.

Results

Based on the FaST-LMM, the fastLmPure function in the R/RcppArmadillo package has been introduced to speed up genome-wide regression scans by a re-weighted least square estimation. When large or highly significant blocks are tested based on EMMAX, the genome-wide haplotype association analysis takes only one to two rounds of genome-wide regression scans. With a genomic dataset of 541,595 SNPs from 513 maize inbred lines, 90,770 haplotype blocks were constructed across the whole genome, and three types of markers (SNPs, haplotype alleles, and haplotype blocks) were genome-widely associated with 17 agronomic traits in maize using the software developed here.

Conclusions

Two SNPs were identified for LNAE, four haplotype alleles for TMAL, LNAE, CD, and DTH, and only three blocks reached the significant level for TMAL, CD, and KNPR. Compared to the R/lm function, the computational time was reduced by ~ 10–15 times.