Whole-genome genotyping on bead arrays

Evidence for Multilevel Chemopreventive Activities of Natural Phenols from Functional Genomic Studies of Curcumin, Resveratrol, Genistein, Quercetin, and Luteolin

Herein, I present an updated and contextualized literature review of functional genomic studies of natural phenols in the context of cancer. I suggest multilevel chemopreventive and anticancer mechanisms of action, which are shared by multiple dietary natural phenols. Specifically, I cite evidence that curcumin and resveratrol have multilevel anti-cancer effects through: (1) inducing either p53-dependent or p53-independent apoptosis in cancer cell lines, (2) acting as potent regulators of expression of oncogenic and anti-oncogenic microRNAs, and (3) inducing complex epigenetic changes that can switch off oncogenes/switch on anti-oncogenes. There is no simple reductionist explanation for anti-cancer effects of curcumin and resveratrol. More generally, multilevel models of chemoprevention are suggested for related natural phenols and flavonoids such as genistein, quercetin, or luteolin.

DNA methylome-wide association study of genetic risk for depression implicates antigen processing and immune responses

BACKGROUND

Depression is a disabling and highly prevalent condition where genetic and epigenetic, such as DNA methylation (DNAm), differences contribute to disease risk. DNA methylation is influenced by genetic variation but the association between polygenic risk of depression and DNA methylation is unknown.

METHODS

We investigated the association between polygenic risk scores (PRS) for depression and DNAm by conducting a methylome-wide association study (MWAS) in Generation Scotland (N = 8898, mean age = 49.8 years) with replication in the Lothian Birth Cohorts of 1921 and 1936 and adults in the Avon Longitudinal Study of Parents and Children (ALSPAC) (Ncombined = 2049, mean age = 79.1, 69.6 and 47.2 years, respectively). We also conducted a replication MWAS in the ALSPAC children (N = 423, mean age = 17.1 years). Gene ontology analysis was conducted for the cytosine-guanine dinucleotide (CpG) probes significantly associated with depression PRS, followed by Mendelian randomisation (MR) analysis to infer the causal relationship between depression and DNAm.

RESULTS

Widespread associations (NCpG = 71, pBonferroni < 0.05, p < 6.3 × 10-8) were found between PRS constructed using genetic risk variants for depression and DNAm in CpG probes that localised to genes involved in immune responses and neural development. The effect sizes for the significant associations were highly correlated between the discovery and replication samples in adults (r = 0.79) and in adolescents (r = 0.82). Gene Ontology analysis showed that significant CpG probes are enriched in immunological processes in the human leukocyte antigen system. Additional MWAS was conducted for each lead genetic risk variant. Over 47.9% of the independent genetic risk variants included in the PRS showed associations with DNAm in CpG probes located in both the same (cis) and distal (trans) locations to the genetic loci (pBonferroni < 0.045). Subsequent MR analysis showed that there are a greater number of causal effects found from DNAm to depression than vice versa (DNAm to depression: pFDR ranged from 0.024 to 7.45 × 10-30; depression to DNAm: pFDR ranged from 0.028 to 0.003).

CONCLUSIONS

PRS for depression, especially those constructed from genome-wide significant genetic risk variants, showed methylome-wide differences associated with immune responses. Findings from MR analysis provided evidence for causal effect of DNAm to depression.

Lifestyle and Genetic Factors Modify Parent-of-Origin Effects on the Human Methylome

BACKGROUND

parent-of-origin effects (POE) play important roles in complex disease and thus understanding their regulation and associated molecular and phenotypic variation are warranted. Previous studies mainly focused on the detection of genomic regions or phenotypes regulated by POE. Understanding whether POE may be modified by environmental or genetic exposures is important for understanding of the source of POE-associated variation, but only a few case studies addressing modifiable POE exist.

METHODS

in order to understand this high order of POE regulation, we screened 101 genetic and environmental factors such as 'predicted mRNA expression levels' of DNA methylation/imprinting machinery genes and environmental exposures. POE-mQTL-modifier interaction models were proposed to test the potential of these factors to modify POE at DNA methylation using data from Generation Scotland: The Scottish Family Health Study(N=2315).

FINDINGS

a set of vulnerable/modifiable POE-CpGs were identified (modifiable-POE-regulated CpGs, N=3). Four factors, 'lifetime smoking status' and 'predicted mRNA expression levels' of TET2, SIRT1 and KDM1A, were found to significantly modify the POE on the three CpGs in both discovery and replication datasets. We further identified plasma protein and health-related phenotypes associated with the methylation level of one of the identified CpGs.

INTERPRETATION

the modifiable POE identified here revealed an important yet indirect path through which genetic background and environmental exposures introduce their effect on DNA methylation, motivating future comprehensive evaluation of the role of these modifiers in complex diseases.

FUNDING

NSFC (81971270),H2020-MSCA-ITN(721815), Wellcome (204979/Z/16/Z,104036/Z/14/Z), MRC (MC_UU_00007/10, MC_PC_U127592696), CSO (CZD/16/6,CZB/4/276, CZB/4/710), SFC (HR03006), EUROSPAN (LSHG-CT-2006-018947), BBSRC (BBS/E/D/30002276), SYSU, Arthritis Research UK, NHLBI, NIH.

Towards next-generation personalization of tacrolimus treatment: a review on advanced diagnostic and therapeutic approaches

The benefit of personalized medicine is that it allows the customization of drug therapy - maximizing efficacy while avoiding side effects. Genetic polymorphisms are one of the major contributors to interindividual variability. Currently, the only gold standard for applying personalized medicine is dose titration. Because of technological advancements, converting genotypic data into an optimum dose has become easier than in earlier years. However, for many medications, determining a personalized dose may be difficult, leading to a trial-and-error method. On the other hand, the technologically oriented pharmaceutical industry has a plethora of smart drug delivery methods that are underutilized in customized medicine. This article elaborates the genetic polymorphisms of tacrolimus as case study, and extensively covers the diagnostic and therapeutic technologies which aid in the delivery of personalized tacrolimus treatment for better clinical outcomes, thereby providing a new strategy for implementing personalized medicine.

Development and validation of a 1 K sika deer (Cervus nippon) SNP Chip

Background

China is the birthplace of the deer family and the country with the most abundant deer resources. However, at present, China’s deer industry faces the problem that pure sika deer and hybrid deer cannot be easily distinguished. Therefore, the development of a SNP identification chip is urgently required.

Results

In this study, 250 sika deer, 206 red deer, 23 first-generation hybrid deer (F1), 20 s-generation hybrid deer (F2), and 20 third-generation hybrid deer (F3) were resequenced. Using the chromosome-level sika deer genome as the reference sequence, mutation detection was performed on all individuals, and a total of 130,306,923 SNP loci were generated. After quality control filtering was performed, the remaining 31,140,900 loci were confirmed. From molecular-level and morphological analyses, the sika deer reference population and the red deer reference population were established. The Fst values of all SNPs in the two reference populations were calculated. According to customized algorithms and strict screening principles, 1000 red deer-specific SNP sites were finally selected for chip design, and 63 hybrid individuals were determined to contain red deer-specific SNP loci. The results showed that the gene content of red deer gradually decreased in subsequent hybrid generations, and this decrease roughly conformed to the law of statistical genetics. Reaction probes were designed according to the screening sites. All candidate sites met the requirements of the Illumina chip scoring system. The average score was 0.99, and the MAF was in the range of 0.3277 to 0.3621. Furthermore, 266 deer (125 sika deer, 39 red deer, 56 F1, 29 F2,17 F3) were randomly selected for 1 K SNP chip verification. The results showed that among the 1000 SNP sites, 995 probes were synthesized, 4 of which could not be typed, while 973 loci were polymorphic. PCA, random forest and ADMIXTURE results showed that the 1 K sika deer SNP chip was able to clearly distinguish sika deer, red deer, and hybrid deer and that this 1 K SNP chip technology may provide technical support for the protection and utilization of pure sika deer species resources.

Conclusion

We successfully developed a low-density identification chip that can quickly and accurately distinguish sika deer from their hybrid offspring, thereby providing technical support for the protection and utilization of pure sika deer germplasm resources.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-021-00994-z.

Design and evaluation of a custom 50K Infinium SNP array for egg-type chickens

With the development of molecular genetics and high-throughput sequencing technology, genotyping arrays consisting of large numbers of SNP have raised great interest in animal and plant research. However, the application of commercial chicken 600K SNP arrays has varied in different populations of egg-type chickens. Moreover, their genotyping cost is too high for large-scale population applications. Herein, we independently developed a custom Illumina 50K BeadChip, named PhenoixChip-I, for egg-type chickens based on SNP from 479 sequenced individuals in 7 lines. We filtered and selected SNP with stringent criteria, such as high polymorphism, genome coverage, design score, and priorities. Finally, a total of 43,681 effective SNP successfully genotyped were included on our custom array. Approximately 14K SNP were previously reported to be associated with important economic traits in egg-type chickens. Subsequently, we verified the applicability and efficiency of the PhenoixChip-I SNP array from many aspects, including evaluating its use scientific research (population structure analysis and genome-wide association study) and the poultry breeding industry (genomic selection). The findings in our study will play a crucial role in accelerating the genetic improvement of egg-type chickens.

