A general approach for haplotype phasing across the full spectrum of relatedness

Polygenetic Risk Scores for Major Psychiatric Disorders Among Schizophrenia Patients, Their First-Degree Relatives, and Healthy Participants

BACKGROUND

The genetic etiology of schizophrenia (SCZ) overlaps with that of other major psychiatric disorders in samples of European ancestry. The present study investigated transethnic polygenetic features shared between Japanese SCZ or their unaffected first-degree relatives and European patients with major psychiatric disorders by conducting polygenic risk score (PRS) analyses.

METHODS

To calculate PRSs for 5 psychiatric disorders (SCZ, bipolar disorder [BIP], major depressive disorder, autism spectrum disorder, and attention-deficit/hyperactivity disorder) and PRSs differentiating SCZ from BIP, we utilized large-scale European genome-wide association study (GWAS) datasets as discovery samples. PRSs derived from these GWASs were calculated for 335 Japanese target participants [SCZ patients, FRs, and healthy controls (HCs)]. We took these PRSs based on GWASs of European psychiatric disorders and investigated their effect on risk in Japanese SCZ patients and unaffected first-degree relatives.

RESULTS

The PRSs obtained from European SCZ and BIP patients were higher in Japanese SCZ patients than in HCs. Furthermore, PRSs differentiating SCZ patients from European BIP patients were higher in Japanese SCZ patients than in HCs. Interestingly, PRSs related to European autism spectrum disorder were lower in Japanese first-degree relatives than in HCs or SCZ patients. The PRSs of autism spectrum disorder were positively correlated with a young onset age of SCZ.

CONCLUSIONS

These findings suggest that polygenic factors related to European SCZ and BIP and the polygenic components differentiating SCZ from BIP can transethnically contribute to SCZ risk in Japanese people. Furthermore, we suggest that reduced levels of an ASD-related genetic factor in unaffected first-degree relatives may help protect against SCZ development.

Genotype Imputation in Genome-Wide Association Studies

Genotype imputation infers missing genotypes in silico using haplotype information from reference samples with genotypes from denser genotyping arrays or sequencing. This approach can confer a number of improvements on genome-wide association studies: it can improve statistical power to detect associations by reducing the number of missing genotypes; it can simplify data harmonization for meta-analyses by improving overlap of genomic variants between differently-genotyped sample sets; and it can increase the overall number and density of genomic variants available for association testing. This article reviews the general concepts behind imputation, describes imputation approaches and methods for various types of genotype data, including family-based data, and identifies web-based resources that can be used in different steps of the imputation process. For practical application, it provides a step-by-step guide to implementation of a two-step imputation process consisting of phasing of the study genotypes and the imputation of reference panel genotypes into the study haplotypes. In addition, this review describes recently developed haplotype reference panel resources and online imputation servers that are capable of remotely and securely implementing an imputation workflow on uploaded genotype array data. © 2019 by John Wiley & Sons, Inc.

Putting RFMix and ADMIXTURE to the test in a complex admixed population

Background

Global and local ancestry inference in admixed human populations can be performed using computational tools implementing distinct algorithms. The development and resulting accuracy of these tools has been tested largely on populations with relatively straightforward admixture histories but little is known about how well they perform in more complex admixture scenarios.

Results

Using simulations, we show that RFMix outperforms ADMIXTURE in determining global ancestry proportions even in a complex 5-way admixed population, in addition to assigning local ancestry with an accuracy of 89%. The ability of RFMix to determine global and local ancestry to a high degree of accuracy, particularly in admixed populations provides the opportunity for more accurate association analyses.

Conclusion

This study highlights the utility of the extension of computational tools to become more compatible to genetically structured populations, as well as the need to expand the sampling of diverse world-wide populations. This is particularly noteworthy as modern-day societies are becoming increasingly genetically complex and some genetic tools and commonly used ancestral populations are less appropriate. Based on these caveats and the results presented here, we suggest that RFMix be used for both global and local ancestry estimation in world-wide complex admixture scenarios particularly when including these estimates in association studies.

CHIT1 at Diagnosis Reflects Long-Term Multiple Sclerosis Disease Activity

OBJECTIVE

Evidence for a role of microglia in the pathogenesis of multiple sclerosis (MS) is growing. We investigated association of microglial markers at time of diagnostic lumbar puncture (LP) with different aspects of disease activity (relapses, disability, magnetic resonance imaging parameters) up to 6 years later in a cohort of 143 patients.

METHODS

In cerebrospinal fluid (CSF), we measured 3 macrophage and microglia-related proteins, chitotriosidase (CHIT1), chitinase-3-like protein 1 (CHI3L1 or YKL-40), and soluble triggering receptor expressed on myeloid cells 2 (sTREM2), as well as a marker of neuronal damage, neurofilament light chain (NfL), using enzyme-linked immunosorbent assay and electrochemiluminescence. We investigated the same microglia-related markers in publicly available RNA expression data from postmortem brain tissue.

RESULTS

CHIT1 levels at diagnostic LP correlated with 2 aspects of long-term disease activity after correction for multiple testing. First, CHIT1 increased with reduced tissue integrity in lesions at a median 3 years later (p = 9.6E-04). Second, CHIT1 reflected disease severity at a median 5 years later (p = 1.2E-04). Together with known clinical covariates, CHIT1 levels explained 12% and 27% of variance in these 2 measures, respectively, and were able to distinguish slow and fast disability progression (area under the curve = 85%). CHIT1 was the best discriminator of chronic active versus chronic inactive lesions and the only marker correlated with NfL (r = 0.3, p = 0.0019). Associations with disease activity were, however, independent of NfL.

INTERPRETATION

CHIT1 CSF levels measured during the diagnostic LP reflect microglial activation early on in MS and can be considered a valuable prognostic biomarker for future disease activity. ANN NEUROL 2020;87:633-645.

PRDM15 Is Associated with Risk of Chronic Obstructive Pulmonary Disease in a Rural Population in Chile

BACKGROUND

Genome-wide association studies (GWAS) have accelerated our understanding of the genetic underpinnings of chronic obstructive pulmonary disease (COPD); however, GWAS populations have typically consisted of European descent, with ∼1% of Latin American ancestry.

OBJECTIVE

To overcome this limitation, we conducted a GWAS in a rural Chilean population with increased COPD risk to investigate genetic variation of COPD risk in this understudied minority population.

METHOD

We carried out a case-control study of 214 COPD patients (defined by the GOLD criteria) and 193 healthy controls in Talca, Chile. DNA was extracted from venous blood and genotyped on the Illumina Global Screening Array (n = 754,159 markers). After exclusion based on Hardy-Weinberg equilibrium (p ≤ 0.001), call rates (<95%), and minor allele frequencies (<0.5%) in controls, 455,564 markers were available for logistic regression.

RESULTS

PRDM15 rs1054761 C allele (p = 2.22 × 10-7) was associated with decreased COPD risk. Three PRDM15 SNPs located on chromosome 21 were significantly associated with COPD risk (p < 10-6). Two of these SNPs, rs1054761 and rs4075967, were located on a noncoding transcript variant region of the gene.

CONCLUSION

PRDM15 overexpression may play a role in the B-cell dysregulation in COPD pathogenesis. To the best of our knowledge, the association between PRDM15 and COPD risk was not previously found in GWAS studies in largely European populations, highlighting the importance of investigating novel variants associated with COPD risk among ethnically diverse populations.

Reconstructed Lost Native American Populations from Eastern Brazil Are Shaped by Differential Jê/Tupi Ancestry

