SNPs in SNCA,MCCC1, DLG2, GBF1 and MBNL2 are associated with Parkinson's disease in southern Chinese population

Numerous single nucleotide polymorphisms (SNPs), which have been identified as susceptibility factors for Parkinson's disease (PD) as per genome‐wide association studies, have not been fully characterized for PD patients in China. This study aimed to replicate the relationship between 12 novel SNPs of 12 genes and PD risk in southern Chinese population. Twelve SNPs of 12 genes were detected in 231 PD patients and 249 controls, using the SNaPshot technique. Meta‐analysis was used to assess heterogeneity of effect sizes between this study and published data. The impact of SNPs on gene expression was investigated by analysing the SNP‐gene association in the expression quantitative trait loci (eQTL) data sets. rs8180209 of SNCA (allele model: P = .047, OR = 0.77; additive model: P = .047, OR = 0.77), rs2270968 of MCCC1 (dominant model: P = .024, OR = 1.52), rs7479949 of DLG2 (recessive model; P = .019, OR = 1.52), rs10748818 of GBF1 (additive model: P < .001, OR = 0.37), and rs4771268 of MBNL2 (recessive model: P = .003, OR = 0.48) were replicated to be significantly associated with the increased risk of PD. Noteworthy, a meta‐analysis of previous studies suggested rs8180209, rs2270968, rs7479949 and rs4771268 were in line with those of our cohort. Our study replicated five novel functional SNPs in SNCA, MCCC1, DLG2, GBF1 and MBNL2 could be associated with increased risk of PD in southern Chinese population.     

Genetic markers that influence feed efficiency phenotypes also affect cattle temperament as measured by flight speed

Flight speed is a predictive indicator of cattle temperament and is associated with feed efficiency phenotypes. Genetic markers associated with both traits may assist with selection of calmer animals with improved economic value. A preliminary genome‐wide association study determined chromosomal regions on BTA9, and 17 were associated with flight speed. The genes quaking (QKI), glutamate receptor, ionotropic, AMPA 2 (GRIA2) and glycine receptor β (GLRB) were identified in these regions as potential functional candidates. Beef steers (n = 1057) were genotyped with SNPs located within and flanking these genes. One SNP located near QKI and one near GRIA2 were nominally associated with flight speed (P ≤ 0.05) although neither was significant after Bonferroni correction. Several studies have shown a correlation between flight speed and feed intake or gain; therefore, we also analyzed SNPs on BTA6:38–39 Mb known to be associated with average daily gain (ADG) and average daily feed intake (ADFI) for association with flight speed. Several SNPs on BTA6 were associated with flight speed (P ≤ 0.005), and three were significant after Bonferroni correction. These results suggest that the genes tested are unlikely to contribute to flight speed variation for our cattle population, but SNPs on BTA6 associated with ADG and ADFI may influence temperament. Use of these markers to select for economically important feed efficiency phenotypes may produce cattle with more desirable temperaments.     

Two missense mutations in exon 9 of caprine PRLR gene were associated with litter size

Guanzhong (n = 321) and Boer (n = 191) goat breeds were used to detect single nucleotide polymorphisms (SNPs) in the coding regions of the prolactin receptor (PRLR) gene by DNA sequencing and PCR‐RFLP. Two SNPs (c.1457G>A and c.1645G>A) were identified that caused amino acid variations p.Ser485Asn and p.Val548Met respectively. Statistical results indicated that the c.1457G>A and c.1645G>A SNPs were significantly associated with litter size in Boer and Guanzhong goat breeds. Further analysis revealed that combined genotype C4 (GGGG) and haplotype G‐G were better than the others for litter size in both goat breeds. These results might contribute to goat genetic resources and breeding.     

Development and evaluation of a genome‐wide Coffee 8.5K SNP array and its application for high‐density genetic mapping and for investigating the origin of Coffea arabica L.

