Family-based association tests for genomewide association scans

Exploring Genetic Factors Associated with Moniezia spp. Tapeworm Resistance in Central Anatolian Merino Sheep via GWAS Approach

Gastrointestinal parasite (GIP) infections pose significant challenges in pasture-based sheep farming, leading to economic losses and welfare concerns. This study aimed to uncover the genetic basis of resistance to Moniezia spp. infections in Central Anatolian Merino (CAM) sheep. Genome-Wide Association Analysis (GWAS) was conducted between Moniezia spp. egg burden and genomic data from 226 CAM lambs. Thirteen significant Single-Nucleotide Polymorphisms (SNPs) were identified, with five surpassing the genome-wide threshold and eight exceeding the chromosome-wide threshold. Functional annotation revealed associations with genes involved in immune function, notably CD79A and MAP3K7. CD79A, integral to B-cell activation and antibody production, plays a key role in the immune response against parasitic infections. Its interaction with helminth-derived proteins modulates B-cell function, highlighting its potential as a therapeutic target. MAP3K7, a central regulator of immune signaling pathways, modulates host responses to helminth infections by influencing NF-κB activity. Additionally, it regulates macrophage function in bacterial infections, showcasing its versatility in mediating immune responses against diverse pathogens. From a practical perspective, the findings of the current research underscore the potential of integrating genomic information into breeding programs to bolster disease resilience in livestock populations for sustainable production purposes. However, further research is needed to elucidate the functional significance of identified SNPs and associated genes. This study underscores the potential of genomic approaches in combating parasitic diseases and promoting sustainable agriculture in sheep production systems.

A Litmus Test for Confounding in Polygenic Scores

Polygenic scores (PGSs) are being rapidly adopted for trait prediction in the clinic and beyond. PGSs are often thought of as capturing the direct genetic effect of one's genotype on their phenotype. However, because PGSs are constructed from population-level associations, they are influenced by factors other than direct genetic effects, including stratification, assortative mating, and dynastic effects ("SAD effects"). Our interpretation and application of PGSs may hinge on the relative impact of SAD effects, since they may often be environmentally or culturally mediated. We developed a method that estimates the proportion of variance in a PGS (in a given sample) that is driven by direct effects, SAD effects, and their covariance. We leverage a comparison of a PGS of interest based on a standard GWAS with a PGS based on a sibling GWAS-which is largely immune to SAD effects-to quantify the relative contribution of each type of effect to variance in the PGS of interest. Our method, Partitioning Genetic Scores Using Siblings (PGSUS, pron. "Pegasus"), breaks down variance components further by axes of genetic ancestry, allowing for a nuanced interpretation of SAD effects. In particular, PGSUS can detect stratification along major axes of ancestry as well as SAD variance that is "isotropic" with respect to axes of ancestry. Applying PGSUS, we found evidence of stratification in PGSs constructed using large meta-analyses of height and educational attainment as well as in a range of PGSs constructed using the UK Biobank. In some instances, a given PGS appears to be stratified along a major axis of ancestry in one prediction sample but not in another (for example, in comparisons of prediction in samples from different countries, or in ancient DNA vs. contemporary samples). Finally, we show that different approaches for adjustment for population structure in GWASs have distinct advantages with respect to mitigation of ancestry-axis-specific and isotropic SAD variance in PGS. Our study illustrates how family-based designs can be combined with standard population-based designs to guide the interpretation and application of genomic predictors.

Unravelling the Genetic Architecture of Serum Biochemical Indicators in Sheep

Serum biochemical indicators serve as vital proxies that reflect the physiological state and functions of different organs. The genetic parameters and molecular mechanisms underlying serum biochemical indicators of sheep (Ovis aries) have not been well understood. Therefore, the aim of the present study was to identify the genetic architecture and genomic loci underlying ten serum biochemical indicators in sheep, including alanine transaminase, aspartate transferase, lactate dehydrogenase, cholesterol, glucose, phosphorus, calcium, creatinine, urea and total protein levels. We implemented genetic parameter estimations and GWASs for each trait in 422 Akkaraman lambs. Overall, low to moderate heritability estimates were found in the range of 0.14-0.55. Additionally, low to high genetic correlations were observed among traits. In total, 23 SNP loci were associated with serum biochemical indicators leading to 19 genes. These were SPTA1, MGST2, CACUL1, IGFBP7, PARD3, PHB1, SLC15A5, TRIM35, RGS6, NUP93, CNTNAP2, SLC7A11, B3GALT5, DPP10, HST2ST1, NRP1, LRP1B, MAP3K9 and ENSOARG00020040484.1, as well as LOC101103187, LOC101117162, LOC105611309 and LOC101118029. To our knowledge, these data provide the first associations between SPTA1 and serum cholesterol and between ENSOARG00020040484.1 and serum glucose. The current findings provide a comprehensive inventory of the relationships between serum biochemical parameters, genetic variants and disease-relevant characteristics. This information may facilitate the identification of therapeutic targets and fluid biomarkers and establish a strong framework for comprehending the pathobiology of complex diseases as well as providing targets for sheep genetic improvement programs.

Genome-wide association analysis uncovers rice blast resistance alleles of Ptr and Pia

A critical step to maximize the usefulness of genome-wide association studies (GWAS) in plant breeding is the identification and validation of candidate genes underlying genetic associations. This is of particular importance in disease resistance breeding where allelic variants of resistance genes often confer resistance to distinct populations, or races, of a pathogen. Here, we perform a genome-wide association analysis of rice blast resistance in 500 genetically diverse rice accessions. To facilitate candidate gene identification, we produce de-novo genome assemblies of ten rice accessions with various rice blast resistance associations. These genome assemblies facilitate the identification and functional validation of novel alleles of the rice blast resistance genes Ptr and Pia. We uncover an allelic series for the unusual Ptr rice blast resistance gene, and additional alleles of the Pia resistance genes RGA4 and RGA5. By linking these associations to three thousand rice genomes we provide a useful tool to inform future rice blast breeding efforts. Our work shows that GWAS in combination with whole-genome sequencing is a powerful tool for gene cloning and to facilitate selection of specific resistance alleles for plant breeding.

Genetic regions determine tolerance to nitrogen deficiency in European elite bread wheats grown under contrasting nitrogen stress scenarios

KEY MESSAGE

Clustering 24 environments in four contrasting nitrogen stress scenarios enabled the detection of genetic regions determining tolerance to nitrogen deficiency in European elite bread wheats. Increasing the nitrogen use efficiency of wheat varieties is an important goal for breeding. However, most genetic studies of wheat grown at different nitrogen levels in the field report significant interactions with the genotype. The chromosomal regions possibly involved in these interactions are largely unknown. The objective of this study was to quantify the response of elite bread wheat cultivars to different nitrogen field stress scenarios and identify genomic regions involved in this response. For this purpose, 212 elite bread wheat varieties were grown in a multi-environment trial at different nitrogen levels. Genomic regions associated with grain yield, protein concentration and grain protein deviation responses to nitrogen deficiency were identified. Environments were clustered according to adjusted means for grain yield, yield components and grain protein concentration. Four nitrogen availability scenarios were identified: optimal condition, moderate early deficiency, severe late deficiency, and severe continuous deficiency. A large range of tolerance to nitrogen deficiency was observed among varieties, which were ranked differently in different nitrogen deficiency scenarios. The well-known negative correlation between grain yield and grain protein concentration also existed between their respective tolerance indices. Interestingly, the tolerance indices for grain yield and grain protein deviation were either null or weakly positive meaning that breeding for the two traits should be less difficult than expected. Twenty-two QTL regions were identified for the tolerance indices. By selecting associated markers, these regions may be selected separately or combined to improve the tolerance to N deficiency within a breeding programme.

Integrated epigenome, whole genome sequence and metabolome analyses identify novel multi-omics pathways in type 2 diabetes: a Middle Eastern study

BACKGROUND

T2D is of high prevalence in the middle east and thus studying its mechanisms is of a significant importance. Using 1026 Qatar BioBank samples, epigenetics, whole genome sequencing and metabolomics were combined to further elucidate the biological mechanisms of T2D in a population with a high prevalence of T2D.