After the colonization of the Americas by Europeans and the consequent Trans-Atlantic Slave Trade, most Native American populations in eastern Brazil disappeared or went through an admixture process that configured a population composed of three main genetic components: the European, the sub-Saharan African, and the Native American. The study of the Native American genetic history is challenged by the lack of availability of genome-wide samples from Native American populations, the technical difficulties to develop ancient DNA studies, and the low proportions of the Native American component in the admixed Brazilian populations (on average 7%). We analyzed genome-wide data of 5,825 individuals from three locations of eastern Brazil: Salvador (North-East), Bambui (South-East), and Pelotas (South) and we reconstructed populations that emulate the Native American groups that were living in the 16th century around the sampling locations. This genetic reconstruction was performed after local ancestry analysis of the admixed Brazilian populations, through the rearrangement of the Native American haplotypes into reconstructed individuals with full Native American ancestry (51 reconstructed individuals in Salvador, 45 in Bambui, and 197 in Pelotas). We compared the reconstructed populations with nonadmixed Native American populations from other regions of Brazil through haplotype-based methods. Our results reveal a population structure shaped by the dichotomy of Tupi-/Jê-speaking ancestry related groups. We also show evidence of a decrease of the diversity of nonadmixed Native American groups after the European contact, in contrast with the reconstructed populations, suggesting a reservoir of the Native American genetic diversity within the admixed Brazilian population.

Analysis of Killer Immunoglobulin-Like Receptor Genes in Colorectal Cancer

Natural killer cells (NK cells) play a major role in the immune response to cancer. An important element of NK target recognition is the binding of human leucocyte antigen (HLA) class I molecules by killer immunoglobulin-like receptors (KIRs). Colorectal carcinoma (CRC) is one of the most common types of inflammation-based cancer. The purpose of the present study was to investigate the presence of KIR genes and HLA class I and II alleles in 1074 CRC patients and 1272 controls. We imputed data from single-nucleotide polymorphism (SNP) Illumina OncoArray to identify associations at HLA (HLA–A, B, C, DPB1, DQA1, DQB1, and DRB1) and KIRs (HIBAG and KIR*IMP, respectively). For association analysis, we used PLINK (v1.9), the PyHLA software, and R version 3.4.0. Only three SNP markers showed suggestive associations (p < 10⁻³; rs16896742, rs28367832, and rs9277952). The frequency of KIR2DS3 was significantly increased in the CRC patients compared to healthy controls (p < 0.005). Our results suggest that the implication of NK cells in CRC may not act through allele combinations in KIR and HLA genes. Much larger studies in ethnically homogeneous populations are needed to rule out the possible role of allelic combinations in KIR and HLA genes in CRC risk.

Somatic mosaicism and common genetic variation contribute to the risk of very-early-onset inflammatory bowel disease

Very-early-onset inflammatory bowel disease (VEO-IBD) is a heterogeneous phenotype associated with a spectrum of rare Mendelian disorders. Here, we perform whole-exome-sequencing and genome-wide genotyping in 145 patients (median age-at-diagnosis of 3.5 years), in whom no Mendelian disorders were clinically suspected. In five patients we detect a primary immunodeficiency or enteropathy, with clinical consequences (XIAP, CYBA, SH2D1A, PCSK1). We also present a case study of a VEO-IBD patient with a mosaic de novo, pathogenic allele in CYBB. The mutation is present in ~70% of phagocytes and sufficient to result in defective bacterial handling but not life-threatening infections. Finally, we show that VEO-IBD patients have, on average, higher IBD polygenic risk scores than population controls (99 patients and 18,780 controls; P < 4 × 10-10), and replicate this finding in an independent cohort of VEO-IBD cases and controls (117 patients and 2,603 controls; P < 5 × 10-10). This discovery indicates that a polygenic component operates in VEO-IBD pathogenesis.

A framework for transcriptome-wide association studies in breast cancer in diverse study populations

BACKGROUND

The relationship between germline genetic variation and breast cancer survival is largely unknown, especially in understudied minority populations who often have poorer survival. Genome-wide association studies (GWAS) have interrogated breast cancer survival but often are underpowered due to subtype heterogeneity and clinical covariates and detect loci in non-coding regions that are difficult to interpret. Transcriptome-wide association studies (TWAS) show increased power in detecting functionally relevant loci by leveraging expression quantitative trait loci (eQTLs) from external reference panels in relevant tissues. However, ancestry- or race-specific reference panels may be needed to draw correct inference in ancestrally diverse cohorts. Such panels for breast cancer are lacking.

RESULTS

We provide a framework for TWAS for breast cancer in diverse populations, using data from the Carolina Breast Cancer Study (CBCS), a population-based cohort that oversampled black women. We perform eQTL analysis for 406 breast cancer-related genes to train race-stratified predictive models of tumor expression from germline genotypes. Using these models, we impute expression in independent data from CBCS and TCGA, accounting for sampling variability in assessing performance. These models are not applicable across race, and their predictive performance varies across tumor subtype. Within CBCS (N = 3,828), at a false discovery-adjusted significance of 0.10 and stratifying for race, we identify associations in black women near AURKA, CAPN13, PIK3CA, and SERPINB5 via TWAS that are underpowered in GWAS.

CONCLUSIONS

We show that carefully implemented and thoroughly validated TWAS is an efficient approach for understanding the genetics underpinning breast cancer outcomes in diverse populations.

Beyond large-effect loci: large-scale GWAS reveals a mixed large-effect and polygenic architecture for age at maturity of Atlantic salmon

BACKGROUND

Understanding genetic architecture is essential for determining how traits will change in response to evolutionary processes such as selection, genetic drift and/or gene flow. In Atlantic salmon, age at maturity is an important life history trait that affects factors such as survival, reproductive success, and growth. Furthermore, age at maturity can seriously impact aquaculture production. Therefore, characterizing the genetic architecture that underlies variation in age at maturity is of key interest.

RESULTS

Here, we refine our understanding of the genetic architecture for age at maturity of male Atlantic salmon using a genome-wide association study of 11,166 males from a single aquaculture strain, using imputed genotypes at 512,397 single nucleotide polymorphisms (SNPs). All individuals were genotyped with a 50K SNP array and imputed to higher density using parents genotyped with a 930K SNP array and pedigree information. We found significant association signals on 28 of 29 chromosomes (P-values: 8.7 × 10^-133-9.8 × 10^-8), including two very strong signals spanning the six6 and vgll3 gene regions on chromosomes 9 and 25, respectively. Furthermore, we identified 116 independent signals that tagged 120 candidate genes with varying effect sizes. Five of the candidate genes found here were previously associated with age at maturity in other vertebrates, including humans.

DISCUSSION

These results reveal a mixed architecture of large-effect loci and a polygenic component that consists of multiple smaller-effect loci, suggesting a more complex genetic architecture of Atlantic salmon age at maturity than previously thought. This more complex architecture will have implications for selection on this key trait in aquaculture and for management of wild salmon populations.

scHaplotyper: haplotype construction and visualization for genetic diagnosis using single cell DNA sequencing data

BACKGROUND

Haplotyping reveals chromosome blocks inherited from parents to in vitro fertilized (IVF) embryos in preimplantation genetic diagnosis (PGD), enabling the observation of the transmission of disease alleles between generations. However, the methods of haplotyping that are suitable for single cells are limited because a whole genome amplification (WGA) process is performed before sequencing or genotyping in PGD, and true haplotype profiles of embryos need to be constructed based on genotypes that can contain many WGA artifacts.

RESULTS

Here, we offer scHaplotyper as a genetic diagnosis tool that reconstructs and visualizes the haplotype profiles of single cells based on the Hidden Markov Model (HMM). scHaplotyper can trace the origin of each haplotype block in the embryo, enabling the detection of carrier status of disease alleles in each embryo. We applied this method in PGD in two families affected with genetic disorders, and the result was the healthy live births of two children in the two families, demonstrating the clinical application of this method.

CONCLUSION

Next generation sequencing (NGS) of preimplantation embryos enable genetic screening for families with genetic disorders, avoiding the birth of affected babies. With the validation and successful clinical application, we showed that scHaplotyper is a convenient and accurate method to screen out embryos. More patients with genetic disorder will benefit from the genetic diagnosis of embryos. The source code of scHaplotyper is available at GitHub repository: https://github.com/yzqheart/scHaplotyper.

Genotype phasing in pedigrees using whole-genome sequence data

