Positional Candidate Genes Research Articles

QTL-mapping in the obese Berlin Fat Mouse identifies additional candidate genes for obesity and fatty liver disease

The Berlin Fat Mouse Inbred line (BFMI) is a model for obesity and the metabolic syndrome. This study aimed to identify genetic variants associated with liver weight, liver triglycerides, and body weight using the obese BFMI sub-line BFMI861-S1. BFMI861-S1 mice are insulin resistant and store ectopic fat in the liver. In generation 10, 58 males and 65 females of the advanced intercross line (AIL) BFMI861-S1xB6N were phenotyped under a standard diet over 20 weeks. QTL analysis was performed after genotyping with the MiniMUGA Genotyping Array. Whole-genome sequencing and gene expression data of the parental lines was used for the prioritization of positional candidate genes. Three QTLs associated with liver weight, body weight, and subcutaneous adipose tissue (scAT) weight were identified. A highly significant QTL on chromosome (Chr) 1 (157–168 Mb) showed an association with liver weight. A QTL for body weight at 20 weeks was found on Chr 3 (34.1–40 Mb) overlapping with a QTL for scAT weight. In a multiple QTL mapping approach, an additional QTL affecting body weight at 16 weeks was identified on Chr 6 (9.5–26.1 Mb). Considering sequence variants and expression differences, Sec16b and Astn1 were prioritized as top positional candidate genes for the liver weight QTL on Chr 1; Met and Ica1 for the body weight QTL on Chr 6. Interestingly, all top candidate genes have previously been linked with metabolic traits. This study shows once more the power of an advanced intercross line for fine mapping. QTL mapping combined with a detailed prioritization approach allowed us to identify additional and plausible candidate genes linked to metabolic traits in the BFMI861-S1xB6N AIL. By reidentifying known candidate genes in a different crossing population the causal link with specific traits is underlined and additional evidence is given for further investigations.

Functional genomics of reproduction in pigs: Are we there yet?

Reproductive failure is the main reason for culling females in swine herds and is both a financial and sustainability issue. Because reproductive traits are complex and lowly to moderately heritable, genomic selection within populations can achieve substantial genetic gain in reproductive efficiency. A better understanding of the physiological components affecting the expression of these traits will facilitate greater understanding of the genes affecting reproductive traits and is necessary to improve and optimize management strategies to maximize reproductive success of gilts and sows. Large-scale genotyping with single-nucleotide polymorphism (SNP) arrays are used for genome-wide association studies (GWAS) and have facilitated identification of positional candidate genes. Transcriptomic data can be used to weight SNP for GWAS and could lead to previously unidentified candidate genes. Resequencing and fine mapping of candidate genes are necessary to identify putative functional variants and some of these have been incorporated into new genotyping arrays. Sequence imputation and genotype by sequence are newer strategies that could reveal novel functional mutations. In this study, these approaches are discussed. Advantages and limitations are highlighted where additional research is needed.

Integration of genome-wide association study and selection signatures reveals genetic determinants for skeletal muscle production traits in an F2 chicken population

Improving the production of broiler chicken meat has been a goal of broiler breeding programs worldwide for many years. However, the genetic architectures of skeletal muscle production traits in chickens have not yet been fully elucidated. In the present study, a total of 519 F2 birds, derived from a cross of Arbor Acres broiler and Baier layer, were re-sequenced (26 F0 individuals were re-sequenced at a 10-fold depth; 519 F2 individuals were re-sequenced at a 3-fold depth) and the coupling of genome-wide association study (GWAS) and selection signatures (FST (fixation index) and θπ (nucleotide diversity)) was carried out to pinpoint the associated loci and genes that contribute to pectoral muscle weight (PMW) and thigh muscle weight (TMW). A total of 7 890 258 single nucleotide polymorphisms (SNPs) remained to be analyzed after quality control and imputation. The integration of GWAS and selection signature analyses revealed that genetic determinants responsible for skeletal muscle production traits were mainly localized on chromosomes 1 (168.95–172.43 Mb) and 4 (74.37–75.23 Mb). A total of 17 positional candidate genes (PCGs) (LRCH1, CDADC1, CAB39L, LOC112531568, LOC112531569, FAM124A, FOXO1, NBEA, GPALPP1, RUBCNL, ARL11, KPNA3, LHFP, GBA3, LOC112532426, KCNIP4, and SLIT2) were identified in these regions. In particular, KPNA3 and FOXO1 were the most promising candidates for meat production in chickens. These findings will help enhance our understanding of the genetic architecture of chicken muscle production traits, and the significant SNPs identified could be promising candidates for integration into practical breeding programs such as genome-wide selection (GS) to improve the meat yield of chickens.