METHODS

An epigenome-wide association study (EWAS) with T2D was performed using the Infinium 850K EPIC array, followed by whole genome-wide sequencing SNP-CpG association analysis (> 5.5 million SNPs) and a methylome-metabolome (CpG-metabolite) analysis of the identified T2D sites.

RESULTS

A total of 66 T2D-CpG associations were identified, including 63 novel sites in pathways of fructose and mannose metabolism, insulin signaling, galactose, starch and sucrose metabolism, and carbohydrate absorption and digestion. Whole genome SNP associations with the 66 CpGs resulted in 688 significant CpG-SNP associations comprising 22 unique CpGs (33% of the 66 CPGs) and included 181 novel pairs or pairs in novel loci. Fourteen of the loci overlapped published GWAS loci for diabetes related traits and were used to identify causal associations of HK1 and PFKFB2 with HbA1c. Methylome-metabolome analysis identified 66 significant CpG-metabolite pairs among which 61 pairs were novel. Using the identified methylome-metabolome associations, methylation QTLs, and metabolic networks, a multi-omics network was constructed which suggested a number of metabolic mechanisms underlying T2D methylated genes. 1-palmitoyl-2-oleoyl-GPE (16:0/18:1) - a triglyceride-associated metabolite, shared a common network with 13 methylated CpGs, including TXNIP, PFKFB2, OCIAD1, and BLCAP. Mannonate - a food component/plant shared a common network with 6 methylated genes, including TXNIP, BLCAP, THBS4 and PEF1, pointing to a common possible cause of methylation in those genes. A subnetwork with alanine, glutamine, urea cycle (citrulline, arginine), and 1-carboxyethylvaline linked to PFKFB2 and TXNIP revealed associations with kidney function, hypertension and triglyceride metabolism. The pathway containing STYXL1-POR was associated with a sphingosine-ceramides subnetwork associated with HDL-C and LDL-C and point to steroid perturbations in T2D.

CONCLUSIONS

This study revealed several novel methylated genes in T2D, with their genomic variants and associated metabolic pathways with several implications for future clinical use of multi-omics associations in disease and for studying therapeutic targets.

Symbiotic Variations among Wheat Genotypes and Detection of Quantitative Trait Loci for Molecular Interaction with Auxin-Producing Azospirillum PGPR

Crop varieties differ in their ability to interact with Plant Growth-Promoting Rhizobacteria (PGPR), but the genetic basis for these differences is unknown. This issue was addressed with the PGPR Azospirillum baldaniorum Sp245, using 187 wheat accessions. We screened the accessions based on the seedling colonization by the PGPR and the expression of the phenylpyruvate decarboxylase gene ppdC (for synthesis of the auxin indole-3-acetic acid), using gusA fusions. Then, the effects of the PGPR on the selected accessions stimulating Sp245 (or not) were compared in soil under stress. Finally, a genome-wide association approach was implemented to identify the quantitative trait loci (QTL) associated with PGPR interaction. Overall, the ancient genotypes were more effective than the modern genotypes for Azospirillum root colonization and ppdC expression. In non-sterile soil, A. baldaniorum Sp245 improved wheat performance for three of the four PGPR-stimulating genotypes and none of the four non-PGPR-stimulating genotypes. The genome-wide association did not identify any region for root colonization but revealed 22 regions spread on 11 wheat chromosomes for ppdC expression and/or ppdC induction rate. This is the first QTL study focusing on molecular interaction with PGPR bacteria. The molecular markers identified provide the possibility to improve the capacity of modern wheat genotypes to interact with Sp245, as well as, potentially, other Azospirillum strains.

Simultaneous selection of multiple important single nucleotide polymorphisms in familial genome wide association studies data

We propose a resampling-based fast variable selection technique for detecting relevant single nucleotide polymorphisms (SNP) in a multi-marker mixed effect model. Due to computational complexity, current practice primarily involves testing the effect of one SNP at a time, commonly termed as 'single SNP association analysis'. Joint modeling of genetic variants within a gene or pathway may have better power to detect associated genetic variants, especially the ones with weak effects. In this paper, we propose a computationally efficient model selection approach-based on the e-values framework-for single SNP detection in families while utilizing information on multiple SNPs simultaneously. To overcome computational bottleneck of traditional model selection methods, our method trains one single model, and utilizes a fast and scalable bootstrap procedure. We illustrate through numerical studies that our proposed method is more effective in detecting SNPs associated with a trait than either single-marker analysis using family data or model selection methods that ignore the familial dependency structure. Further, we perform gene-level analysis in Minnesota Center for Twin and Family Research (MCTFR) dataset using our method to detect several SNPs using this that have been implicated to be associated with alcohol consumption.

Genome-Wide Scan of Wool Production Traits in Akkaraman Sheep

The objective of this study was to uncover the genetic background of wool quality, a production trait, by estimating genomic heritability and implementing GWAS in Akkaraman sheep. The wool characteristics measured included fibre diameter (FD) and staple length (SL) at the age of 8 months and yearling fibre diameter (YFD), yearling staple length (YSL) and yearling greasy fleece weight (YGFW) at 18 months of age. Animals were genotyped using the Axiom 50 K Ovine Genotyping Array. Maximum likelihood estimations of a linear mixed model (LMM) were used to estimate genomic heritability, where GWAS was conducted following a score test of each trait. Genomic heritability estimates for the traits ranged between 0.22 and 0.63, indicating that phenotypes have a moderate range of heritability. One genome- and six chromosome-wide significant SNPs were associated with the wool traits in Akkaraman lambs. Accordingly, TRIM2, MND1, TLR2, RNF175, CEP290, TMTC3, RERE, SLC45A1, SOX2, MORN1, SKI, FAAP20, PRKCZ, GABRD, CFAP74, CALML6 and TMEM52 genes as well as nine uncharacterized regions (LOC101118971, LOC105609137, LOC105603067, LOC101122892, LOC106991694, LOC106991467, LOC106991455, LOC105616534 and LOC105609719) were defined as plausible candidates. The findings of this study shed light on the genetics of wool quality and yield for the Akkaraman breed and suggests targets for breeders during systematic breeding programmes.

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.

Estimates of genomic heritability and genome-wide association studies for blood parameters in Akkaraman sheep

The aim of this study was to estimate genomic heritability and the impact that genetic backgrounds have on blood parameters in Akkaraman sheep by conducting genome-wide association studies and regional heritability mapping analysis. Genomic heritability estimates for blood parameters ranged from 0.00 to 0.55, indicating that measured phenotypes have a low to moderate heritability. A total of 7 genome- and 13 chromosome-wide significant SNPs were associated with phenotypic changes in 15 blood parameters tested. Accordingly, SCN7A, SCN9A, MYADM-like, CCDC67, ITGA9, MGAT5, SLC19A1, AMPH, NTRK2, MSRA, SLC35F3, SIRT6, CREB3L3, and NAV3 genes as well as three undefined regions (LOC101117887, LOC106991526 and LOC105608461) were suggested as candidates. Most of the identified genes were involved in basic biological processes that are essential to immune system function and cellular growth; specific functions include cellular transport, histone deacetylation, cell differentiation, erythropoiesis, and endocytosis. The top significant SNP for HCT, MCH, and MCHC was found within a genomic region mainly populated by the MYADM-like gene family. This region was previously suggested to be under historical selection pressure in many sheep breeds from various parts of the world. These results have implications on animal breeding program studies due to the effect that the genetic background has on blood parameters, which underlying many productive and wellness related traits.

Genetic Parameters and Genomic Regions Underlying Growth and Linear Type Traits in Akkaraman Sheep

