Accounting for multiple comparisons in a genome-wide association study (GWAS)

Genome-wide association study uncovers new genetic loci and candidate genes underlying seed chilling-germination in maize

As one of the major crops, maize (Zea mays L.) is mainly distributed in tropical and temperate regions. However, with the changes of the environments, chilling stress has become a significantly abiotic stress affecting seed germination and thus the reproductive and biomass accumulation of maize. Herein, we investigated five seed germination-related phenotypes among 300 inbred lines under low-temperature condition (10 °C). By combining 43,943 single nucleotide polymorphisms (SNPs), a total of 15 significant (P < 2.03 ×  10^-6) SNPs were identified to correlate with seed germination under cold stress based on the FarmCPU model in GWAS, among which three loci were repeatedly associated with multiple traits. Ten gene models were closely linked to these three variations, among which Zm00001d010454, Zm00001d010458, Zm00001d010459, and Zm00001d050021 were further verified by candidate gene association study and expression pattern analysis. Importantly, these candidate genes were previously reported to involve plant tolerance to chilling stress and other abiotic stress. Our findings contribute to the understanding of the genetic and molecular mechanisms underlying chilling germination in maize.

Genomic Regions Related to White/Black Tail Feather Color in Dwarf Chickens Identified Using a Genome-Wide Association Study

Although the genetic foundation of chicken body feather color has been extensively explored, that of tail feather color remains poorly understood. In the present study, we used a synthetic chicken dwarf line (DW), derived from hybrids bred between a black tail chicken breed, Rhode Island Red (RIR), and a white tail breed, dwarf layer (DL), to investigate the genetic rules associated white/black tail color. Even though the body feathers are predominantly red, the DW line still comprises individuals with black or white tails after more than 10 generations of self-crossing and selection for the body feather color. We first performed four crosses using the DW chickens, including black-tailed males to females, reciprocal crosses between the black and white, and white males to females to elucidate the inheritance pattern of the white/black tail. We also performed a genome-wide association (GWA) analysis to determine the candidate genomic regions underlying the tail feather color using black tail chickens from the RIR and DW lines and white individuals from the DW line. In the crossing experiment, we found that (i) the white/black tail feather color is independent of body feather color; (ii) the phenotype is a simple autosomal trait; and (iii) the white is dominant to the black in the DW line. The GWA results showed that seven single-nucleotide polymorphisms (SNPs) on chromosome 24 were significantly correlated with tail feather color. The significant region (3.97-4.26 Mb) comprises nine known genes (NECTIN1, THY1, gga-mir-1466, USP2, C1QTNF5, RNF26, MCAM, CBL, and CCDC153) and five anonymous genes. This study revealed that the white/black tail feather trait is autosome-linked in DW chickens. Fourteen genes were found in the significant ~0.29 Mb genomic region, and some, especially MCAM, are suggested to play critical roles in the determination of white/black tail feather color. Our research is the first study on the genetics underlying tail feather color and could help further the understanding of feather pigmentation in chickens.

Correlation of [18F]florbetaben textural features and age of onset of Alzheimer's disease: a principal components analysis approach

BACKGROUND

When Alzheimer's disease (AD) is occurring at an early onset before 65 years old, its clinical course is generally more aggressive than in the case of a late onset. We aim at identifying [[Formula: see text]F]florbetaben PET biomarkers sensitive to differences between early-onset Alzheimer's disease (EOAD) and late-onset Alzheimer's disease (LOAD). We conducted [[Formula: see text]F]florbetaben PET/CT scans of 43 newly diagnosed AD subjects. We calculated 93 textural parameters for each of the 83 Hammers areas. We identified 41 independent principal components for each brain region, and we studied their Spearman correlation with the age of AD onset, by taking into account multiple comparison corrections. Finally, we calculated the probability that EOAD and LOAD patients have different amyloid-[Formula: see text] ([Formula: see text]) deposition by comparing the mean and the variance of the significant principal components obtained in the two groups with a 2-tailed Student's t-test.

RESULTS

We found that four principal components exhibit a significant correlation at a 95% confidence level with the age of onset in the left lateral part of the anterior temporal lobe, the right anterior orbital gyrus of the frontal lobe, the right lateral orbital gyrus of the frontal lobe and the left anterior part of the superior temporal gyrus. The data are consistent with the hypothesis that EOAD patients have a significantly different [[Formula: see text]F]florbetaben uptake than LOAD patients in those four brain regions.

CONCLUSIONS

Early-onset AD implies a very irregular pattern of [Formula: see text] deposition. The authors suggest that the identified textural features can be used as quantitative biomarkers for the diagnosis and characterization of EOAD patients.

Analysis of Brugada syndrome loci reveals that fine-mapping clustered GWAS hits enhances the annotation of disease-relevant variants

Genome-wide association studies (GWASs) are instrumental in identifying loci harboring common single-nucleotide variants (SNVs) that affect human traits and diseases. GWAS hits emerge in clusters, but the focus is often on the most significant hit in each trait- or disease-associated locus. The remaining hits represent SNVs in linkage disequilibrium (LD) and are considered redundant and thus frequently marginally reported or exploited. Here, we interrogate the value of integrating the full set of GWAS hits in a locus repeatedly associated with cardiac conduction traits and arrhythmia, SCN5A-SCN10A. Our analysis reveals 5 common 7-SNV haplotypes (Hap1-5) with 2 combinations associated with life-threatening arrhythmia-Brugada syndrome (the risk Hap1/1 and protective Hap2/3 genotypes). Hap1 and Hap2 share 3 SNVs; thus, this analysis suggests that assuming redundancy among clustered GWAS hits can lead to confounding disease-risk associations and supports the need to deconstruct GWAS data in the context of haplotype composition.

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.

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.

Multi-Parental Populations Suitable for Identifying Sources of Resistance to Powdery Mildew in Winter Wheat

Wheat (Triticum aestivum L.) is one of the world's staple food crops and one of the most devastating foliar diseases attacking wheat is powdery mildew (PM). In Denmark only a few specific fungicides are available for controlling PM and the use of resistant cultivars is often recommended. In this study, two Chinese wheat landraces and two synthetic hexaploid wheat lines were used as donors for creating four multi-parental populations with a total of 717 individual lines to identify new PM resistance genetic variants. These lines and the nine parental lines (including the elite cultivars used to create the populations) were genotyped using a 20 K Illumina SNP chip, which resulted in 8,902 segregating single nucleotide polymorphisms for assessment of the population structure and whole genome association study. The largest genetic difference among the lines was between the donors and the elite cultivars, the second largest genetic difference was between the different donors; a difference that was also reflected in differences between the four multi-parental populations. The 726 lines were phenotyped for PM resistance in 2017 and 2018. A high PM disease pressure was observed in both seasons, with severities ranging from 0 to >50%. Whole genome association studies for genetic variation in PM resistance in the populations revealed significant markers mapped to either chromosome 2A, B, or D in each of the four populations. However, linkage disequilibrium between these putative quantitative trait loci (QTL) were all above 0.80, probably representing a single QTL. A combined analysis of all the populations confirmed this result and the most associated marker explained 42% of the variation in PM resistance. This study gives both knowledge about the resistance as well as molecular tools and plant material that can be utilised in marker-assisted selection. Additionally, the four populations produced in this study are highly suitable for association studies of other traits than PM resistance.

Genome-wide association studies for mottled eggs in chickens using a high-density single-nucleotide polymorphism array