Coffee species such as Coffea canephora P. (Robusta) and C. arabica L. (Arabica) are important cash crops in tropical regions around the world. C. arabica is an allotetraploid (2n = 4x = 44) originating from a hybridization event of the two diploid species C. canephora and C. eugenioides (2n = 2x = 22). Interestingly, these progenitor species harbour a greater level of genetic variability and are an important source of genes to broaden the narrow Arabica genetic base. Here, we describe the development, evaluation and use of a single‐nucleotide polymorphism (SNP) array for coffee trees. A total of 8580 unique and informative SNPs were selected from C. canephora and C. arabica sequencing data, with 40% of the SNP located in annotated genes. In particular, this array contains 227 markers associated to 149 genes and traits of agronomic importance. Among these, 7065 SNPs (~82.3%) were scorable and evenly distributed over the genome with a mean distance of 54.4 Kb between markers. With this array, we improved the Robusta high‐density genetic map by adding 1307 SNP markers, whereas 945 SNPs were found segregating in the Arabica mapping progeny. A panel of C. canephora accessions was successfully discriminated and over 70% of the SNP markers were transferable across the three species. Furthermore, the canephora‐derived subgenome of C. arabica was shown to be more closely related to C. canephora accessions from northern Uganda than to other current populations. These validated SNP markers and high‐density genetic maps will be useful to molecular genetics and for innovative approaches in coffee breeding.     

Association between putative functional variants in the PSMB9 gene and risk of melanoma – re‐analysis of published melanoma genome‐wide association studies

To mine possibly hidden causal single‐nucleotide polymorphisms (SNPs) of melanoma, we investigated the association of SNPs in 76 M/G1 transition genes with melanoma risk using our published genome‐wide association study (GWAS) data set with 1804 melanoma cases and 1026 cancer‐free controls. We found multiple SNPs with P < 0.01 and performed validation studies for 18 putative functional SNPs in PSMB9 in two other GWAS data sets. Two SNPs (rs1351383 and rs2127675) were associated with melanoma risk in the GenoMEL data set (P = 0.013 and 0.004, respectively), but failed in validation using the Australian data set. Genotype–phenotype analysis revealed these two SNPs were significantly correlated with mRNA expression level of PSMB9. Further experiments revealed that SNP rs2071480, which is in high LD with rs1351383 and rs2127675, may have a weak effect on the promoter activity of PSMB9. Taken together, our data suggested that functional variants in PSMB9 may contribute to melanoma susceptibility.     

Genome‐wide detection of genetic loci and candidate genes for teat number and body conformation traits at birth in Chinese Sushan pigs

Body conformation at birth and teat number are economically important traits in the pig industry, as these traits are usually explored to evaluate the growth and reproductive potential of piglets. To detect genetic loci and candidate genes for these traits, we performed a GWAS on 269 pigs from a recently developed Chinese breed (Sushan) using 38  128 informative SNPs on the Affymetrix Porcine SNP 55K Array. In total, we detected one genome‐wide significant (P = 1.31e‐6) SNP for teat number on chromosome X and 15 chromosome‐wide significant SNPs for teat number, body weight, body length, chest circumference and cannon circumference at birth on chromosomes 1, 3, 4, 6, 7, 9, 10, 13, 14, 15, 17 and 18. The most significant SNP had an additive effect of 0.74 × total teat number, explaining 20% of phenotypic variance. Five significant SNPs resided in the previously reported quantitative trait loci for these traits and seven significant SNPs had a pleiotropic effect on multiple traits. Intriguingly, 12 of the genes nearest to the significant SNPs are functionally related to body conformation and teat number traits, including SPRED2, MKX, TMSB4X and ESR1. GO analysis revealed that candidate genes proximal to the significant SNPs were enriched in the G‐protein coupled receptor and steroid hormone‐mediated signaling pathway. Our findings shed light on the genetic basis of the measured traits and provide molecular markers especially for the genetic improvement of teat number in Sushan and related pigs.     