In the current study, the genetic architecture of growth and linear type traits were investigated in Akkaraman sheep. Estimations of genomic heritability, genetic correlations, and phenotypic correlations were implemented for 17 growth and linear type traits of 473 Akkaraman lambs by the univariate and multivariate analysis of animal mixed models. Correspondingly, moderate heritability estimates, as well as high and positive genetic/phenotypic correlations were found between growth and type traits. On the other hand, 2 genome-wide and 19 chromosome-wide significant single nucleotide polymorphisms were found to be associated with the traits as a result of animal mixed model-based genome-wide association analyses. Accordingly, we propose several genes located on different chromosomes (e.g., PRDM2, PTGDR, PTPRG, KCND2, ZNF260, CPE, GRID2, SCD5, SPIDR, ZNF407, HCN3, TMEM50A, FKBP1A, TLE4, SP1, SLC44A1, and MYOM3) as putative quantitative trait loci for the 22 growth and linear type traits studied. In our study, specific genes (e.g., TLE4, PTGDR, and SCD5) were found common between the traits studied, suggesting an interplay between the genetic backgrounds of these traits. The fact that four of the proposed genes (TLE4, MYOM3, SLC44A1, and TMEM50A) are located on sheep chromosome 2 confirms the importance of these genomic regions for growth and morphological structure in sheep. The results of our study are therefore of great importance for the development of efficient selection indices and marker-assisted selection programs, as well as for the understanding of the genetic architecture of growth and linear traits in sheep.

Benchmarking statistical methods for analyzing parent-child dyads in genetic association studies

Genetic association studies of child health outcomes often employ family-based study designs. One of the most popular family-based designs is the case-parent trio design that considers the smallest possible nuclear family consisting of two parents and their affected child. This trio design is particularly advantageous for studying relatively rare disorders because it is less prone to type 1 error inflation due to population stratification compared to population-based study designs (e.g., case-control studies). However, obtaining genetic data from both parents is difficult, from a practical perspective, and many large studies predominantly measure genetic variants in mother-child dyads. While some statistical methods for analyzing parent-child dyad data (most commonly involving mother-child pairs) exist, it is not clear if they provide the same advantage as trio methods in protecting against population stratification, or if a specific dyad design (e.g., case-mother dyads vs. case-mother/control-mother dyads) is more advantageous. In this article, we review existing statistical methods for analyzing genome-wide marker data on dyads and perform extensive simulation experiments to benchmark their type I errors and statistical power under different scenarios. We extend our evaluation to existing methods for analyzing a combination of case-parent trios and dyads together. We apply these methods on genotyped and imputed data from multiethnic mother-child pairs only, case-parent trios only or combinations of both dyads and trios from the Gene, Environment Association Studies consortium (GENEVA), where each family was ascertained through a child affected by nonsyndromic cleft lip with or without cleft palate. Results from the GENEVA study corroborate the findings from our simulation experiments. Finally, we provide recommendations for using statistical genetic association methods for dyads.

A robust biostatistical method leverages informative but uncertainly determined qPCR data for biomarker detection, early diagnosis, and treatment

As a common medium-throughput technique, qPCR (quantitative real-time polymerase chain reaction) is widely used to measure levels of nucleic acids. In addition to accurate and complete data, experimenters have unavoidably observed some incomplete and uncertainly determined qPCR data because of intrinsically low overall amounts of biological materials, such as nucleic acids present in biofluids. When there are samples with uncertainly determined qPCR data, some investigators apply the statistical complete-case method by excluding the subset of samples with uncertainly determined data from analysis (CO), while others simply choose not to analyze (CNA) these datasets altogether. To include as many observations as possible in analysis for interesting differential changes between groups, some investigators set incomplete observations equal to the maximum quality qPCR cycle (MC), such as 32 and 40. Although straightforward, these methods may decrease the sample size, skew the data distribution, and compromise statistical power and research reproducibility across replicate qPCR studies. To overcome the shortcomings of the existing, commonly-used qPCR data analysis methods and to join the efforts in advancing statistical analysis in rigorous preclinical research, we propose a robust nonparametric statistical cycle-to-threshold method (CTOT) to analyze incomplete qPCR data for two-group comparisons. CTOT incorporates important characteristics of qPCR data and time-to-event statistical methodology, resulting in a novel analytical method for qPCR data that is built around good quality data from all subjects, certainly determined or not. Considering the benchmark full data (BFD), we compared the abilities of CTOT, CO, MC, and CNA statistical methods to detect interesting differential changes between groups with informative but uncertainly determined qPCR data. Our simulations and applications show that CTOT improves the power of detecting and confirming differential changes in many situations over the three commonly used methods without excess type I errors. The robust nonparametric statistical method of CTOT helps leverage qPCR technology and increase the power to detect differential changes that may assist decision making with respect to biomarker detection and early diagnosis, with the goal of improving the management of patient healthcare.

Controlling for background genetic effects using polygenic scores improves the power of genome-wide association studies

Ongoing increases in the size of human genotype and phenotype collections offer the promise of improved understanding of the genetics of complex diseases. In addition to the biological insights that can be gained from the nature of the variants that contribute to the genetic component of complex trait variability, these data bring forward the prospect of predicting complex traits and the risk of complex genetic diseases from genotype data. Here we show that advances in phenotype prediction can be applied to improve the power of genome-wide association studies. We demonstrate a simple and efficient method to model genetic background effects using polygenic scores derived from SNPs that are not on the same chromosome as the target SNP. Using simulated and real data we found that this can result in a substantial increase in the number of variants passing genome-wide significance thresholds. This increase in power to detect trait-associated variants also translates into an increase in the accuracy with which the resulting polygenic score predicts the phenotype from genotype data. Our results suggest that advances in methods for phenotype prediction can be exploited to improve the control of background genetic effects, leading to more accurate GWAS results and further improvements in phenotype prediction.

Genome-wide association studies of preweaning growth and in vivo carcass composition traits in Esme sheep

Sheep are considered as a major contributor of global food security. Moreover, sheep preweaning growth traits as well as in vivo carcass composition traits such as ultrasonic measurements of Longissimus dorsi muscle depth (UMD) and back-fat thickness (UFD) are crucially important indicators of meat yield and hot carcass composition. Despite their relative importance for productivity and profitability of a sheep production system, detected QTL for these traits are quite scarce. Therefore, we implemented GWAS for these traits using animal mixed model-based association approach provided by GenABEL in Esme sheep. Three genome-wide and 14 individual chromosome-wide associated SNPs were discovered. As a result, ESRP1, LOC105613082, ZNF641, DUSP5, TEAD1, SMOX, PTPRT, RALYL, POM121C, PHIP, LOC101106051, ZIM3, PEG3, TRPC7, FBXL4, LOC105610397, LOC105616489 and DNAAF2 were suggested as candidates. Some of the discovered genes and involved pathways were already annotated to contribute growth and development in various species including human, mice and cattle. All in all, the results of this study are expected to strongly contribute to shed a light on the underlying molecular mechanisms behind growth and carcass composition traits, with potential implications on studies aiming faster genetic improvement, targeted low-resolution SNP panel designs and genome-editing studies.

A further survey of the quantitative trait loci affecting swine body size and carcass traits in five related pig populations

Breeding for good meat quality performance while maintaining large body size and desirable carcass traits has been the major challenge for modern swine selective breeding. To address this goal, in the present work we studied five related populations produced by two commercial breeds (Berkshire and Duroc) and two Chinese breeds (Licha black pig and Lulai black pig). A single-trait GWAS performed on 20 body size and carcass traits using a self-developed China Chip-1 porcine SNP50K BeadChip identified 11 genome-wide significant QTL on nine chromosomes and 22 suggestive QTL on 15 chromosomes. For the 11 genome-wide significant QTL, eight were detected in at least two populations, and the rest were population-specific and only mapped in Shanxia black pig. Most of the genome-wide significant QTL were pleiotropic; for example, the QTL around 75.65 Mb on SSC4 was associated with four traits at genome-wide significance level. After screening the genes within 50 kb of the top SNP for each genome-wide significant QTL, NR6A1 and VRTN were chosen as candidate genes for vertebrae number; PLAG1 and BMP2 were identified as candidate genes for body size; and MC4R was the strong candidate gene for body weight. The four genes have been reported as candidates for thoracic vertebrae number, lumbar vertebrae number, carcass length and body weight respectively in previous studies. The effects of VRTN on thoracic vertebrae number, carcass length and body length have been verified in Shanxia black pig. Therefore, the VRTN genotype could be used in gene-assisted selection, and this could accelerate genetic improvement of body size and carcass traits in Shanxia black pig.