Mottled eggs in layer chickens are gaining increasing attention because of the economic impact on the egg industry caused by the reduced sale value of commodity eggs. However, the genetic architecture underlying mottled eggs is not well understood. The genetic architecture underlying the mottled egg trait was investigated using genome-wide association studies (GWAS) by high-density arrays, using a total of 407 pink eggs and 799 blue eggs from an F₂ resource population generated by crossing Dongxiang Blue-shelled and White Leghorn chickens. The mottled egg score in blue eggs was found to be higher than that in pink eggs. The single-nucleotide polymorphism heritability of mottled egg at laying day and storage for 7 days was 0.18 and 0.20, respectively. Bivariate GWAS provided 29 significant loci, mainly located on GGA2, GGA3, GGA8, GGA10, GGA15, GGA17, and GGA23, affecting mottled egg on laying day. Candidate genes RIMS2, SLC25A32, RIMBP2, VPS13B, and RGS3 were obtained for mottled eggshell by bivariate GWAS and gene annotation. Our findings provide new insights into the genetic architecture of mottled egg in hens, and demonstrate that a genomic selection method would be profitable for breeding out the mottled egg trait.

Identification of Novel Genomic Associations and Gene Candidates for Grain Starch Content in Sorghum

Starch accumulated in the endosperm of cereal grains as reserve energy for germination serves as a staple in human and animal nutrition. Unraveling genetic control for starch metabolism is important for breeding grains with high starch content. In this study, we used a sorghum association panel with 389 individuals and 141,557 single nucleotide polymorphisms (SNPs) to fit linear mixed models (LMM) for identifying genomic regions and potential candidate genes associated with starch content. Three associated genomic regions, one in chromosome (chr) 1 and two novel associations in chr-8, were identified using combination of LMM and Bayesian sparse LMM. All significant SNPs were located within protein coding genes, with SNPs ∼ 52 Mb of chr-8 encoding a Casperian strip membrane protein (CASP)-like protein (Sobic.008G111500) and a heat shock protein (HSP) 90 (Sobic.008G111600) that were highly expressed in reproductive tissues including within the embryo and endosperm. The HSP90 is a potential hub gene with gene network of 75 high-confidence first interactors that is enriched for five biochemical pathways including protein processing. The first interactors of HSP90 also showed high transcript abundance in reproductive tissues. The candidates of this study are likely involved in intricate metabolic pathways and represent candidate gene targets for source-sink activities and drought and heat stress tolerance during grain filling.

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

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

Genome-Wide Association Studies for the Concentration of Albumin in Colostrum and Serum in Chinese Holstein

Albumin can be of particular benefit in fighting infections for newborn calves due to its anti-inflammatory and anti-oxidative stress properties. To identify the candidate genes related to the concentration of albumin in colostrum and serum, we collected the colostrum and blood samples from 572 Chinese Holstein cows within 24 h after calving and measured the concentration of albumin in the colostrum and serum using the ELISA methods. The cows were genotyped with GeneSeek 150 K chips (containing 140,668 single nucleotide polymorphisms; SNPs). After quality control, we performed GWASs via GCTA software with 91,620 SNPs and 563 cows. Consequently, 9 and 7 genome-wide significant SNPs (false discovery rate (FDR) at 1%) were identified. Correspondingly, 42 and 206 functional genes that contained or were approximate to (±1 Mbp) the significant SNPs were acquired. Integrating the biological process of these genes and the reported QTLs for immune and inflammation traits in cattle, 3 and 12 genes were identified as candidates for the concentration of colostrum and serum albumin, respectively; these are RUNX1, CBR1, OTULIN,CDK6, SHARPIN, CYC1, EXOSC4, PARP10, NRBP2, GFUS, PYCR3, EEF1D, GSDMD, PYCR2 and CXCL12. Our findings provide important information for revealing the genetic mechanism behind albumin concentration and for molecular breeding of disease-resistance traits in dairy cattle.

The genetic basis for the inverse relationship between rheumatoid arthritis and schizophrenia

Introduction

Rheumatoid arthritis is a common autoimmune disease and schizophrenia is a relatively common and debilitating neurological disorder. There are several common features between rheumatoid arthritis and schizophrenia. The inverse relationship between rheumatoid arthritis and schizophrenia has been replicated in several studies. Despite evidence for an inverse epidemiological relationship and negative correlations for risk between rheumatoid arthritis and schizophrenia, there are no biological data that directly support this inverse relationship.

Materials and Methods’

We meta‐analyzed the genome‐wide association studies to investigate the shared association loci between rheumatoid arthritis and schizophrenia at the genome‐wide scale. Rheumatoid arthritis‐ and schizophrenia‐associated loci in most recent genome‐wide association studies of rheumatoid arthritis and schizophrenia were tested. Genetic risk score analysis was also conducted to investigate the collective contribution of schizophrenia risk loci to rheumatoid arthritis risk.

Results

Rheumatoid arthritis and schizophrenia meta‐genome‐wide association study showed a significant peak at the major histocompatibility complex locus on chromosome 6 in both rheumatoid arthritis‐schizophrenia meta‐genome‐wide association study and inverted meta‐genome‐wide association study datasets. Testing rheumatoid arthritis‐ and schizophrenia‐associated loci outside the human leukocyte antigen region showed no association with both rheumatoid arthritis and schizophrenia at a genome‐wide level of significance. Weighted genetic risk scores showed no evidence for a statistically significant association between rheumatoid arthritis and schizophrenia.

Conclusion

The finding of our study is consistent with the role of the major histocompatibility complex locus in the genetic correlation between rheumatoid arthritis and schizophrenia, and suggests that either schizophrenia has an autoimmune basis and/or rheumatoid arthritis has an active neurological component.

Using imputed whole-genome sequence variants to uncover candidate mutations and genes affecting milking speed and temperament in Holstein cattle

Milking speed (MS) and temperament (MT) are 2 workability traits of great importance in dairy cattle production and breeding. This is mainly due to an increased intensification of the worldwide production systems and greater adoption of precision technologies with less human-cattle interaction. Both MS and MT are heritable traits and thus, genomic selection is a promising tool to expedite their genetic progress. However, the genetic architecture and biological mechanisms underlying the phenotypic expression of these traits remain underexplored. In this study, we investigated the association of >5.7 million imputed whole-genome sequence variants with MT and MS in 4,381 and 4,219 North American Holstein cattle, respectively. The statistical analyses were performed using a mixed linear model fitting a polygenic effect. We detected 40 and 35 significant SNPs independently associated with MT and MS, respectively, which were distributed across 26 chromosomes. Eight candidate genes (GRIN3A, KCNJ3, BOSTAUV1R417, BOSTAUV1R419, MAP2K5, KCTD3, GAP43, and LSAMP) were suggested to play an important role in MT as they are involved in biologically relevant pathways, such as glutamatergic synapse, vomeronasal receptor and oxytocin signaling. Within their coding and upstream sequences, we used an independent data set to further detect or validate significantly differentiated SNP between cattle breeds with known differences in MT. There were fewer candidate genes potentially implicated in MS, but immunity-related genes (e.g., BOLA-NC1 and LOC512672), also identified in other populations, were validated in this study. The significant SNP and novel candidate genes identified contribute to a better understanding of the biological mechanisms underlying both traits in dairy cattle. This information will also be useful for the optimization of prediction of genomic breeding values by giving greater weights to SNP located in the genomic regions identified.

Genome-Wide Association Study and Pathway Analysis for Heterophil/Lymphocyte (H/L) Ratio in Chicken

Disease control and prevention have been critical factors in the dramatic growth of the poultry industry. Disease resistance in chickens can be improved through genetic selection for immunocompetence. The heterophil/lymphocyte ratio (H/L) in the blood reflects the immune system status of chickens. Our objective was to conduct a genome-wide association study (GWAS) and pathway analysis to identify possible biological mechanisms involved in H/L traits. In this study, GWAS for H/L was performed in 1317 Cobb broilers to identify significant single-nucleotide polymorphisms (SNPs) associated with H/L. Eight SNPs (p < 1/8068) reached a significant level of association. The significant SNP on GGA 19 (chicken chromosome 19) was in the gene for complement C1q binding protein (C1QBP). The wild-type and mutant individuals showed significant differences in H/L at five identified SNPs (p < 0.05). According to the results of pathway analysis, nine associated pathways (p < 0.05) were identified. By combining GWAS with pathway analysis, we found that all SNPs after QC explained 12.4% of the phenotypic variation in H/L, and 52 SNPs associated with H/L explained as much as 9.7% of the phenotypic variation in H/L. Our findings contribute to understanding of the genetic regulation of H/L and provide theoretical support.