In Vivo Evaluation of the Cidec R46S Polymorphism for an Obesity QTL, tabw2a

ObjectivesA quantitative trait locus (QTL), tabw2a, was linked to obesity in TALLYHO/Jng (TH) mice, a polygenic model of obesity and type 2 diabetes. In combination with fine mapping information, we used whole genome sequencing data of TH mice to search positional candidate genes for tabw2a. From this search, the cell death-inducing DFFA-like effector c (Cidec) gene, containing a SNP identified as deleterious (likely to affect protein function) by both SIFT and PROVEAN, was raised as a potential candidate for tabw2a. CIDEC is a lipid droplet protein involved in the regulation of cellular lipid droplet size and lipid storage in adipocytes. There was one nucleotide substitution in the coding sequence of Cidec in TH mice, 136 C > A, compared to C57BL/6J (B6) reference sequence, that results in an amino acid substitution R46S (Arginine to Serine). The aim of this study was to evaluate the Cidec R46S polymorphism in vivo for the tabw2a QTL.MethodsWe generated knock-in mice where the TH allele (S46) is exchanged for the B6 allele (R46) at the Cidec gene on a B6 background. Targeting, ES clone selection, blastocyst injection, and breeding of ES cell chimeras were performed by Transgenics Core facility at Pennington Biomedical Research Center (Baton Rouge, LA). We bred the chimeras to B6 mice for 10 generation and collected heterozygous mice with germline transmission of the knock-in allele of Cidec S46. We then interbred the heterozygous mice and obtained homozygous mice for the knock-in allele of Cidec S46. Male and female Cidec R46S knock-in and B6 mice were weaned onto chow and high fat (HF) diets (D12266B, Research Diets) at 4 weeks of age and maintained on the same diets. At 20 weeks of age, body composition including fat and lean mass was measured using EchoMRI-100 whole body composition analyzer (Houston, TX).ResultsOn chow, no genotype effects were shown on adiposity between knock-in and B6 mice for both sexes. However, on HF diets, for both males and females, knock-in homozygous mice had significantly larger body fat mass than B6, although their body weights were comparable.ConclusionsWe conclude that the Cidec gene is a likely candidate for the obesity effect of tabw2a, and that Cidec S46 variant contributes to the obesity susceptibility.Funding SourcesAmerican Heart Association and National Institutes of Health.

Genetic variants associated with two major bovine milk fatty acids offer opportunities to breed for altered milk fat composition