Phasing is the process of inferring haplotypes from genotype data. Efficient algorithms and associated software for accurate phasing in pedigrees are needed, especially for populations lacking reference panels of sequenced individuals. We present a novel method for phasing genotypes from whole-genome sequence data in pedigrees, called PULSAR (Phasing Using Lineage Specific Alleles/Rare variants). The method is based on the property that alleles specific to a single founding chromosome within a pedigree are highly informative for identifying haplotypes that are shared identical by descent. Simulation studies are used to assess the performance of PULSAR with various pedigree sizes and structures, and the effect of genotyping errors and the presence of nonsequenced individuals is investigated. In pedigrees with complete sequencing and realistic genotyping error rates, PULSAR correctly phases >99.9% of heterozygous genotypes, excluding sites at which all individuals are heterozygous, and does so with a switch error rate frequently below 10^-4. PULSAR is highly accurate, capable of genotype error correction and imputation, and computationally competitive with alternative phasing software applicable to pedigrees. Our method has the significant advantage of not requiring reference panels that are essential for other population-based phasing algorithms. A software implementation of PULSAR is freely available.

Lack of association of genetic variants for diabetic retinopathy in Taiwanese patients with diabetic nephropathy

OBJECTIVE

Diabetic nephropathy (DN) and diabetic retinopathy (DR) comprise major microvascular complications of diabetes that occur with a high concordance rate in patients and are considered to potentially share pathogeneses. In this case-control study, we sought to investigate whether DR-related single nucleotide polymorphisms (SNPs) exert pleiotropic effects on renal function outcomes among patients with diabetes.

RESEARCH DESIGN AND METHODS

A total of 33 DR-related SNPs were identified by replicating published SNPs and via a genome-wide association study. Furthermore, we assessed the cumulative effects by creating a weighted genetic risk score and evaluated the discriminatory and prediction ability of these genetic variants using DN cases according to estimated glomerular filtration rate (eGFR) status along with a cohort with early renal functional decline (ERFD).

RESULTS

Multivariate logistic regression models revealed that the DR-related SNPs afforded no individual or cumulative genetic effect on the nephropathy risk, eGFR status or ERFD outcome among patients with type two diabetes in Taiwan.

CONCLUSION

Our findings indicate that larger studies would be necessary to clearly ascertain the effects of individual genetic variants and further investigation is also required to identify other genetic pathways underlying DN.

Investigating genetic links between grapheme-colour synaesthesia and neuropsychiatric traits

Synaesthesia is a neurological phenomenon affecting perception, where triggering stimuli (e.g. letters and numbers) elicit unusual secondary sensory experiences (e.g. colours). Family-based studies point to a role for genetic factors in the development of this trait. However, the contributions of common genomic variation to synaesthesia have not yet been investigated. Here, we present the SynGenes cohort, the largest genotyped collection of unrelated people with grapheme-colour synaesthesia (n = 723). Synaesthesia has been associated with a range of other neuropsychological traits, including enhanced memory and mental imagery, as well as greater sensory sensitivity. Motivated by the prior literature on putative trait overlaps, we investigated polygenic scores derived from published genome-wide scans of schizophrenia and autism spectrum disorder (ASD), comparing our SynGenes cohort to 2181 non-synaesthetic controls. We found a very slight association between schizophrenia polygenic scores and synaesthesia (Nagelkerke's R2 = 0.0047, empirical p = 0.0027) and no significant association for scores related to ASD (Nagelkerke's R2 = 0.00092, empirical p = 0.54) or body mass index (R2 = 0.00058, empirical p = 0.60), included as a negative control. As sample sizes for studying common genomic variation continue to increase, genetic investigations of the kind reported here may yield novel insights into the shared biology between synaesthesia and other traits, to complement findings from neuropsychology and brain imaging. This article is part of a discussion meeting issue 'Bridging senses: novel insights from synaesthesia'.

De novo diploid genome assembly for genome-wide structural variant detection

Detection of structural variants (SVs) on the basis of read alignment to a reference genome remains a difficult problem. De novo assembly, traditionally used to generate reference genomes, offers an alternative for SV detection. However, it has not been applied broadly to human genomes because of fundamental limitations of short-fragment approaches and high cost of long-read technologies. We here show that 10× linked-read sequencing supports accurate SV detection. We examined variants in six de novo 10× assemblies with diverse experimental parameters from two commonly used human cell lines: NA12878 and NA24385. The assemblies are effective for detecting mid-size SVs, which were discovered by simple pairwise alignment of the assemblies’ contigs to the reference (hg38). Our study also shows that the base-pair level SV breakpoint accuracy is high, with a majority of SVs having precisely correct sizes and breakpoints. Setting the ancestral state of SV loci by comparing to ape orthologs allows inference of the actual molecular mechanism (insertion or deletion) causing the mutation. In about half of cases, the mechanism is the opposite of the reference-based call. We uncover 214 SVs that may have been maintained as polymorphisms in the human lineage since before our divergence from chimp. Overall, we show that de novo assembly of 10× linked-read data can achieve cost-effective SV detection for personal genomes.

Population history and genetic adaptation of the Fulani nomads: inferences from genome-wide data and the lactase persistence trait

Background

Human population history in the Holocene was profoundly impacted by changes in lifestyle following the invention and adoption of food-production practices. These changes triggered significant increases in population sizes and expansions over large distances. Here we investigate the population history of the Fulani, a pastoral population extending throughout the African Sahel/Savannah belt.

Results

Based on genome-wide analyses we propose that ancestors of the Fulani population experienced admixture between a West African group and a group carrying both European and North African ancestries. This admixture was likely coupled with newly adopted herding practices, as it resulted in signatures of genetic adaptation in contemporary Fulani genomes, including the control element of the LCT gene enabling carriers to digest lactose throughout their lives. The lactase persistence (LP) trait in the Fulani is conferred by the presence of the allele T-13910, which is also present at high frequencies in Europe. We establish that the T-13910 LP allele in Fulani individuals analysed in this study lies on a European haplotype background thus excluding parallel convergent evolution. We furthermore directly link the T-13910 haplotype with the Lactase Persistence phenotype through a Genome Wide Association study (GWAS) and identify another genomic region in the vicinity of the SPRY2 gene associated with glycaemic measurements after lactose intake.

Conclusions

Our findings suggest that Eurasian admixture and the European LP allele was introduced into the Fulani through contact with a North African population/s. We furthermore confirm the link between the lactose digestion phenotype in the Fulani to the MCM6/LCT locus by reporting the first GWAS of the lactase persistence trait. We also explored other signals of recent adaptation in the Fulani and identified additional candidates for selection to adapt to herding life-styles.

Comparative genetic architectures of schizophrenia in East Asian and European populations

Schizophrenia is a debilitating psychiatric disorder with approximately 1% lifetime risk globally. Large-scale schizophrenia genetic studies have reported primarily on European ancestry samples, potentially missing important biological insights. Here, we report the largest study to date of East Asian participants (22,778 schizophrenia cases and 35,362 controls), identifying 21 genome-wide-significant associations in 19 genetic loci. Common genetic variants that confer risk for schizophrenia have highly similar effects between East Asian and European ancestries (genetic correlation = 0.98 ± 0.03), indicating that the genetic basis of schizophrenia and its biology are broadly shared across populations. A fixed-effect meta-analysis including individuals from East Asian and European ancestries identified 208 significant associations in 176 genetic loci (53 novel). Trans-ancestry fine-mapping reduced the sets of candidate causal variants in 44 loci. Polygenic risk scores had reduced performance when transferred across ancestries, highlighting the importance of including sufficient samples of major ancestral groups to ensure their generalizability across populations.

Stable longitudinal associations of family income with children's hippocampal volume and memory persist after controlling for polygenic scores of educational attainment