Myocardial Infarction Associates With a Distinct Pericoronary Adipose Tissue Radiomic Phenotype: A Prospective Case-Control Study

OBJECTIVES

This study sought to determine whether coronary computed tomography angiography (CCTA)-based radiomic analysis of pericoronary adipose tissue (PCAT) could distinguish patients with acute myocardial infarction (MI) from patients with stable or no coronary artery disease (CAD).

BACKGROUND

Imaging of PCAT with CCTA enables detection of coronary inflammation. Radiomics involves extracting quantitative features from medical images to create big data and identify novel imaging biomarkers.

METHODS

In a prospective case-control study, 60 patients with acute MI underwent CCTA within 48 h of admission, before invasive angiography. These subjects were matched to patients with stable CAD (n = 60) and controls with no CAD (n = 60) by age, sex, risk factors, medications, and CT tube voltage. PCAT was segmented around the proximal right coronary artery (RCA) in all patients and around culprit and nonculprit lesions in patients with MI. PCAT segmentations were analyzed using Radiomics Image Analysis software.

RESULTS

Of 1,103 calculated radiomic parameters, 20.3% differed significantly between MI patients and controls, and 16.5% differed between patients with MI and stable CAD (critical p < 0.0006); whereas none differed between patients with stable CAD and controls. On cluster analysis, the most significant radiomic parameters were texture or geometry based. At 6 months post-MI, there was no significant change in the PCAT radiomic profile around the proximal RCA or nonculprit lesions. Using machine learning (XGBoost), a model integrating clinical features (risk factors, serum lipids, high-sensitivity C-reactive protein), PCAT attenuation, and radiomic parameters provided superior discrimination of acute MI (area under the receiver operator characteristic curve [AUC]: 0.87) compared with a model with clinical features and PCAT attenuation (AUC: 0.77; p = 0.001) or clinical features alone (AUC: 0.76; p < 0.001).

CONCLUSIONS

Patients with acute MI have a distinct PCAT radiomic phenotype compared with patients with stable or no CAD. Using machine learning, a radiomics-based model outperforms a PCAT attenuation-based model in accurately identifying patients with MI.

Thickness of the cerebral cortex shows positive association with blood levels of triacylglycerols carrying 18-carbon fatty acids

Perturbations in fatty acid (FA) metabolism as well as thinning of the cerebral cortex have been associated with cognitive decline in the elderly. Predominant FAs in the brain are docosahexaenoic acid (DHA) and arachidonic acid (ARA). Approximately 2-8% of esterified DHA and 3-5% of esterified ARA in the brain are replaced daily. DHA and ARA are derivatives of 18-carbon essential FAs, α-linolenic acid and linoleic acid, that must be imported into the brain from the circulation. In blood, FAs are primarily transported in triacylglycerols (TAGs) from which they can be released at the blood-brain-barrier and transported inside the brain. We show that circulating levels of TAGs carrying 18-carbon FAs are positively associated with cortical thickness in middle-aged adults. These associations are stronger in cortical regions with higher expression of genes regulating long-chain FA metabolism and cellular membranes, and cortical thickness in the same regions may be related to cognitive performance.

Combining Random Forests and a Signal Detection Method Leads to the Robust Detection of Genotype-Phenotype Associations

Genome wide association studies (GWAS) are a well established methodology to identify genomic variants and genes that are responsible for traits of interest in all branches of the life sciences. Despite the long time this methodology has had to mature the reliable detection of genotype-phenotype associations is still a challenge for many quantitative traits mainly because of the large number of genomic loci with weak individual effects on the trait under investigation. Thus, it can be hypothesized that many genomic variants that have a small, however real, effect remain unnoticed in many GWAS approaches. Here, we propose a two-step procedure to address this problem. In a first step, cubic splines are fitted to the test statistic values and genomic regions with spline-peaks that are higher than expected by chance are considered as quantitative trait loci (QTL). Then the SNPs in these QTLs are prioritized with respect to the strength of their association with the phenotype using a Random Forests approach. As a case study, we apply our procedure to real data sets and find trustworthy numbers of, partially novel, genomic variants and genes involved in various egg quality traits.

Find the Needle in the Haystack, Then Find It Again: Replication and Validation in the 'Omics Era

Advancements in high-throughput technologies have made it feasible to study thousands of biological pathways simultaneously for a holistic assessment of health and disease risk via ‘omics platforms. A major challenge in ‘omics research revolves around the reproducibility of findings—a feat that hinges upon balancing false-positive associations with generalizability. Given the foundational role of reproducibility in scientific inference, replication and validation of ‘omics findings are cornerstones of this effort. In this narrative review, we define key terms relevant to replication and validation, present issues surrounding each concept with historical and contemporary examples from genomics (the most well-established and upstream ‘omics), discuss special issues and unique considerations for replication and validation in metabolomics (an emerging field and most downstream ‘omics for which best practices remain yet to be established), and make suggestions for future research leveraging multiple ‘omics datasets.

Genomics models in radiotherapy: From mechanistic to machine learning

Machine learning (ML) provides a broad framework for addressing high-dimensional prediction problems in classification and regression. While ML is often applied for imaging problems in medical physics, there are many efforts to apply these principles to biological data toward questions of radiation biology. Here, we provide a review of radiogenomics modeling frameworks and efforts toward genomically guided radiotherapy. We first discuss medical oncology efforts to develop precision biomarkers. We next discuss similar efforts to create clinical assays for normal tissue or tumor radiosensitivity. We then discuss modeling frameworks for radiosensitivity and the evolution of ML to create predictive models for radiogenomics.

Identification of nodulation-related genes in Medicago truncatula using genome-wide association studies and co-expression networks

Genome-wide association studies (GWAS) have proven to be a valuable approach for identifying genetic intervals associated with phenotypic variation in Medicago truncatula. These intervals can vary in size, depending on the historical local recombination. Typically, significant intervals span numerous gene models, limiting the ability to resolve high-confidence candidate genes underlying the trait of interest. Additional genomic data, including gene co-expression networks, can be combined with the genetic mapping information to successfully identify candidate genes. Co-expression network analysis provides information about the functional relationships of each gene through its similarity of expression patterns to other well-defined clusters of genes. In this study, we integrated data from GWAS and co-expression networks to pinpoint candidate genes that may be associated with nodule-related phenotypes in M. truncatula. We further investigated a subset of these genes and confirmed that several had existing evidence linking them nodulation, including MEDTR2G101090 (PEN3-like), a previously validated gene associated with nodule number.

Genetic Architecture of Early Vigor Traits in Wild Soybean

A worldwide food shortage has been projected as a result of the current increase in global population and climate change. In order to provide sufficient food to feed more people, we must develop crops that can produce higher yields. Plant early vigor traits, early growth rate (EGR), early plant height (EPH), inter-node length, and node count are important traits that are related to crop yield. Glycine soja, the wild counterpart to cultivated soybean, Glycine max, harbors much higher genetic diversity and can grow in diverse environments. It can also cross easily with cultivated soybean. Thus, it holds a great potential in developing soybean cultivars with beneficial agronomic traits. In this study, we used 225 wild soybean accessions originally from diverse environments across its geographic distribution in East Asia. We quantified the natural variation of several early vigor traits, investigated the relationships among them, and dissected the genetic basis of these traits by applying a Genome-Wide Association Study (GWAS) with genome-wide single nucleotide polymorphism (SNP) data. Our results showed positive correlation between all early vigor traits studied. A total of 12 SNPs significantly associated with EPH were identified with 4 shared with EGR. We also identified two candidate genes, Glyma.07G055800.1 and Glyma.07G055900.1, playing important roles in influencing trait variation in both EGR and EPH in G. soja.