SNP-based heritability and genetic architecture of cranial cruciate ligament rupture in Labrador Retrievers

Cranial cruciate ligament rupture (CCLR) is one of the leading causes of pelvic limb lameness in dogs. About 6% of Labrador Retrievers suffer from this orthopedic problem. The aim of this study was to determine the heritability of CCLR in this breed using SNP array genotyping data. DNA samples were collected from CCLR-affected dogs (n = 190) and unaffected dogs over the age of 8 years (n = 143). All 333 dogs were genotyped directly or imputed up to approximately 710k SNPs on the Affymetrix Axiom CanineHD SNP array. Heritability of CCLR was calculated using multiple methodologies, including linear mixed models, Bayesian models and a model that incorporates LD. The covariates of sex and sterilization status were added to each analysis to assess their impact. Across the algorithms of these models, heritability ranged from 0.550 to 0.886, depending on covariate inclusion. The relatively high heritability for this disease indicates that a substantial genetic component contributes to CCLR in the Labrador Retriever.     

Comparative genomics analysis in Prunoideae to identify biologically relevant polymorphisms

Prunus is an economically important genus with a wide range of physiological and biological variability. Using the peach genome as a reference, sequencing reads from four almond accessions and one sweet cherry cultivar were used for comparative analysis of these three Prunus species. Reference mapping enabled the identification of many biological relevant polymorphisms within the individuals. Examining the depth of the polymorphisms and the overall scaffold coverage, we identified many potentially interesting regions including hundreds of small scaffolds with no coverage from any individual. Non‐sense mutations account for about 70 000 of the 13 million identified single nucleotide polymorphisms (SNPs). Blast2GO analyses on these non‐sense SNPs revealed several interesting results. First, non‐sense SNPs were not evenly distributed across all gene ontology terms. Specifically, in comparison with peach, sweet cherry is found to have non‐sense SNPs in two 1‐aminocyclopropane‐1‐carboxylate synthase (ACS) genes and two 1‐aminocyclopropane‐1‐carboxylate oxidase (ACO) genes. These polymorphisms may be at the root of the nonclimacteric ripening of sweet cherry. A set of candidate genes associated with bitterness in almond were identified by comparing sweet and bitter almond sequences. To the best of our knowledge, this is the first report in plants of non‐sense SNP abundance in a genus being linked to specific GO terms.     

Identification of genetic associations of ECHS1 gene with milk fatty acid traits in dairy cattle

Our previous genome‐wide association study identified 83 genome‐wide significant SNPs and 20 novel promising candidate genes for milk fatty acids in Chinese Holstein. Among them, the enoyl‐CoA hydratase, short chain 1 (ECHS1) and enoyl‐CoA hydratase and 3‐hydroxyacyl CoA dehydrogenase (EHHADH) genes were located near two SNPs and one SNP respectively, and they play important roles in fatty acid metabolism pathways. We herein validated whether the two genes have genetic effects on milk fatty acid traits in dairy cattle. By re‐sequencing the full‐length coding region, partially adjacent introns and 3000 bp up/downstream flanking sequences, we identified 12 SNPs in ECHS1: two in exons, four in the 3′ flanking region and six in introns. The g.25858322C>T SNP results in an amino acid replacement from leucine to phenylalanine and changes the secondary structure of the ECHS1 protein, and single‐locus association analysis showed that it was significantly associated with three milk fatty acids (P = 0.0002–0.0013). The remaining 11 SNPs were found to be significantly associated with at least one milk fatty acid (P = <0.0001–0.0040). Also, we found that two haplotype blocks, consisting of nine and two SNPs respectively, were significantly associated with eight milk fatty acids (P = <0.0001–0.0125). However, none of polymorphisms was observed in the EHHADH gene. In conclusion, our findings are the first to indicate that the ECHS1 gene has a significant genetic impact on long‐chain unsaturated and medium‐chain saturated fatty acid traits in dairy cattle, although the biological mechanism is still undetermined and requires further in‐depth validation.     