A genome-wide association study, supported by a new chromosome-level genome assembly, suggests sox2 as a main driver of the undifferentiatiated ZZ/ZW sex determination of turbot (Scophthalmus maximus)

BACKGROUND

Understanding sex determination (SD) across taxa is a major challenge for evolutionary biology. The new genomic tools are paving the way to identify genomic features underlying SD in fish, a group frequently showing limited sex chromosome differentiation and high SD evolutionary turnover. Turbot (Scophthalmus maximus) is a commercially important flatfish with an undifferentiated ZW/ZZ SD system and remarkable sexual dimorphism. Here we describe a new long-read turbot genome assembly used to disentangle the genetic architecture of turbot SD by combining genomics and classical genetics approaches.

RESULTS

The new turbot genome assembly consists of 145 contigs (N50 = 22.9 Mb), 27 of them representing >95% of its estimated genome size. A genome wide association study (GWAS) identified a ~ 6.8 Mb region on chromosome 12 associated with sex in 69.4% of the 36 families analyzed. The highest associated markers flanked sox2, the only gene in the region showing differential expression between sexes before gonad differentiation. A single SNP showed consistent differences between Z and W chromosomes. The analysis of a broad sample of families suggested the presence of additional genetic and/or environmental factors on turbot SD.

CONCLUSIONS

The new chromosome-level turbot genome assembly, one of the most contiguous fish assemblies to date, facilitated the identification of sox2 as a consistent candidate gene putatively driving SD in this species. This chromosome SD system barely showed any signs of differentiation, and other factors beyond the main QTL seem to control SD in a certain proportion of families.

Estimating the effective sample size in association studies of quantitative traits

The effective sample size (ESS) is a metric used to summarize in a single term the amount of correlation in a sample. It is of particular interest when predicting the statistical power of genome-wide association studies (GWAS) based on linear mixed models. Here, we introduce an analytical form of the ESS for mixed-model GWAS of quantitative traits and relate it to empirical estimators recently proposed. Using our framework, we derived approximations of the ESS for analyses of related and unrelated samples and for both marginal genetic and gene-environment interaction tests. We conducted simulations to validate our approximations and to provide a quantitative perspective on the statistical power of various scenarios, including power loss due to family relatedness and power gains due to conditioning on the polygenic signal. Our analyses also demonstrate that the power of gene-environment interaction GWAS in related individuals strongly depends on the family structure and exposure distribution. Finally, we performed a series of mixed-model GWAS on data from the UK Biobank and confirmed the simulation results. We notably found that the expected power drop due to family relatedness in the UK Biobank is negligible.

Investigating mechanisms underlying genetic resistance to Salmon Rickettsial Syndrome in Atlantic salmon using RNA sequencing

Background

Salmon Rickettsial Syndrome (SRS), caused by Piscirickettsia salmonis, is one of the primary causes of morbidity and mortality in Atlantic salmon aquaculture, particularly in Chile. Host resistance is a heritable trait, and functional genomic studies have highlighted genes and pathways important in the response of salmon to the bacteria. However, the functional mechanisms underpinning genetic resistance are not yet well understood. In the current study, a large population of salmon pre-smolts were challenged with P. salmonis, with mortality levels recorded and samples taken for genotyping. In parallel, head kidney and liver samples were taken from animals of the same population with high and low genomic breeding values for resistance, and used for RNA-Sequencing to compare their transcriptome profile both pre and post infection.

Results

A significant and moderate heritability (h² = 0.43) was shown for the trait of binary survival. Genome-wide association analyses using 38 K imputed SNP genotypes across 2265 animals highlighted that resistance is a polygenic trait. Several thousand genes were identified as differentially expressed between controls and infected samples, and enriched pathways related to the host immune response were highlighted. In addition, several networks with significant correlation with SRS resistance breeding values were identified, suggesting their involvement in mediating genetic resistance. These included apoptosis, cytoskeletal organisation, and the inflammasome.

Conclusions

While resistance to SRS is a polygenic trait, this study has highlighted several relevant networks and genes that are likely to play a role in mediating genetic resistance. These genes may be future targets for functional studies, including genome editing, to further elucidate their role underpinning genetic variation in host resistance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07443-2.

Four pleiotropic loci associated with fat mass and lean mass

BACKGROUND

Fat mass and lean mass are two biggest components of body mass. Both fat mass and lean mass are under strong genetic determinants and are correlated.

METHODS

We performed a bivariate genome-wide association meta-analysis of (lean adjusted) leg fat mass and (fat adjusted) leg lean mass in 12,517 subjects from 6 samples, and followed by in silico replication in large-scale UK biobank cohort sample (N = 370 097).

RESULTS

We identified four loci that were significant at the genome-wide significance (GWS, α = 5.0 × 10-8) level at the discovery meta-analysis, and successfully replicated in the replication sample: 2q36.3 (rs1024137, pdiscovery = 3.32 × 10-8, preplication = 4.07 × 10-13), 5q13.1 (rs4976033, pdiscovery = 1.93 × 10-9, preplication = 6.35 × 10-7), 12q24.31 (rs4765528, pdiscovery = 7.19 × 10-12, preplication = 1.88 × 10-11) and 18q21.32 (rs371326986, pdiscovery = 9.04 × 10-9, preplication = 2.35 × 10-95). The above four pleiotropic loci may play a pleiotropic role for fat mass and lean mass development.

CONCLUSIONS

Our findings further enhance the understanding of the genetic association between fat mass and lean mass and provide a new theoretical basis for their understanding.

Statistically efficient association analysis of quantitative traits with haplotypes and untyped SNPs in family studies

BACKGROUND

Associations between haplotypes and quantitative traits provide valuable information about the genetic basis of complex human diseases. Haplotypes also provide an effective way to deal with untyped SNPs. Two major challenges arise in haplotype-based association analysis of family data. First, haplotypes may not be inferred with certainty from genotype data. Second, the trait values within a family tend to be correlated because of common genetic and environmental factors.

RESULTS

To address these challenges, we present an efficient likelihood-based approach to analyzing associations of quantitative traits with haplotypes or untyped SNPs. This approach properly accounts for within-family trait correlations and can handle general pedigrees with arbitrary patterns of missing genotypes. We characterize the genetic effects on the quantitative trait by a linear regression model with random effects and develop efficient likelihood-based inference procedures. Extensive simulation studies are conducted to examine the performance of the proposed methods. An application to family data from the Childhood Asthma Management Program Ancillary Genetic Study is provided. A computer program is freely available.

CONCLUSIONS

Results from extensive simulation studies show that the proposed methods for testing the haplotype effects on quantitative traits have correct type I error rates and are more powerful than some existing methods.

Genome-wide association studies in Samoans give insight into the genetic architecture of fasting serum lipid levels

The current understanding of the genetic architecture of lipids has largely come from genome-wide association studies (GWAS). To date, few GWAS have examined the genetic architecture of lipids in Polynesians, and none have in Samoans, whose unique population history, including many population bottlenecks, may provide insight into the biological foundations of variation in lipid levels. Here we performed a GWAS of four fasting serum lipid levels: total cholesterol (TC), high-density lipoprotein (HDL), low-density lipoprotein (LDL), and triglycerides (TG) in a sample of 2849 Samoans, with validation genotyping for associations in a replication cohort comprising 1798 Samoans and American Samoans. We identified multiple genome-wide significant associations (P < 5 × 10^-8) previously seen in other populations-APOA1 with TG, CETP with HDL, and APOE with TC and LDL-and several suggestive associations (P < 1 × 10^-5), including an association of variants downstream of MGAT1 and RAB21 with HDL. However, we observed different association signals for variants near APOE than what has been previously reported in non-Polynesian populations. The association with several known lipid loci combined with the newly identified associations with variants near MGAT1 and RAB21 suggest that while some of the genetic architecture of lipids is shared between Samoans and other populations, part of the genetic architecture may be Polynesian-specific.