Identification of Age-Specific and Common Key Regulatory Mechanisms Governing Eggshell Strength in Chicken Using Random Forests

In today's chicken egg industry, maintaining the strength of eggshells in longer laying cycles is pivotal for improving the persistency of egg laying. Eggshell development and mineralization underlie a complex regulatory interplay of various proteins and signaling cascades involving multiple organ systems. Understanding the regulatory mechanisms influencing this dynamic trait over time is imperative, yet scarce. To investigate the temporal changes in the signaling cascades, we considered eggshell strength at two different time points during the egg production cycle and studied the genotype-phenotype associations by employing the Random Forests algorithm on chicken genotypic data. For the analysis of corresponding genes, we adopted a well established systems biology approach to delineate gene regulatory pathways and master regulators underlying this important trait. Our results indicate that, while some of the master regulators (Slc22a1 and Sox11) and pathways are common at different laying stages of chicken, others (e.g., Scn11a, St8sia2, or the TGF- β pathway) represent age-specific functions. Overall, our results provide: (i) significant insights into age-specific and common molecular mechanisms underlying the regulation of eggshell strength; and (ii) new breeding targets to improve the eggshell quality during the later stages of the chicken production cycle.

Biological Machine Learning Combined with Campylobacter Population Genomics Reveals Virulence Gene Allelic Variants Cause Disease

Highly dimensional data generated from bacterial whole-genome sequencing is providing an unprecedented scale of information that requires an appropriate statistical analysis framework to infer biological function from populations of genomes. The application of genome-wide association study (GWAS) methods is an appropriate framework for bacterial population genome analysis that yields a list of candidate genes associated with a phenotype, but it provides an unranked measure of importance. Here, we validated a novel framework to define infection mechanism using the combination of GWAS, machine learning, and bacterial population genomics that ranked allelic variants that accurately identified disease. This approach parsed a dataset of 1.2 million single nucleotide polymorphisms (SNPs) and indels that resulted in an importance ranked list of associated alleles of porA in Campylobacter jejuni using spatiotemporal analysis over 30 years. We validated this approach using previously proven laboratory experimental alleles from an in vivo guinea pig abortion model. This framework, termed µPathML, defined intestinal and extraintestinal groups that have differential allelic porA variants that cause abortion. Divergent variants containing indels that defeated automated annotation were rescued using biological context and knowledge that resulted in defining rare, divergent variants that were maintained in the population over two continents and 30 years. This study defines the capability of machine learning coupled with GWAS and population genomics to simultaneously identify and rank alleles to define their role in infectious disease mechanisms.

The emerging landscape of health research based on biobanks linked to electronic health records: Existing resources, statistical challenges, and potential opportunities

Biobanks linked to electronic health records provide rich resources for health-related research. With improvements in administrative and informatics infrastructure, the availability and utility of data from biobanks have dramatically increased. In this paper, we first aim to characterize the current landscape of available biobanks and to describe specific biobanks, including their place of origin, size, and data types. The development and accessibility of large-scale biorepositories provide the opportunity to accelerate agnostic searches, expedite discoveries, and conduct hypothesis-generating studies of disease-treatment, disease-exposure, and disease-gene associations. Rather than designing and implementing a single study focused on a few targeted hypotheses, researchers can potentially use biobanks' existing resources to answer an expanded selection of exploratory questions as quickly as they can analyze them. However, there are many obvious and subtle challenges with the design and analysis of biobank-based studies. Our second aim is to discuss statistical issues related to biobank research such as study design, sampling strategy, phenotype identification, and missing data. We focus our discussion on biobanks that are linked to electronic health records. Some of the analytic issues are illustrated using data from the Michigan Genomics Initiative and UK Biobank, two biobanks with two different recruitment mechanisms. We summarize the current body of literature for addressing these challenges and discuss some standing open problems. This work complements and extends recent reviews about biobank-based research and serves as a resource catalog with analytical and practical guidance for statisticians, epidemiologists, and other medical researchers pursuing research using biobanks.

A Geometry-Based Multiple Testing Correction for Contingency Tables by Truncated Normal Distribution

Inference procedure is a critical step of experimental researches to draw scientific conclusions especially in multiple testing. The false positive rate increases unless the unadjusted marginal p-values are corrected. Therefore, a multiple testing correction is necessary to adjust the p-values based on the number of tests to control type I error. We propose a multiple testing correction of MAX-test for a contingency table, where multiple χ2-tests are applied based on a truncated normal distribution (TND) estimation method by Botev. The table and tests are defined geometrically by contour hyperplanes in the degrees of freedom (df) dimensional space. A linear algebraic method called spherization transforms the shape of the space, defined by the contour hyperplanes of the distribution of tables sharing the same marginal counts. So, the stochastic distributions of these tables are transformed into a standard multivariate normal distribution in df-dimensional space. Geometrically, the p-value is defined by a convex polytope consisted of truncating hyperplanes of test’s contour lines in df-dimensional space. The TND approach of the Botev method was used to estimate the corrected p. Finally, the features of our approach were extracted using a real GWAS data.

Genetic Landscape of Nonobstructive Azoospermia and New Perspectives for the Clinic

Nonobstructive azoospermia (NOA) represents the most severe expression of male infertility, involving around 1% of the male population and 10% of infertile men. This condition is characterised by the inability of the testis to produce sperm cells, and it is considered to have an important genetic component. During the last two decades, different genetic anomalies, including microdeletions of the Y chromosome, karyotype defects, and missense mutations in genes involved in the reproductive function, have been described as the primary cause of NOA in many infertile men. However, these alterations only explain around 25% of azoospermic cases, with the remaining patients showing an idiopathic origin. Recent studies clearly suggest that the so-called idiopathic NOA has a complex aetiology with a polygenic inheritance, which may alter the spermatogenic process. Although we are far from a complete understanding of the molecular mechanisms underlying NOA, the use of the new technologies for genetic analysis has enabled a considerable increase in knowledge during the last years. In this review, we will provide a comprehensive and updated overview of the genetic basis of NOA, with a special focus on the possible application of the recent insights in clinical practice.

Biomarkers in patients with heart failure and central sleep apnoea: findings from the SERVE-HF trial

Aims

The Treatment of Sleep‐Disordered Breathing with Predominant Central Sleep Apnoea by Adaptive Servo Ventilation in Patients with Heart Failure trial investigated the effects of adaptive servo‐ventilation (ASV) (vs. control) on outcomes of 1325 patients with heart failure and reduced ejection fraction (HFrEF) and central sleep apnoea (CSA). The primary outcome (a composite of all‐cause death or unplanned HF hospitalization) did not differ between the two groups. However, all‐cause and cardiovascular (CV) mortality were higher in the ASV group. Circulating biomarkers may help in better ascertain patients' risk, and this is the first study applying a large set of circulating biomarkers in patients with both HFrEF and CSA.

Methods and results

Circulating protein‐biomarkers (n = 276) ontologically involved in CV pathways, were studied in 749 (57% of the trial population) patients (biomarker substudy), to investigate their association with the study outcomes (primary outcome, CV death and all‐cause death). The mean age was 69 ± 10 years, and > 90% were male. The groups (ASV vs. control and biomarker substudy vs. no biomarker) were well balanced. The “best” clinical prognostic model included male sex, systolic blood pressure < 120 mmHg, diabetes, loop diuretic, cardiac device, 6‐min walking test distance, and N‐terminal pro BNP as the strongest prognosticators. On top of the “best” clinical prognostic model, the biomarkers that significantly improved both the discrimination (c‐index) and the net reclassification index (NRI) of the model were soluble suppression of tumorigenicity 2 for the primary outcome; neurogenic locus notch homolog protein 3 (Notch‐3) for CV‐death and all‐cause death; and growth differentiation factor 15 (GDF‐15) for all‐cause death only.