Genetic signatures of natural selection in response to air pollution in red spruce (Picea rubens,Pinaceae)

One of the most important drivers of local adaptation for forest trees is climate. Coupled to these patterns, however, are human‐induced disturbances through habitat modification and pollution. The confounded effects of climate and disturbance have rarely been investigated with regard to selective pressure on forest trees. Here, we have developed and used a population genetic approach to search for signals of selection within a set of 36 candidate genes chosen for their putative effects on adaptation to climate and human‐induced air pollution within five populations of red spruce (Picea rubens Sarg.), distributed across its natural range and air pollution gradient in eastern North America. Specifically, we used F_ST outlier and environmental correlation analyses to highlight a set of seven single nucleotide polymorphisms (SNPs) that were overly correlated with climate and levels of sulphate pollution after correcting for the confounding effects of population history. Use of three age cohorts within each population allowed the effects of climate and pollution to be separated temporally, as climate‐related SNPs (n = 7) showed the strongest signals in the oldest cohort, while pollution‐related SNPs (n = 3) showed the strongest signals in the youngest cohorts. These results highlight the usefulness of population genetic scans for the identification of putatively nonneutral evolution within genomes of nonmodel forest tree species, but also highlight the need for the development and application of robust methodologies to deal with the inherent multivariate nature of the genetic and ecological data used in these types of analyses.     

Genome‐wide analyses suggest parallel selection for universal traits may eclipse local environmental selection in a highly mobile carnivore

Ecological and environmental heterogeneity can produce genetic differentiation in highly mobile species. Accordingly, local adaptation may be expected across comparatively short distances in the presence of marked environmental gradients. Within the European continent, wolves (Canis lupus) exhibit distinct north–south population differentiation. We investigated more than 67‐K single nucleotide polymorphism (SNP) loci for signatures of local adaptation in 59 unrelated wolves from four previously identified population clusters (northcentral Europe n = 32, Carpathian Mountains n = 7, Dinaric‐Balkan n = 9, Ukrainian Steppe n = 11). Our analyses combined identification of outlier loci with findings from genome‐wide association study of individual genomic profiles and 12 environmental variables. We identified 353 candidate SNP loci. We examined the SNP position and neighboring megabase (1 Mb, one million bases) regions in the dog (C. lupus familiaris) genome for genes potentially under selection, including homologue genes in other vertebrates. These regions included functional genes for, for example, temperature regulation that may indicate local adaptation and genes controlling for functions universally important for wolves, including olfaction, hearing, vision, and cognitive functions. We also observed strong outliers not associated with any of the investigated variables, which could suggest selective pressures associated with other unmeasured environmental variables and/or demographic factors. These patterns are further supported by the examination of spatial distributions of the SNPs associated with universally important traits, which typically show marked differences in allele frequencies among population clusters. Accordingly, parallel selection for features important to all wolves may eclipse local environmental selection and implies long‐term separation among population clusters.     

A targeted genotyping approach enhances identification of variants in taste receptor and appetite/reward genes of potential functional importance for obesity‐related porcine traits

Taste receptors (TASRs) and appetite and reward (AR) mechanisms influence eating behaviour, which in turn affects food intake and risk of obesity. In a previous study, we used next generation sequencing to identify potentially functional mutations in TASR and AR genes and found indications for genetic associations between identified variants and growth and fat deposition in a subgroup of animals (n = 38) from the UNIK resource pig population. This population was created for studying obesity and obesity‐related diseases. In the present study we validated results from our previous study by investigating genetic associations between 24 selected single nucleotide variants in TASR and AR gene variants and 35 phenotypes describing obesity and metabolism in the entire UNIK population (n = 564). Fifteen variants showed significant association with specific obesity‐related phenotypes after Bonferroni correction. Six of the 15 genes, namely SIM1, FOS, TAS2R4, TAS2R9, MCHR2 and LEPR, showed good correlation between known biological function and associated phenotype. We verified a genetic association between potentially functional variants in TASR/AR genes and growth/obesity and conclude that the combination of identification of potentially functional variants by next generation sequencing followed by targeted genotyping and association studies is a powerful and cost‐effective approach for increasing the power of genetic association studies.     