Impact of Toll-like Receptors 2(TLR2) and TLR 4 Gene Variations on HCV Susceptibility, Response to Treatment and Development of Hepatocellular Carcinoma in Cirrhotic HCV Patients

Background and Aims: Genetic polymorphisms of Toll-like receptors (TLRs) have been proposed to affect susceptibility to HCV infection and progression to end-stage liver disease. This study was conducted to clarify the association of SNPS of TLR2 and TLR4 with clinical outcome of hepatitis C, response to treatment and development of HCC.Methods: The current study examined 3295 individuals from 725 families that were categorized into groups comprising chronic HCV (CH), spontaneous viral clearance (SC) and control subjects. Treated patients were classified into responders (RT) and non-responders (NRT). In addition, patients with liver cirrhotic (LC), and hepatocellular carcinoma (HCC) were also included. All subjects were genotyped for five single nucleotide polymorphisms (SNPs) of TLR2 and four SNPs of TLR4 and their haplotypes using allelic discrimination real-time PCR.Results: Results demonstrated strong association with allele A of rs13105517 of TLR2 and allele C of rs10116253 of TLR4 with CH in comparison to SC group. However, The peak of risk of HCC was observed with allele C of rs3804099 of TLR2 and C allele of rs10116253 TLR4 (p < 0.001).A strong association was found with allele T of rs1816702 of TLR2 and allele A of rs5030728 of TLR4 in non responder group in comparison to responders (p < 0.001). Haplotypes CAGT of TLR4 and ATAC of TLR2 showed significant association with CH and HCC groups in comparison to other groups.Conclusions: This study shows an association of minor alleles of TLR2 and TLR4 with outcome of HCV infection, response to therapy and development of HCC in cirrhotic patients.

Genomic regions influencing intramuscular fat in divergently selected rabbit lines

Intramuscular fat (IMF) is one of the main meat quality traits for breeding programmes in livestock species. The main objective of this study was to identify genomic regions associated with IMF content comparing two rabbit populations divergently selected for this trait, and to generate a list of putative candidate genes. Animals were genotyped using the Affymetrix Axiom OrcunSNP Array (200k). After quality control, the data involved 477 animals and 93 540 SNPs. Two methods were used in this research: single marker regressions with the data adjusted by genomic relatedness, and a Bayesian multiple marker regression. Associated genomic regions were located on the rabbit chromosomes (OCU) OCU1, OCU8 and OCU13. The highest value for the percentage of the genomic variance explained by a genomic region was found in two consecutive genomic windows on OCU8 (7.34%). Genes in the associated regions of OCU1 and OCU8 presented biological functions related to the control of adipose cell function, lipid binding, transportation and localisation (APOLD1, PLBD1, PDE6H, GPRC5D and GPRC5A) and lipid metabolic processes (MTMR2). The EWSR1 gene, underlying the OCU13 region, is linked to the development of brown adipocytes. The findings suggest that there is a large component of polygenic effect behind the differences in IMF content in these two lines, as the variance explained by most of the windows was low. The genomic regions of OCU1, OCU8 and OCU13 revealed novel candidate genes. Further studies would be needed to validate the associations and explore their possible application in selection programmes.

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.

Venous blood gas and chemistry components are moderately heritable in commercial white egg-laying hens under acute or chronic heat exposure

Heat stress has a large negative impact on poultry around the world in both intensive and small-scale production systems. Better understanding of genetic factors contributing to response to high ambient temperatures would provide a basis to develop strategies for alleviating negative impacts of heat on poultry production. The objective of this work was to characterize the genetic control (heritability estimate and quantitative trait loci (QTL)) of blood chemistry components before and after exposure to acute and chronic high ambient temperature in a commercial egg laying line Hy-Line W-36 female parent line mature hens were exposed to 4 wk of daily cyclic heat exposure. Blood was collected pre-heat, on the first day of heat, and 2 and 4 wk post heat initiation and analyzed immediately using an i-STAT® hand-held blood analyzer. Thirteen blood components were quantified at the 4 time points: pH, pCO2, pO2, HCO3, TCO2, sO2, iCa, Na, K, base excess, glucose, "hematocrit" (estimated from blood electrical conductivity, BEC), and "hemoglobin" (calculated from BEC). Heritabilities were estimated using genomic relationship information obtained from 600k SNP chip data. All 13 parameters exhibited a significant change after 5 h of heat exposure and most did not return to pre-heat levels throughout the duration of the study. Eight parameters (base excess, glucose, hemoglobin, HCO3, hematocrit, K, pCO2, TCO2) had heritability estimates differing from zero at one or more time points (0.21 to 0.45). The traits with significant heritability would be good candidates for use as biomarkers in a selection program if they are correlated with traits of economic importance that are more difficult to measure. QTL were identified for nine of the traits at one or more time point. These nine traits, however, did not have significant heritability estimates suggesting that while some QTL have been identified their effects are generally small.

Genome-wide association analyses reveal polygenic genomic architecture underlying divergent shell morphology in Spanish Littorina saxatilis ecotypes

Gene flow between diverging populations experiencing dissimilar ecological conditions can theoretically constrain adaptive evolution. To minimize the effect of gene flow, alleles underlying traits essential for local adaptation are predicted to be located in linked genome regions with reduced recombination. Local reduction in gene flow caused by selection is expected to produce elevated divergence in these regions. The highly divergent crab-adapted and wave-adapted ecotypes of the marine snail Littorina saxatilis present a model system to test these predictions. We used genome-wide association (GWA) analysis of geometric morphometric shell traits associated with microgeographic divergence between the two L. saxatilis ecotypes within three separate sampling sites. A total of 477 snails that had individual geometric morphometric data and individual genotypes at 4,066 single nucleotide polymorphisms (SNPs) were analyzed using GWA methods that corrected for population structure among the three sites. This approach allowed dissection of the genomic architecture of shell shape divergence between ecotypes across a wide geographic range, spanning two glacial lineages. GWA revealed 216 quantitative trait loci (QTL) with shell size or shape differences between ecotypes, with most loci explaining a small proportion of phenotypic variation. We found that QTL were evenly distributed across 17 linkage groups, and exhibited elevated interchromosomal linkage, suggesting a genome-wide response to divergent selection on shell shape between the two ecotypes. Shell shape trait-associated loci showed partial overlap with previously identified outlier loci under divergent selection between the two ecotypes, supporting the hypothesis of diversifying selection on these genomic regions. These results suggest that divergence in shell shape between the crab-adapted and wave-adapted ecotypes is produced predominantly by a polygenic genomic architecture with positive linkage disequilibrium among loci of small effect.

Using genetic data to strengthen causal inference in observational research

Causal inference is essential across the biomedical, behavioural and social sciences.By progressing from confounded statistical associations to evidence of causal relationships, causal inference can reveal complex pathways underlying traits and diseases and help to prioritize targets for intervention. Recent progress in genetic epidemiology - including statistical innovation, massive genotyped data sets and novel computational tools for deep data mining - has fostered the intense development of methods exploiting genetic data and relatedness to strengthen causal inference in observational research. In this Review, we describe how such genetically informed methods differ in their rationale, applicability and inherent limitations and outline how they should be integrated in the future to offer a rich causal inference toolbox.

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.

Genetic analysis of production, physiological, and egg quality traits in heat-challenged commercial white egg-laying hens using 600k SNP array data

Background

Heat stress negatively affects the welfare and production of chickens. High ambient temperature is considered one of the most ubiquitous abiotic environmental challenges to laying hens around the world. In this study, we recorded several production traits, feed intake, body weight, digestibility, and egg quality of 400 commercial white egg-laying hens before and during a 4-week heat treatment. For the phenotypes that had estimated heritabilities (using 600k SNP chip data) higher than 0, SNP associations were tested using the same 600k genotype data.

Results