Despite common notion that the correlation of socioeconomic status with child cognitive performance may be driven by both environmentally- and genetically-mediated transactional pathways, there is a lack of longitudinal and genetically informed research that examines these postulated associations. The present study addresses whether family income predicts associative memory growth and hippocampal development in middle childhood and tests whether these associations persist when controlling for DNA-based polygenic scores of educational attainment. Participants were 142 6-to-7-year-old children, of which 127 returned when they were 8-to-9 years old. Longitudinal analyses indicated that the association of family income with children's memory performance and hippocampal volume remained stable over this age range and did not predict change. On average, children from economically disadvantaged background showed lower memory performance and had a smaller hippocampal volume. There was no evidence to suggest that differences in memory performance were mediated by differences in hippocampal volume. Further exploratory results suggested that the relationship of income with hippocampal volume and memory in middle childhood is not primarily driven by genetic variance captured by polygenic scores of educational attainment, despite the fact that polygenic scores significantly predicted family income.

Local ancestry at APOE modifies Alzheimer's disease risk in Caribbean Hispanics

INTRODUCTION

Although the relationship between APOE and Alzheimer's disease (AD) is well established in populations of European descent, the effects of APOE and ancestry on AD risk in diverse populations is not well understood.

METHODS

Logistic mixed model regression and survival analyses were performed in a sample of 3067 Caribbean Hispanics and 3028 individuals of European descent to assess the effects of APOE genotype, local ancestry, and genome-wide ancestry on AD risk and age at onset.

RESULTS

Among the Caribbean Hispanics, individuals with African-derived ancestry at APOE had 39% lower odds of AD than individuals with European-derived APOE, after adjusting for APOE genotype, age, and genome-wide ancestry. While APOE E2 and E4 effects on AD risk and age at onset were significant in the Caribbean Hispanics, they were substantially attenuated compared with those in European ancestry individuals.

DISCUSSION

These results suggest that additional genetic variation in the APOE region influences AD risk beyond APOE E2/E3/E4.

Screening Human Embryos for Polygenic Traits Has Limited Utility

The increasing proportion of variance in human complex traits explained by polygenic scores, along with progress in preimplantation genetic diagnosis, suggests the possibility of screening embryos for traits such as height or cognitive ability. However, the expected outcomes of embryo screening are unclear, which undermines discussion of associated ethical concerns. Here, we use theory, simulations, and real data to evaluate the potential gain of embryo screening, defined as the difference in trait value between the top-scoring embryo and the average embryo. The gain increases very slowly with the number of embryos but more rapidly with the variance explained by the score. Given current technology, the average gain due to screening would be ≈2.5 cm for height and ≈2.5 IQ points for cognitive ability. These mean values are accompanied by wide prediction intervals, and indeed, in large nuclear families, the majority of children top-scoring for height are not the tallest.

Novel genetic loci affecting facial shape variation in humans

The human face represents a combined set of highly heritable phenotypes, but knowledge on its genetic architecture remains limited, despite the relevance for various fields. A series of genome-wide association studies on 78 facial shape phenotypes quantified from 3-dimensional facial images of 10,115 Europeans identified 24 genetic loci reaching study-wide suggestive association (p < 5 × 10-8), among which 17 were previously unreported. A follow-up multi-ethnic study in additional 7917 individuals confirmed 10 loci including six unreported ones (padjusted < 2.1 × 10-3). A global map of derived polygenic face scores assembled facial features in major continental groups consistent with anthropological knowledge. Analyses of epigenomic datasets from cranial neural crest cells revealed abundant cis-regulatory activities at the face-associated genetic loci. Luciferase reporter assays in neural crest progenitor cells highlighted enhancer activities of several face-associated DNA variants. These results substantially advance our understanding of the genetic basis underlying human facial variation and provide candidates for future in-vivo functional studies.

Kinpute: using identity by descent to improve genotype imputation

MOTIVATION

Genotype imputation, though generally accurate, often results in many genotypes being poorly imputed, particularly in studies where the individuals are not well represented by standard reference panels. When individuals in the study share regions of the genome identical by descent (IBD), it is possible to use this information in combination with a study-specific reference panel (SSRP) to improve the imputation results. Kinpute uses IBD information-due to recent, familial relatedness or distant, unknown ancestors-in conjunction with the output from linkage disequilibrium (LD) based imputation methods to compute more accurate genotype probabilities. Kinpute uses a novel method for IBD imputation, which works even in the absence of a pedigree, and results in substantially improved imputation quality.

RESULTS

Given initial estimates of average IBD between subjects in the study sample, Kinpute uses a novel algorithm to select an optimal set of individuals to sequence and use as an SSRP. Kinpute is designed to use as input both this SSRP and the genotype probabilities output from other LD-based imputation software, and uses a new method to combine the LD imputed genotype probabilities with IBD configurations to substantially improve imputation. We tested Kinpute on a human population isolate where 98 individuals have been sequenced. In half of this sample, whose sequence data was masked, we used Impute2 to perform LD-based imputation and Kinpute was used to obtain higher accuracy genotype probabilities. Measures of imputation accuracy improved significantly, particularly for those genotypes that Impute2 imputed with low certainty.

AVAILABILITY AND IMPLEMENTATION

Kinpute is an open-source and freely available C++ software package that can be downloaded from.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Assessment of human diploid genome assembly with 10x Linked-Reads data

BACKGROUND

Producing cost-effective haplotype-resolved personal genomes remains challenging. 10x Linked-Read sequencing, with its high base quality and long-range information, has been demonstrated to facilitate de novo assembly of human genomes and variant detection. In this study, we investigate in depth how the parameter space of 10x library preparation and sequencing affects assembly quality, on the basis of both simulated and real libraries.

RESULTS

We prepared and sequenced eight 10x libraries with a diverse set of parameters from standard cell lines NA12878 and NA24385 and performed whole-genome assembly on the data. We also developed the simulator LRTK-SIM to follow the workflow of 10x data generation and produce realistic simulated Linked-Read data sets. We found that assembly quality could be improved by increasing the total sequencing coverage (C) and keeping physical coverage of DNA fragments (CF) or read coverage per fragment (CR) within broad ranges. The optimal physical coverage was between 332× and 823× and assembly quality worsened if it increased to >1,000× for a given C. Long DNA fragments could significantly extend phase blocks but decreased contig contiguity. The optimal length-weighted fragment length (W${\mu _{FL}}$) was ∼50-150 kb. When broadly optimal parameters were used for library preparation and sequencing, ∼80% of the genome was assembled in a diploid state.

CONCLUSIONS

The Linked-Read libraries we generated and the parameter space we identified provide theoretical considerations and practical guidelines for personal genome assemblies based on 10x Linked-Read sequencing.

Uganda Genome Resource Enables Insights into Population History and Genomic Discovery in Africa

Genomic studies in African populations provide unique opportunities to understand disease etiology, human diversity, and population history. In the largest study of its kind, comprising genome-wide data from 6,400 individuals and whole-genome sequences from 1,978 individuals from rural Uganda, we find evidence of geographically correlated fine-scale population substructure. Historically, the ancestry of modern Ugandans was best represented by a mixture of ancient East African pastoralists. We demonstrate the value of the largest sequence panel from Africa to date as an imputation resource. Examining 34 cardiometabolic traits, we show systematic differences in trait heritability between European and African populations, probably reflecting the differential impact of genes and environment. In a multi-trait pan-African GWAS of up to 14,126 individuals, we identify novel loci associated with anthropometric, hematological, lipid, and glycemic traits. We find that several functionally important signals are driven by Africa-specific variants, highlighting the value of studying diverse populations across the region.

Transethnic meta-analysis of metabolic syndrome in a multiethnic study

Genome-wide association studies (GWAS) have been used to establish thousands of genetic associations across numerous phenotypes. To improve the power of GWAS and generalize associations across ethnic groups, transethnic meta-analysis methods are used to combine the results of several GWAS from diverse ancestries. The goal of this study is to identify genetic associations for eight quantitative metabolic syndrome (MetS) traits through a meta-analysis across four ethnic groups. Traits were measured in the GENetics of Noninsulin dependent Diabetes Mellitus (GENNID) Study which consists of African-American (families = 73, individuals = 288), European-American (families = 79, individuals = 519), Japanese-American (families = 17, individuals = 132), and Mexican-American (families = 113, individuals = 610) samples. Genome-wide association results from these four ethnic groups were combined using four meta-analysis methods: fixed effects, random effects, TransMeta, and MR-MEGA. We provide an empirical comparison of the four meta-analysis methods from the GENNID results, discuss which types of loci (characterized by allelic heterogeneity) appear to be better detected by each of the four meta-analysis methods in the GENNID Study, and validate our results using previous genetic discoveries. We specifically compare the two transethnic methods, TransMeta and MR-MEGA, and discuss how each transethnic method's framework relates to the types of loci best detected by each method.