Development and preliminary evaluation of a genomewide single nucleotide polymorphisms resource generated by RAD‐seq for the small yellow croaker (Larimichthys polyactis)

Recent advances in high‐throughput sequencing technologies have offered the possibility to generate genomewide sequence data to delineate previously unidentified genetic structure, obtain more accurate estimates of demographic parameters and to evaluate potential adaptive divergence. Here, we identified 27 556 single nucleotide polymorphisms for the small yellow croaker (Larimichthys polyactis) using restriction‐site‐associated DNA (RAD) sequencing of 24 individuals from two populations. Significant sources of genetic variation were identified, with an average nucleotide diversity (π) of 0.00105 ± 0.000425 across individuals, and long‐term effective population size was thus estimated to range between 26 172 and 261 716. According to the results, no differentiation between the two populations was detected based on the SNP data set of top quality score per contig or neutral loci. However, the two analysed populations were highly differentiated based on SNP data set of both top F_ST value per contig and the outlier SNPs. Moreover, local adaptation was highlighted by an F_ST‐based outlier tests implemented in LOSITAN and a total of 538 potentially locally selected SNPs were identified. blast2go annotation of contigs containing the outlier SNPs yielded hits for 37 (66%) of 56 significant blastx matches. Candidate genes for local adaptation constituted a wide array of biological functions, including cellular response to oxidative stress, actin filament binding, ion transmembrane transport and synapse assembly. The generated SNP resources in this study provided a valuable tool for future population genetics and genomics studies of L. polyactis.     

Diversity and population structure of northern switchgrass as revealed through exome capture sequencing

Panicum virgatum L. (switchgrass) is a polyploid, perennial grass species that is native to North America, and is being developed as a future biofuel feedstock crop. Switchgrass is present primarily in two ecotypes: a northern upland ecotype, composed of tetraploid and octoploid accessions, and a southern lowland ecotype, composed of primarily tetraploid accessions. We employed high‐coverage exome capture sequencing (~2.4 Tb) to genotype 537 individuals from 45 upland and 21 lowland populations. From these data, we identified ~27 million single‐nucleotide polymorphisms (SNPs), of which 1 590 653 high‐confidence SNPs were used in downstream analyses of diversity within and between the populations. From the 66 populations, we identified five primary population groups within the upland and lowland ecotypes, a result that was further supported through genetic distance analysis. We identified conserved, ecotype‐restricted, non‐synonymous SNPs that are predicted to affect the protein function of CONSTANS (CO) and EARLY HEADING DATE 1 (EHD1), key genes involved in flowering, which may contribute to the phenotypic differences between the two ecotypes. We also identified, relative to the near‐reference Kanlow population, 17 228 genes present in more copies than in the reference genome (up‐CNVs), 112 630 genes present in fewer copies than in the reference genome (down‐CNVs) and 14 430 presence/absence variants (PAVs), affecting a total of 9979 genes, including two upland‐specific CNV clusters. In total, 45 719 genes were affected by an SNP, CNV, or PAV across the panel, providing a firm foundation to identify functional variation associated with phenotypic traits of interest for biofuel feedstock production.     

Single nucleotide polymorphism discovery in albacore and Atlantic bluefin tuna provides insights into worldwide population structure