Conclusions

We studied 276 circulating biomarkers in patients with HFrEF and central sleep apnoea; of these biomarkers, three added significant prognostic information on top of the best clinical model: soluble suppression of tumorigenicity 2 (primary outcome), Notch‐3 (CV and all‐cause death), and GDF‐15 (all‐cause death).

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

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

Genetic Architecture and Genomic Prediction of Cooking Time in Common Bean (Phaseolus vulgaris L.)

Cooking time of the common bean is an important trait for consumer preference, with implications for nutrition, health, and environment. For efficient germplasm improvement, breeders need more information on the genetics to identify fast cooking sources with good agronomic properties and molecular breeding tools. In this study, we investigated a broad genetic variation among tropical germplasm from both Andean and Mesoamerican genepools. Four populations were evaluated for cooking time (CKT), water absorption capacity (WAC), and seed weight (SdW): a bi-parental RIL population (DxG), an eight-parental Mesoamerican MAGIC population, an Andean (VEF), and a Mesoamerican (MIP) breeding line panel. A total of 922 lines were evaluated in this study. Significant genetic variation was found in all populations with high heritabilities, ranging from 0.64 to 0.89 for CKT. CKT was related to the color of the seed coat, with the white colored seeds being the ones that cooked the fastest. Marker trait associations were investigated by QTL analysis and GWAS, resulting in the identification of 10 QTL. In populations with Andean germplasm, an inverse correlation of CKT and WAC, and also a QTL on Pv03 that inversely controls CKT and WAC (CKT3.2/WAC3.1) were observed. WAC7.1 was found in both Mesoamerican populations. QTL only explained a small part of the variance, and phenotypic distributions support a more quantitative mode of inheritance. For this reason, we evaluated how genomic prediction (GP) models can capture the genetic variation. GP accuracies for CKT varied, ranging from good results for the MAGIC population (0.55) to lower accuracies in the MIP panel (0.22). The phenotypic characterization of parental material will allow for the cooking time trait to be implemented in the active germplasm improvement programs. Molecular breeding tools can be developed to employ marker-assisted selection or genomic selection, which looks to be a promising tool in some populations to increase the efficiency of breeding activities.

Genetic dissection of canine hip dysplasia phenotypes and osteoarthritis reveals three novel loci

Background

Hip dysplasia and osteoarthritis continue to be prevalent problems in veterinary and human medicine. Canine hip dysplasia is particularly problematic as it massively affects several large-sized breeds and can cause a severe impairment of the quality of life. In Finland, the complex condition is categorized to five classes from normal to severe dysplasia, but the categorization includes several sub-traits: congruity of the joint, Norberg angle, subluxation degree of the joint, shape and depth of the acetabulum, and osteoarthritis. Hip dysplasia and osteoarthritis have been proposed to have separate genetic etiologies.

Results

Using Fédération Cynologique Internationale -standardized ventrodorsal radiographs, German shepherds were rigorously phenotyped for osteoarthritis, and for joint incongruity by Norberg angle and femoral head center position in relation to dorsal acetabular edge. The affected dogs were categorized into mild, moderate and severe dysplastic phenotypes using official hip scores. Three different genome-wide significant loci were uncovered. The strongest candidate genes for hip joint incongruity were noggin (NOG), a bone and joint developmental gene on chromosome 9, and nanos C2HC-type zinc finger 1 (NANOS1), a regulator of matrix metalloproteinase 14 (MMP14) on chromosome 28. Osteoarthritis mapped to a long intergenic region on chromosome 1, between genes encoding for NADPH oxidase 3 (NOX3), an intriguing candidate for articular cartilage degradation, and AT-rich interactive domain 1B (ARID1B) that has been previously linked to joint laxity.

Conclusions

Our findings highlight the complexity of canine hip dysplasia phenotypes. In particular, the results of this study point to the potential involvement of specific and partially distinct loci and genes or pathways in the development of incongruity, mild dysplasia, moderate-to-severe dysplasia and osteoarthritis of canine hip joints. Further studies should unravel the unique and common mechanisms for the various sub-traits.

Landscape of Loci and Candidate Genes for Muscle Fatty Acid Composition in Pigs Revealed by Multiple Population Association Analysis

Genome wide association analyses in diverse populations can identify complex trait loci that are specifically present in one population or shared across multiple populations, which help to better understand the genetic architecture of complex traits in a broader genetic context. In this study, we conducted genome-wide association studies and meta-analysis for 38 fatty acid composition traits with 12–19 million imputed genome sequence SNPs in 2446 pigs from six populations, encompassing White Duroc × Erhualian F₂, Sutai, Duroc-Landrace-Yorkshire (DLY) three-way cross, Laiwu, Erhualian, and Bamaxiang pigs that were originally genotyped with 60 K or 1.4 million single nucleotide polymorphism (SNP) chips. The analyses uncovered 285 lead SNPs (P < 5 × 10^-8), among which 78 locate more than 1 Mb to the lead chip SNPs were considered as novel, largely augmented the landscape of loci for porcine muscle fatty acid composition. Meta-analysis enhanced the association significance at loci near FADS2, ABCD2, ELOVL5, ELOVL6, ELOVL7, SCD, and THRSP genes, suggesting possible existence of population shared mutations underlying these loci. Further haplotype analysis at SCD loci identified a shared 3.7 kb haplotype in F₂, Sutai and DLY pigs showing consistent effects of decreasing C18:0 contents in the three populations. In contrast, at FASN loci, we found an Erhualian specific haplotype explaining the population specific association signals in Erhualian pigs. This study refines our understanding on landscape of loci and candidate genes for fatty acid composition traits of pigs.

Massive computational identification of somatic variants in exonic splicing enhancers using The Cancer Genome Atlas

Owing to the development of next-generation sequencing (NGS) technologies, a large number of somatic variants have been identified in various types of cancer. However, the functional significance of most somatic variants remains unknown. Somatic variants that occur in exonic splicing enhancer (ESE) regions are thought to prevent serine and arginine-rich (SR) proteins from binding to ESE sequence motifs, which leads to exon skipping. We computationally identified somatic variants in ESEs by compiling numerous open-access datasets from The Cancer Genome Atlas (TCGA). Using somatic variants and RNA-seq data from 9635 patients across 32 TCGA projects, we identified 646 ESE-disrupting variants. The false positive rate of our method, estimated using a permutation test, was approximately 1%. Of these ESE-disrupting variants, approximately 71% were located in the binding motifs of four classical SR proteins. ESE-disrupting variants occurred in proportion to the number of somatic variants, but not necessarily in the specific genes associated with the biological processes of cancer. Existing bioinformatics tools could not predict the pathogenicity of ESE-disrupting variants identified in this study, although these variants could cause exon skipping. We demonstrated that ESE-disrupting nonsense variants tended to escape nonsense-mediated decay surveillance. Using integrated analyses of open access data, we could specifically identify ESE-disrupting variants. We have generated a powerful tool, which can handle datasets without normal samples or raw data, and thus contribute to reducing variants of uncertain significance because our statistical approach only uses the exon-junction read counts from the tumor samples.

Regulatory variants: from detection to predicting impact