Interval breast cancer is associated with other types of tumors

Breast cancer (BC) patients diagnosed between two screenings (interval cancers) are more likely than screen-detected patients to carry rare deleterious mutations in cancer genes potentially leading to increased risk for other non-breast cancer (non-BC) tumors. In this study, we include 14,846 women diagnosed with BC of which 1,772 are interval and 13,074 screen-detected. Compared to women with screen-detected cancers, interval breast cancer patients are more likely to have a non-BC tumor before (Odds ratio (OR): 1.43 [1.19–1.70], P = 9.4 x 10⁻⁵) and after (OR: 1.28 [1.14–1.44], P = 4.70 x 10⁻⁵) breast cancer diagnosis, are more likely to report a family history of non-BC tumors and have a lower genetic risk score based on common variants for non-BC tumors. In conclusion, interval breast cancer is associated with other tumors and common cancer variants are unlikely to be responsible for this association. These findings could have implications for future screening and prevention programs.

Interval cancer patients are more likely to carry rare gene mutations than screen-detected breast cancer patients. Here, the authors report that interval cancer patients are more likely cancer survivors and are at a greater risk of developing other non-breast tumors.

Inter-individual genomic heterogeneity within European population isolates

A number of studies carried out since the early ‘70s has investigated the effects of isolation on genetic variation within and among human populations in diverse geographical contexts. However, no extensive analysis has been carried out on the heterogeneity among genomes within isolated populations. This issue is worth exploring since events of recent admixture and/or subdivision could potentially disrupt the genetic homogeneity which is to be expected when isolation is prolonged and constant over time. Here, we analyze literature data relative to 87,815 autosomal single-nucleotide polymorphisms, which were obtained from a total of 28 European populations. Our results challenge the traditional paradigm of population isolates as structured as genetically (and genomically) uniform entities. In fact, focusing on the distribution of variance of intra-population diversity measures across individuals, we show that the inter-individual heterogeneity of isolated populations is at least comparable to the open ones. More in particular, three small and highly inbred isolates (Sappada, Sauris and Timau in Northeastern Italy) were found to be characterized by levels of inter-individual heterogeneity largely exceeding that of all other populations, possibly due to relatively recent events of genetic introgression. Finally, we propose a way to monitor the effects of inter-individual heterogeneity in disease-gene association studies.

International meta-analysis of PTSD genome-wide association studies identifies sex- and ancestry-specific genetic risk loci

The risk of posttraumatic stress disorder (PTSD) following trauma is heritable, but robust common variants have yet to be identified. In a multi-ethnic cohort including over 30,000 PTSD cases and 170,000 controls we conduct a genome-wide association study of PTSD. We demonstrate SNP-based heritability estimates of 5-20%, varying by sex. Three genome-wide significant loci are identified, 2 in European and 1 in African-ancestry analyses. Analyses stratified by sex implicate 3 additional loci in men. Along with other novel genes and non-coding RNAs, a Parkinson's disease gene involved in dopamine regulation, PARK2, is associated with PTSD. Finally, we demonstrate that polygenic risk for PTSD is significantly predictive of re-experiencing symptoms in the Million Veteran Program dataset, although specific loci did not replicate. These results demonstrate the role of genetic variation in the biology of risk for PTSD and highlight the necessity of conducting sex-stratified analyses and expanding GWAS beyond European ancestry populations.

A genome-wide association and replication study of blood pressure in Ugandan early adolescents

BACKGROUND

Genetic association studies of blood pressure (BP) have mostly been conducted in non-African populations. Using the Entebbe Mother and Baby Study (EMaBS), we aimed to identify genetic variants associated with BP among Ugandan adolescents.

METHODS

Systolic and diastolic BP were measured among 10- and 11-year olds. Whole-genome genotype data were generated using Illumina omni 2.5M arrays and untyped variants were imputed. Genome-wide association study (GWAS) was conducted using linear mixed model regression to account for population structure. Linear regression analysis was used to assess whether variants previously associated with BP (p < 5.0 × 10-8 ) in published BP GWASs were replicated in our study.

RESULTS

Of the 14 million variants analyzed among 815 adolescents, none reached genome-wide significance (p < 5.0×10-8 ) for association with systolic or diastolic BP. The most strongly associated variants were rs181430167 (p = 6.8 × 10-7 ) for systolic BP and rs12991132 (p = 4.0 × 10-7 ) for diastolic BP. Thirty-three (17 single nucleotide polymorphisms (SNPs) for systolic BP, 15 SNPs for diastolic BP and one SNP for both) of 330 variants previously identified as associated with BP were replicated in this study, but none remained significant after accounting for multiple testing.

CONCLUSION

Variants showing suggestive associations are worthy of future investigation. Replication results suggest that variants influencing adolescent BP may overlap somewhat with those already established in previous studies, largely based on adults in Western settings.

Noninvasive prenatal diagnosis of β-thalassemia by relative haplotype dosage without analyzing proband

BACKGROUND

β-thalassemia is one of the most common monogenic diseases in the world. Southeast China is a highly infected area affected by four β-thalassemia mutation types (HBB:c.-78A>G, HBB:c.52A>T, HBB:c.126_129delCTTT, and HBB:c.316-197C>T). Relative haplotype dosage (RHDO), a haplotype-based approach, has shown promise as an application for noninvasive prenatal diagnosis (NIPD); however, additional family members (such as the proband) are required for haplotype construction. The abovementioned circumstances make RHDO-based NIPD cost prohibitive; additionally, the genetic information of the proband is not always available. Thus, it is necessary to find a practical method to solve these problems.

METHODS

Targeted sequencing was applied to sequence parental genomic DNA and cell-free fetal DNA (cffDNA). Parental haplotypes were constructed with the SHAPEIT software based on the 1000 Genomes Project (1000G) Phase 3 v5 Southern Han Chinese (CHS) haplotype dataset. Single-nucleotide polymorphisms (SNPs) in the target region were called and classified, and the fetal mutation inheritance status was deduced using the RHDO method.

RESULTS

Construction of the parental haplotypes and detection of the inherited parental mutations were successfully achieved in five families, despite a suspected recombination event. The status of the affected fetuses is consistent with the results of traditional reverse dot blot (RDB) diagnosis.

CONCLUSION

This research introduced SHAPEIT into the classical RHDO workflow and proved that it is applicable to construct parental haplotypes without information from other family members.

Distribution of local ancestry and evidence of adaptation in admixed populations

Admixed American populations have different global proportions of European, Sub-Saharan African, and Native-American ancestry. However, individuals who display the same global ancestry could exhibit remarkable differences in the distribution of local ancestry blocks. We studied for the first time the distribution of local ancestry across the genome of 264 Brazilian admixed individuals, ascertained within the scope of the Brazilian Initiative on Precision Medicine. We found a decreased proportion of European ancestry together with an excess of Native-American ancestry on chromosome 8p23.1 and showed that this is due to haplotypes created by chromosomal inversion events. Furthermore, Brazilian non-inverted haplotypes were more similar to Native-American haplotypes than to European haplotypes, in contrast to what was found in other American admixed populations. We also identified signals of recent positive selection on chromosome 8p23.1, and one gene within this locus, PPP1R3B, is related to glycogenesis and has been associated with an increased risk of type 2 diabetes and obesity. These findings point to a selection event after admixture, which is still not entirely understood in recent admixture events.

European Roma groups show complex West Eurasian admixture footprints and a common South Asian genetic origin