Typing and determination of SNP functional gene based on highly selective and signal-amplified fluorescence double-probe with the help of ExoIII nuclease and magnetic bead

We have developed a fluorescence double-probe detection system with signal amplification for simple typing and determination of single nucleotide polymorphism (SNP) functional gene based on non-sequence dependence of ExoIII nuclease on dsDNA and rapid separation of magnetic bead. Matched detected gene can cyclically release abundant fluorescence-labeled ssDNA from the probe and the corresponding measured fluorescence signal is amplified up to 6063 times. In this case, the probe cannot release the measured fluorescence signal for the point mutation gene and then the corresponding measured signal is inhibited. According to signal amplification and inhabitation of the probe, we proposed both an accurate genotyping approach with strong specificity and a sensitive determination approach with high selectivity for SNP functional gene. For qualitative genotyping, there are obvious genotype-based differences of measured fluorescence phenotypes among three kinds of the samples of the investigated SNP. The quantitative determinations of its wild-type gene and mutant gene have all a good linearity in the range from 0.5 to 500 pmol/L with the correlation coefficients R² of 0.9940 and 0.9911, and a high sensitivity with the detection limits of 0.11 and 0.20 pmol/L, respectively. Compared to the usual single-probe detection system, the developed double-probe system can achieve not only accurate genotyping but also the sensitive gene determination. Meanwhile, it is also a simple and reliable method for both quantitative and qualitative analysis of functional gene.

An actionable KCNH2 Long QT Syndrome variant detected by sequence and haplotype analysis in a population research cohort

The Viking Health Study Shetland is a population-based research cohort of 2,122 volunteer participants with ancestry from the Shetland Isles in northern Scotland. The high kinship and detailed phenotype data support a range of approaches for associating rare genetic variants, enriched in this isolate population, with quantitative traits and diseases. As an exemplar, the c.1750G > A; p.Gly584Ser variant within the coding sequence of the KCNH2 gene implicated in Long QT Syndrome (LQTS), which occurred once in 500 whole genome sequences from this population, was investigated. Targeted sequencing of the KCNH2 gene in family members of the initial participant confirmed the presence of the sequence variant and identified two further members of the same family pedigree who shared the variant. Investigation of these three related participants for whom single nucleotide polymorphism (SNP) array genotypes were available allowed a unique shared haplotype of 1.22 Mb to be defined around this locus. Searching across the full cohort for this haplotype uncovered two additional apparently unrelated individuals with no known genealogical connection to the original kindred. All five participants with the defined haplotype were shown to share the rare variant by targeted Sanger sequencing. If this result were verified in a healthcare setting, it would be considered clinically actionable, and has been actioned in relatives ascertained independently through clinical presentation. The General Practitioners of four study participants with the rare variant were alerted to the research findings by letters outlining the phenotype (prolonged electrocardiographic QTc interval). A lack of detectable haplotype sharing between c.1750G > A; p.Gly584Ser chromosomes from previously reported individuals from Finland and those in this study from Shetland suggests that this mutation has arisen more than once in human history. This study showcases the potential value of isolate population-based research resources for genomic medicine. It also illustrates some challenges around communication of actionable findings in research participants in this context.

Parent of origin genetic effects on methylation in humans are common and influence complex trait variation

Parent-of-origin effects (POE) exist when there is differential expression of alleles inherited from the two parents. A genome-wide scan for POE on DNA methylation at 639,238 CpGs in 5,101 individuals identifies 733 independent methylation CpGs potentially influenced by POE at a false discovery rate ≤ 0.05 of which 331 had not previously been identified. Cis and trans methylation quantitative trait loci (mQTL) regulate methylation variation through POE at 54% (399/733) of the identified POE-influenced CpGs. The combined results provide strong evidence for previously unidentified POE-influenced CpGs at 171 independent loci. Methylation variation at 14 of the POE-influenced CpGs is associated with multiple metabolic traits. A phenome-wide association analysis using the POE mQTL SNPs identifies a previously unidentified imprinted locus associated with waist circumference. These results provide a high resolution population-level map for POE on DNA methylation sites, their local and distant regulators and potential consequences for complex traits.

Genome-wide by environment interaction studies of depressive symptoms and psychosocial stress in UK Biobank and Generation Scotland

Stress is associated with poorer physical and mental health. To improve our understanding of this link, we performed genome-wide association studies (GWAS) of depressive symptoms and genome-wide by environment interaction studies (GWEIS) of depressive symptoms and stressful life events (SLE) in two UK population-based cohorts (Generation Scotland and UK Biobank). No SNP was individually significant in either GWAS, but gene-based tests identified six genes associated with depressive symptoms in UK Biobank (DCC, ACSS3, DRD2, STAG1, FOXP2 and KYNU; p < 2.77 × 10-6). Two SNPs with genome-wide significant GxE effects were identified by GWEIS in Generation Scotland: rs12789145 (53-kb downstream PIWIL4; p = 4.95 × 10-9; total SLE) and rs17070072 (intronic to ZCCHC2; p = 1.46 × 10-8; dependent SLE). A third locus upstream CYLC2 (rs12000047 and rs12005200, p < 2.00 × 10-8; dependent SLE) when the joint effect of the SNP main and GxE effects was considered. GWEIS gene-based tests identified: MTNR1B with GxE effect with dependent SLE in Generation Scotland; and PHF2 with the joint effect in UK Biobank (p < 2.77 × 10-6). Polygenic risk scores (PRSs) analyses incorporating GxE effects improved the prediction of depressive symptom scores, when using weights derived from either the UK Biobank GWAS of depressive symptoms (p = 0.01) or the PGC GWAS of major depressive disorder (p = 5.91 × 10-3). Using an independent sample, PRS derived using GWEIS GxE effects provided evidence of shared aetiologies between depressive symptoms and schizotypal personality, heart disease and COPD. Further such studies are required and may result in improved treatments for depression and other stress-related conditions.

A catalog of annotated high-confidence SNPs from exome capture and sequencing reveals highly polymorphic genes in Norway spruce (Picea abies)

BACKGROUND

Norway spruce [Picea abies (L.) Karst.] is ecologically and economically one of the most important conifer worldwide. Our main goal was to develop a large catalog of annotated high confidence gene SNPs that should sustain the development of genomic tools for the conservation of natural and domesticated genetic diversity resources, and hasten tree breeding efforts in this species.

RESULTS

Targeted sequencing was achieved by capturing P. abies exome with probes previously designed from the sequenced transcriptome of white spruce (Picea glauca (Moench) Voss). Capture efficiency was high (74.5%) given a high level of exome conservation between the two species. Using stringent criteria, we delimited a set of 61,771 high-confidence SNPs across 13,543 genes. To validate SNPs, a high-throughput genotyping array was developed for a subset of 5571 predicted SNPs representing as many different gene loci, and was used to genotype over 1000 trees. The estimated true positive rate of the resource was 84.2%, which was comparable with the genotyping success rate obtained for P. abies control SNPs recycled from previous genotyping efforts. We also analyzed SNP abundance across various gene functional categories. Several GO terms and gene families involved in stress response were found over-represented in highly polymorphic genes.

CONCLUSION

The annotated high-confidence SNP catalog developed herein represents a valuable genomic resource, being representative of over 13 K genes distributed across the P. abies genome. This resource should serve a variety of population genomics and breeding applications in Norway spruce.

High-density genetic mapping of a major QTL for resistance to multiple races of loose smut in a tetraploid wheat cross

Loose smut, caused by Ustilago tritici (Pers.) Rostr., is a systemic disease of tetraploid durum wheat (Triticum turgidum L.). Loose smut can be economically controlled by growing resistant varieties, making it important to find and deploy new sources of resistance. Blackbird, a variety of T. turgidum L. subsp. carthlicum (Nevski) A. Love & D. Love, carries a high level of resistance to loose smut. Blackbird was crossed with the loose smut susceptible durum cultivar Strongfield to produce a doubled haploid (DH) mapping population. The parents and progenies were inoculated with U. tritici races T26, T32 and T33 individually and as a mixture at Swift Current, Canada in 2011 and 2012 and loose smut incidence (LSI) was assessed. Genotyping of the DH population and parents using an Infinium iSelect 90K single nucleotide polymorphism (SNP) array identified 12,952 polymorphic SNPs. The SNPs and 426 SSRs (previously genotyped in the same population) were mapped to 16 linkage groups spanning 3008.4 cM at an average inter-marker space of 0.2 cM in a high-density genetic map. Composite interval mapping analysis revealed three significant quantitative trait loci (QTL) for loose smut resistance on chromosomes 3A, 6B and 7A. The loose smut resistance QTL on 6B (QUt.spa-6B.2) and 7A (QUt.spa-7A.2) were derived from Blackbird. Strongfield contributed the minor QTL on 3A (QUt.spa-3A.2). The resistance on 6B was a stable major QTL effective against all individual races and the mixture of the three races; it explained up to 74% of the phenotypic variation. This study is the first attempt in durum wheat to identify and map loose smut resistance QTL using a high-density genetic map. The QTL QUt.spa-6B.2 would be an effective source for breeding resistance to multiple races of the loose smut pathogen because it provides near-complete broad resistance to the predominant virulence on the Canadian prairies.