BackgroundAlthough bovine milk is regarded as healthy and nutritious, its high content of saturated fatty acids (FA) may be harmful to cardiovascular health. Palmitic acid (C16:0) is the predominant saturated FA in milk with adverse health effects that could be countered by substituting it with higher levels of unsaturated FA, such as oleic acid (C18:1cis-9). In this work, we performed genome-wide association analyses for milk fatty acids predicted from FTIR spectroscopy data using 1811 Norwegian Red cattle genotyped and imputed to a high-density 777k single nucleotide polymorphism (SNP)-array. In a follow-up analysis, we used imputed whole-genome sequence data to detect genetic variants that are involved in FTIR-predicted levels of C16:0 and C18:1cis-9 and explore the transcript profile and protein level of candidate genes.ResultsGenome-wise significant associations were detected for C16:0 on Bos taurus (BTA) autosomes 11, 16 and 27, and for C18:1cis-9 on BTA5, 13 and 19. Closer examination of a significant locus on BTA11 identified the PAEP gene, which encodes the milk protein β-lactoglobulin, as a particularly attractive positional candidate gene. At this locus, we discovered a tightly linked cluster of genetic variants in coding and regulatory sequences that have opposing effects on the levels of C16:0 and C18:1cis-9. The favourable haplotype, linked to reduced levels of C16:0 and increased levels of C18:1cis-9 was also associated with a marked reduction in PAEP expression and β-lactoglobulin protein levels. β-lactoglobulin is the most abundant whey protein in milk and lower levels are associated with important dairy production parameters such as improved cheese yield.ConclusionsThe genetic variants detected in this study may be used in breeding to produce milk with an improved FA health-profile and enhanced cheese-making properties.

Genotype by environment interactions for chronic wasting disease in farmed US white-tailed deer.

Despite implementation of enhanced management practices, chronic wasting disease in US white-tailed deer (Odocoileus virginianus) continues to expand geographically. Herein, we perform the largest genome-wide association analysis to date for chronic wasting disease (n = 412 chronic wasting disease-positive; n = 758 chronic wasting disease-nondetect) using a custom Affymetrix Axiom single-nucleotide polymorphism array (n = 121,010 single-nucleotide polymorphisms), and confirm that differential susceptibility to chronic wasting disease is a highly heritable () polygenic trait in farmed US white-tailed deer, but with greater trait complexity than previously appreciated. We also confirm PRNP codon 96 (G96S) as having the largest-effects on risk (P ≤ 3.19E-08; phenotypic variance explained ≥ 0.025) across 3 US regions (Northeast, Midwest, South). However, 20 chronic wasting disease-positive white-tailed deer possessing codon 96SS genotypes were also observed, including one that was lymph node and obex positive. Beyond PRNP, we also detected 23 significant single-nucleotide polymorphisms (P-value ≤ 5E-05) implicating ≥24 positional candidate genes; many of which have been directly implicated in Parkinson’s, Alzheimer’s and prion diseases. Genotype-by-environment interaction genome-wide association analysis revealed a single-nucleotide polymorphism in the lysosomal enzyme gene ARSB as having the most significant regional heterogeneity of effects on chronic wasting disease (P ≤ 3.20E-06); with increasing copy number of the minor allele increasing susceptibility to chronic wasting disease in the Northeast and Midwest; but with opposite effects in the South. In addition to ARSB, 38 significant genotype-by-environment single-nucleotide polymorphisms (P-value ≤ 5E-05) were also detected, thereby implicating ≥ 36 positional candidate genes; the majority of which have also been associated with aspects of Parkinson’s, Alzheimer’s, and prion diseases.

Identifying pleiotropic variants and candidate genes for fertility and reproduction traits in Holstein cattle via association studies based on imputed whole-genome sequence genotypes