The Roma population is the largest transnational ethnic minority in Europe, characterized by a linguistic, cultural and historical heterogeneity. Comparative linguistics and genetic studies have placed the origin of European Roma in the Northwest of India. After their migration across Persia, they entered into the Balkan Peninsula, from where they spread into Europe, arriving in the Iberian Peninsula in the 15th century. Their particular demographic history has genetic implications linked to rare and common diseases. However, the South Asian source of the proto-Roma remains still untargeted and the West Eurasian Roma component has not been yet deeply characterized. Here, in order to describe both the South Asian and West Eurasian ancestries, we analyze previously published genome-wide data of 152 European Roma and 34 new Iberian Roma samples at a fine-scale and haplotype-based level, with special focus on the Iberian Roma genetic substructure. Our results suggest that the putative origin of the proto-Roma involves a Punjabi group with low levels of West Eurasian ancestry. In addition, we have identified a complex West Eurasian component (around 65%) in the Roma, as a result of the admixture events occurred with non-proto-Roma populations between 1270–1580. Particularly, we have detected the Balkan genetic footprint in all European Roma, and the Baltic and Iberian components in the Northern and Western Roma groups, respectively. Finally, our results show genetic substructure within the Iberian Roma, with different levels of West Eurasian admixture, as a result of the complex historical events occurred in the Peninsula.

Human demographic processes and admixture events leave traceable footprints in the genomes of the populations and they can modulate the genetic architecture of complex diseases. Here, we aim to study the Roma people, an admixed population with a particular demographic history recognized as the largest ethnic minority in Europe. Previous studies suggest that they originated in South Asia 1,500 years ago and followed a diaspora towards Europe with extensive admixture with non-Roma West Eurasian groups. However, the genetic components of the Roma have not been deeply characterized. Our study reveals a common South Asian origin of all European Roma, closely related to a Punjabi group from Northwestern India. Through fine-scale haplotype-based methods, we describe a complex West Eurasian genetic component in the Roma groups, identifying a common Balkan ancestry and country-specific admixture footprints consistent with the dispersion through Europe. Our findings provide new insights into the demographic history and recent admixture events that have shaped the genetic composition of European Roma groups and could enable a better genetic characterization of complex disease in this population.

A comprehensive assessment of demographic, environmental, and host genetic associations with gut microbiome diversity in healthy individuals

BACKGROUND

The gut microbiome is an important determinant of human health. Its composition has been shown to be influenced by multiple environmental factors and likely by host genetic variation. In the framework of the Milieu Intérieur Consortium, a total of 1000 healthy individuals of western European ancestry, with a 1:1 sex ratio and evenly stratified across five decades of life (age 20-69), were recruited. We generated 16S ribosomal RNA profiles from stool samples for 858 participants. We investigated genetic and non-genetic factors that contribute to individual differences in fecal microbiome composition.

RESULTS

Among 110 demographic, clinical, and environmental factors, 11 were identified as significantly correlated with α-diversity, ß-diversity, or abundance of specific microbial communities in multivariable models. Age and blood alanine aminotransferase levels showed the strongest associations with microbiome diversity. In total, all non-genetic factors explained 16.4% of the variance. We then searched for associations between > 5 million single nucleotide polymorphisms and the same indicators of fecal microbiome diversity, including the significant non-genetic factors as covariates. No genome-wide significant associations were identified after correction for multiple testing. A small fraction of previously reported associations between human genetic variants and specific taxa could be replicated in our cohort, while no replication was observed for any of the diversity metrics.

CONCLUSION

In a well-characterized cohort of healthy individuals, we identified several non-genetic variables associated with fecal microbiome diversity. In contrast, host genetics only had a negligible influence. Demographic and environmental factors are thus the main contributors to fecal microbiome composition in healthy individuals.

TRIAL REGISTRATION

ClinicalTrials.gov identifier NCT01699893.

Early progression to active tuberculosis is a highly heritable trait driven by 3q23 in Peruvians

Of the 1.8 billion people worldwide infected with Mycobacterium tuberculosis, 5-15% will develop active tuberculosis (TB). Approximately half will progress to active TB within the first 18 months after infection, presumably because they fail to mount an effective initial immune response. Here, in a genome-wide genetic study of early TB progression, we genotype 4002 active TB cases and their household contacts in Peru. We quantify genetic heritability ([Formula: see text]) of early TB progression to be 21.2% (standard error 0.08). This suggests TB progression has a strong genetic basis, and is comparable to traits with well-established genetic bases. We identify a novel association between early TB progression and variants located in a putative enhancer region on chromosome 3q23 (rs73226617, OR = 1.18; P = 3.93 × 10-8). With in silico and in vitro analyses we identify rs73226617 or rs148722713 as the likely functional variant and ATP1B3 as a potential causal target gene with monocyte specific function.

Genotyping by low-coverage whole-genome sequencing in intercross pedigrees from outbred founders: a cost-efficient approach

BACKGROUND

Experimental intercrosses between outbred founder populations are powerful resources for mapping loci that contribute to complex traits i.e. quantitative trait loci (QTL). Here, we present an approach and its accompanying software for high-resolution reconstruction of founder mosaic genotypes in the intercross offspring from such populations using whole-genome high-coverage sequence data on founder individuals (~ 30×) and very low-coverage sequence data on intercross individuals (< 0.5×). Sets of founder-line informative markers were selected for each full-sib family and used to infer the founder mosaic genotypes of the intercross individuals. The application of this approach and the quality of the estimated genome-wide genotypes are illustrated in a large F₂ pedigree between two divergently selected lines of chickens.

RESULTS

We describe how we obtained whole-genome genotype data for hundreds of individuals in a cost- and time-efficient manner by using a Tn5-based library preparation protocol and an imputation algorithm that was optimized for this application. In total, 7.6 million markers segregated in this pedigree and, within each full-sib family, between 10.0 and 13.7% of these were fully informative, i.e. fixed for alternative alleles in the founders from the divergent lines, and were used for reconstruction of the offspring mosaic genotypes. The genotypes that were estimated based on the low-coverage sequence data were highly consistent (> 95% agreement) with those obtained using individual single nucleotide polymorphism (SNP) genotyping. The estimated resolution of the inferred recombination breakpoints was relatively high, with 50% of them being defined on regions shorter than 10 kb.

CONCLUSIONS

A method and software for inferring founder mosaic genotypes in intercross offspring from low-coverage whole-genome sequencing in pedigrees from heterozygous founders are described. They provide high-quality, high-resolution genotypes in a time- and cost-efficient manner. The software is freely available at https://github.com/CarlborgGenomics/Stripes .

Assessment of Novel Genome-Wide Significant Gene Loci and Lesion Growth in Geographic Atrophy Secondary to Age-Related Macular Degeneration

Importance

Age-related macular degeneration (AMD) is a common threat to vision loss in individuals older than 50 years. While neovascular complications in AMD are treatable, there is currently no therapy for geographic atrophy secondary to AMD. Geographic atrophy lesion progression over time shows considerable interindividual variability, but little is known about prognostic factors.

Objective

To elucidate the contribution of common genetic variants to geographic atrophy lesion growth.

Design, Setting, and Participants

This pooled analysis combined 4 independent studies: the Fundus Autofluorescence Imaging in Age-Related Macular Degeneration (FAM) study, the Directional Spread in Geographic Atrophy (DSGA) study, the Age-Related Eye Disease Study (AREDS), and the Geographic Atrophy Treatment Evaluation (GATE) study. Each provided data for geographic atrophy lesion growth in specific designs. Patients with geographic atrophy secondary to AMD were recruited to these studies. Genotypes were retrieved through the database of Genotypes and Phenotypes (for AREDS) or generated at the Cologne Center for Genomics (for FAM, DSGA, and GATE).

Main Outcomes

The correlation between square root-transformed geographic atrophy growth rate and 7 596 219 genetic variants passing quality control was estimated using linear regression. The calculations were adjusted for known factors influencing geographic atrophy growth, such as the presence of bilateral geographic atrophy as well as the number of lesion spots and follow-up times.

Main Outcomes and Measures

Slopes per allele, 95% CIs, and P values of genetic variants correlated with geographic atrophy lesion growth.

Results