Genetic and environmental contributions to psychological resilience and coping

Background: Twin studies indicate that genetic and environmental factors contribute to both psychological resilience and coping style, but estimates of their relative molecular and shared environmental contributions are limited. The degree of overlap in the genetic architectures of these traits is also unclear.

Methods: Using data from a large population- and family-based cohort Generation Scotland (N = 8,734), we estimated the genetic and shared environmental variance components for resilience, task-, emotion-, and avoidance-oriented coping style in a linear mixed model (LMM). Bivariate LMM analyses were used to estimate the genetic correlations between these traits. Resilience and coping style were measured using the Brief Resilience Scale and Coping Inventory for Stressful Situations, respectively.

Results: The greatest proportion of the phenotypic variance in resilience remained unexplained, although significant contributions from common genetic variants and family-shared environment were found. Both task- and avoidance-oriented coping had significant contributions from common genetic variants, sibling- and couple-shared environments, variance in emotion-oriented coping was attributable to common genetic variants, family- and couple-shared environments. The estimated correlation between resilience and emotion-oriented coping was high for both common-variant-associated genetic effects (r _G = -0.79, se = 0.19), and for the additional genetic effects from the pedigree (r _K = -0.94, se = 0.30). Genetic correlations between resilience and task- and avoidance-oriented coping did not meet statistical significance.

Conclusions: Both genetics and shared environmental effects were major contributing factors to coping style, whilst the variance in resilience remains largely unexplained. Strong genetic overlap between resilience and emotion-oriented coping suggests a relationship whereby genetic factors that increase negative emotionality also lead to decreased resilience. We suggest that genome-wide family-based studies of resilience and coping may help to elucidate tractable methodologies to identify genetic architectures and modifiable environmental risk factors to protect against psychiatric illness, although further work with larger sample sizes is needed.

Do regional brain volumes and major depressive disorder share genetic architecture? A study of Generation Scotland (n=19 762), UK Biobank (n=24 048) and the English Longitudinal Study of Ageing (n=5766)

Major depressive disorder (MDD) is a heritable and highly debilitating condition. It is commonly associated with subcortical volumetric abnormalities, the most replicated of these being reduced hippocampal volume. Using the most recent published data from Enhancing Neuroimaging Genetics through Meta-analysis (ENIGMA) consortium's genome-wide association study of regional brain volume, we sought to test whether there is shared genetic architecture between seven subcortical brain volumes and intracranial volume (ICV) and MDD. We explored this using linkage disequilibrium score regression, polygenic risk scoring (PRS) techniques, Mendelian randomisation (MR) analysis and BUHMBOX. Utilising summary statistics from ENIGMA and Psychiatric Genomics Consortium, we demonstrated that hippocampal volume was positively genetically correlated with MDD (r_G=0.46, P=0.02), although this did not survive multiple comparison testing. None of the other six brain regions studied were genetically correlated and amygdala volume heritability was too low for analysis. Using PRS analysis, no regional volumetric PRS demonstrated a significant association with MDD or recurrent MDD. MR analysis in hippocampal volume and MDD identified no causal association, however, BUHMBOX analysis identified genetic subgrouping in GS:SFHS MDD cases only (P=0.00281). In this study, we provide some evidence that hippocampal volume and MDD may share genetic architecture in a subgroup of individuals, albeit the genetic correlation did not survive multiple testing correction and genetic subgroup heterogeneity was not replicated. In contrast, we found no evidence to support a shared genetic architecture between MDD and other regional subcortical volumes or ICV.

The Peach v2.0 release: high-resolution linkage mapping and deep resequencing improve chromosome-scale assembly and contiguity

BACKGROUND

The availability of the peach genome sequence has fostered relevant research in peach and related Prunus species enabling the identification of genes underlying important horticultural traits as well as the development of advanced tools for genetic and genomic analyses. The first release of the peach genome (Peach v1.0) represented a high-quality WGS (Whole Genome Shotgun) chromosome-scale assembly with high contiguity (contig L50 214.2 kb), large portions of mapped sequences (96%) and high base accuracy (99.96%). The aim of this work was to improve the quality of the first assembly by increasing the portion of mapped and oriented sequences, correcting misassemblies and improving the contiguity and base accuracy using high-throughput linkage mapping and deep resequencing approaches.

RESULTS

Four linkage maps with 3,576 molecular markers were used to improve the portion of mapped and oriented sequences (from 96.0% and 85.6% of Peach v1.0 to 99.2% and 98.2% of v2.0, respectively) and enabled a more detailed identification of discernible misassemblies (10.4 Mb in total). The deep resequencing approach fixed 859 homozygous SNPs (Single Nucleotide Polymorphisms) and 1347 homozygous indels. Moreover, the assembled NGS contigs enabled the closing of 212 gaps with an improvement in the contig L50 of 19.2%.

CONCLUSIONS

The improved high quality peach genome assembly (Peach v2.0) represents a valuable tool for the analysis of the genetic diversity, domestication, and as a vehicle for genetic improvement of peach and related Prunus species. Moreover, the important phylogenetic position of peach and the absence of recent whole genome duplication (WGD) events make peach a pivotal species for comparative genomics studies aiming at elucidating plant speciation and diversification processes.

The Peach v2.0 release: high-resolution linkage mapping and deep resequencing improve chromosome-scale assembly and contiguity

BACKGROUND

The availability of the peach genome sequence has fostered relevant research in peach and related Prunus species enabling the identification of genes underlying important horticultural traits as well as the development of advanced tools for genetic and genomic analyses. The first release of the peach genome (Peach v1.0) represented a high-quality WGS (Whole Genome Shotgun) chromosome-scale assembly with high contiguity (contig L50 214.2 kb), large portions of mapped sequences (96%) and high base accuracy (99.96%). The aim of this work was to improve the quality of the first assembly by increasing the portion of mapped and oriented sequences, correcting misassemblies and improving the contiguity and base accuracy using high-throughput linkage mapping and deep resequencing approaches.

RESULTS

Four linkage maps with 3,576 molecular markers were used to improve the portion of mapped and oriented sequences (from 96.0% and 85.6% of Peach v1.0 to 99.2% and 98.2% of v2.0, respectively) and enabled a more detailed identification of discernible misassemblies (10.4 Mb in total). The deep resequencing approach fixed 859 homozygous SNPs (Single Nucleotide Polymorphisms) and 1347 homozygous indels. Moreover, the assembled NGS contigs enabled the closing of 212 gaps with an improvement in the contig L50 of 19.2%.

CONCLUSIONS

The improved high quality peach genome assembly (Peach v2.0) represents a valuable tool for the analysis of the genetic diversity, domestication, and as a vehicle for genetic improvement of peach and related Prunus species. Moreover, the important phylogenetic position of peach and the absence of recent whole genome duplication (WGD) events make peach a pivotal species for comparative genomics studies aiming at elucidating plant speciation and diversification processes.

Exploration of haplotype research consortium imputation for genome-wide association studies in 20,032 Generation Scotland participants

Background

The Generation Scotland: Scottish Family Health Study (GS:SFHS) is a family-based population cohort with DNA, biological samples, socio-demographic, psychological and clinical data from approximately 24,000 adult volunteers across Scotland. Although data collection was cross-sectional, GS:SFHS became a prospective cohort due to of the ability to link to routine Electronic Health Record (EHR) data. Over 20,000 participants were selected for genotyping using a large genome-wide array.

Methods

GS:SFHS was analysed using genome-wide association studies (GWAS) to test the effects of a large spectrum of variants, imputed using the Haplotype Research Consortium (HRC) dataset, on medically relevant traits measured directly or obtained from EHRs. The HRC dataset is the largest available haplotype reference panel for imputation of variants in populations of European ancestry and allows investigation of variants with low minor allele frequencies within the entire GS:SFHS genotyped cohort.

Results

Genome-wide associations were run on 20,032 individuals using both genotyped and HRC imputed data. We present results for a range of well-studied quantitative traits obtained from clinic visits and for serum urate measures obtained from data linkage to EHRs collected by the Scottish National Health Service. Results replicated known associations and additionally reveal novel findings, mainly with rare variants, validating the use of the HRC imputation panel. For example, we identified two new associations with fasting glucose at variants near to Y_RNA and WDR4 and four new associations with heart rate at SNPs within CSMD1 and ASPH, upstream of HTR1F and between PROKR2 and GPCPD1. All were driven by rare variants (minor allele frequencies in the range of 0.08–1%). Proof of principle for use of EHRs was verification of the highly significant association of urate levels with the well-established urate transporter SLC2A9.

Conclusions

GS:SFHS provides genetic data on over 20,000 participants alongside a range of phenotypes as well as linkage to National Health Service laboratory and clinical records. We have shown that the combination of deeper genotype imputation and extended phenotype availability make GS:SFHS an attractive resource to carry out association studies to gain insight into the genetic architecture of complex traits.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-017-0414-4) contains supplementary material, which is available to authorized users.

A Combined Pathway and Regional Heritability Analysis Indicates NETRIN1 Pathway Is Associated With Major Depressive Disorder