BackgroundGenetic progress for fertility and reproduction traits in dairy cattle has been limited due to the low heritability of most indicator traits. Moreover, most of the quantitative trait loci (QTL) and candidate genes associated with these traits remain unknown. In this study, we used 5.6 million imputed DNA sequence variants (single nucleotide polymorphisms, SNPs) for genome-wide association studies (GWAS) of 18 fertility and reproduction traits in Holstein cattle. Aiming to identify pleiotropic variants and increase detection power, multiple-trait analyses were performed using a method to efficiently combine the estimated SNP effects of single-trait GWAS based on a chi-square statistic.ResultsThere were 87, 72, and 84 significant SNPs identified for heifer, cow, and sire traits, respectively, which showed a wide and distinct distribution across the genome, suggesting that they have relatively distinct polygenic nature. The biological functions of immune response and fatty acid metabolism were significantly enriched for the 184 and 124 positional candidate genes identified for heifer and cow traits, respectively. No known biological function was significantly enriched for the 147 positional candidate genes found for sire traits. The most important chromosomes that had three or more significant QTL identified are BTA22 and BTA23 for heifer traits, BTA8 and BTA17 for cow traits, and BTA4, BTA7, BTA17, BTA22, BTA25, and BTA28 for sire traits. Several novel and biologically important positional candidate genes were strongly suggested for heifer (SOD2, WTAP, DLEC1, PFKFB4, TRIM27, HECW1, DNAH17, and ADAM3A), cow (ANXA1, PCSK5, SPESP1, and JMJD1C), and sire (ELMO1, CFAP70, SOX30, DGCR8, SEPTIN14, PAPOLB, JMJD1C, and NELL2) traits.ConclusionsThese findings contribute to better understand the underlying biological mechanisms of fertility and reproduction traits measured in heifers, cows, and sires, which may contribute to improve genomic evaluation for these traits in dairy cattle.

Genome-wide association mapping of ethanol sensitivity in the Diversity Outbred mouse population.

A strong predictor for the development of alcohol use disorder (AUD) is altered sensitivity to the intoxicating effects of alcohol. Individual differences in the initial sensitivity to alcohol are controlled in part by genetic factors. Mice offer a powerful tool to elucidate the genetic basis of behavioral and physiological traits relevant to AUD, but conventional experimental crosses have only been able to identify large chromosomal regions rather than specific genes. Genetically diverse, highly recombinant mouse populations make it possible to observe a wider range of phenotypic variation, offer greater mapping precision, and thus increase the potential for efficient gene identification. We have taken advantage of the Diversity Outbred (DO) mouse population to identify and precisely map quantitative trait loci (QTL) associated with ethanol sensitivity. We phenotyped 798male J:DO mice for three measures of ethanol sensitivity: ataxia, hypothermia, and loss of the righting response. We used high-density MegaMUGA and GigaMUGA to obtain genotypes ranging from 77,808 to 143,259 SNPs. We also performed RNA sequencing in striatum to map expression QTLs and identify gene expression-trait correlations. We then applied a systems genetic strategy to identify narrow QTLs and construct the network of correlations that exists between DNA sequence, gene expression values, and ethanol-related phenotypes to prioritize our list of positional candidate genes. We observed large amounts of phenotypic variation with the DO population and identified suggestive and significant QTLs associated with ethanol sensitivity on chromosomes 1, 2, and 16. The implicated regions were narrow (4.5-6.9 Mb in size) and each QTL explained ~4-5% of the variance. Our results can be used to identify alleles that contribute to AUD in humans, elucidate causative biological mechanisms, or assist in the development of novel therapeutic interventions.

Identification of Key Candidate Genes for Beak Length Phenotype by Whole-Genome Resequencing in Geese.

The domestic goose is an important economic animal in agriculture and its beak, a trait with high heritability, plays an important role in promoting food intake and defending against attacks. In this study, we sequenced 772 420-day-old Xingguo gray geese (XGG) using a low-depth (~1 ×) whole-genome resequencing strategy. We detected 12,490,912 single nucleotide polymorphisms (SNPs) using the standard GATK and imputed with STITCH. We then performed a genome-wide association study on the beak length trait in XGG. The results indicated that 57 SNPs reached genome-wide significance levels for the beak length trait and were assigned to seven genes, including TAPT1, DHX15, CCDC149, LGI2, SEPSECS, ANAPC4, and Slc34a2. The different genotypes of the most significant SNP (top SNP), which was located upstream of LGI2 and explained 7.24% of the phenotypic variation in beak length, showed significant differences in beak length. Priority-based significance analysis concluded that CCDC149, LGI2, and SEPSECS genes in the most significant quantitative trait locus interval were the most plausible positional and functional candidate genes for beak length development in the XGG population. These findings not only enhance our understanding of the genetic mechanism of the beak length phenotype in geese, but also lay the foundation for further studies to facilitate the genetic selection of traits in geese.