Seventeen phenotypes had heritability estimates higher than 0, including measurements at various time points for feed intake, feed efficiency, body weight, albumen weight, egg quality expressed in Haugh units, egg mass, and also for change in egg mass from prior to heat exposure to various time points during the 4-week heat treatment. Quantitative trait loci (QTL) were identified for 10 of these 17 phenotypes. Some of the phenotypes shared QTL including Haugh units before heat exposure and after 4 weeks of heat treatment.

Conclusions

Estimated heritabilities differed from 0 for 17 traits, which indicates that they are under genetic control and that there is potential for improving these traits through selective breeding. The association of different QTL with the same phenotypes before heat exposure and during heat treatment indicates that genomic control of traits under heat stress is distinct from that under thermoneutral conditions. This study contributes to the knowledge on the genomic control of response to heat stress in laying hens.

Electronic supplementary material

The online version of this article (10.1186/s12711-019-0474-6) contains supplementary material, which is available to authorized users.

Genome-Wide Association Studies Based on Equine Joint Angle Measurements Reveal New QTL Affecting the Conformation of Horses

The evaluation of conformation traits is an important part of selection for breeding stallions and mares. Some of these judged conformation traits involve joint angles that are associated with performance, health, and longevity. To improve our understanding of the genetic background of joint angles in horses, we have objectively measured the angles of the poll, elbow, carpal, fetlock (front and hind), hip, stifle, and hock joints based on one photograph of each of the 300 Franches-Montagnes (FM) and 224 Lipizzan (LIP) horses. After quality control, genome-wide association studies (GWASs) for these traits were performed on 495 horses, using 374,070 genome-wide single nucleotide polymorphisms (SNPs) in a mixed-effect model. We identified two significant quantitative trait loci (QTL) for the poll angle on ECA28 (p = 1.36 × 10⁻⁷), 50 kb downstream of the ALX1 gene, involved in cranial morphology, and for the elbow joint on ECA29 (p = 1.69 × 10⁻⁷), 49 kb downstream of the RSU1 gene, and 75 kb upstream of the PTER gene. Both genes are associated with bone mineral density in humans. Furthermore, we identified other suggestive QTL associated with the stifle joint on ECA8 (p = 3.10 × 10⁻⁷); the poll on ECA1 (p = 6.83 × 10⁻⁷); the fetlock joint of the hind limb on ECA27 (p = 5.42 × 10⁻⁷); and the carpal joint angle on ECA3 (p = 6.24 × 10⁻⁷), ECA4 (p = 6.07 × 10⁻⁷), and ECA7 (p = 8.83 × 10⁻⁷). The application of angular measurements in genetic studies may increase our understanding of the underlying genetic effects of important traits in equine breeding.

Mapping of QTL for chicken body weight, carcass composition, and meat quality traits in a slow-growing line

Slow-growing chicken lines are valuable genetic resources for the development of well-perceived alternative free-range production. While there is no constraint on increasing growth rate, breeding programs have to evolve in order to include new traits improving the positioning of such lines in the growing market for parts and processed products. In this study, we used dense genotyping to fine map QTL for chicken growth, body composition, and meat quality traits in view of developing new tools for selection of a slow-growing line. The dataset included a total of 836 birds (10 sires, 87 dams, 739 descendants) and 40,203 SNP. QTL for the 15 traits analyzed were detected by 3 different methods, i.e., linkage and linkage disequilibrium haplotype-based analysis (LDLA), family-based single marker association (FASTA), and Bayesian multi-marker regression (Bayes Cπ). After filtering for QTL redundancy, we found 16, 16, and 9 QTL when using the FASTA, LDLA, and Bayes Cπ methods, respectively, with a threshold of 2.49 × 10-5 for FASTA and LDLA, and a Bayes factor of 150 for the Bayes Cπ analysis. They comprised 17 QTL for body weight, 9 QTL for body composition, and 15 QTL for breast meat quality or behavior at slaughter. The 3 methods agreed in the detection of highly significant QTL such as that detected on GGA24 for body weight at 3, 6, and 9 wk, and the 2 QTL detected on GGA17 and GGA18 for breast meat yield. Several significant QTL were also detected for the different components of breast meat quality. This study provided new locations for investigation in order to improve our understanding of the genetic architecture of growth, carcass composition, and meat quality in the chicken and to develop molecular tools for the selection of these traits in a slow-growing line.

OPENMENDEL: a cooperative programming project for statistical genetics