BACKGROUND

Genome-wide association studies (GWASs) of major depressive disorder (MDD) have identified few significant associations. Testing the aggregation of genetic variants, in particular biological pathways, may be more powerful. Regional heritability analysis can be used to detect genomic regions that contribute to disease risk.

METHODS

We integrated pathway analysis and multilevel regional heritability analyses in a pipeline designed to identify MDD-associated pathways. The pipeline was applied to two independent GWAS samples [Generation Scotland: The Scottish Family Health Study (GS:SFHS, N = 6455) and Psychiatric Genomics Consortium (PGC:MDD) (N = 18,759)]. A polygenic risk score (PRS) composed of single nucleotide polymorphisms from the pathway most consistently associated with MDD was created, and its accuracy to predict MDD, using area under the curve, logistic regression, and linear mixed model analyses, was tested.

RESULTS

In GS:SFHS, four pathways were significantly associated with MDD, and two of these explained a significant amount of pathway-level regional heritability. In PGC:MDD, one pathway was significantly associated with MDD. Pathway-level regional heritability was significant in this pathway in one subset of PGC:MDD. For both samples the regional heritabilities were further localized to the gene and subregion levels. The NETRIN1 signaling pathway showed the most consistent association with MDD across the two samples. PRSs from this pathway showed competitive predictive accuracy compared with the whole-genome PRSs when using area under the curve statistics, logistic regression, and linear mixed model.

CONCLUSIONS

These post-GWAS analyses highlight the value of combining multiple methods on multiple GWAS data for the identification of risk pathways for MDD. The NETRIN1 signaling pathway is identified as a candidate pathway for MDD and should be explored in further large population studies.

Shared Genetics and Couple-Associated Environment Are Major Contributors to the Risk of Both Clinical and Self-Declared Depression

Background

Both genetic and environmental factors contribute to risk of depression, but estimates of their relative contributions are limited. Commonalities between clinically-assessed major depressive disorder (MDD) and self-declared depression (SDD) are also unclear.

Methods

Using data from a large Scottish family-based cohort (GS:SFHS, N = 19,994), we estimated the genetic and environmental variance components for MDD and SDD. The components representing the genetic effect associated with genome-wide common genetic variants (SNP heritability), the additional pedigree-associated genetic effect and non-genetic effects associated with common environments were estimated in a linear mixed model (LMM).

Findings

Both MDD and SDD had significant contributions from components representing the effect from common genetic variants, the additional genetic effect associated with the pedigree and the common environmental effect shared by couples. The estimate of correlation between SDD and MDD was high (r = 1.00, se = 0.20) for common-variant-associated genetic effect and lower for the additional genetic effect from the pedigree (r = 0.57, se = 0.08) and the couple-shared environmental effect (r = 0.53, se = 0.22).

Interpretation

Both genetics and couple-shared environmental effects were major factors influencing liability to depression. SDD may provide a scalable alternative to MDD in studies seeking to identify common risk variants. Rarer variants and environmental effects may however differ substantially according to different definitions of depression.

•

Shared genetics and couple-associated environment are major contributors to the risk of depression.

•

The common-variant-associated genetic correlation between MDD and SDD was very high (r = 1.00, se = 0.20).

•

Lower correlations were detected for pedigree-associated genetics and couple-shared environmental components.

It is important to understand the differences between clinical depression (MDD) and self-declared depression (SDD) for which there has been recent genetics progress. We found major contributions from genetics and couple-shared environment to both traits. There is a very high correlation between the traits associated with common genetic variants but a lower correlation for other (probably rarer) genetic variation. MDD and SDD also likely differ in their shared environmental risk factors. Thus for studies of common genetic variation, SDD is a potential alternative to MDD. In clinical practice and research, the spousal depression status should also be considered a risk indicator.

Dissection of major depressive disorder using polygenic risk scores for schizophrenia in two independent cohorts

Major depressive disorder (MDD) is known for its substantial clinical and suspected causal heterogeneity. It is characterized by low mood, psychomotor slowing and increased levels of the personality trait neuroticism; factors also associated with schizophrenia (SCZ). It is possible that some cases of MDD may have a substantial genetic loading for SCZ. The presence of SCZ-like MDD subgroups would be indicated by an interaction between MDD status and polygenic risk of SCZ on cognitive, personality and mood measures. Here, we hypothesized that higher SCZ polygenic risk would define larger MDD case-control differences in cognitive ability, and smaller differences in distress and neuroticism. Polygenic risk scores (PRSs) for SCZ and their association with cognitive variables, neuroticism, mood and psychological distress were estimated in a large population-based cohort (Generation Scotland: Scottish Family Health Study, GS:SFHS). The individuals were divided into those with, and without, depression (n=2587 and n=16 764, respectively) to test for the interactions between MDD status and schizophrenia risk. Replication was sought in UK Biobank (UKB; n=6049 and n=27 476 cases and controls, respectively). In both the cohorts, we found significant interactions between SCZ-PRS and MDD status for measures of psychological distress (βGS=-0.04, PGS=0.014 and βUKB=-0.09, PUKB⩽0.001 for GS:SFHS and UKB, respectively) and neuroticism (βGS=-0.04, PGS=0.002 and βUKB=-0.06, PUKB=0.023). In both the cohorts, there was a reduction of case-control differences on a background of higher genetic risk of SCZ. These findings suggest that depression on a background of high genetic risk for SCZ may show attenuated associations with distress and neuroticism. This may represent a causally distinct form of MDD more closely related to SCZ.

Assessing the genetic overlap between BMI and cognitive function

Obesity and low cognitive function are associated with multiple adverse health outcomes across the life course. They have a small phenotypic correlation (r=-0.11; high body mass index (BMI)-low cognitive function), but whether they have a shared genetic aetiology is unknown. We investigated the phenotypic and genetic correlations between the traits using data from 6815 unrelated, genotyped members of Generation Scotland, an ethnically homogeneous cohort from five sites across Scotland. Genetic correlations were estimated using the following: same-sample bivariate genome-wide complex trait analysis (GCTA)-GREML; independent samples bivariate GCTA-GREML using Generation Scotland for cognitive data and four other samples (n=20 806) for BMI; and bivariate LDSC analysis using the largest genome-wide association study (GWAS) summary data on cognitive function (n=48 462) and BMI (n=339 224) to date. The GWAS summary data were also used to create polygenic scores for the two traits, with within- and cross-trait prediction taking place in the independent Generation Scotland cohort. A large genetic correlation of -0.51 (s.e. 0.15) was observed using the same-sample GCTA-GREML approach compared with -0.10 (s.e. 0.08) from the independent-samples GCTA-GREML approach and -0.22 (s.e. 0.03) from the bivariate LDSC analysis. A genetic profile score using cognition-specific genetic variants accounts for 0.08% (P=0.020) of the variance in BMI and a genetic profile score using BMI-specific variants accounts for 0.42% (P=1.9 × 10(-7)) of the variance in cognitive function. Seven common genetic variants are significantly associated with both traits at P<5 × 10(-5), which is significantly more than expected by chance (P=0.007). All these results suggest there are shared genetic contributions to BMI and cognitive function.

Genetics of Coronary Artery Disease in Taiwan: A Cardiometabochip Study by the Taichi Consortium

By means of a combination of genome-wide and follow-up studies, recent large-scale association studies of populations of European descent have now identified over 46 loci associated with coronary artery disease (CAD). As part of the TAICHI Consortium, we have collected and genotyped 8556 subjects from Taiwan, comprising 5423 controls and 3133 cases with coronary artery disease, for 9087 CAD SNPs using the CardioMetaboChip. We applied penalized logistic regression to ascertain the top SNPs that contribute together to CAD susceptibility in Taiwan. We observed that the 9p21 locus contributes to CAD at the level of genome-wide significance (rs1537372, with the presence of C, the major allele, the effect estimate is -0.216, standard error 0.033, p value 5.8x10-10). In contrast to a previous report, we propose that the 9p21 locus is a single genetic contribution to CAD in Taiwan because: 1) the penalized logistic regression and the follow-up conditional analysis suggested that rs1537372 accounts for all of the CAD association in 9p21, and 2) the high linkage disequilibrium observed for all associated SNPs in 9p21. We also observed evidence for the following loci at a false discovery rate >5%: SH2B3, ADAMTS7, PHACTR1, GGCX, HTRA1, COL4A1, and LARP6-LRRC49. We also took advantage of the fact that penalized methods are an efficient approach to search for gene-by-gene interactions, and observed that two-way interactions between the PHACTR1 and ADAMTS7 loci and between the SH2B3 and COL4A1 loci contribute to CAD risk. Both the similarities and differences between the significance of these loci when compared with significance of loci in studies of populations of European descent underscore the fact that further genetic association of studies in additional populations will provide clues to identify the genetic architecture of CAD across all populations worldwide.

Polygenic risk for alcohol dependence associates with alcohol consumption, cognitive function and social deprivation in a population-based cohort