Congenital idiopathic megaesophagus in the German shepherd dog is a sex-differentiated trait and is associated with an intronic variable number tandem repeat in Melanin-Concentrating Hormone Receptor 2.

Congenital idiopathic megaesophagus (CIM) is a gastrointestinal (GI) motility disorder of dogs in which reduced peristaltic activity and dilation of the esophagus prevent the normal transport of food into the stomach. Affected puppies regurgitate meals and water, fail to thrive, and experience complications such as aspiration pneumonia that may necessitate euthanasia. The German shepherd dog (GSD) has the highest disease incidence, indicative of a genetic predisposition. Here, we discover that male GSDs are twice as likely to be affected as females and show that the sex bias is independent of body size. We propose that female endogenous factors (e.g., estrogen) are protective via their role in promoting relaxation of the sphincter between the esophagus and stomach, facilitating food passage. A genome-wide association study for CIM revealed an association on canine chromosome 12 (P-val = 3.12x10-13), with the lead SNPs located upstream or within Melanin-Concentrating Hormone Receptor 2 (MCHR2), a compelling positional candidate gene having a role in appetite, weight, and GI motility. Within the first intron of MCHR2, we identified a 33 bp variable number tandem repeat (VNTR) containing a consensus binding sequence for the T-box family of transcription factors. Across dogs and wolves, the major allele includes two copies of the repeat, whereas the predominant alleles in GSDs have one or three copies. The single-copy allele is strongly associated with CIM (P-val = 1.32x10-17), with homozygosity for this allele posing the most significant risk. Our findings suggest that the number of T-box protein binding motifs may correlate with MCHR2 expression and that an imbalance of melanin-concentrating hormone plays a role in CIM. We describe herein the first genetic factors identified in CIM: sex and a major locus on chromosome 12, which together predict disease state in the GSD with greater than 75% accuracy.

Genome-wide association study using haplotype libraries and repeated-measures model to identify candidate genomic regions for stillbirth in Holstein cattle

Reduced fertility is one of the main causes of economic losses on dairy farms, resulting in economic losses estimated at $938 per stillbirth case in Holstein herds. The identification of genomic regions associated with stillbirth could help to develop better management and breeding strategies aimed to reduce the frequency of undesirable gestation outcomes. Here, 10,570 cows and 50,541 birth records were used to perform a haplotype-based GWAS. A total of 41 significantly associated pseudo-SNPs (haplotypes within haplotype blocks converted to a binary classification) were identified after Bonferroni adjustment for multiple tests. A total of 117 positional candidate genes were annotated within or close (in a 200-kb interval) to significant pseudo-SNPs (haplotype blocks). The guilt-by-association functional prioritization identified 31 potential functional candidate genes for reproductive performance out of the 117 positional candidate genes annotated. These genes play crucial roles in biological processes associated with pregnancy persistence, fetus development, immune response, among others. These results helped us to better understand the genetic basis of stillbirth in dairy cattle and may be useful for the prediction of stillbirth in Holstein cattle, helping to reduce the related economic losses caused by this phenotype.

Novel Positional and Biological Candidate Gene for Grip Strength in Older Adults: The Long Life Family Study