Variants within non-coding genomic regions can greatly affect disease. In recent years, increasing focus has been given to these variants, and how they can alter regulatory elements, such as enhancers, transcription factor binding sites and DNA methylation regions. Such variants can be considered regulatory variants. Concurrently, much effort has been put into establishing international consortia to undertake large projects aimed at discovering regulatory elements in different tissues, cell lines and organisms, and probing the effects of genetic variants on regulation by measuring gene expression. Here, we describe methods and techniques for discovering disease-associated non-coding variants using sequencing technologies. We then explain the computational procedures that can be used for annotating these variants using the information from the aforementioned projects, and prediction of their putative effects, including potential pathogenicity, based on rule-based and machine learning approaches. We provide the details of techniques to validate these predictions, by mapping chromatin-chromatin and chromatin-protein interactions, and introduce Clustered Regularly Interspaced Short Palindromic Repeats-Associated Protein 9 (CRISPR-Cas9) technology, which has already been used in this field and is likely to have a big impact on its future evolution. We also give examples of regulatory variants associated with multiple complex diseases. This review is aimed at bioinformaticians interested in the characterization of regulatory variants, molecular biologists and geneticists interested in understanding more about the nature and potential role of such variants from a functional point of views, and clinicians who may wish to learn about variants in non-coding genomic regions associated with a given disease and find out what to do next to uncover how they impact on the underlying mechanisms.

Data Integration in Poplar: 'Omics Layers and Integration Strategies

Populus trichocarpa is an important biofuel feedstock that has been the target of extensive research and is emerging as a model organism for plants, especially woody perennials. This research has generated several large ‘omics datasets. However, only few studies in Populus have attempted to integrate various data types. This review will summarize various ‘omics data layers, focusing on their application in Populus species. Subsequently, network and signal processing techniques for the integration and analysis of these data types will be discussed, with particular reference to examples in Populus.

Rediscover and Refine QTLs for Pig Scrotal Hernia by Increasing a Specially Designed F3 Population and Using Whole-Genome Sequence Imputation Technology

Pig scrotal hernia is one of the most common congenital defects triggered by both genetic and environmental factors, leading to severe economic loss as well as poor animal welfare in the pig industry. Identification and implementation of genomic regions controlling scrotal hernia in breeding is of great appeal to reduce incidences of hernia in pig production. The aim of this study was to identify such regions or molecular markers affecting scrotal hernia in pigs. First of all, we summarized and analyzed the results of some international teams on scrotal hernia and designed a specially population which contains 246 male individuals. We then performed genome-wide association study (GWAS) in this specially designed population using two scenarios, i.e., the target panel data before and after imputation, which contain 42,365 SNPs and 18,756,672 SNPs, respectively. In addition, a series of methods including genetic differentiation analysis, linkage disequilibrium and linkage analysis (LDLA), and haplotype sharing analysis were appropriate to provide for further analysis to identify the potential gene underlying the QTL. The GWAS in this report detected a highly significant region affecting scrotal hernia within a 24.8Mb region (114.1-138.9Mb) on SSC8. And the result of genetic differentiation analysis also showed a strong genetic differentiation signal between 116.1 and 132.7Mb on SSC8. In addition, the QTL interval was refined to 2.99Mb by combining LDLA and genetic differentiation analysis. Finally, two susceptibility haplotypes were identified through haplotype sharing analysis, with one potential causal gene in it. Our study provided deeper insights into the genetic architecture of pig scrotal hernia and contributed to further fine-mapping and characterize haplotype and gene that influence scrotal hernia in pigs.

Scanning of Genetic Variants and Genetic Mapping of Phenotypic Traits in Gilthead Sea Bream Through ddRAD Sequencing

Gilthead sea bream (Sparus aurata) is a teleost of considerable economic importance in Southern European aquaculture. The aquaculture industry shows a growing interest in the application of genetic methods that can locate phenotype-genotype associations with high economic impact. Through selective breeding, the aquaculture industry can exploit this information to maximize the financial yield. Here, we present a Genome Wide Association Study (GWAS) of 112 samples belonging to seven different sea bream families collected from a Greek commercial aquaculture company. Through double digest Random Amplified DNA (ddRAD) Sequencing, we generated a per-sample genetic profile consisting of 2,258 high-quality Single Nucleotide Polymorphisms (SNPs). These profiles were tested for association with four phenotypes of major financial importance: Fat, Weight, Tag Weight, and the Length to Width ratio. We applied two methods of association analysis. The first is the typical single-SNP to phenotype test, and the second is a feature selection (FS) method through two novel algorithms that are employed for the first time in aquaculture genomics and produce groups with multiple SNPs associated to a phenotype. In total, we identified 9 single SNPs and 6 groups of SNPs associated with weight-related phenotypes (Weight and Tag Weight), 2 groups associated with Fat, and 16 groups associated with the Length to Width ratio. Six identified loci (Chr4:23265532, Chr6:12617755, Chr:8:11613979, Chr13:1098152, Chr15:3260819, and Chr22:14483563) were present in genes associated with growth in other teleosts or even mammals, such as semaphorin-3A and neurotrophin-3. These loci are strong candidates for future studies that will help us unveil the genetic mechanisms underlying growth and improve the sea bream aquaculture productivity by providing genomic anchors for selection programs.

Whole-genome re-sequencing association study for direct genetic effects and social genetic effects of six growth traits in Large White pigs

Socially affected traits are affected by direct genetic effects (DGE) and social genetic effects (SGE). DGE and SGE of an individual directly quantify the genetic influence of its own phenotypes and the phenotypes of other individuals, respectively. In the current study, a total of 3,276 Large White pigs from different pens were used, and each pen contained 10 piglets. DGE and SGE were estimated for six socially affected traits, and then a GWAS was conducted to identify SNPs associated with DGE and SGE. Based on the whole-genome re-sequencing, 40 Large White pigs were genotyped and 10,501,384 high quality SNPs were retained for single-locus and multi-locus GWAS. For single-locus GWAS, a total of 54 SNPs associated with DGE and 33 SNPs with SGE exceeded the threshold (P < 5.00E-07) were detected for six growth traits. Of these, 22 SNPs with pleiotropic effects were shared by DGE and SGE. For multi-locus GWAS, a total of 72 and 110 putative QTNs were detected for DGE and SGE, respectively. Of these, 5 SNPs with pleiotropic effects were shared by DGE and SGE. It is noteworthy that 2 SNPs (SSC8: 16438396 for DGE and SSC17: 9697454 for SGE) were detected in single-locus and multi-locus GWAS. Furthermore, 15 positional candidate genes shared by SGE and DGE were identified because of their roles in behaviour, health and disease. Identification of genetic variants and candidate genes for DGE and SGE for socially affected traits will provide a new insight to understand the genetic architecture of socially affected traits in pigs.

Novel protective and risk loci in hip dysplasia in German Shepherds

Canine hip dysplasia is a common, non-congenital, complex and hereditary disorder. It can inflict severe pain via secondary osteoarthritis and lead to euthanasia. An analogous disorder exists in humans. The genetic background of hip dysplasia in both species has remained ambiguous despite rigorous studies. We aimed to investigate the genetic causes of this disorder in one of the high-risk breeds, the German Shepherd. We performed genetic analyses with carefully phenotyped case-control cohorts comprising 525 German Shepherds. In our genome-wide association studies we identified four suggestive loci on chromosomes 1 and 9. Targeted resequencing of the two loci on chromosome 9 from 24 affected and 24 control German Shepherds revealed deletions of variable sizes in a putative enhancer element of the NOG gene. NOG encodes for noggin, a well-described bone morphogenetic protein inhibitor affecting multiple developmental processes, including joint development. The deletion was associated with the healthy controls and mildly dysplastic dogs suggesting a protective role against canine hip dysplasia. Two enhancer variants displayed a decreased activity in a dual luciferase reporter assay. Our study identifies novel loci and candidate genes for canine hip dysplasia, with potential regulatory variants in the NOG gene. Further research is warranted to elucidate how the identified variants affect the expression of noggin in canine hips, and what the potential effects of the other identified loci are.

Using imputed whole-genome sequence data to improve the accuracy of genomic prediction for parasite resistance in Australian sheep

Background