Alcohol dependence is frequently co‐morbid with cognitive impairment. The relationship between these traits is complex as cognitive dysfunction may arise as a consequence of heavy drinking or exist prior to the onset of dependence. In the present study, we tested the genetic overlap between cognitive abilities and alcohol dependence using polygenic risk scores (PGRS). We created two independent PGRS derived from two recent genome‐wide association studies (GWAS) of alcohol dependence (SAGE GWAS: n = 2750; Yale‐Penn GWAS: n = 2377) in a population‐based cohort, Generation Scotland: Scottish Family Health Study (GS:SFHS) (n = 9863). Data on alcohol consumption and four tests of cognitive function [Mill Hill Vocabulary (MHV), digit symbol coding, phonemic verbal fluency (VF) and logical memory] were available. PGRS for alcohol dependence were negatively associated with two measures of cognitive function: MHV (SAGE: P = 0.009, β = −0.027; Yale‐Penn: P = 0.001, β = −0.034) and VF (SAGE: P = 0.0008, β = −0.036; Yale‐Penn: P = 0.00005, β = −0.044). VF remained robustly associated after adjustment for education and social deprivation; however, the association with MHV was substantially attenuated. Shared genetic variants may account for some of the phenotypic association between cognitive ability and alcohol dependence. A significant negative association between PGRS and social deprivation was found (SAGE: P = 5.2 × 10⁻⁷, β = −0.054; Yale‐Penn: P = 0.000012, β = −0.047). Individuals living in socially deprived regions were found to carry more alcohol dependence risk alleles which may contribute to the increased prevalence of problem drinking in regions of deprivation. Future work to identify genes which affect both cognitive impairment and alcohol dependence will help elucidate biological processes common to both disorders.

Common polygenic risk for autism spectrum disorder (ASD) is associated with cognitive ability in the general population

Cognitive impairment is common among individuals diagnosed with autism spectrum disorder (ASD) and attention-deficit hyperactivity disorder (ADHD). It has been suggested that some aspects of intelligence are preserved or even superior in people with ASD compared with controls, but consistent evidence is lacking. Few studies have examined the genetic overlap between cognitive ability and ASD/ADHD. The aim of this study was to examine the polygenic overlap between ASD/ADHD and cognitive ability in individuals from the general population. Polygenic risk for ADHD and ASD was calculated from genome-wide association studies of ASD and ADHD conducted by the Psychiatric Genetics Consortium. Risk scores were created in three independent cohorts: Generation Scotland Scottish Family Health Study (GS:SFHS) (n=9863), the Lothian Birth Cohorts 1936 and 1921 (n=1522), and the Brisbane Adolescent Twin Sample (BATS) (n=921). We report that polygenic risk for ASD is positively correlated with general cognitive ability (beta=0.07, P=6 × 10(-7), r(2)=0.003), logical memory and verbal intelligence in GS:SFHS. This was replicated in BATS as a positive association with full-scale intelligent quotient (IQ) (beta=0.07, P=0.03, r(2)=0.005). We did not find consistent evidence that polygenic risk for ADHD was associated with cognitive function; however, a negative correlation with IQ at age 11 years (beta=-0.08, Z=-3.3, P=0.001) was observed in the Lothian Birth Cohorts. These findings are in individuals from the general population, suggesting that the relationship between genetic risk for ASD and intelligence is partly independent of clinical state. These data suggest that common genetic variation relevant for ASD influences general cognitive ability.

Polygenic risk for coronary artery disease is associated with cognitive ability in older adults

Background: Coronary artery disease (CAD) is associated with cognitive decrements and risk of later dementia, but it is not known if shared genetic factors underlie this association. We tested whether polygenic risk for CAD was associated with cognitive ability in community-dwelling cohorts of middle-aged and older adults.

Methods: Individuals from Generation Scotland: Scottish Family Health Study (GS:SFHS, N = 9865) and from the Lothian Birth Cohorts of 1921 (LBC1921, N  = 517) and 1936 (LBC1936, N  = 1005) provided cognitive data and genome-wide genotype data. Polygenic risk profile scores for CAD were calculated for all of the cohorts using the largest available genome-wide association studies (GWAS) data set, the CARDIoGRAM consortium (22 233 cases and 64 762 controls). Polygenic risk profile scores for CAD were then tested for their association with cognitive abilities in the presence and absence of manifest cardiovascular disease.

Results: A meta-analysis of all three cohorts showed a negative association between CAD polygenic risk and fluid cognitive ability (β = −0.022, P  = 0.016), verbal intelligence (β = −0.024, P  = 0.011) and memory (β = −0.021, P  = 0.028).

Conclusions: Increased polygenic risk for CAD is associated with lower cognitive ability in older adults. Common genetic variants may underlie some of the association between age-related cognitive decrements and the risk for CAD.

Development of two major resources for pea genomics: the GenoPea 13.2K SNP Array and a high-density, high-resolution consensus genetic map

Single nucleotide polymorphism (SNP) arrays represent important genotyping tools for innovative strategies in both basic research and applied breeding. Pea is an important food, feed and sustainable crop with a large (about 4.45 Gbp) but not yet available genome sequence. In the present study, 12 pea recombinant inbred line populations were genotyped using the newly developed GenoPea 13.2K SNP Array. Individual and consensus genetic maps were built providing insights into the structure and organization of the pea genome. Largely collinear genetic maps of 3918-8503 SNPs were obtained from all mapping populations, and only two of these exhibited putative chromosomal rearrangement signatures. Similar distortion patterns in different populations were noted. A total of 12 802 transcript-derived SNP markers placed on a 15 079-marker high-density, high-resolution consensus map allowed the identification of ohnologue-rich regions within the pea genome and the localization of local duplicates. Dense syntenic networks with sequenced legume genomes were further established, paving the way for the identification of the molecular bases of important agronomic traits segregating in the mapping populations. The information gained on the structure and organization of the genome from this research will undoubtedly contribute to the understanding of the evolution of the pea genome and to its assembly. The GenoPea 13.2K SNP Array and individual and consensus genetic maps are valuable genomic tools for plant scientists to strengthen pea as a model for genetics and physiology and enhance breeding.

Development of highly reliable in silico SNP resource and genotyping assay from exome capture and sequencing: an example from black spruce (Picea mariana)

Picea mariana is a widely distributed boreal conifer across Canada and the subject of advanced breeding programmes for which population genomics and genomic selection approaches are being developed. Targeted sequencing was achieved after capturing P. mariana exome with probes designed from the sequenced transcriptome of Picea glauca, a distant relative. A high capture efficiency of 75.9% was reached although spruce has a complex and large genome including gene sequences interspersed by some long introns. The results confirmed the relevance of using probes from congeneric species to perform successfully interspecific exome capture in the genus Picea. A bioinformatics pipeline was developed including stringent criteria that helped detect a set of 97,075 highly reliable in silico SNPs. These SNPs were distributed across 14,909 genes. Part of an Infinium iSelect array was used to estimate the rate of true positives by validating 4267 of the predicted in silico SNPs by genotyping trees from P. mariana populations. The true positive rate was 96.2% for in silico SNPs, compared to a genotyping success rate of 96.7% for a set 1115 P. mariana control SNPs recycled from previous genotyping arrays. These results indicate the high success rate of the genotyping array and the relevance of the selection criteria used to delineate the new P. mariana in silico SNP resource. Furthermore, in silico SNPs were generally of medium to high frequency in natural populations, thus providing high informative value for future population genomics applications.

Mining microsatellites in the peach genome: development of new long-core SSR markers for genetic analyses in five Prunus species

A wide inventory of molecular markers is nowadays available for individual fingerprinting. Microsatellites, or simple sequence repeats (SSRs), play a relevant role due to their relatively ease of use, their abundance in the plant genomes, and their co-dominant nature, together with the availability of primer sequences in many important agricultural crops. Microsatellites with long-core motifs are more easily scored and were adopted long ago in human genetics but they were developed only in few crops, and Prunus species are not among them. In the present work the peach whole-genome sequence was used to select 216 SSRs containing long-core motifs with tri-, tetra- and penta-nucleotide repeats. Microsatellite primer pairs were designed and tested for polymorphism in the five diploid Prunus species of economic relevance (almond, apricot, Japanese plum, peach and sweet cherry). A set of 26 microsatellite markers covering all the eight chromosomes, was also selected and used in the molecular characterization, population genetics and structure analyses of a representative sample of the five diploid Prunus species, assessing their transportability and effectiveness. The combined probability of identity between two random individuals for the whole set of 26 SSRs was quite low, ranging from 2.30 × 10⁻⁷ in peach to 9.48 × 10⁻¹⁰ in almond, confirming the usefulness of the proposed set for fingerprinting analyses in Prunus species.

Major depressive disorder and current psychological distress moderate the effect of polygenic risk for obesity on body mass index