Grip strength is a robust indicator of overall health, is moderately heritable and predicts longevity in older adults. Using genome-wide linkage analysis, we identified a novel locus on chromosome 18p linked to grip strength in 4534 individuals from 582 families (age 70.0 ± 15.8, range 24–110 years; 54% women). DNA sequencing was completed to identify single nucleotide variants (SNVs) in the 3.44 – 4.04 mega-basepair region on chromosome 18p. Using the sequencing data, we performed association analyses between the 7312 SNVs in the region and grip strength in families exhibiting evidence for linkage. Models were adjusted for age, age2, sex, height, field center and population substructure. There were 23 families (263 individuals) that contributed to the linkage peak (cumulative logarithm of the odds [LOD] score = 12.4). Six families (112 individuals) accounted for most of the linkage signal (LOD = 6.4). In these 6 families, we found highly significant associations between SNVs in the Disks Large-associated Protein 1 (DLGAP1) gene and grip strength (lead SNV: β= -0.75kg ± 0.15, p-value= 4.3*10-6). Correcting for the top SNV in DLGAP1 reduces the LOD by 72% in these families. Further, the effect allele frequency is much higher in these 6 families (39.7%) compared with both the NHLBI’s Trans-OMICs for Precision Medicine (23.5%) and 1000Genomes (28.0%) references panels. The DLGAP1 gene plays an important role in post-synaptic density of neurons; thus, it is a novel positional and biological candidate gene for follow-up studies aimed at uncovering genetic determinants of muscle strength.

Novel Positional Candidate Gene for Gait Speed Change in Older Adults: The Long Life Family Study

Older adults living alone are at higher risk of mortality, morbidity and healthcare utilization. As more older adults live alone, Emergency Department (ED) admissions could rapidly increase, particularly among those with multimorbidity. We studied the association of living alone on ED admissions among older adults with multimorbidity. We used data from 16,785 older adults of the population-based Singapore Chinese Health Study (mean age: 73 years, range: 61-96 years) who were interviewed in 2014-2016 for living arrangements and medical history. Participants were followed-up for one year on ED admission outcomes (number of admissions, inpatient days and hospitalization costs). We used multivariable logistic regression to study the association between living alone and ED admission, and ran two-part models (probit & generalised linear model) to estimate the association of living alone on inpatient days and hospitalization cost. We found that compared to living with others, living alone was associated with a higher odds of ED admissions [Odds Ratio (OR) 1.28, 95% Confidence Interval (CI) 1.08-1.51)], longer inpatient days (+0.61, 95% CI 0.25-0.97) and higher hospitalization costs (+322 USD, 95% CI 54-591). Compared to those living with others without multimorbidity, living alone with multimorbidity was associated with higher odds of ED admission (OR 1.64 95% CI 1.33-2.03), longer inpatient days (+0.73, 95% CI 0.29-1.17) and higher hospitalization costs (+567 USD, 95% CI 230-906). In conclusion, living alone is associated with higher odds of ED admission, longer inpatient days and higher hospitalization costs among older adults, particularly among those with multimorbidity.

Novel Positional Candidate Gene for Gait Speed Change in Older Adults: The Long Life Family Study

Research has shown cognitive ability in older age is associated with activity engagement, but little is known about what psychological mechanisms are linking the two constructs. This study investigates an emotional pathway, in which affective states mediate the temporal associations between momentary working memory and momentary activities in older age. We examined data from 153 healthy older adults aged 65 to 91 who completed a smartphone-based ambulatory assessment survey seven times a day over 15 days. In each assessment point, participants reported their momentary activities (e.g., social activities, mentally stimulating activities) and affective states (i.e., positive affect, negative affect) and took a working memory task. Initial results suggest that during an approximate time period of six hours (i.e., across three assessment points), working memory performance influences subsequent likelihood of social activity engagement. Moreover, positive affect mediates this temporal association. Results will be discussed in the context of cognitive aging research.

Candidate genes for congenital malformations in pigs

Candidate genes for congenital malformations in pigs

Identification of four novel QTL linked to the metabolic syndrome in the Berlin Fat Mouse