A total of 935 patients (mean [SD] age, 74.7 [7.8] years; 547 female participants [59.0%]) were included. Two gene loci with conservative genome-wide significance were identified. Each minor allele of the genome-wide associated variants increased the geographic atrophy growth rate by a mean of about 15% or 0.05 mm per year. Gene prioritization within each locus suggests the protein arginine methyltransferase 6 gene (PRMT6; chromosome 1; slope, 0.046 [95% CI, 0.026-0.066]; P = 4.09 × 10-8) and the lanosterol synthase gene (LSS; chromosome 21; slope, 0.105 [95% CI, 0.068-0.143]; P = 4.07 × 10-7) as the most likely progression-associated genes.

Conclusions and Relevance

These data provide further insight into the genetic architecture of geographic atrophy lesion growth. Geographic atrophy is a clinical outcome with a high medical need for effective therapy. The genes PRMT6 and LSS are promising candidates for future studies aimed at understanding functional aspects of geographic atrophy progression and also for designing novel and targeted treatment options.

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.

Genome Analysis Reveals Genetic Admixture and Signature of Selection for Productivity and Environmental Traits in Iraqi Cattle

The Near East cattle are adapted to different agro-ecological zones including desert areas, mountains habitats, and humid regions along the Tigris and Euphrates rivers system. The region was one of the earliest and most significant areas of cattle husbandry. Currently, four main breeds of Iraqi cattle are recognized. Among these, the Jenoubi is found in the southern more humid part of Iraq, while the Rustaqi is found in the middle and drier region of the country. Despite their importance, Iraqi cattle have up to now been poorly characterized at the genome level. Here, we report at a genome-wide level the diversity and signature of positive selection in these two breeds. Thirty-five unrelated Jenoubi cattle, sampled in the Maysan and Basra regions, and 60 Rustaqi cattle, from around Baghdad and Babylon, were genotyped using the Illumina Bovine HD BeadChip (700K). Genetic population structure and diversity level were studied using principal component analysis (PCA), expected heterozygosity (He), observed heterozygosity (Ho), and admixture. Signatures of selection were studied using extended haplotype homozygosity (EHH) (iHS and Rsb) and inter-population Wright's Fst. The results of PCA and admixture analysis, including European taurine, Asian indicine, African indicine, and taurine indicate that the two breeds are crossbreed zebu × taurine, with more zebu background in Jenoubi cattle compared with Rustaqi. The Rustaqi has the greatest mean heterozygosity (He = 0.37) among all breeds. iHS and Rsb signatures of selection analyses identify 68 candidate genes under positive selection in the two Iraqi breeds, while Fst analysis identifies 220 candidate genes including genes related to the innate and acquired immunity responses, different environmental selection pressures (e.g., tick resistance and heat stress), and genes of commercial interest (e.g., marbling score).

Integrative Genomics Reveals the Genetics and Evolution of the Honey Bee's Social Immune System

Social organisms combat pathogens through individual innate immune responses or through social immunity—behaviors among individuals that limit pathogen transmission within groups. Although we have a relatively detailed understanding of the genetics and evolution of the innate immune system of animals, we know little about social immunity. Addressing this knowledge gap is crucial for understanding how life-history traits influence immunity, and identifying if trade-offs exist between innate and social immunity. Hygienic behavior in the Western honey bee, Apis mellifera, provides an excellent model for investigating the genetics and evolution of social immunity in animals. This heritable, colony-level behavior is performed by nurse bees when they detect and remove infected or dead brood from the colony. We sequenced 125 haploid genomes from two artificially selected highly hygienic populations and a baseline unselected population. Genomic contrasts allowed us to identify a minimum of 73 genes tentatively associated with hygienic behavior. Many genes were within previously discovered QTLs associated with hygienic behavior and were predictive of hygienic behavior within the unselected population. These genes were often involved in neuronal development and sensory perception in solitary insects. We found that genes associated with hygienic behavior have evidence of positive selection within honey bees (Apis), supporting the hypothesis that social immunity contributes to fitness. Our results indicate that genes influencing developmental neurobiology and behavior in solitary insects may have been co-opted to give rise to a novel and adaptive social immune phenotype in honey bees.

Association of BUD13-ZNF259-APOA5-APOA1-SIK3 cluster polymorphism in 11q23.3 and structure of APOA5 with increased plasma triglyceride levels in a Korean population

In this association study on chromosome 11, the data from 12,537 Korean individuals within the Health Examinee (HEXA) and the Korea Association Resource (KARE) projects were analysed to identify genetic loci correlating with increased and decreased plasma triglyceride (TG) levels. We identified a locus in chromosomal region 11q23.3 that harbours genes BUD13, ZNF259, APOA5, APOA1, and SIK3, which may be associated with plasma TG levels. In this locus, 13 relevant single-nucleotide polymorphisms (SNPs) were found: rs184616707, rs118175510, rs60954647, rs79408961, and rs180373 (near BUD13); rs11604424 (in ZNF259); rs2075291, rs651821, and rs7123666 (in or near APOA5); rs525028 (near APOA1), and rs645258, rs10160754, and rs142395187 (in or near SIK3). All 13 SNPs satisfied the genome-wide significance level (P < 5.0 × 10^-8) in both meta-analysis and conditional analysis. Haplotype analysis of six SNPs (rs79408961, rs180373, rs2075291, rs651821, rs525028, and rs10160754) that were selected based on the β coefficient and conditional P values, revealed nine common haplotypes (with frequency 0.02-0.34) associated with both increased and reduced TG levels. Furthermore, to shed light on possible structural implications, we modelled and simulated the G185C variant of APOA5 (corresponding to rs2075291), which showed the strongest association. Molecular dynamics simulation results showed that this polymorphic variant of APOA5 has a different hydrogen bond network, increased average distance between chains, and an ability to form distinct clusters. Owing to the orientation of cysteine, the possibility of disulphide bond formation with other proteins is evident. In summary, our association and modelling analyses provided evidence that genetic variations in chromosomal region 11q23.3 are associated with elevated TG levels.

New insight into human sweet taste: a genome-wide association study of the perception and intake of sweet substances

BACKGROUND

Individual differences in human perception of sweetness are partly due to genetics; however, which genes are associated with the perception and the consumption of sweet substances remains unclear.

OBJECTIVE

The aim of this study was to verify previous reported associations within genes involved in the peripheral receptor systems (i.e., TAS1R2, TAS1R3, and GNAT3) and reveal novel loci.

METHODS

We performed genome-wide association scans (GWASs) of the perceived intensity of 2 sugars (glucose and fructose) and 2 high-potency sweeteners (neohesperidin dihydrochalcone and aspartame) in an Australian adolescent twin sample (n = 1757), and the perceived intensity and sweetness and the liking of sucrose in a US adult twin sample (n = 686). We further performed GWASs of the intake of total sugars (i.e., total grams of all dietary mono- and disaccharides per day) and sweets (i.e., handfuls of candies per day) in the UK Biobank sample (n = ≤174,424 white-British individuals). All participants from the 3 independent samples were of European ancestry.

RESULTS

We found a strong association between the intake of total sugars and the single nucleotide polymorphism rs11642841 within the FTO gene on chromosome 16 (P = 3.8 × 10-8) and many suggestive associations (P < 1.0 × 10-5) for each of the sweet perception and intake phenotypes. We showed genetic evidence for the involvement of the brain in both sweet taste perception and sugar intake. There was limited support for the associations with TAS1R2, TAS1R3, and GNAT3 in all 3 European samples.

CONCLUSIONS

Our findings indicate that genes additional to those involved in the peripheral receptor system are also associated with the sweet taste perception and intake of sweet-tasting foods. The functional potency of the genetic variants within TAS1R2, TAS1R3, and GNAT3 may be different between ethnic groups and this warrants further investigations.

Palaeo-Eskimo genetic ancestry and the peopling of Chukotka and North America