Major depressive disorder (MDD) and obesity are frequently co-morbid and this correlation is partly due to genetic factors. Although specific genetic risk variants are associated with body mass index (BMI) and with larger effect sizes in depressed individuals, the genetic overlap and interaction with depression has not been addressed using whole-genome data. Polygenic profile scores for MDD and BMI were created in 13,921 members of Generation Scotland: the Scottish Family Health Study and tested for their association with BMI, MDD, neuroticism and scores on the General Health Questionnaire (GHQ) (current psychological distress). The association between BMI polygenic profile scores and BMI was tested fitting GHQ, neuroticism or MDD status as an interaction term to test for a moderating effect of mood disorder. BMI polygenic profile scores were not associated with lifetime MDD status or neuroticism although a significant positive association with GHQ scores was found (P = 0.0001, β = 0.034, r(2) = 0.001). Polygenic risk for MDD was not associated with BMI. A significant interaction between BMI polygenic profile scores and MDD (P = 0.0003, β = 0.064), GHQ (P = 0.0005, β = 0.027) and neuroticism (P = 0.003, β = 0.023) was found when BMI was the dependent variable. The effect of BMI-increasing alleles was greater in those with MDD, high neuroticism or current psychological distress. MDD, neuroticism and current psychological distress amplify the effect of BMI polygenic profile scores on BMI. Depressed individuals with a greater polygenic load for obesity are at greater risk of becoming obese than control individuals.

On-demand and location selective particle assembly via electrophoretic deposition for fabricating structures with particle-to-particle precision

Programmable positioning of 2 μm polystyrene (PS) beads with single particle precision and location selective, "on-demand", particle deposition was demonstrated by utilizing patterned electrodes and electrophoretic deposition (EPD). An electrode with differently sized hole patterns, from 0.5 to 5 μm, was used to illustrate the discriminatory particle deposition events based on the voltage and particle-to-hole size ratio. With decreasing patterned hole size, a larger electric field was required for a particle deposition event to occur in that hole. For the 5 μm hole, particle deposition began to occur at 10 V/cm where as an electric field of 15 V/cm was required for particles to begin depositing in the 2 μm holes. The likelihood of particle depositions continued to increase for smaller sized holes as the electric field increased. Eventually, a monolayer of particles began to form at approximately 20 V/cm. In essence, a voltage threshold was found for each hole pattern of different sizes, allowing fine adjustments in pattern hole size and voltage to control when a particle deposition event took place, even with the patterns on the same electrode. This phenomenon opens a route toward controlled, multimaterial deposition and assembly onto substrates without repatterning of the electrode or complicated surface modification of the particles. An analytical approach using the theories for electrophoresis and dielectrophoresis found the former to be the dominating force for depositing a particle into a patterned hole. Ebeam lithography was used to pattern spherical holes in precise configurations onto electrode surfaces, where each hole accompanied a polystyrene (PS) particle placement and attachment during EPD. The versatility of e-beam lithography was utilized to create arbitrary pattern configurations to fabricate particle assemblies of limitless configurations, enabling fabrication of unique materials assemblies and interfaces.

Estimation of linkage disequilibrium and interspecific gene flow in Ficedula flycatchers by a newly developed 50k single-nucleotide polymorphism array

With the access to draft genome sequence assemblies and whole-genome resequencing data from population samples, molecular ecology studies will be able to take truly genome-wide approaches. This now applies to an avian model system in ecological and evolutionary research: Old World flycatchers of the genus Ficedula, for which we recently obtained a 1.1 Gb collared flycatcher genome assembly and identified 13 million single-nucleotide polymorphism (SNP)s in population resequencing of this species and its sister species, pied flycatcher. Here, we developed a custom 50K Illumina iSelect flycatcher SNP array with markers covering 30 autosomes and the Z chromosome. Using a number of selection criteria for inclusion in the array, both genotyping success rate and polymorphism information content (mean marker heterozygosity = 0.41) were high. We used the array to assess linkage disequilibrium (LD) and hybridization in flycatchers. Linkage disequilibrium declined quickly to the background level at an average distance of 17 kb, but the extent of LD varied markedly within the genome and was more than 10-fold higher in ‘genomic islands’ of differentiation than in the rest of the genome. Genetic ancestry analysis identified 33 F₁ hybrids but no later-generation hybrids from sympatric populations of collared flycatchers and pied flycatchers, contradicting earlier reports of backcrosses identified from much fewer number of markers. With an estimated divergence time as recently as <1 Ma, this suggests strong selection against F₁ hybrids and unusually rapid evolution of reproductive incompatibility in an avian system.

Molecular genetic contributions to socioeconomic status and intelligence

Education, socioeconomic status, and intelligence are commonly used as predictors of health outcomes, social environment, and mortality. Education and socioeconomic status are typically viewed as environmental variables although both correlate with intelligence, which has a substantial genetic basis. Using data from 6815 unrelated subjects from the Generation Scotland study, we examined the genetic contributions to these variables and their genetic correlations. Subjects underwent genome-wide testing for common single nucleotide polymorphisms (SNPs). DNA-derived heritability estimates and genetic correlations were calculated using the 'Genome-wide Complex Trait Analyses' (GCTA) procedures. 21% of the variation in education, 18% of the variation in socioeconomic status, and 29% of the variation in general cognitive ability was explained by variation in common SNPs (SEs ~ 5%). The SNP-based genetic correlations of education and socioeconomic status with general intelligence were 0.95 (SE 0.13) and 0.26 (0.16), respectively. There are genetic contributions to intelligence and education with near-complete overlap between common additive SNP effects on these traits (genetic correlation ~ 1). Genetic influences on socioeconomic status are also associated with the genetic foundations of intelligence. The results are also compatible with substantial environmental contributions to socioeconomic status.

Multiple sporadic colorectal cancers display a unique methylation phenotype

Epigenetics are thought to play a major role in the carcinogenesis of multiple sporadic colorectal cancers (CRC). Previous studies have suggested concordant DNA hypermethylation between tumor pairs. However, only a few methylation markers have been analyzed. This study was aimed at describing the epigenetic signature of multiple CRC using a genome-scale DNA methylation profiling. We analyzed 12 patients with synchronous CRC and 29 age-, sex-, and tumor location-paired patients with solitary tumors from the EPICOLON II cohort. DNA methylation profiling was performed using the Illumina Infinium HM27 DNA methylation assay. The most significant results were validated by Methylight. Tumors samples were also analyzed for the CpG Island Methylator Phenotype (CIMP); KRAS and BRAF mutations and mismatch repair deficiency status. Functional annotation clustering was performed. We identified 102 CpG sites that showed significant DNA hypermethylation in multiple tumors with respect to the solitary counterparts (difference in β value ≥0.1). Methylight assays validated the results for 4 selected genes (p = 0.0002). Eight out of 12(66.6%) multiple tumors were classified as CIMP-high, as compared to 5 out of 29(17.2%) solitary tumors (p = 0.004). Interestingly, 76 out of the 102 (74.5%) hypermethylated CpG sites found in multiple tumors were also seen in CIMP-high tumors. Functional analysis of hypermethylated genes found in multiple tumors showed enrichment of genes involved in different tumorigenic functions. In conclusion, multiple CRC are associated with a distinct methylation phenotype, with a close association between tumor multiplicity and CIMP-high. Our results may be important to unravel the underlying mechanism of tumor multiplicity.

A high-density SNP genotyping array for rice biology and molecular breeding

A high-density single nucleotide polymorphism (SNP) array is critically important for geneticists and molecular breeders. With the accumulation of huge amounts of genomic re-sequencing data and available technologies for accurate SNP detection, it is possible to design high-density and high-quality rice SNP arrays. Here we report the development of a high-density rice SNP array and its utility. SNP probes were designed by screening more than 10 000 000 SNP loci extracted from the re-sequencing data of 801 rice varieties and an array named RiceSNP50 was produced on the Illumina Infinium platform. The array contained 51 478 evenly distributed markers, 68% of which were within genic regions. Several hundred rice plants with parent/F1 relationships were used to generate a high-quality cluster file for accurate SNP calling. Application tests showed that this array had high genotyping accuracy, and could be used for different objectives. For example, a core collection of elite rice varieties was clustered with fine resolution. Genome-wide association studies (GWAS) analysis correctly identified a characterized QTL. Further, this array was successfully used for variety verification and trait introgression. As an accurate high-throughput genotyping tool, RiceSNP50 will play an important role in both functional genomics studies and molecular breeding.

Common genetic variants explain the majority of the correlation between height and intelligence: the generation Scotland study

Greater height and higher intelligence test scores are predictors of better health outcomes. Here, we used molecular (single-nucleotide polymorphism) data to estimate the genetic correlation between height and general intelligence (g) in 6,815 unrelated subjects (median age 57, IQR 49–63) from the Generation Scotland: Scottish Family Health Study cohort. The phenotypic correlation between height and g was 0.16 (SE 0.01). The genetic correlation between height and g was 0.28 (SE 0.09) with a bivariate heritability estimate of 0.71. Understanding the molecular basis of the correlation between height and intelligence may help explain any shared role in determining health outcomes. This study identified a modest genetic correlation between height and intelligence with the majority of the phenotypic correlation being explained by shared genetic influences.

Design and characterization of a 52K SNP chip for goats