BackgroundThe Berlin Fat Mouse Inbred line (BFMI) is a model for obesity and the metabolic syndrome. This study aimed to identify genetic variants associated with impaired glucose metabolism using the obese lines BFMI861-S1 and BFMI861-S2, which are genetically closely related, but differ in several traits. BFMI861-S1 is insulin resistant and stores ectopic fat in the liver, whereas BFMI861-S2 is insulin sensitive.MethodsIn generation 10, 397 males of an advanced intercross line (AIL) BFMI861-S1 × BFMI861-S2 were challenged with a high-fat, high-carbohydrate diet and phenotyped over 25 weeks. QTL-analysis was performed after selective genotyping of 200 mice using the GigaMUGA Genotyping Array. Additional 197 males were genotyped for 7 top SNPs in QTL regions. For the prioritization of positional candidate genes whole genome sequencing and gene expression data of the parental lines were used.ResultsOverlapping QTL for gonadal adipose tissue weight and blood glucose concentration were detected on chromosome (Chr) 3 (95.8–100.1 Mb), and for gonadal adipose tissue weight, liver weight, and blood glucose concentration on Chr 17 (9.5–26.1 Mb). Causal modeling suggested for Chr 3-QTL direct effects on adipose tissue weight, but indirect effects on blood glucose concentration. Direct effects on adipose tissue weight, liver weight, and blood glucose concentration were suggested for Chr 17-QTL. Prioritized positional candidate genes for the identified QTL were Notch2 and Fmo5 (Chr 3) and Plg and Acat2 (Chr 17). Two additional QTL were detected for gonadal adipose tissue weight on Chr 15 (67.9–74.6 Mb) and for body weight on Chr 16 (3.9–21.4 Mb).ConclusionsQTL mapping together with a detailed prioritization approach allowed us to identify candidate genes associated with traits of the metabolic syndrome. In addition, we provided evidence for direct and indirect genetic effects on blood glucose concentration in the insulin-resistant mouse line BFMI861-S1.

Discovering of Genomic Variations Associated to Growth Traits by GWAS in Braunvieh Cattle.

A genome-wide association study (GWAS) was performed to elucidate genetic architecture of growth traits in Braunvieh cattle. Methods: The study included 300 genotyped animals by the GeneSeek® Genomic Profiler Bovine LDv.4 panel; after quality control, 22,734 SNP and 276 animals were maintained in the analysis. The examined phenotypic data considered birth (BW), weaning (WW), and yearling weights. The association analysis was performed using the principal components method via the egscore function of the GenABEL version 1.8-0 package in the R environment. The marker rs133262280 located in BTA 22 was associated with BW, and two SNPs were associated with WW, rs43668789 (BTA 11) and rs136155567 (BTA 27). New QTL associated with these liveweight traits and four positional and functional candidate genes potentially involved in variations of the analyzed traits were identified. The most important genes in these genomic regions were MCM2 (minichromosome maintenance complex component 2), TPRA1 (transmembrane protein adipocyte associated 1), GALM (galactose mutarotase), and NRG1 (neuregulin 1), related to embryonic cleavage, bone and tissue growth, cell adhesion, and organic development. This study is the first to present a GWAS conducted in Braunvieh cattle in Mexico providing evidence for genetic architecture of assessed growth traits. Further specific analysis of found associated genes and regions will clarify its contribution to the genetic basis of growth-related traits.

Genetic differentiation of mainland-island sheep of Greece: Implications for identifying candidate genes for long-term local adaptation.