The optimal management of the commercially important, but mostly over‐exploited, pelagic tunas, albacore (Thunnus alalunga Bonn., 1788) and Atlantic bluefin tuna (BFT; Thunnus thynnus L., 1758), requires a better understanding of population structure than has been provided by previous molecular methods. Despite numerous studies of both species, their population structures remain controversial. This study reports the development of single nucleotide polymorphisms (SNPs) in albacore and BFT and the application of these SNPs to survey genetic variability across the geographic ranges of these tunas. A total of 616 SNPs were discovered in 35 albacore tuna by comparing sequences of 54 nuclear DNA fragments. A panel of 53 SNPs yielded F_ST values ranging from 0.0 to 0.050 between samples after genotyping 460 albacore collected throughout the distribution of this species. No significant heterogeneity was detected within oceans, but between‐ocean comparisons (Atlantic, Pacific and Indian oceans along with Mediterranean Sea) were significant. Additionally, a 17‐SNP panel was developed in Atlantic BFT by cross‐species amplification in 107 fish. This limited number of SNPs discriminated between samples from the two major spawning areas of Atlantic BFT (F_ST = 0.116). The SNP markers developed in this study can be used to genotype large numbers of fish without the need for standardizing alleles among laboratories.     

Polymorphisms in the AKT1 and AKT2 genes and oesophageal squamous cell carcinoma risk in an Eastern Chinese population

Ethnic Han Chinese are at high risk of developing oesophageal squamous cell carcinoma (ESCC). Aberrant activation of the AKT signalling pathway is involved in many cancers, including ESCC. Some single nucleotide polymorphisms (SNPs) in genes involved in this pathway may contribute to ESCC susceptibility. We selected five potentially functional SNPs in AKT1 (rs2494750, rs2494752 and rs10138277) and AKT2 (rs7254617 and rs2304186) genes and investigated their associations with ESCC risk in 1117 ESCC cases and 1096 controls in an Eastern Chinese population. None of individual SNPs exhibited an association with ESCC risk. However, the combined analysis of three AKT1 SNPs suggested that individuals carrying one of AKT1 variant genotypes had a decreased ESCC risk [adjusted odds ratio (OR) = 0.60, 95% CI = 0.42–0.87]. Further stratified analysis found that AKT1 rs2294750 SNP was associated with significantly decreased ESCC risk among women (adjusted OR = 0.63, 95% CI = 0.43–0.94) and non‐drinkers (OR = 0.79, 95% CI = 0.64–0.99). Similar protective effects on women (adjusted OR = 0.56, 95% CI = 0.37–0.83) and non‐drinker (adjusted OR = 0.75, 95% CI = 0.60–0.94) were also observed for the combined genotypes of AKT1 SNPs. Consistently, logistic regression analysis indicated significant gene–gene interactions among three AKT1 SNPs (P < 0.015). A three‐AKT1 SNP haplotype (C‐A‐C) showed a significant association with a decreased ESCC risk (adjusted OR = 0.70, 95% CI = 0.52–0.94). Multifactor dimensionality reduction analysis confirmed a high‐order gene–environment interaction in ESCC risk. Overall, we found that three AKT1 SNPs might confer protection against ESCC risk; nevertheless, these effects may be dependent on other risk factors. Our results provided evidence of important gene–environment interplay in ESCC carcinogenesis.     

Genetic variation at KIT locus may predispose to melanoma

As loss of KIT frequently occurs in melanoma progression, we hypothesized that KIT is implicated in predisposition to melanoma (MM). Thus, we sequenced the KIT coding region in 112 familial MM cases and 143 matched controls and genotyped tag single‐nucleotide polymorphisms (SNPs) in two cohorts of melanoma patients and matched controls. Five rare KIT substitutions, all predicted possibly or probably deleterious, were identified in five patients, but none in controls [RR = 2.26 (1.26–2.26)]. Expressed in melanocyte lines, three substitutions inhibited KIT signaling. Comparison with exomes database (7020 alleles) confirmed a significant excess of rare deleterious KIT substitutions in patients. Additionally, a common SNP, rs2237028, was associated with MM risk, and 6 KIT variants were associated with nevus count. Our data strongly suggest that rare KIT substitutions predispose to melanoma and that common variants at KIT locus may also impact nevus count and melanoma risk.