The success of Genome Wide Association Studies in the discovery of sequence variation linked to complex traits in humans has increased interest in high throughput SNP genotyping assays in livestock species. Primary goals are QTL detection and genomic selection. The purpose here was design of a 50–60,000 SNP chip for goats. The success of a moderate density SNP assay depends on reliable bioinformatic SNP detection procedures, the technological success rate of the SNP design, even spacing of SNPs on the genome and selection of Minor Allele Frequencies (MAF) suitable to use in diverse breeds. Through the federation of three SNP discovery projects consolidated as the International Goat Genome Consortium, we have identified approximately twelve million high quality SNP variants in the goat genome stored in a database together with their biological and technical characteristics. These SNPs were identified within and between six breeds (meat, milk and mixed): Alpine, Boer, Creole, Katjang, Saanen and Savanna, comprising a total of 97 animals. Whole genome and Reduced Representation Library sequences were aligned on >10 kb scaffolds of the de novo goat genome assembly. The 60,000 selected SNPs, evenly spaced on the goat genome, were submitted for oligo manufacturing (Illumina, Inc) and published in dbSNP along with flanking sequences and map position on goat assemblies (i.e. scaffolds and pseudo-chromosomes), sheep genome V2 and cattle UMD3.1 assembly. Ten breeds were then used to validate the SNP content and 52,295 loci could be successfully genotyped and used to generate a final cluster file. The combined strategy of using mainly whole genome Next Generation Sequencing and mapping on a contig genome assembly, complemented with Illumina design tools proved to be efficient in producing this GoatSNP50 chip. Advances in use of molecular markers are expected to accelerate goat genomic studies in coming years.

Tomato breeding in the genomics era: insights from a SNP array

BACKGROUND

The major bottle neck in genetic and linkage studies in tomato has been the lack of a sufficient number of molecular markers. This has radically changed with the application of next generation sequencing and high throughput genotyping. A set of 6000 SNPs was identified and 5528 of them were used to evaluate tomato germplasm at the level of species, varieties and segregating populations.

RESULTS

From the 5528 SNPs, 1980 originated from 454-sequencing, 3495 from Illumina Solexa sequencing and 53 were additional known markers. Genotyping different tomato samples allowed the evaluation of the level of heterozygosity and introgressions among commercial varieties. Cherry tomatoes were especially different from round/beefs in chromosomes 4, 5 and 12. We were able to identify a set of 750 unique markers distinguishing S. lycopersicum 'Moneymaker' from all its distantly related wild relatives. Clustering and neighbour joining analysis among varieties and species showed expected grouping patterns, with S. pimpinellifolium as the most closely related to commercial tomatoes earlier results.

CONCLUSIONS

Our results show that a SNP search in only a few breeding lines already provides generally applicable markers in tomato and its wild relatives. It also shows that the Illumina bead array generated data are highly reproducible. Our SNPs can roughly be divided in two categories: SNPs of which both forms are present in the wild relatives and in domesticated tomatoes (originating from common ancestors) and SNPs unique for the domesticated tomato (originating from after the domestication event). The SNPs can be used for genotyping, identification of varieties, comparison of genetic and physical linkage maps and to confirm (phylogenetic) relations. In the SNPs used for the array there is hardly any overlap with the SolCAP array and it is strongly recommended to combine both SNP sets and to select a core collection of robust SNPs completely covering the entire tomato genome.

Microarray in parasitic infections

Modern biology and genomic sciences are rooted in parasitic disease research. Genome sequencing efforts have provided a wealth of new biological information that promises to have a major impact on our understanding of parasites. Microarrays provide one of the major high-throughput platforms by which this information can be exploited in the laboratory. Many excellent reviews and technique articles have recently been published on applying microarrays to organisms for which fully annotated genomes are at hand. However, many parasitologists work on organisms whose genomes have been only partially sequenced. This review is mainly focused on how to use microarray in these situations.

A SNP resource for Douglas-fir: de novo transcriptome assembly and SNP detection and validation

BACKGROUND

Douglas-fir (Pseudotsuga menziesii), one of the most economically and ecologically important tree species in the world, also has one of the largest tree breeding programs. Although the coastal and interior varieties of Douglas-fir (vars. menziesii and glauca) are native to North America, the coastal variety is also widely planted for timber production in Europe, New Zealand, Australia, and Chile. Our main goal was to develop a SNP resource large enough to facilitate genomic selection in Douglas-fir breeding programs. To accomplish this, we developed a 454-based reference transcriptome for coastal Douglas-fir, annotated and evaluated the quality of the reference, identified putative SNPs, and then validated a sample of those SNPs using the Illumina Infinium genotyping platform.

RESULTS

We assembled a reference transcriptome consisting of 25,002 isogroups (unique gene models) and 102,623 singletons from 2.76 million 454 and Sanger cDNA sequences from coastal Douglas-fir. We identified 278,979 unique SNPs by mapping the 454 and Sanger sequences to the reference, and by mapping four datasets of Illumina cDNA sequences from multiple seed sources, genotypes, and tissues. The Illumina datasets represented coastal Douglas-fir (64.00 and 13.41 million reads), interior Douglas-fir (80.45 million reads), and a Yakima population similar to interior Douglas-fir (8.99 million reads). We assayed 8067 SNPs on 260 trees using an Illumina Infinium SNP genotyping array. Of these SNPs, 5847 (72.5%) were called successfully and were polymorphic.

CONCLUSIONS

Based on our validation efficiency, our SNP database may contain as many as ~200,000 true SNPs, and as many as ~69,000 SNPs that could be genotyped at ~20,000 gene loci using an Infinium II array-more SNPs than are needed to use genomic selection in tree breeding programs. Ultimately, these genomic resources will enhance Douglas-fir breeding and allow us to better understand landscape-scale patterns of genetic variation and potential responses to climate change.

A SNP resource for Douglas-fir: de novo transcriptome assembly and SNP detection and validation

BACKGROUND

Douglas-fir (Pseudotsuga menziesii), one of the most economically and ecologically important tree species in the world, also has one of the largest tree breeding programs. Although the coastal and interior varieties of Douglas-fir (vars. menziesii and glauca) are native to North America, the coastal variety is also widely planted for timber production in Europe, New Zealand, Australia, and Chile. Our main goal was to develop a SNP resource large enough to facilitate genomic selection in Douglas-fir breeding programs. To accomplish this, we developed a 454-based reference transcriptome for coastal Douglas-fir, annotated and evaluated the quality of the reference, identified putative SNPs, and then validated a sample of those SNPs using the Illumina Infinium genotyping platform.

RESULTS

We assembled a reference transcriptome consisting of 25,002 isogroups (unique gene models) and 102,623 singletons from 2.76 million 454 and Sanger cDNA sequences from coastal Douglas-fir. We identified 278,979 unique SNPs by mapping the 454 and Sanger sequences to the reference, and by mapping four datasets of Illumina cDNA sequences from multiple seed sources, genotypes, and tissues. The Illumina datasets represented coastal Douglas-fir (64.00 and 13.41 million reads), interior Douglas-fir (80.45 million reads), and a Yakima population similar to interior Douglas-fir (8.99 million reads). We assayed 8067 SNPs on 260 trees using an Illumina Infinium SNP genotyping array. Of these SNPs, 5847 (72.5%) were called successfully and were polymorphic.

CONCLUSIONS

Based on our validation efficiency, our SNP database may contain as many as ~200,000 true SNPs, and as many as ~69,000 SNPs that could be genotyped at ~20,000 gene loci using an Infinium II array-more SNPs than are needed to use genomic selection in tree breeding programs. Ultimately, these genomic resources will enhance Douglas-fir breeding and allow us to better understand landscape-scale patterns of genetic variation and potential responses to climate change.

Curcumin modulates DNA methylation in colorectal cancer cells

Aim

Recent evidence suggests that several dietary polyphenols may exert their chemopreventive effect through epigenetic modifications. Curcumin is one of the most widely studied dietary chemopreventive agents for colon cancer prevention, however, its effects on epigenetic alterations, particularly DNA methylation, remain unclear. Using systematic genome-wide approaches, we aimed to elucidate the effect of curcumin on DNA methylation alterations in colorectal cancer cells.

Materials and Methods

To evaluate the effect of curcumin on DNA methylation, three CRC cell lines, HCT116, HT29 and RKO, were treated with curcumin. 5-aza-2′-deoxycytidine (5-aza-CdR) and trichostatin A treated cells were used as positive and negative controls for DNA methylation changes, respectively. Methylation status of LINE-1 repeat elements, DNA promoter methylation microarrays and gene expression arrays were used to assess global methylation and gene expression changes. Validation was performed using independent microarrays, quantitative bisulfite pyrosequencing, and qPCR.

Results

As expected, genome-wide methylation microarrays revealed significant DNA hypomethylation in 5-aza-CdR-treated cells (mean β-values of 0.12), however, non-significant changes in mean β-values were observed in curcumin-treated cells. In comparison to mock-treated cells, curcumin-induced DNA methylation alterations occurred in a time-dependent manner. In contrast to the generalized, non-specific global hypomethylation observed with 5-aza-CdR, curcumin treatment resulted in methylation changes at selected, partially-methylated loci, instead of fully-methylated CpG sites. DNA methylation alterations were supported by corresponding changes in gene expression at both up- and down-regulated genes in various CRC cell lines.

Conclusions

Our data provide previously unrecognized evidence for curcumin-mediated DNA methylation alterations as a potential mechanism of colon cancer chemoprevention. In contrast to non-specific global hypomethylation induced by 5-aza-CdR, curcumin-induced methylation changes occurred only in a subset of partially-methylated genes, which provides additional mechanistic insights into the potent chemopreventive effect of this dietary nutraceutical.