In Greece, a number of local sheep breeds are raised in a wide range of ecological niches across the country. These breeds can be used for the identification of genetic variants that contribute to local adaptation. To this end, 50k genotypes of 90 local sheep from mainland Greece (Epirus, n = 35 and Peloponnesus, n = 55) were used, as well as 147 genotypes of sheep from insular Greece (Skyros, n = 21), Lemnos, n = 36 and Lesvos, n = 90). Principal components and phylogenetic analysis along with admixture and spatial point patterns analyses suggested genetic differentiation of ‘mainland-island’ populations. Genome scans for signatures of selection and genome-wide association analysis (GWAS) pointed to one highly differentiating marker on OAR4 (FST = 0.39, FLK = 21.93, FDR p-value = 0.10) that also displayed genome wide significance (FDR p-value = 0.002) during GWAS. A total number of 6 positional candidate genes (LOC106990429, ZNF804B, TEX47, STEAP4, SRI and ADAM22) were identified within 500 kb flanking regions around the significant marker. In addition, two QTLs related to fat tail deposition are reported in genomic regions 800 kb downstream the significant marker. Based on gene ontology analysis and literature evidence, the identified candidate genes possess biological functions relevant to local adaptation that worth further investigation.

Whole-genome association study searching for QTL for Aeromonas salmonicida resistance in rainbow trout

Aeromonas salmonicida subsp. salmonicida, the causative agent of furunculosis, has extensive negative effects on wild and farmed salmonids worldwide. Vaccination induces some protection under certain conditions but disease outbreaks occur even in vaccinated fish. Therefore, alternative disease control approaches are required to ensure the sustainable expansion of rainbow trout aquaculture. Selective breeding can be applied to enhance host resistance to pathogens. The present work used genome-wide association study (GWAS) to identify quantitative trait loci (QTL) associated with A. salmonicida resistance in rainbow trout. A total 798 rainbow trout exposed to A. salmonicida by bath challenge revealed 614 susceptible and 138 resistant fish. Genotyping was conducted using the 57 K single nucleotide polymorphism (SNP) array and the GWAS was performed for survival and time to death phenotypes. We identified a QTL on chromosome 16 and located positional candidate genes in the proximity of the most significant SNPs. In addition, samples from exposed fish were examined for expression of 24 immune-relevant genes indicating a systematic immune response to the infection. The present work demonstrated that resistance to A. salmonicida is moderately heritable with oligogenic architecture. These result will be useful for the future breeding programs for improving the natural resistance of rainbow trout against furunculosis.

Genome-wide association studies for canine hip dysplasia in single and multiple populations \u2013 implications and potential novel risk loci

BackgroundAssociation mapping studies of quantitative trait loci (QTL) for canine hip dysplasia (CHD) can contribute to the understanding of the genetic background of this common and debilitating disease and might contribute to its genetic improvement. The power of association studies for CHD is limited by relatively small sample numbers for CHD records within countries, suggesting potential benefits of joining data across countries. However, this is complicated due to the use of different scoring systems across countries. In this study, we incorporated routinely assessed CHD records and genotype data of German Shepherd dogs from two countries (UK and Sweden) to perform genome-wide association studies (GWAS) within populations using different variations of CHD phenotypes. As phenotypes, dogs were either classified into cases and controls based on the Fédération Cynologique Internationale (FCI) five-level grading of the worst hip or the FCI grade was treated as an ordinal trait. In a subsequent meta-analysis, we added publicly available data from a Finnish population and performed the GWAS across all populations. Genetic associations for the CHD phenotypes were evaluated in a linear mixed model using 62,089 SNPs.ResultsMultiple SNPs with genome-wide significant and suggestive associations were detected in single-population GWAS and the meta-analysis. Few of these SNPs overlapped between populations or between single-population GWAS and the meta-analysis, suggesting that many CHD-related QTL are population-specific. More significant or suggestive SNPs were identified when FCI grades were used as phenotypes in comparison to the case-control approach. MED13 (Chr 9) and PLEKHA7 (Chr 21) emerged as novel positional candidate genes associated with hip dysplasia.ConclusionsOur findings confirm the complex genetic nature of hip dysplasia in dogs, with multiple loci associated with the trait, most of which are population-specific. Routinely assessed CHD information collected across countries provide an opportunity to increase sample sizes and statistical power for association studies. While the lack of standardisation of CHD assessment schemes across countries poses a challenge, we showed that conversion of traits can be utilised to overcome this obstacle.