This study aimed at (1) comparing the accuracies of genomic prediction for parasite resistance in sheep based on whole-genome sequence (WGS) data to those based on 50k and high-density (HD) single nucleotide polymorphism (SNP) panels; (2) investigating whether the use of variants within quantitative trait loci (QTL) regions that were selected from regional heritability mapping (RHM) in an independent dataset improved the accuracy more than variants selected from genome-wide association studies (GWAS); and (3) comparing the prediction accuracies between variants selected from WGS data to variants selected from the HD SNP panel.

Results

The accuracy of genomic prediction improved marginally from 0.16 ± 0.02 and 0.18 ± 0.01 when using all the variants from 50k and HD genotypes, respectively, to 0.19 ± 0.01 when using all the variants from WGS data. Fitting a GRM from the selected variants alongside a GRM from the 50k SNP genotypes improved the prediction accuracy substantially compared to fitting the 50k SNP genotypes alone. The gain in prediction accuracy was slightly more pronounced when variants were selected from WGS data compared to when variants were selected from the HD panel. When sequence variants that passed the GWAS \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$- log_{10} (p\,value)$$\end{document}-log10(pvalue) threshold of 3 across the entire genome were selected, the prediction accuracy improved by 5% (up to 0.21 ± 0.01), whereas when selection was limited to sequence variants that passed the same GWAS \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$- log_{10} (p\,value)$$\end{document}-log10(pvalue) threshold of 3 in regions identified by RHM, the accuracy improved by 9% (up to 0.25 ± 0.01).

Conclusions

Our results show that through careful selection of sequence variants from the QTL regions, the accuracy of genomic prediction for parasite resistance in sheep can be improved. These findings have important implications for genomic prediction in sheep.

Rare genetic variation and balanced polymorphisms are important for survival in global change conditions

Standing genetic variation is important for population persistence in extreme environmental conditions. While some species may have the capacity to adapt to predicted average future global change conditions, the ability to survive extreme events is largely unknown. We used single-generation selection experiments on hundreds of thousands of Strongylocentrotus purpuratus sea urchin larvae generated from wild-caught adults to identify adaptive genetic variation responsive to moderate (pH 8.0) and extreme (pH 7.5) low-pH conditions. Sequencing genomic DNA from pools of larvae, we identified consistent changes in allele frequencies across replicate cultures for each pH condition and observed increased linkage disequilibrium around selected loci, revealing selection on recombined standing genetic variation. We found that loci responding uniquely to either selection regime were at low starting allele frequencies while variants that responded to both pH conditions (11.6% of selected variants) started at high frequencies. Loci under selection performed functions related to energetics, pH tolerance, cell growth and actin/cytoskeleton dynamics. These results highlight that persistence in future conditions will require two classes of genetic variation: common, pH-responsive variants maintained by balancing selection in a heterogeneous environment, and rare variants, particularly for extreme conditions, that must be maintained by large population sizes.

Association of candidate genes with drought tolerance traits in zoysiagrass germplasm

Drought stress negatively influences the growth and physiology of perennial grasses. The objective of this study was to identify associations of candidate genes with drought tolerance traits in 96 zoysiagrass (Zoysia Willd.) accessions. Germplasm varied largely in leaf wilting, canopy and air temperature difference (CAD), leaf water content (LWC), chlorophyll fluorescence (Fv/Fm), leaf dry weight (LDW), stolon dry weight (SDW), rhizome dry weight (RZW), and root dry weight (RDW) under drought stress across the two experiments in 2014 and 2015 in a greenhouse. The population exhibited three subgroups based on molecular marker analysis and had minimum relative kinship. Associations between single nucleotide polymorphisms (SNPs) in BADH encoding betaine aldehyde dehydrogenase, DREB1 encoding DREB-like protein 1, Ndhf encoding NADH dehydrogenase subunit F, CAT encoding catalase, and VP1 encoding H⁺-pyrophosphatase were analyzed with trait under drought stress (D) and relative values compared to the control (R). Twenty-seven mark and trait associations were detected in year 2014, 2015, and a two-year combination across four genes, including 13 associations in 7 SNP loci in BADH, 9 associations in 5 SNP loci in DREB1, 3 associations in one SNP locus in Ndhf, and 2 associations in one SNP locus in CAT. Of them, one SNP in BADH was associated with D-RDW or D-SDW, three SNPs in DREB1 were associated with D-RZW, D-RDW, R-LWC, and D-CAD, and one SNP in CAT was associated with D-SDW. Nucleotide changes in these SNP loci caused non-synonymous amino acid substitutions. The results indicated that allelic diversity in genes involved in antioxidant metabolism, osmotic homeostasis, and dehydration responsive transcription factor could contribute to growth and physiological variations in zoysiagrass under drought stress.

Genetic influence on scar height and pliability after burn injury in individuals of European ancestry: A prospective cohort study

After similar extent of injury there is considerable variability in scarring between individuals, in part due to genetic factors. This study aimed to identify genetic variants associated with scar height and pliability after burn injury. An exome-wide array association study and gene pathway analysis were performed on a prospective cohort of 665 patients treated for burn injury. Outcomes were scar height (SH) and scar pliability (SP) sub-scores of the modified Vancouver Scar Scale (mVSS). DNA was genotyped using the Infinium^® HumanCoreExome-24 BeadChip. Associations between genetic variants (single nucleotide polymorphisms) and SH and SP were estimated using an additive genetic model adjusting for age, sex, number of surgical procedures and % total body surface area of burn in subjects of European ancestry. No individual genetic variants achieved the cut-off threshold of significance. Gene regions were analysed for spatially correlated single nucleotide polymorphisms and significant regions identified using comb-p software. This gene list was subject to gene pathway analysis to find which biological process terms were over-represented. Using this approach biological processes related to the nervous system and cell adhesion were the predominant gene pathways associated with both SH and SP. This study suggests genes associated with innervation may be important in scar fibrosis. Further studies using similar and larger datasets will be essential to validate these findings.

Genome-Wide Association Studies Reveal Genomic Regions Associated With the Response of Wheat (Triticum aestivum L.) to Mycorrhizae Under Drought Stress Conditions

In the majority of wheat growing areas worldwide, the incidence of drought stress has increased significantly resulting in a negative impact on plant development and grain yield. Arbuscular mycorrhizal symbiosis is known to improve drought stress tolerance of wheat. However, quantitative trait loci (QTL) involved in the response to drought stress conditions in the presence of mycorrhizae are largely unknown. Therefore, a diverse set consisting of 94 bread wheat genotypes was phenotyped under drought stress and well watered conditions in the presence and absence of mycorrhizae. Grain yield and yield components, drought stress related traits as well as response to mycorrhizae were assessed. In parallel, wheat accessions were genotyped by using the 90k iSelect chip, resulting in a set of 15511 polymorphic and mapped SNP markers, which were used for genome-wide association studies (GWAS). In general, drought stress tolerance of wheat was significantly increased in the presence of mycorrhizae compared to drought stress tolerance in the absence of mycorrhizae. However, genotypes differed in their response to mycorrhizae under drought stress conditions. Several QTL regions on different chromosomes were detected associated with grain yield and yield components under drought stress conditions. Furthermore, two genome regions on chromosomes 3D and 7D were found to be significantly associated with the response to mycorrhizae under drought stress conditions. Overall, the results reveal that inoculation of wheat with mycorrhizal fungi significantly improves drought stress tolerance and that QTL regions associated with the response to mycorrhizae under drought stress conditions exist in wheat. Further research is necessary to validate detected QTL regions. However, this study may be the starting point for the identification of candidate genes associated with drought stress tolerance and response to mycorrhizae under drought stress conditions. Maybe in future, these initial results will help to contribute to use mycorrhizal fungi effectively in agriculture and combine new approaches i.e., use of genotypic variation in response to mycorrhizae under drought stress conditions with existing drought tolerance breeding programs to develop new drought stress tolerant genotypes.