Genome wide SNP discovery in flax through next generation sequencing of reduced representation libraries

BACKGROUND

Flax (Linum usitatissimum L.) is a significant fibre and oilseed crop. Current flax molecular markers, including isozymes, RAPDs, AFLPs and SSRs are of limited use in the construction of high density linkage maps and for association mapping applications due to factors such as low reproducibility, intense labour requirements and/or limited numbers. We report here on the use of a reduced representation library strategy combined with next generation Illumina sequencing for rapid and large scale discovery of SNPs in eight flax genotypes. SNP discovery was performed through in silico analysis of the sequencing data against the whole genome shotgun sequence assembly of flax genotype CDC Bethune. Genotyping-by-sequencing of an F6-derived recombinant inbred line population provided validation of the SNPs.

RESULTS

Reduced representation libraries of eight flax genotypes were sequenced on the Illumina sequencing platform resulting in sequence coverage ranging from 4.33 to 15.64X (genome equivalents). Depending on the relatedness of the genotypes and the number and length of the reads, between 78% and 93% of the reads mapped onto the CDC Bethune whole genome shotgun sequence assembly. A total of 55,465 SNPs were discovered with the largest number of SNPs belonging to the genotypes with the highest mapping coverage percentage. Approximately 84% of the SNPs discovered were identified in a single genotype, 13% were shared between any two genotypes and the remaining 3% in three or more. Nearly a quarter of the SNPs were found in genic regions. A total of 4,706 out of 4,863 SNPs discovered in Macbeth were validated using genotyping-by-sequencing of 96 F6 individuals from a recombinant inbred line population derived from a cross between CDC Bethune and Macbeth, corresponding to a validation rate of 96.8%.

CONCLUSIONS

Next generation sequencing of reduced representation libraries was successfully implemented for genome-wide SNP discovery from flax. The genotyping-by-sequencing approach proved to be efficient for validation. The SNP resources generated in this work will assist in generating high density maps of flax and facilitate QTL discovery, marker-assisted selection, phylogenetic analyses, association mapping and anchoring of the whole genome shotgun sequence.

Identifying breast cancer risk loci by global differential allele-specific expression (DASE) analysis in mammary epithelial transcriptome

BACKGROUND

The significant mortality associated with breast cancer (BCa) suggests a need to improve current research strategies to identify new genes that predispose women to breast cancer. Differential allele-specific expression (DASE) has been shown to contribute to phenotypic variables in humans and recently to the pathogenesis of cancer. We previously reported that nonsense-mediated mRNA decay (NMD) could lead to DASE of BRCA1/2, which is associated with elevated susceptibility to breast cancer. In addition to truncation mutations, multiple genetic and epigenetic factors can contribute to DASE, and we propose that DASE is a functional index for cis-acting regulatory variants and pathogenic mutations, and that global analysis of DASE in breast cancer precursor tissues can be used to identify novel causative alleles for breast cancer susceptibility.

RESULTS

To test our hypothesis, we employed the Illumina(®) Omni1-Quad BeadChip in paired genomic DNA (gDNA) and double-stranded cDNA (ds-cDNA) samples prepared from eight BCa patient-derived normal mammary epithelial lines (HMEC). We filtered original array data according to heterozygous genotype calls and calculated DASE values using the Log ratio of cDNA allele intensity, which was normalized to the corresponding gDNA. We developed two statistical methods, SNP- and gene-based approaches, which allowed us to identify a list of 60 candidate DASE loci (DASE ≥ 2.00, P ≤ 0.01, FDR ≤ 0.05) by both methods. Ingenuity Pathway Analysis of DASE loci revealed one major breast cancer-relevant interaction network, which includes two known cancer causative genes, ZNF331 (DASE = 2.31, P = 0.0018, FDR = 0.040) and USP6 (DASE = 4.80, P = 0.0013, FDR = 0.013), and a breast cancer causative gene, DMBT1 (DASE=2.03, P = 0.0017, FDR = 0.014). Sequence analysis of a 5' RACE product of DMBT1 demonstrated that rs2981745, a putative breast cancer risk locus, appears to be one of the causal variants leading to DASE in DMBT1.

CONCLUSIONS

Our study demonstrated for the first time that global DASE analysis is a powerful new approach to identify breast cancer risk allele(s).

Development and evaluation of a 9K SNP array for peach by internationally coordinated SNP detection and validation in breeding germplasm

Although a large number of single nucleotide polymorphism (SNP) markers covering the entire genome are needed to enable molecular breeding efforts such as genome wide association studies, fine mapping, genomic selection and marker-assisted selection in peach [Prunus persica (L.) Batsch] and related Prunus species, only a limited number of genetic markers, including simple sequence repeats (SSRs), have been available to date. To address this need, an international consortium (The International Peach SNP Consortium; IPSC) has pursued a coordinated effort to perform genome-scale SNP discovery in peach using next generation sequencing platforms to develop and characterize a high-throughput Illumina Infinium® SNP genotyping array platform. We performed whole genome re-sequencing of 56 peach breeding accessions using the Illumina and Roche/454 sequencing technologies. Polymorphism detection algorithms identified a total of 1,022,354 SNPs. Validation with the Illumina GoldenGate® assay was performed on a subset of the predicted SNPs, verifying ∼75% of genic (exonic and intronic) SNPs, whereas only about a third of intergenic SNPs were verified. Conservative filtering was applied to arrive at a set of 8,144 SNPs that were included on the IPSC peach SNP array v1, distributed over all eight peach chromosomes with an average spacing of 26.7 kb between SNPs. Use of this platform to screen a total of 709 accessions of peach in two separate evaluation panels identified a total of 6,869 (84.3%) polymorphic SNPs.

The almost 7,000 SNPs verified as polymorphic through extensive empirical evaluation represent an excellent source of markers for future studies in genetic relatedness, genetic mapping, and dissecting the genetic architecture of complex agricultural traits. The IPSC peach SNP array v1 is commercially available and we expect that it will be used worldwide for genetic studies in peach and related stone fruit and nut species.

DNA methylation profiles in precancerous tissue and cancers: carcinogenetic risk estimation and prognostication based on DNA methylation status

Alterations in DNA methylation, which are associated with DNA methyltransferase abnormalities and result in silencing of tumor-related genes and chromosomal instability, are involved even in precancerous changes in various organs. DNA methylation alterations also account for the histological heterogeneity and clinicopathological diversity of human cancers. Therefore, we have analyzed DNA methylation on a genome-wide scale in clinical tissue samples. Our approach using the bacterial artificial chromosome array-based methylated CpG island amplification method has revealed that DNA methylation alterations correlated with the future development of more malignant cancers are already accumulated at the precancerous stage in the kidney, liver and urinary tract. DNA methylation profiles at precancerous stages are basically inherited by the corresponding cancers developing in individual patients. Such DNA methylation alterations may confer vulnerability to further genetic and epigenetic alterations, generate more malignant cancers, and thus determine patient outcome. On the basis of bacterial artificial chromosome array-based methylated CpG island amplification data, indicators for carcinogenetic risk estimation have been established using liver tissue specimens from patients with hepatitis virus infection, chronic hepatitis and liver cirrhosis or histologically normal urothelia, and for prognostication using biopsy or surgically resected specimens from patients with renal cell carcinoma, hepatocellular carcinoma and urothelial carcinoma. Such genome-wide DNA methylation profiling has now firmly established the clinical relevance of translational epigenetics.

The relationship of haplotype in lactotransferrin and its expression levels in Chinese Han ovarian cancer

Chromosomal DNA sequence polymorphisms may contribute to individuality, confer risk for diseases, and most commonly are used as genetic markers in association study. The iron-binding protein lactoferrin inhibits bacterial growth by sequestering essential iron and also exhibits antitumor, anti-inflammatory, and immunoregulatory activities. The gene coding for lactotransferrin (LTF) is polymorphic, with the occurrence of several common alleles in the general population. This genetically determined variation can affect LTF functions. In this study, we determined the distribution of LTF gene polymorphisms (rs1126477, rs1126478, rs2073495, and rs9110) in the Chinese Han population and investigated whether these polymorphisms were associated with increased risk of ovarian carcinoma in the Chinese. It was found that the rs1126477 was correlated significantly with ovarian cancer. The frequency of A allele of rs1126477 was significantly higher in 700 ovarian cancer patients compared with that in the control group of 700 cases (P< 0.01, χ(2)= 6.79). The frequency of AA genotype was significantly higher in ovarian cancer patients compared with that in the control group (P< 0.05, χ(2)= 6.49). AA genotype is the risk factor of ovarian cancer. The odds ratio (OR) was 2.24 and the 95% confidence interval (CI) was 1.08-4.59, respectively. The 'A-G-C-C' haplotype constructed with rs1126477, rs1126478, rs2073495, and rs9110 was the risk factor to be ovarian cancer. The expression of LTF gene was lower in individuals with 'A-G-C-C' haplotype compared with that in individuals without 'A-G-C-C' haplotype. These findings suggested that rs1126477 could play important roles in ovarian carcinoma physiological processes in the Chinese.