Much of the American Arctic was first settled 5,000 years ago, by groups of people known as Palaeo-Eskimos. They were subsequently joined and largely displaced around 1,000 years ago by ancestors of the present-day Inuit and Yup'ik1-3. The genetic relationship between Palaeo-Eskimos and Native American, Inuit, Yup'ik and Aleut populations remains uncertain4-6. Here we present genomic data for 48 ancient individuals from Chukotka, East Siberia, the Aleutian Islands, Alaska, and the Canadian Arctic. We co-analyse these data with data from present-day Alaskan Iñupiat and West Siberian populations and published genomes. Using methods based on rare-allele and haplotype sharing, as well as established techniques4,7-9, we show that Palaeo-Eskimo-related ancestry is ubiquitous among people who speak Na-Dene and Eskimo-Aleut languages. We develop a comprehensive model for the Holocene peopling events of Chukotka and North America, and show that Na-Dene-speaking peoples, people of the Aleutian Islands, and Yup'ik and Inuit across the Arctic region all share ancestry from a single Palaeo-Eskimo-related Siberian source.

Signatures of selection in Charolais beef cattle identified by genome-wide analysis

Charolais cattle are one of the most important breeds for meat production worldwide; in México, its selection is mainly made by live weight traits. One strategy for mapping important genomic regions that might influence productive traits is the identification of signatures of selection. This type of genomic features contains loci with extended linkage disequilibrium (LD) and homozygosity patterns that are commonly associated with sites of quantitative trait locus (QTL). Therefore, the objective of this study was to identify the signatures of selection in Charolais cattle genotyped with the GeneSeek Genomic Profiler Bovine HD panel consisting of 77 K single nucleotide polymorphisms (SNPs). A total 61,311 SNPs and 819 samples were used for the analysis. Identification of signatures of selection was carried out using the integrated haplotype score (iHS) methodology implemented in the rehh R package. The top ten SNPs with the highest piHS values were located on BTA 4, 5, 6 and 14. By identifying markers in LD with top ten SNPs, the candidate regions defined were mapped to 52.8-59.3 Mb on BTA 4; 67.5-69.3 on BTA 5; 39.5-41.0 Mb on BTA 6; and 26.4-29.6 Mb on BTA 14. The comparison of these candidate regions with the bovine QTLdb effectively confirmed the association (p < 0.05) with QTL related to growth traits and other important productive traits. The genomic regions identified in this study indicated selection for growth traits on the Charolais population via the conservation of haplotypes on various chromosomes. These genomic regions and their associated genes could serve as the basis for haplotype association studies and for the identification of causal genes related to growth traits.

Cardiac Troponin T and Troponin I in the General Population

BACKGROUND

There is great interest in widening the use of high-sensitivity cardiac troponins for population cardiovascular disease (CVD) and heart failure screening. However, it is not clear whether cardiac troponin T (cTnT) and troponin I (cTnI) are equivalent measures of risk in this setting. We aimed to compare and contrast (1) the association of cTnT and cTnI with CVD and non-CVD outcomes, and (2) their determinants in a genome-wide association study.

METHODS

High-sensitivity cTnT and cTnI were measured in serum from 19 501 individuals in Generation Scotland Scottish Family Health Study. Median follow-up was 7.8 years (quartile 1 to quartile 3, 7.1-9.2). Associations of each troponin with a composite CVD outcome (1177 events), CVD death (n=266), non-CVD death (n=374), and heart failure (n=216) were determined by using Cox models. A genome-wide association study was conducted using a standard approach developed for the cohort.

RESULTS

Both cTnI and cTnT were strongly associated with CVD risk in unadjusted models. After adjusting for classical risk factors, the hazard ratio for a 1 SD increase in log transformed troponin was 1.24 (95% CI, 1.17-1.32) and 1.11 (1.04-1.19) for cTnI and cTnT, respectively; ratio of hazard ratios 1.12 (1.04-1.21). cTnI, but not cTnT, was associated with myocardial infarction and coronary heart disease. Both cTnI and cTnT had strong associations with CVD death and heart failure. By contrast, cTnT, but not cTnI, was associated with non-CVD death; ratio of hazard ratios 0.77 (0.67-0.88). We identified 5 loci (53 individual single-nucleotide polymorphisms) that had genome-wide significant associations with cTnI, and a different set of 4 loci (4 single-nucleotide polymorphisms) for cTnT.

CONCLUSIONS

The upstream genetic causes of low-grade elevations in cTnI and cTnT appear distinct, and their associations with outcomes also differ. Elevations in cTnI are more strongly associated with some CVD outcomes, whereas cTnT is more strongly associated with the risk of non-CVD death. These findings help inform the selection of an optimal troponin assay for future clinical care and research in this setting.

Genetic components of human pain sensitivity: a protocol for a genome-wide association study of experimental pain in healthy volunteers

INTRODUCTION

Pain constitutes a major component of the global burden of diseases. Recent studies suggest a strong genetic contribution to pain susceptibility and severity. Whereas most of the available evidence relies on candidate gene association or linkage studies, research on the genetic basis of pain sensitivity using genome-wide association studies (GWAS) is still in its infancy. This protocol describes a proposed GWAS on genetic contributions to baseline pain sensitivity and nociceptive sensitisation in a sample of unrelated healthy individuals of mixed Latin American ancestry.

METHODS AND ANALYSIS

A GWAS on genetic contributions to pain sensitivity in the naïve state and following nociceptive sensitisation will be conducted in unrelated healthy individuals of mixed ancestry. Mechanical and thermal pain sensitivity will be evaluated with a battery of quantitative sensory tests evaluating pain thresholds. In addition, variation in mechanical and thermal sensitisation following topical application of mustard oil to the skin will be evaluated.

ETHICS AND DISSEMINATION

This study received ethical approval from the University College London research ethics committee (3352/001) and from the bioethics committee of the Odontology Faculty at the University of Antioquia (CONCEPTO 01-2013). Findings will be disseminated to commissioners, clinicians and service users via papers and presentations at international conferences.

Genome-wide association study of anti-Müllerian hormone levels in pre-menopausal women of late reproductive age and relationship with genetic determinants of reproductive lifespan

Anti-Müllerian hormone (AMH) is required for sexual differentiation in the fetus, and in adult females AMH is produced by growing ovarian follicles. Consequently, AMH levels are correlated with ovarian reserve, declining towards menopause when the oocyte pool is exhausted. A previous genome-wide association study identified three genetic variants in and around the AMH gene that explained 25% of variation in AMH levels in adolescent males but did not identify any genetic associations reaching genome-wide significance in adolescent females. To explore the role of genetic variation in determining AMH levels in women of late reproductive age, we carried out a genome-wide meta-analysis in 3344 pre-menopausal women from five cohorts (median age 44-48 years at blood draw). A single genetic variant, rs16991615, previously associated with age at menopause, reached genome-wide significance at P = 3.48 × 10-10, with a per allele difference in age-adjusted inverse normal AMH of 0.26 standard deviations (SD) (95% confidence interval (CI) [0.18,0.34]). We investigated whether genetic determinants of female reproductive lifespan were more generally associated with pre-menopausal AMH levels. Genetically-predicted age at menarche had no robust association but genetically-predicted age at menopause was associated with lower AMH levels by 0.18 SD (95% CI [0.14,0.21]) in age-adjusted inverse normal AMH per one-year earlier age at menopause. Our findings provide genetic support for the well-established use of AMH as a marker of ovarian reserve.

Genome-wide sequence analyses of ethnic populations across Russia

The Russian Federation is the largest and one of the most ethnically diverse countries in the world, however no centralized reference database of genetic variation exists to date. Such data are crucial for medical genetics and essential for studying population history. The Genome Russia Project aims at filling this gap by performing whole genome sequencing and analysis of peoples of the Russian Federation. Here we report the characterization of genome-wide variation of 264 healthy adults, including 60 newly sequenced samples. People of Russia carry known and novel genetic variants of adaptive, clinical and functional consequence that in many cases show allele frequency divergence from neighboring populations. Population genetics analyses revealed six phylogeographic partitions among indigenous ethnicities corresponding to their geographic locales. This study presents a characterization of population-specific genomic variation in Russia with results important for medical genetics and for understanding the dynamic population history of the world's largest country.