Genome wide association study on feed conversion ratio using imputed sequence data in chickens

Objective

Feed consumption contributes a large percentage for total production costs in the poultry industry. Detecting genes associated with feeding traits will be of benefit to improve our understanding of the molecular determinants for feed efficiency. The objective of this study was to identify candidate genes associated with feed conversion ratio (FCR) via genome-wide association study (GWAS) using sequence data imputed from single nucleotide polymorphism (SNP) panel in a Chinese indigenous chicken population.

Methods

A total of 435 Chinese indigenous chickens were phenotyped for FCR and were genotyped using a 600K SNP genotyping array. Twenty-four birds were selected for sequencing, and the 600K SNP panel data were imputed to whole sequence data with the 24 birds as the reference. The GWAS were performed with GEMMA software.

Results

After quality control, 8,626,020 SNPs were used for sequence based GWAS, in which ten significant genomic regions were detected to be associated with FCR. Ten candidate genes, ubiquitin specific peptidase 44, leukotriene A4 hydrolase, ETS transcription factor, R-spondin 2, inhibitor of apoptosis protein 3, sosondowah ankyrin repeat domain family member D, calmodulin regulated spectrin associated protein family member 2, zinc finger and BTB domain containing 41, potassium sodium-activated channel subfamily T member 2, and member of RAS oncogene family were annotated. Several of them were within or near the reported FCR quantitative trait loci, and others were newly reported.

Conclusion

Results from this study provide valuable prior information on chicken genomic breeding programs, and potentially improve our understanding of the molecular mechanism for feeding traits.

Imputation-Based Whole-Genome Sequence Association Study Reveals Constant and Novel Loci for Hematological Traits in a Large-Scale Swine F2 Resource Population

The whole-genome sequences of progenies with low-density single-nucleotide polymorphism (SNP) genotypes can be imputed with high accuracy based on the deep-coverage sequences of key ancestors. With this imputation technology, a more powerful genome-wide association study (GWAS) can be carried out using imputed whole-genome variants and the phenotypes of interest to overcome the shortcomings of low-power detection and the large confidence interval derived from low-density SNP markers in classic association studies. In this study, 19 ancestors of a large-scale swine F₂ White Duroc × Erhualian population were deeply sequenced for their genome with an average coverage of 25×. Considering 98 pigs from 10 different breeds with high-quality deep sequenced genomes, we imputed the whole genomic variants of 1020 F₂ pigs genotyped by the PorcineSNP60 BeadChip with high accuracy and obtained 14,851,440 sequence variants after quality control. Based on this, 87 novel quantitative traits loci (QTLs) for 18 hematological traits at three different physiological stages of the F₂ pigs were identified, among which most of the novel QTLs have been repeated in two of the three stages. Literature mining pinpointed that the FGF14 and LCLAT1 genes at SSC11 and SSC3 may affect the MCH at day 240 and MCV at day 18, respectively. The present study shows that combining high-quality imputed genomic variants and correlated phenomic traits into GWAS can improve the capability to detect QTL considerably. The large number of different QTLs for hematological traits identified at multiple growth stages implies the complexity and time specificity of these traits.

Linkage mapping combined with association analysis reveals QTL and candidate genes for three husk traits in maize

Key message Combined linkage and association mapping analyses facilitate the emphasis on the candidate genes putatively involved in maize husk growth. The maize (Zea mays L.) husk consists of multiple leafy layers and plays important roles in protecting the ear from pathogen infection and in preventing grain dehydration. Although husk morphology varies widely among different maize inbred lines, the genetic basis of such variation is poorly understood. In this study, we used three maize recombinant inbred line (RIL) populations to dissect the genetic basis of three husk traits: i.e., husk length (HL), husk width (HW), and the number of husk layers (HN). Three husk traits in all three RIL populations showed wide phenotypic variation and high heritability. The HL showed stronger correlations with ear traits than did HW and HN. A total of 21 quantitative trait loci (QTL) were identified for the three traits in three RIL populations, and some of them were commonly observed for the same trait in different populations. The proportions of total phenotypic variation explained by QTL in three RIL populations were 31.8, 35.3, and 44.5% for HL, HW, and HN, respectively. The highest proportions of phenotypic variation explained by a single QTL were 14.7% for HL in the By815/K22 RIL population (BYK), 13.5% for HW in the By815/DE3 RIL population (BYD), and 19.4% for HN in the BYD population. A combined analysis of linkage mapping with a previous genome-wide association study revealed five candidate genes related to husk morphology situated within three QTL loci. These five genes were related to metabolism, gene expression regulation, and signal transduction.

Genome-Scale Association Study of Abnormal Scale Pattern in Yellow River Carp Identified Previously Known Causative Gene in European Mirror Carp

Common carp (Cyprinus carpio) is one of the most widely studied fish species due to its great economic value and strong environmental adaptability. Scattered scale, a typical phenotype of the mirror carp that is derived from Europe, has never been observed in the Yellow River carp previously. We recently identified approximately one fourth of the F1 progenies displaying scattered scale in a full-sib Yellow River carp family in our breeding program, despite both parents that showed wild type with normal scale patterns. This family provides us unique materials to investigate the genetic basis underlying the abnormal scale mutant in Yellow River carp population. Genome-wide association study (GWAS) and association mapping were performed based on genome-wide single nucleotide polymorphisms (SNP) genotyped with common carp 250 K SNP genotyping array in 82 samples of the Yellow River carp family. We identified a 1.4 Mb genome region that was significantly associated with abnormal scattered scale patterns. We further identified a deletion mutation in fibroblast growth factor receptor 1 a1 (fgfr1a1) gene within this genome region. Amplification and sequencing analysis of this gene revealed a 311-bp deletion in intron 10 and exon 11, which proved that fgfr1a1 could be the causal gene responsible for abnormal scattered scale in the Yellow River carp family. Since similar fragment mutation with 306-bp and 310-bp deletions had been previously reported as causal mutation of scattered scale patterns in the mirror carp, we speculate that either the deletion mutation was introduced from Europe-derived mirror carp or the deletion independently occurred in the mutation hotspot in fgfr1a1 gene. The results provided insights into the genetic basis of scale pattern mutant in Yellow River carp population, which would help us to eliminate the recessive allele of the abnormal scale patterns in Yellow River carp population by molecular marker-assisted breeding.

The Influence of Genetic Variation on Social Disposition, Romantic Relationships and Social Networks: a Replication Study

OBJECTIVES

Sociality is underpinned by a variety of neurochemicals. We previously showed, in a large healthy Caucasian sample, that genes for different neurochemicals are typically associated with differing social domains (disposition, romantic relationships and networks). Here we seek to confirm the validity of these findings by asking whether they replicate in other population samples.

METHODS

We test for associations between the same 24 Single Nucleotide Polymorphisms (SNPs) and measures of sociality as previously, in two smaller independent samples: Caucasian individuals with histories of mental illness (subclinical sample) (N = 140), and non-Caucasian individuals (N = 66). We also combined the relevant SNPs and social measures into 18 distinct neurochemical/social domain categories to examine the distribution of significant associations across these.

RESULTS

In the subclinical Caucasian sample, we confirm previous associations between (i) specific oxytocin and dopamine receptor gene SNPs and sexual attitudes and behavior, and (ii) two SNPs associated with dopamine receptor 2 and feelings of inclusion in the local community. In the non-Caucasian sample, we replicate the previous association between an oxytocin receptor SNP and anxious attachment. More generally, chi-squared tests indicated that the distribution of significant associations for each neurochemical across the three social domains did not differ significantly between the original sample and either of the new samples, except for oxytocin in the non-Caucasian sample.

CONCLUSIONS

These results corroborate both the SNP-specific and broader neurochemical associations with particular facets of sociality in two new populations, thereby confirming the validity of the previous findings.