Statistical methods for genome-wide association studies (GWAS) continue to improve. However, the increasing volume and variety of genetic and genomic data make computational speed and ease of data manipulation mandatory in future software. In our view, a collaborative effort of statistical geneticists is required to develop open source software targeted to genetic epidemiology. Our attempt to meet this need is called the OPENMENDEL project (https://openmendel.github.io). It aims to (1) enable interactive and reproducible analyses with informative intermediate results, (2) scale to big data analytics, (3) embrace parallel and distributed computing, (4) adapt to rapid hardware evolution, (5) allow cloud computing, (6) allow integration of varied genetic data types, and (7) foster easy communication between clinicians, geneticists, statisticians, and computer scientists. This article reviews and makes recommendations to the genetic epidemiology community in the context of the OPENMENDEL project.

Family-Based Genome-Wide Association Study of South Indian Pedigrees Supports WNT7B as a Central Corneal Thickness Locus

Purpose

To identify genetic risk factors contributing to central corneal thickness (CCT) in individuals from South India, a population with a high prevalence of ocular disorders.

Methods

One hundred ninety-five individuals from 15 large South Indian pedigrees were genotyped using the Omni2.5 bead array. Family-based association for CCT was conducted using the score test in MERLIN.

Results

Genome-wide association study (GWAS) identified strongest association for single nucleotide polymorphisms (SNPs) in the first intron of WNT7B and CCT (top SNP rs9330813; β = −0.57, 95% confidence interval [CI]: −0.78 to −0.36; P = 1.7 × 10⁻⁷). We further investigated rs9330813 in a Latino cohort and four independent European cohorts. A meta-analysis of these data sets demonstrated statistically significant association between rs9330813 and CCT (β = −3.94, 95% CI: −5.23 to −2.66; P = 1.7 × 10⁻⁹). WNT7B SNPs located in the same genomic region that includes rs9330813 have previously been associated with CCT in Latinos but with other ocular quantitative traits related to myopia (corneal curvature and axial length) in a Japanese population (rs10453441 and rs200329677). To evaluate the specificity of the observed WNT7B association with CCT in the South Indian families, we completed an ocular phenome-wide association study (PheWAS) for the top WNT7B SNPs using 45 ocular traits measured in these same families including corneal curvature and axial length. The ocular PheWAS results indicate that in the South Indian families WNT7B SNPs are primarily associated with CCT.

Conclusions

The results indicate robust evidence for association between WNT7B SNPs and CCT in South Indian pedigrees, and suggest that WNT7B SNPs can have population-specific effects on ocular quantitative traits.

Revisit Population-based and Family-based Genotype Imputation

Genome-Wide Association (GWA) with population-based imputation (PBI) has been successful in identifying common variants associated with complex diseases; however, much heritability remains to be explained and low frequency variants (LFV) may contribute. To identify LFV, a study of unrelated individuals may no longer be as efficient as a family study, where rare population variants can be frequent in families. Family-based imputation (FBI) provides an opportunity to evaluate LFV. To compare the performance of PBI and FBI, we conducted extensive simulations, generating genotypes using SeqSIMLA from various reference panels for families. We masked genotype information for variants unavailable in Framingham 550 K GWA genotype data in less informative subjects selected by GIGI-Pick. We implemented IMPUTE2 with duoHMM in SHAPEIT (Impute2_duoHMM) for PBI, MERLIN and GIGI for FBI and PedBLIMP for a hybrid approach. In general, FBI in both MERLIN and GIGI outperformed other approaches with imputation accuracy greater than 0.99 for the squared correlation and imputation quality scores (IQS) especially for LFV, although imputation accuracy from MERLIN depends on pedigree splitting for larger families. PBI performed worst with the exception of good imputation accuracy for common variants when a closely ancestry matched reference is used. In summary, linkage disequilibrium (LD) information from large available genotype resources provides good imputation for common variants with well-selected reference panels without requiring densely sequenced data in family members, while imputation of LFV with FBI benefits more from information on inheritance patterns within families yielding better imputation.

Discovery and Functional Annotation of Quantitative Trait Loci Affecting Resistance to Sea Lice in Atlantic Salmon

Sea lice (Caligus rogercresseyi) are ectoparasitic copepods which have a large negative economic and welfare impact in Atlantic salmon (Salmo salar) aquaculture, particularly in Chile. A multi-faceted prevention and control strategy is required to tackle lice, and selective breeding contributes via cumulative improvement of host resistance to the parasite. While host resistance has been shown to be heritable, little is yet known about the individual loci that contribute to this resistance, the potential underlying genes, and their mechanisms of action. In this study we took a multifaceted approach to identify and characterize quantitative trait loci (QTL) affecting host resistance in a population of 2,688 Caligus-challenged Atlantic salmon post-smolts from a commercial breeding program. We used low and medium density genotyping with imputation to collect genome-wide SNP marker data for all animals. Moderate heritability estimates of 0.28 and 0.24 were obtained for lice density (as a measure of host resistance) and growth during infestation, respectively. Three QTL explaining between 7 and 13% of the genetic variation in resistance to sea lice (as represented by the traits of lice density) were detected on chromosomes 3, 18, and 21. Characterisation of these QTL regions was undertaken using RNA sequencing and pooled whole genome sequencing data. This resulted in the identification of a shortlist of potential underlying causative genes, and candidate functional mutations for further study. For example, candidates within the chromosome 3 QTL include a putative premature stop mutation in TOB1 (an anti-proliferative transcription factor involved in T cell regulation) and an uncharacterized protein which showed significant differential allelic expression (implying the existence of a cis-acting regulatory mutation). While host resistance to sea lice is polygenic in nature, the results of this study highlight significant QTL regions together explaining between 7 and 13 % of the heritability of the trait. Future investigation of these QTL may enable improved knowledge of the functional mechanisms of host resistance to sea lice, and incorporation of functional variants to improve genomic selection accuracy.

Association score testing for rare variants and binary traits in family data with shared controls

Genome-wide association studies have been an important approach used to localize trait loci, with primary focus on common variants. The multiple rare variant-common disease hypothesis may explain the missing heritability remaining after accounting for identified common variants. Advances of sequencing technologies with their decreasing costs, coupled with methodological advances in the context of association studies in large samples, now make the study of rare variants at a genome-wide scale feasible. The resurgence of family-based association designs because of their advantage in studying rare variants has also stimulated more methods development, mainly based on linear mixed models (LMMs). Other tests such as score tests can have advantages over the LMMs, but to date have mainly been proposed for single-marker association tests. In this article, we extend several score tests (χcorrected2, WQLS, and SKAT) to the multiple variant association framework. We evaluate and compare their statistical performances relative with the LMM. Moreover, we show that three tests can be cast as the difference between marker allele frequencies (AFs) estimated in each of the group of affected and unaffected subjects. We show that these tests are flexible, as they can be based on related, unrelated or both related and unrelated subjects. They also make feasible an increasingly common design that only sequences a subset of affected subjects (related or unrelated) and uses for comparison publicly available AFs estimated in a group of healthy subjects. Finally, we show the great impact of linkage disequilibrium on the performance of all these tests.

Genome scan for the possibility of identifying candidate resistance genes for goat lentiviral infections in the Italian Garfagnina goat breed

Small ruminant lentiviruses (SRLVs) are a heterogeneous group of viruses of sheep, goat, and wild ruminants responsible of lifelong persistent infection leading to a multisystem chronic disease. Increased evidences indicate that host genetic factors could influence the individual SRLV resistance. The present study was conducted on the Garfagnina goat breed, an Italian goat population registered on the Tuscan regional repertory of genetic resources at risk of extinction. Forty-eight adult goats belonging to a single flock were studied. SRLV diagnosis was achieved by serological tests and 21 serologically positive animals were identified. All animals were genotyped with the Illumina GoatSNP60 BeadChip and a genome-wide scan was then performed on the individual marker genotypes, in an attempt to identify genomic regions associated with the infection. One SNP was found significant (P < 5 × 10^-5) on CHR 18 at 62,360,918 bp. The SNP was an intron of the zinc finger protein 331 (ZNF331) protein. In the region 1 Mb upstream the significant SNP, the NLRP12 (NLR family pyrin domain containing 12), the PRKCG (protein kinase C gamma), and the CACNG7 (calcium voltage-gated channel auxiliary subunit gamma 7) were found.

Genomic Selection for Growth Traits in Pacific Oyster (Crassostrea gigas): Potential of Low-Density Marker Panels for Breeding Value Prediction

Pacific oysters are a key aquaculture species globally, and genetic improvement via selective breeding is a major target. Genomic selection has the potential to expedite genetic gain for key target traits of a breeding program, but has not yet been evaluated in oyster. The recent development of SNP arrays for Pacific oyster (Crassostrea gigas) raises the opportunity to test genomic selection strategies for polygenic traits. In this study, a population of 820 oysters (comprising 23 full-sibling families) were genotyped using a medium density SNP array (23 K informative SNPs), and the genetic architecture of growth-related traits [shell height (SH), shell length (SL), and wet weight (WW)] was evaluated. Heritability was estimated to be moderate for the three traits (0.26 ± 0.06 for SH, 0.23 ± 0.06 for SL and 0.35 ± 0.05 for WW), and results of a GWAS indicated that the underlying genetic architecture was polygenic. Genomic prediction approaches were used to estimate breeding values for growth, and compared to pedigree based approaches. The accuracy of the genomic prediction models (GBLUP) outperformed the traditional pedigree approach (PBLUP) by ∼25% for SL and WW, and ∼30% for SH. Further, reduction in SNP marker density had little impact on prediction accuracy, even when density was reduced to a few hundred SNPs. These results suggest that the use of genomic selection in oyster breeding could offer benefits for the selection of breeding candidates to improve complex economic traits at relatively modest cost.

Coping with family structure in genome-wide association studies: a comparative evaluation

In this paper, a fully statistical investigation of the control of family structure as random effects is analyzed and discussed, using both the genome-wide association studies (GWAS) data and simulated data. Three modeling strategies are proposed and the analysis results suggest the hybrid use of results from all possible models should be combined in practice.

Rare variant association test in family-based sequencing studies

The objective of this article is to introduce valid and robust methods for the analysis of rare variants for family-based exome chips, whole-exome sequencing or whole-genome sequencing data. Family-based designs provide unique opportunities to detect genetic variants that complement studies of unrelated individuals. Currently, limited methods and software tools have been developed to assist family-based association studies with rare variants, especially for analyzing binary traits. In this article, we address this gap by extending existing burden and kernel-based gene set association tests for population data to related samples, with a particular emphasis on binary phenotypes. The proposed approach blends the strengths of kernel machine methods and generalized estimating equations. Importantly, the efficient generalized kernel score test can be applied as a mega-analysis framework to combine studies with different designs. We illustrate the application of the proposed method using data from an exome sequencing study of autism. Methods discussed in this article are implemented in an R package 'gskat', which is available on CRAN and GitHub.

Genetic Analysis of a Commercial Egg Laying Line Challenged With Newcastle Disease Virus

In low income countries, chickens play a vital role in daily life. They provide a critical source of protein through egg production and meat. Newcastle disease, caused by avian paramyxovirus type 1, has been ranked as the most devastating disease for scavenging chickens in Africa and Asia. High mortality among flocks infected with velogenic strains leads to a devastating loss of dietary protein and buying power for rural households. Improving the genetic resistance of chickens to Newcastle Disease virus (NDV), in addition to vaccination, is a practical target for improvement of poultry production in low income countries. Because response to NDV has a component of genetic control, it can be influenced through selective breeding. Adding genomic information to a breeding program can increase the amount of genetic progress per generation. In this study, we challenged a commercial egg-laying line with a lentogenic strain of NDV, measured phenotypic responses, collected genotypes, and associated genotypes with phenotypes. Collected phenotypes included viral load at 2 and 6 days post-infection (dpi), antibody levels pre-challenge and 10 dpi, and growth rates pre- and post-challenge. Six suggestive QTL associated with response to NDV and/or growth were identified, including novel and known QTL confirming previously reported associations with related traits. Additionally, previous RNA-seq analysis provided support for several of the genes located in or near the identified QTL. Considering the trend of negative genetic correlation between antibody and Newcastle Disease tolerance (growth under disease) and estimates of moderate to high heritability, we provide evidence that these NDV response traits can be influenced through selective breeding. Producing chickens that perform favorably in challenging environments will ultimately increase the supply of quality protein for human consumption.

Discordant association of the CREBRF rs373863828 A allele with increased BMI and protection from type 2 diabetes in Māori and Pacific (Polynesian) people living in Aotearoa/New Zealand

AIMS/HYPOTHESIS

The A (minor) allele of CREBRF rs373863828 has been associated with increased BMI and reduced risk of type 2 diabetes in the Samoan populations of Samoa and American Samoa. Our aim was to test rs373863828 for associations with BMI and the odds of type 2 diabetes, gout and chronic kidney disease (CKD) in Māori and Pacific (Polynesian) people living in Aotearoa/New Zealand.

METHODS

Linear and logistic regression models were used to analyse the association of the A allele of CREBRF rs373863828 with BMI, log-transformed BMI, waist circumference, type 2 diabetes, gout and CKD in 2286 adults. The primary analyses were adjusted for age, sex, the first four genome-wide principal components and (where appropriate) BMI, waist circumference and type 2 diabetes. The primary analysis was conducted in ancestrally defined groups and association effects were combined using meta-analysis.

RESULTS

For the A allele of rs373863828, the effect size was 0.038 (95% CI 0.022, 0.055, p = 4.8 × 10-6) for log-transformed BMI, with OR 0.59 (95% CI 0.47, 0.73, p = 1.9 × 10-6) for type 2 diabetes. There was no evidence for an association of genotype with variance in BMI (p = 0.13), and nor was there evidence for associations with serum urate (β = 0.012 mmol/l, pcorrected = 0.10), gout (OR 1.00, p = 0.98) or CKD (OR 0.91, p = 0.59).

CONCLUSIONS/INTERPRETATION

Our results in New Zealand Polynesian adults replicate, with very similar effect sizes, the association of the A allele of rs373863828 with higher BMI but lower odds of type 2 diabetes among Samoan adults living in Samoa and American Samoa.

Genome-wide association study for bone strength in laying hens

Bone fracture in egg laying hens is a growing welfare and economic concern in the industry. Although environmental conditions and management (especially nutrition) can exacerbate it, the primary cause of bone weakness and the resulting fractures is believed to have a genetic basis. To test this hypothesis, we performed a genome-wide association study to identify the loci associated with bone strength in laying hens. Genotype and phenotype data were obtained from 752 laying hens belonging to the same pure line population. These hens were genotyped for 580,961 SNPs, with 232,021 SNPs remaining after quality control. Each of the SNPs were tested for association with tibial breaking strength using the family-based score test for association. A total of 52 SNPs across chromosomes 1, 3, 8, and 16 were significantly associated with tibial breaking strength with the genome-wide significance threshold set as a corrected P value of 10e-5. Based on the local linkage disequilibrium around the significant SNPs, 5 distinct and novel QTLs were identified on chromosomes 1 (2 QTLs), 3 (1 QTL), 8 (1 QTL) and 16 (1 QTL). The strongest association was detected within the QTL region on chromosome 8, with the most significant SNP having a corrected P value of 4e-7. A number of candidate genes were identified within the QTL regions, including the BRD2 gene that is required for normal bone physiology. Bone-related pathways involving some of the genes were also identified including chloride channel activity, which regulates bone reabsorption, and intermediate filament organization, which plays a role in the regulation of bone mass. Our result supports previous studies that suggest that bone strength is highly regulated by genetics. It is therefore possible to reduce bone fractures in laying hens through genetic selection and ultimately improve hen welfare.

Genome-wide association meta-analysis highlights light-induced signaling as a driver for refractive error

Refractive errors, including myopia, are the most frequent eye disorders worldwide and an increasingly common cause of blindness. This genome-wide association meta-analysis in 160,420 participants and replication in 95,505 participants increased the number of established independent signals from 37 to 161 and showed high genetic correlation between Europeans and Asians (>0.78). Expression experiments and comprehensive in silico analyses identified retinal cell physiology and light processing as prominent mechanisms, and also identified functional contributions to refractive-error development in all cell types of the neurosensory retina, retinal pigment epithelium, vascular endothelium and extracellular matrix. Newly identified genes implicate novel mechanisms such as rod-and-cone bipolar synaptic neurotransmission, anterior-segment morphology and angiogenesis. Thirty-one loci resided in or near regions transcribing small RNAs, thus suggesting a role for post-transcriptional regulation. Our results support the notion that refractive errors are caused by a light-dependent retina-to-sclera signaling cascade and delineate potential pathobiological molecular drivers.

Genome-wide association study identifies a novel locus for cannabis dependence

Despite moderate heritability, only one study has identified genome-wide significant loci for cannabis-related phenotypes. We conducted meta-analyses of genome-wide association study data on 2080 cannabis-dependent cases and 6435 cannabis-exposed controls of European descent. A cluster of correlated single-nucleotide polymorphisms (SNPs) in a novel region on chromosome 10 was genome-wide significant (lowest P=1.3E-8). Among the SNPs, rs1409568 showed enrichment for H3K4me1 and H3K427ac marks, suggesting its role as an enhancer in addiction-relevant brain regions, such as the dorsolateral prefrontal cortex and the angular and cingulate gyri. This SNP is also predicted to modify binding scores for several transcription factors. We found modest evidence for replication for rs1409568 in an independent cohort of African American (896 cases and 1591 controls; P=0.03) but not European American (EA; 781 cases and 1905 controls) participants. The combined meta-analysis (3757 cases and 9931 controls) indicated trend-level significance for rs1409568 (P=2.85E-7). No genome-wide significant loci emerged for cannabis dependence criterion count (n=8050). There was also evidence that the minor allele of rs1409568 was associated with a 2.1% increase in right hippocampal volume in an independent sample of 430 EA college students (fwe-P=0.008). The identification and characterization of genome-wide significant loci for cannabis dependence is among the first steps toward understanding the biological contributions to the etiology of this psychiatric disorder, which appears to be rising in some developed nations.

A Genome-Wide Association Study for Host Resistance to Ostreid Herpesvirus in Pacific Oysters (Crassostrea gigas)

Ostreid herpesvirus (OsHV) can cause mass mortality events in Pacific oyster aquaculture. While various factors impact on the severity of outbreaks, it is clear that genetic resistance of the host is an important determinant of mortality levels. This raises the possibility of selective breeding strategies to improve the genetic resistance of farmed oyster stocks, thereby contributing to disease control. Traditional selective breeding can be augmented by use of genetic markers, either via marker-assisted or genomic selection. The aim of the current study was to investigate the genetic architecture of resistance to OsHV in Pacific oyster, to identify genomic regions containing putative resistance genes, and to inform the use of genomics to enhance efforts to breed for resistance. To achieve this, a population of ∼1,000 juvenile oysters were experimentally challenged with a virulent form of OsHV, with samples taken from mortalities and survivors for genotyping and qPCR measurement of viral load. The samples were genotyped using a recently-developed SNP array, and the genotype data were used to reconstruct the pedigree. Using these pedigree and genotype data, the first high density linkage map was constructed for Pacific oyster, containing 20,353 SNPs mapped to the ten pairs of chromosomes. Genetic parameters for resistance to OsHV were estimated, indicating a significant but low heritability for the binary trait of survival and also for viral load measures (h2 0.12 – 0.25). A genome-wide association study highlighted a region of linkage group 6 containing a significant QTL affecting host resistance. These results are an important step toward identification of genes underlying resistance to OsHV in oyster, and a step toward applying genomic data to enhance selective breeding for disease resistance in oyster aquaculture.