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

The Pop-Gen Pipeline Platform: A Software Platform for Population Genomic Analyses

The Pop-Gen Pipeline Platform (PPP) is a software platform for population genomic analyses. The PPP was designed as a collection of scripts that facilitate common population genomic workflows in a consistent and standardized Python environment. Functions were developed to encompass entire workflows, including input preparation, file format conversion, various population genomic analyses, and output generation. The platform has also been developed with reproducibility and extensibility of analyses in mind. The PPP is an open-source package that is available for download and use at https://ppp.readthedocs.io/en/latest/PPP_pages/install.html.

Genetic architecture underlying nascent speciation - The evolution of Eurasian pigs under domestication

Speciation is a process whereby the evolution of reproductive barriers leads to isolated species. Although many studies have addressed large-effect genetic footprints in the advanced stages of speciation, the genetics of reproductive isolation in nascent stage of speciation remains unclear. Here, we show that pig domestication offers an interesting model for studying the early stages of speciation in great details. Pig breeds have not evolved the large X-effect of hybrid incompatibility commonly observed between “good species.” Instead, deleterious epistatic interactions among multiple autosomal loci are common. These weak Dobzhansky–Muller incompatibilities confer partial hybrid inviability with sex biases in crosses between European and East Asian domestic pigs. The genomic incompatibility is enriched in pathways for angiogenesis, androgen receptor signaling and immunity, with an observation of many highly differentiated cis-regulatory variants. Our study suggests that partial hybrid inviability caused by pervasive but weak interactions among autosomal loci may be a hallmark of nascent speciation in mammals.

Common Variants Associated With OSMR Expression Contribute to Carotid Plaque Vulnerability, but Not to Cardiovascular Disease in Humans

Background and Aims: Oncostatin M (OSM) signaling is implicated in atherosclerosis, however the mechanism remains unclear. We investigated the impact of common genetic variants in OSM and its receptors, OSMR and LIFR, on overall plaque vulnerability, plaque phenotype, intraplaque OSMR and LIFR expression, coronary artery calcification burden and cardiovascular disease susceptibility. Methods and Results: We queried Genotype-Tissue Expression data and found that rs13168867 (C allele) was associated with decreased OSMR expression and that rs10491509 (A allele) was associated with increased LIFR expression in arterial tissues. No variant was significantly associated with OSM expression. We associated these two variants with plaque characteristics from 1,443 genotyped carotid endarterectomy patients in the Athero-Express Biobank Study. After correction for multiple testing, rs13168867 was significantly associated with an increased overall plaque vulnerability (β = 0.118 ± s.e. = 0.040, p = 3.00 × 10-3, C allele). Looking at individual plaque characteristics, rs13168867 showed strongest associations with intraplaque fat (β = 0.248 ± s.e. = 0.088, p = 4.66 × 10-3, C allele) and collagen content (β = -0.259 ± s.e. = 0.095, p = 6.22 × 10-3, C allele), but these associations were not significant after correction for multiple testing. rs13168867 was not associated with intraplaque OSMR expression. Neither was intraplaque OSMR expression associated with plaque vulnerability and no known OSMR eQTLs were associated with coronary artery calcification burden, or cardiovascular disease susceptibility. No associations were found for rs10491509 in the LIFR locus. Conclusions: Our study suggests that rs1316887 in the OSMR locus is associated with increased plaque vulnerability, but not with coronary calcification or cardiovascular disease risk. It remains unclear through which precise biological mechanisms OSM signaling exerts its effects on plaque morphology. However, the OSM-OSMR/LIFR pathway is unlikely to be causally involved in lifetime cardiovascular disease susceptibility.

Toward a fine-scale population health monitoring system

Understanding population health disparities is an essential component of equitable precision health efforts. Epidemiology research often relies on definitions of race and ethnicity, but these population labels may not adequately capture disease burdens and environmental factors impacting specific sub-populations. Here, we propose a framework for repurposing data from electronic health records (EHRs) in concert with genomic data to explore the demographic ties that can impact disease burdens. Using data from a diverse biobank in New York City, we identified 17 communities sharing recent genetic ancestry. We observed 1,177 health outcomes that were statistically associated with a specific group and demonstrated significant differences in the segregation of genetic variants contributing to Mendelian diseases. We also demonstrated that fine-scale population structure can impact the prediction of complex disease risk within groups. This work reinforces the utility of linking genomic data to EHRs and provides a framework toward fine-scale monitoring of population health.

A Likelihood Approach for Uncovering Selective Sweep Signatures from Haplotype Data

Selective sweeps are frequent and varied signatures in the genomes of natural populations, and detecting them is consequently important in understanding mechanisms of adaptation by natural selection. Following a selective sweep, haplotypic diversity surrounding the site under selection decreases, and this deviation from the background pattern of variation can be applied to identify sweeps. Multiple methods exist to locate selective sweeps in the genome from haplotype data, but none leverages the power of a model-based approach to make their inference. Here, we propose a likelihood ratio test statistic T to probe whole-genome polymorphism data sets for selective sweep signatures. Our framework uses a simple but powerful model of haplotype frequency spectrum distortion to find sweeps and additionally make an inference on the number of presently sweeping haplotypes in a population. We found that the T statistic is suitable for detecting both hard and soft sweeps across a variety of demographic models, selection strengths, and ages of the beneficial allele. Accordingly, we applied the T statistic to variant calls from European and sub-Saharan African human populations, yielding primarily literature-supported candidates, including LCT, RSPH3, and ZNF211 in CEU, SYT1, RGS18, and NNT in YRI, and HLA genes in both populations. We also searched for sweep signatures in Drosophila melanogaster, finding expected candidates at Ace, Uhg1, and Pimet. Finally, we provide open-source software to compute the T statistic and the inferred number of presently sweeping haplotypes from whole-genome data.

Haplotype-resolved diverse human genomes and integrated analysis of structural variation

Long-read and strand-specific sequencing technologies together facilitate the de novo assembly of high-quality haplotype-resolved human genomes without parent-child trio data. We present 64 assembled haplotypes from 32 diverse human genomes. These highly contiguous haplotype assemblies (average minimum contig length needed to cover 50% of the genome: 26 million base pairs) integrate all forms of genetic variation, even across complex loci. We identified 107,590 structural variants (SVs), of which 68% were not discovered with short-read sequencing, and 278 SV hotspots (spanning megabases of gene-rich sequence). We characterized 130 of the most active mobile element source elements and found that 63% of all SVs arise through homology-mediated mechanisms. This resource enables reliable graph-based genotyping from short reads of up to 50,340 SVs, resulting in the identification of 1526 expression quantitative trait loci as well as SV candidates for adaptive selection within the human population.

Obesity and ADHD: Exploring the role of body composition, BMI polygenic risk score, and reward system genes

The association between obesity and attention-deficit hyperactivity disorder (ADHD) has been extensively reported in the literature. However, the potential mechanisms underlying this association are not completely understood. This study aimed to evaluate the association between body composition and ADHD and explore the possible genetic mechanisms involved. We used data from the 1982 Pelotas (Brazil) Birth Cohort at age 30-year follow-up (N = 3630). We first used logistic regression analysis to test whether body mass index (BMI), fat mass (FM), and fat-free mass (FFM) were associated with ADHD. We further tested the association between BMI polygenic risk score (BMI-PRS) and ADHD and the role of the genes upregulated in the reward system using a gene-set association approach. BMI (odds ratio [OR] = 1.05; 95% confidence interval [CI], 1.00-1.09; p = 0.038) and FM (OR = 1.04; 95% CI, 1.00-1.07; p = 0.043) were associated with ADHD. The BMI-PRS was associated with ADHD (using p-value threshold (P_T) = 0.4; OR = 1.65; 95% CI, 1.02-2.65) at a nominal level. In gene-set analysis, the reward system genes were associated with BMI in subjects with a high BMI-PRS score, considering P_T = 0.4 (p = 0.014). The results suggest that BMI genetic components, especially those genes related to the reward system, may be involved in this association.

Exploring a Region on Chromosome 8p23.1 Displaying Positive Selection Signals in Brazilian Admixed Populations: Additional Insights Into Predisposition to Obesity and Related Disorders

We recently reported a deviation of local ancestry on the chromosome (ch) 8p23.1, which led to positive selection signals in a Brazilian population sample. The deviation suggested that the genetic variability of candidate genes located on ch 8p23.1 may have been evolutionarily advantageous in the early stages of the admixture process. In the present work, we aim to extend the previous work by studying additional Brazilian admixed individuals and examining DNA sequencing data from the ch 8p23.1 candidate region. Thus, we inferred the local ancestry of 125 exomes from individuals born in five towns within the Southeast region of Brazil (São Paulo, Campinas, Barretos, and Ribeirão Preto located in the state of São Paulo and Belo Horizonte, the capital of the state of Minas Gerais), and compared to data from two public Brazilian reference genomic databases, BIPMed and ABraOM, and with information from the 1000 Genomes Project phase 3 and gnomAD databases. Our results revealed that ancestry is similar among individuals born in the five Brazilian towns assessed; however, an increased proportion of sub-Saharan African ancestry was observed in individuals from Belo Horizonte. In addition, individuals from the five towns considered, as well as those from the ABRAOM dataset, had the same overrepresentation of Native-American ancestry on the ch 8p23.1 locus that was previously reported for the BIPMed reference sample. Sequencing analysis of ch 8p23.1 revealed the presence of 442 non-synonymous variants, including frameshift, inframe deletion, start loss, stop gain, stop loss, and splicing site variants, which occurred in 24 genes. Among these genes, 13 were associated with obesity, type II diabetes, lipid levels, and waist circumference (PRAG1, MFHAS1, PPP1R3B, TNKS, MSRA, PRSS55, RP1L1, PINX1, MTMR9, FAM167A, BLK, GATA4, and CTSB). These results strengthen the hypothesis that a set of variants located on ch 8p23.1 that result from positive selection during early admixture events may influence obesity-related disease predisposition in admixed individuals of the Brazilian population. Furthermore, we present evidence that the exploration of local ancestry deviation in admixed individuals may provide information with the potential to be translated into health care improvement.

Population History and Gene Divergence in Native Mexicans Inferred from 76 Human Exomes

Native American genetic variation remains underrepresented in most catalogs of human genome sequencing data. Previous genotyping efforts have revealed that Mexico’s Indigenous population is highly differentiated and substructured, thus potentially harboring higher proportions of private genetic variants of functional and biomedical relevance. Here we have targeted the coding fraction of the genome and characterized its full site frequency spectrum by sequencing 76 exomes from five Indigenous populations across Mexico. Using diffusion approximations, we modeled the demographic history of Indigenous populations from Mexico with northern and southern ethnic groups splitting 7.2 KYA and subsequently diverging locally 6.5 and 5.7 KYA, respectively. Selection scans for positive selection revealed BCL2L13 and KBTBD8 genes as potential candidates for adaptive evolution in Rarámuris and Triquis, respectively. BCL2L13 is highly expressed in skeletal muscle and could be related to physical endurance, a well-known phenotype of the northern Mexico Rarámuri. The KBTBD8 gene has been associated with idiopathic short stature and we found it to be highly differentiated in Triqui, a southern Indigenous group from Oaxaca whose height is extremely low compared to other Native populations.

The Counteracting Effects of Demography on Functional Genomic Variation: The Roma Paradigm

Demographic history plays a major role in shaping the distribution of genomic variation. Yet the interaction between different demographic forces and their effects in the genomes is not fully resolved in human populations. Here, we focus on the Roma population, the largest transnational ethnic minority in Europe. They have a South Asian origin and their demographic history is characterized by recent dispersals, multiple founder events, and extensive gene flow from non-Roma groups. Through the analyses of new high-coverage whole exome sequences and genome-wide array data for 89 Iberian Roma individuals together with forward simulations, we show that founder effects have reduced their genetic diversity and proportion of rare variants, gene flow has counteracted the increase in mutational load, runs of homozygosity show ancestry-specific patterns of accumulation of deleterious homozygotes, and selection signals primarily derive from preadmixture adaptation in the Roma population sources. The present study shows how two demographic forces, bottlenecks and admixture, act in opposite directions and have long-term balancing effects on the Roma genomes. Understanding how demography and gene flow shape the genome of an admixed population provides an opportunity to elucidate how genomic variation is modeled in human populations.

MOSTWAS: Multi-Omic Strategies for Transcriptome-Wide Association Studies

Traditional predictive models for transcriptome-wide association studies (TWAS) consider only single nucleotide polymorphisms (SNPs) local to genes of interest and perform parameter shrinkage with a regularization process. These approaches ignore the effect of distal-SNPs or other molecular effects underlying the SNP-gene association. Here, we outline multi-omics strategies for transcriptome imputation from germline genetics to allow more powerful testing of gene-trait associations by prioritizing distal-SNPs to the gene of interest. In one extension, we identify mediating biomarkers (CpG sites, microRNAs, and transcription factors) highly associated with gene expression and train predictive models for these mediators using their local SNPs. Imputed values for mediators are then incorporated into the final predictive model of gene expression, along with local SNPs. In the second extension, we assess distal-eQTLs (SNPs associated with genes not in a local window around it) for their mediation effect through mediating biomarkers local to these distal-eSNPs. Distal-eSNPs with large indirect mediation effects are then included in the transcriptomic prediction model with the local SNPs around the gene of interest. Using simulations and real data from ROS/MAP brain tissue and TCGA breast tumors, we show considerable gains of percent variance explained (1-2% additive increase) of gene expression and TWAS power to detect gene-trait associations. This integrative approach to transcriptome-wide imputation and association studies aids in identifying the complex interactions underlying genetic regulation within a tissue and important risk genes for various traits and disorders.

A case-only study to identify genetic modifiers of breast cancer risk for BRCA1/BRCA2 mutation carriers

Breast cancer (BC) risk for BRCA1 and BRCA2 mutation carriers varies by genetic and familial factors. About 50 common variants have been shown to modify BC risk for mutation carriers. All but three, were identified in general population studies. Other mutation carrier-specific susceptibility variants may exist but studies of mutation carriers have so far been underpowered. We conduct a novel case-only genome-wide association study comparing genotype frequencies between 60,212 general population BC cases and 13,007 cases with BRCA1 or BRCA2 mutations. We identify robust novel associations for 2 variants with BC for BRCA1 and 3 for BRCA2 mutation carriers, P < 10-8, at 5 loci, which are not associated with risk in the general population. They include rs60882887 at 11p11.2 where MADD, SP11 and EIF1, genes previously implicated in BC biology, are predicted as potential targets. These findings will contribute towards customising BC polygenic risk scores for BRCA1 and BRCA2 mutation carriers.

Genetic Effects on Transcriptome Profiles in Colon Epithelium Provide Functional Insights for Genetic Risk Loci

BACKGROUND & AIMS

The association of genetic variation with tissue-specific gene expression and alternative splicing guides functional characterization of complex trait-associated loci and may suggest novel genes implicated in disease. Here, our aims were as follows: (1) to generate reference profiles of colon mucosa gene expression and alternative splicing and compare them across colon subsites (ascending, transverse, and descending), (2) to identify expression and splicing quantitative trait loci (QTLs), (3) to find traits for which identified QTLs contribute to single-nucleotide polymorphism (SNP)-based heritability, (4) to propose candidate effector genes, and (5) to provide a web-based visualization resource.

METHODS

We collected colonic mucosal biopsy specimens from 485 healthy adults and performed bulk RNA sequencing. We performed genome-wide SNP genotyping from blood leukocytes. Statistical approaches and bioinformatics software were used for QTL identification and downstream analyses.

RESULTS

We provided a complete quantification of gene expression and alternative splicing across colon subsites and described their differences. We identified thousands of expression and splicing QTLs and defined their enrichment at genome-wide regulatory regions. We found that part of the SNP-based heritability of diseases affecting colon tissue, such as colorectal cancer and inflammatory bowel disease, but also of diseases affecting other tissues, such as psychiatric conditions, can be explained by the identified QTLs. We provided candidate effector genes for multiple phenotypes. Finally, we provided the Colon Transcriptome Explorer web application.

CONCLUSIONS

We provide a large characterization of gene expression and splicing across colon subsites. Our findings provide greater etiologic insight into complex traits and diseases influenced by transcriptomic changes in colon tissue.

A GWAS in Latin Americans identifies novel face shape loci, implicating VPS13B and a Denisovan introgressed region in facial variation

To characterize the genetic basis of facial features in Latin Americans, we performed a genome-wide association study (GWAS) of more than 6000 individuals using 59 landmark-based measurements from two-dimensional profile photographs and ~9,000,000 genotyped or imputed single-nucleotide polymorphisms. We detected significant association of 32 traits with at least 1 (and up to 6) of 32 different genomic regions, more than doubling the number of robustly associated face morphology loci reported until now (from 11 to 23). These GWAS hits are strongly enriched in regulatory sequences active specifically during craniofacial development. The associated region in 1p12 includes a tract of archaic adaptive introgression, with a Denisovan haplotype common in Native Americans affecting particularly lip thickness. Among the nine previously unidentified face morphology loci we identified is the VPS13B gene region, and we show that variants in this region also affect midfacial morphology in mice.

Inferring recombination patterns in African populations

Although several high-resolution recombination maps exist for European-descent populations, the recombination landscape of African populations remains relatively understudied. Given that there is high genetic divergence among groups in Africa, it is possible that recombination hotspots also diverge significantly. Both limitations and opportunities exist for developing recombination maps for these populations. In this review, we discuss various recombination inference methods, and the strengths and weaknesses of these methods in analyzing recombination in African-descent populations. Furthermore, we provide a decision tree and recommendations for which inference method to use in various research contexts. Establishing an appropriate methodology for recombination rate inference in a particular study will improve the accuracy of various downstream analyses including but not limited to local ancestry inference, haplotype phasing, fine-mapping of GWAS loci and genome assemblies.

Exome-wide scan identifies significant association of rs4788084 in IL27 promoter with increase in hepatic fat content among Indians

BACKGROUND

Non-alcoholic fatty liver disease (NAFLD) is a global epidemic that often progresses to liver cirrhosis and hepatocellular carcinoma. In contrast to most world populations where NAFLD is mostly prevalent among obese, NAFLD among Indians and generally among South and South-East Asians is unique and highly prevalent among individuals who are lean. Genetics of NAFLD in Indian populations is understudied. In this study, we have used an exome-wide approach to identify genetic determinants of hepatic fat content (HFC) in India.

METHODS

HFC was measured in 244 participants using Proton magnetic resonance spectroscopy (H1-MRS). Quantitative trait loci (QTL) mapping was done exome-wide, to identify SNPs associated with HFC. The effects of the interaction between adiposity and QTLs on HFC were studied using a regression model. Association of the significant loci with disease severity was studied in 146 NAFLD patients among 244 participants, who underwent liver biopsy.

RESULTS

Our study identified 4 significantly associated SNPs (rs738409 and rs2281135 (PNPLA3), rs3761472 (SAMM50), rs17513722 (FAM161A) and rs4788084), with HFC after adjusting for the effects of covariates (p-value < 0.0005). rs738409, rs2281135 (PNPLA3), and rs3761472 (SAMM50) were associated with hepatocyte ballooning, lobular and portal inflammation and non-alcoholic steatohepatitis (NASH) (p-value < 0.05). rs4788048 is an eQTL for IL27 and SULT1A2 genes, both of which are highly expressed in healthy livers and are likely to be involved in NAFLD pathogenesis.

CONCLUSIONS

Our study identified the novel association of rs4788084 with HFC, which regulates the expression of IL-27, an immune regulatory gene. We further showed that adiposity affected the HFC, irrespective of the genetic predisposition.

High-resolution population-specific recombination rates and their effect on phasing and genotype imputation

Previous research has shown that using population-specific reference panels has a significant effect on downstream population genomic analyses like haplotype phasing, genotype imputation, and association, especially in the context of population isolates. Here, we developed a high-resolution recombination rate mapping at 10 and 50 kb scale using high-coverage (20–30×) whole-genome sequenced data of 55 family trios from Finland and compared it to recombination rates of non-Finnish Europeans (NFE). We tested the downstream effects of the population-specific recombination rates in statistical phasing and genotype imputation in Finns as compared to the same analyses performed by using the NFE-based recombination rates. We found that Finnish recombination rates have a moderately high correlation (Spearman’s ρ = 0.67–0.79) with NFE, although on average (across all autosomal chromosomes), Finnish rates (2.268 ± 0.4209 cM/Mb) are 12–14% lower than NFE (2.641 ± 0.5032 cM/Mb). Finnish recombination map was found to have no significant effect in haplotype phasing accuracy (switch error rates ~2%) and average imputation concordance rates (97–98% for common, 92–96% for low frequency and 78–90% for rare variants). Our results suggest that haplotype phasing and genotype imputation mostly depend on population-specific contexts like appropriate reference panels and their sample size, but not on population-specific recombination maps. Even though recombination rate estimates had some differences between the Finnish and NFE populations, haplotyping and imputation had not been noticeably affected by the recombination map used. Therefore, the currently available HapMap recombination maps seem robust for population-specific phasing and imputation pipelines, even in the context of relatively isolated populations like Finland.

Alternative Applications of Genotyping Array Data Using Multivariant Methods

One of the forerunners that pioneered the revolution of high-throughput genomic technologies is the genotyping microarray technology, which can genotype millions of single-nucleotide variants simultaneously. Owing to apparent benefits, such as high speed, low cost, and high throughput, the genotyping array has gained lasting applications in genome-wide association studies (GWAS) and thus accumulated an enormous amount of data. Empowered by continuous manufactural upgrades and analytical innovation, unconventional applications of genotyping array data have emerged to address more diverse genetic problems, holding promise of boosting genetic research into human diseases through the re-mining of the rich accumulated data. Here, we review several unconventional genotyping array analysis techniques that have been built on the idea of large-scale multivariant analysis and provide empirical application examples. These unconventional outcomes of genotyping arrays include polygenic score, runs of homozygosity (ROH)/heterozygosity ratio, distant pedigree computation, and mitochondrial DNA (mtDNA) copy number inference.

Genome-wide association study in a Korean population identifies six novel susceptibility loci for rheumatoid arthritis

OBJECTIVE

Genome-wide association studies (GWAS) in rheumatoid arthritis (RA) have discovered over 100 RA loci, explaining patient-relevant RA pathogenesis but showing a large fraction of missing heritability. As a continuous effort, we conducted GWAS in a large Korean RA case-control population.

METHODS

We newly generated genome-wide variant data in two independent Korean cohorts comprising 4068 RA cases and 36 487 controls, followed by a whole-genome imputation and a meta-analysis of the disease association results in the two cohorts. By integrating publicly available omics data with the GWAS results, a series of bioinformatic analyses were conducted to prioritise the RA-risk genes in RA loci and to dissect biological mechanisms underlying disease associations.

RESULTS

We identified six new RA-risk loci (SLAMF6, CXCL13, SWAP70, NFKBIA, ZFP36L1 and LINC00158) with pmeta<5×10-8 and consistent disease effect sizes in the two cohorts. A total of 122 genes were prioritised from the 6 novel and 13 replicated RA loci based on physical distance, regulatory variants and chromatin interaction. Bioinformatics analyses highlighted potentially RA-relevant tissues (including immune tissues, lung and small intestine) with tissue-specific expression of RA-associated genes and suggested the immune-related gene sets (such as CD40 pathway, IL-21-mediated pathway and citrullination) and the risk-allele sharing with other diseases.

CONCLUSION

This study identified six new RA-associated loci that contributed to better understanding of the genetic aetiology and biology in RA.

High-depth African genomes inform human migration and health

The African continent is regarded as the cradle of modern humans and African genomes contain more genetic variation than those from any other continent, yet only a fraction of the genetic diversity among African individuals has been surveyed¹. Here we performed whole-genome sequencing analyses of 426 individuals-comprising 50 ethnolinguistic groups, including previously unsampled populations-to explore the breadth of genomic diversity across Africa. We uncovered more than 3 million previously undescribed variants, most of which were found among individuals from newly sampled ethnolinguistic groups, as well as 62 previously unreported loci that are under strong selection, which were predominantly found in genes that are involved in viral immunity, DNA repair and metabolism. We observed complex patterns of ancestral admixture and putative-damaging and novel variation, both within and between populations, alongside evidence that Zambia was a likely intermediate site along the routes of expansion of Bantu-speaking populations. Pathogenic variants in genes that are currently characterized as medically relevant were uncommon-but in other genes, variants denoted as 'likely pathogenic' in the ClinVar database were commonly observed. Collectively, these findings refine our current understanding of continental migration, identify gene flow and the response to human disease as strong drivers of genome-level population variation, and underscore the scientific imperative for a broader characterization of the genomic diversity of African individuals to understand human ancestry and improve health.

Molecular Diversification of the Seminal Fluid Proteome in a Recently Diverged Passerine Species Pair

Seminal fluid proteins (SFPs) mediate an array of postmating reproductive processes that influence fertilization and fertility. As such, it is widely held that SFPs may contribute to postmating, prezygotic reproductive barriers between closely related taxa. We investigated seminal fluid (SF) diversification in a recently diverged passerine species pair (Passer domesticus and Passer hispaniolensis) using a combination of proteomic and comparative evolutionary genomic approaches. First, we characterized and compared the SF proteome of the two species, revealing consistencies with known aspects of SFP biology and function in other taxa, including the presence and diversification of proteins involved in immunity and sperm maturation. Second, using whole-genome resequencing data, we assessed patterns of genomic differentiation between house and Spanish sparrows. These analyses detected divergent selection on immunity-related SF genes and positive selective sweeps in regions containing a number of SF genes that also exhibited protein abundance diversification between species. Finally, we analyzed the molecular evolution of SFPs across 11 passerine species and found a significantly higher rate of positive selection in SFPs compared with the rest of the genome, as well as significant enrichments for functional pathways related to immunity in the set of positively selected SF genes. Our results suggest that selection on immunity pathways is an important determinant of passerine SF composition and evolution. Assessing the role of immunity genes in speciation in other recently diverged taxa should be prioritized given the potential role for immunity-related proteins in reproductive incompatibilities in Passer sparrows.

Causal inference in genetic trio studies

We introduce a method to draw causal inferences-inferences immune to all possible confounding-from genetic data that include parents and offspring. Causal conclusions are possible with these data because the natural randomness in meiosis can be viewed as a high-dimensional randomized experiment. We make this observation actionable by developing a conditional independence test that identifies regions of the genome containing distinct causal variants. The proposed digital twin test compares an observed offspring to carefully constructed synthetic offspring from the same parents to determine statistical significance, and it can leverage any black-box multivariate model and additional nontrio genetic data to increase power. Crucially, our inferences are based only on a well-established mathematical model of recombination and make no assumptions about the relationship between the genotypes and phenotypes. We compare our method to the widely used transmission disequilibrium test and demonstrate enhanced power and localization.

Genetic variation associated with childhood and adult stature and risk of MYCN-amplified neuroblastoma

Background

Neuroblastoma is the most common pediatric solid tumor. MYCN‐amplification is an important negative prognostic indicator and inherited genetic contributions to risk are incompletely understood. Genetic determinants of stature increase risk of several adult and childhood cancers, but have not been studied in neuroblastoma despite elevated neuroblastoma incidence in children with congenital overgrowth syndromes.

Methods

We investigated the association between genetic determinants of height and neuroblastoma risk in 1538 neuroblastoma cases, stratified by MYCN‐amplification status, and compared to 3390 European‐ancestry controls using polygenic scores for birth length (five variants), childhood height (six variants), and adult height (413 variants). We further examined the UK Biobank to evaluate the association of known neuroblastoma risk loci and stature.

Results

An increase in the polygenic score for childhood stature, corresponding to a ~0.5 cm increase in pre‐pubertal height, was associated with greater risk of MYCN‐amplified neuroblastoma (OR = 1.14, P = .047). An increase in the polygenic score for adult stature, corresponding to a ~1.7 cm increase in adult height attainment, was associated with decreased risk of MYCN‐amplified neuroblastoma (OR = 0.87, P = .047). These associations persisted in case‐case analyses comparing MYCN‐amplified to MYCN‐unamplified neuroblastoma. No polygenic height scores were associated with MYCN‐unamplified neuroblastoma risk. Previously identified genome‐wide association study hits for neuroblastoma (N = 10) were significantly enriched for association with both childhood (P = 4.0 × 10⁻³) and adult height (P = 8.9 × 10⁻³) in >250 000 UK Biobank study participants.

Conclusions

Genetic propensity to taller childhood height and shorter adult height were associated with MYCN‐amplified neuroblastoma risk, suggesting that biological pathways affecting growth trajectories and pubertal timing may contribute to MYCN‐amplified neuroblastoma etiology.

Transcriptomic signatures across human tissues identify functional rare genetic variation

Rare genetic variants are abundant across the human genome, and identifying their function and phenotypic impact is a major challenge. Measuring aberrant gene expression has aided in identifying functional, large-effect rare variants (RVs). Here, we expanded detection of genetically driven transcriptome abnormalities by analyzing gene expression, allele-specific expression, and alternative splicing from multitissue RNA-sequencing data, and demonstrate that each signal informs unique classes of RVs. We developed Watershed, a probabilistic model that integrates multiple genomic and transcriptomic signals to predict variant function, validated these predictions in additional cohorts and through experimental assays, and used them to assess RVs in the UK Biobank, the Million Veterans Program, and the Jackson Heart Study. Our results link thousands of RVs to diverse molecular effects and provide evidence to associate RVs affecting the transcriptome with human traits.

Polygenic risk score for the prediction of breast cancer is related to lesser terminal duct lobular unit involution of the breast

Terminal duct lobular units (TDLUs) are the predominant anatomical structures where breast cancers originate. Having lesser degrees of age-related TDLU involution, measured as higher TDLUs counts or more epithelial TDLU substructures (acini), is related to increased breast cancer risk among women with benign breast disease (BBD). We evaluated whether a recently developed polygenic risk score (PRS) based on 313-common variants for breast cancer prediction is related to TDLU involution in the background, normal breast tissue, as this could provide mechanistic clues on the genetic predisposition to breast cancer. Among 1398 women without breast cancer, higher values of the PRS were significantly associated with higher TDLU counts (P = 0.004), but not with acini counts (P = 0.808), in histologically normal tissue samples from donors and diagnostic BBD biopsies. Mediation analysis indicated that TDLU counts may explain a modest proportion (≤10%) of the association of the 313-variant PRS with breast cancer risk. These findings suggest that TDLU involution might be an intermediate step in the association between common genetic variation and breast cancer risk.

Low-frequency variation near common germline susceptibility loci are associated with risk of Ewing sarcoma

BACKGROUND

Ewing sarcoma (EwS) is a rare, aggressive solid tumor of childhood, adolescence and young adulthood associated with pathognomonic EWSR1-ETS fusion oncoproteins altering transcriptional regulation. Genome-wide association studies (GWAS) have identified 6 common germline susceptibility loci but have not investigated low-frequency inherited variants with minor allele frequencies below 5% due to limited genotyped cases of this rare tumor.

METHODS

We investigated the contribution of rare and low-frequency variation to EwS susceptibility in the largest EwS genome-wide association study to date (733 EwS cases and 1,346 unaffected controls of European ancestry).

RESULTS

We identified two low-frequency variants, rs112837127 and rs2296730, on chromosome 20 that were associated with EwS risk (OR = 0.186 and 2.038, respectively; P-value < 5×10-8) and located near previously reported common susceptibility loci. After adjusting for the most associated common variant at the locus, only rs112837127 remained a statistically significant independent signal (OR = 0.200, P-value = 5.84×10-8).

CONCLUSIONS

These findings suggest rare variation residing on common haplotypes are important contributors to EwS risk.

IMPACT

Motivate future targeted sequencing studies for a comprehensive evaluation of low-frequency and rare variation around common EwS susceptibility loci.

A comparison of humans and baboons suggests germline mutation rates do not track cell divisions

In humans, most germline mutations are inherited from the father. This observation has been widely interpreted as reflecting the replication errors that accrue during spermatogenesis. If so, the male bias in mutation should be substantially lower in a closely related species with similar rates of spermatogonial stem cell divisions but a shorter mean age of reproduction. To test this hypothesis, we resequenced two 3-4 generation nuclear families (totaling 29 individuals) of olive baboons (Papio anubis), who reproduce at approximately 10 years of age on average, and analyzed the data in parallel with three 3-generation human pedigrees (26 individuals). We estimated a mutation rate per generation in baboons of 0.57×10-8 per base pair, approximately half that of humans. Strikingly, however, the degree of male bias in germline mutations is approximately 4:1, similar to that of humans-indeed, a similar male bias is seen across mammals that reproduce months, years, or decades after birth. These results mirror the finding in humans that the male mutation bias is stable with parental ages and cast further doubt on the assumption that germline mutations track cell divisions. Our mutation rate estimates for baboons raise a further puzzle, suggesting a divergence time between apes and Old World monkeys of 65 million years, too old to be consistent with the fossil record; reconciling them now requires not only a slowdown of the mutation rate per generation in humans but also in baboons.

AlphaFamImpute: high-accuracy imputation in full-sib families from genotype-by-sequencing data

SUMMARY

AlphaFamImpute is an imputation package for calling, phasing and imputing genome-wide genotypes in outbred full-sib families from single nucleotide polymorphism (SNP) array and genotype-by-sequencing (GBS) data. GBS data are increasingly being used to genotype individuals, especially when SNP arrays do not exist for a population of interest. Low-coverage GBS produces data with a large number of missing or incorrect naïve genotype calls, which can be improved by identifying shared haplotype segments between full-sib individuals. Here, we present AlphaFamImpute, an algorithm specifically designed to exploit the genetic structure of full-sib families. It performs imputation using a two-step approach. In the first step, it phases and imputes parental genotypes based on the segregation states of their offspring (i.e. which pair of parental haplotypes the offspring inherited). In the second step, it phases and imputes the offspring genotypes by detecting which haplotype segments the offspring inherited from their parents. With a series of simulations, we find that AlphaFamImpute obtains high-accuracy genotypes, even when the parents are not genotyped and individuals are sequenced at <1x coverage.

AVAILABILITY AND IMPLEMENTATION

AlphaFamImpute is available as a Python package from the AlphaGenes website http://www.AlphaGenes.roslin.ed.ac.uk/AlphaFamImpute.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Leveraging genome-wide data to investigate differences between opioid use vs. opioid dependence in 41,176 individuals from the Psychiatric Genomics Consortium

To provide insights into the biology of opioid dependence (OD) and opioid use (i.e., exposure, OE), we completed a genome-wide analysis comparing 4503 OD cases, 4173 opioid-exposed controls, and 32,500 opioid-unexposed controls, including participants of European and African descent (EUR and AFR, respectively). Among the variants identified, rs9291211 was associated with OE (exposed vs. unexposed controls; EUR z = -5.39, p = 7.2 × 10-8). This variant regulates the transcriptomic profiles of SLC30A9 and BEND4 in multiple brain tissues and was previously associated with depression, alcohol consumption, and neuroticism. A phenome-wide scan of rs9291211 in the UK Biobank (N > 360,000) found association of this variant with propensity to use dietary supplements (p = 1.68 × 10-8). With respect to the same OE phenotype in the gene-based analysis, we identified SDCCAG8 (EUR + AFR z = 4.69, p = 10-6), which was previously associated with educational attainment, risk-taking behaviors, and schizophrenia. In addition, rs201123820 showed a genome-wide significant difference between OD cases and unexposed controls (AFR z = 5.55, p = 2.9 × 10-8) and a significant association with musculoskeletal disorders in the UK Biobank (p = 4.88 × 10-7). A polygenic risk score (PRS) based on a GWAS of risk-tolerance (n = 466,571) was positively associated with OD (OD vs. unexposed controls, p = 8.1 × 10-5; OD cases vs. exposed controls, p = 0.054) and OE (exposed vs. unexposed controls, p = 3.6 × 10-5). A PRS based on a GWAS of neuroticism (n = 390,278) was positively associated with OD (OD vs. unexposed controls, p = 3.2 × 10-5; OD vs. exposed controls, p = 0.002) but not with OE (p = 0.67). Our analyses highlight the difference between dependence and exposure and the importance of considering the definition of controls in studies of addiction.

A pooled genome-wide association study identifies pancreatic cancer susceptibility loci on chromosome 19p12 and 19p13.3 in the full-Jewish population

Jews are estimated to be at increased risk of pancreatic cancer compared to non-Jews, but their observed 50-80% excess risk is not explained by known non-genetic or genetic risk factors. We conducted a GWAS in a case-control sample of American Jews, largely Ashkenazi, including 406 pancreatic cancer patients and 2332 controls, identified in the dbGaP, PanScan I/II, PanC4 and GERA data sets. We then examined resulting SNPs with P < 10^-7 in an expanded sample set, of 539 full- or part-Jewish pancreatic cancer patients and 4117 full- or part-Jewish controls from the same data sets. Jewish ancestries were genetically determined using seeded FastPCA. Among the full Jews, a novel genome-wide significant association was detected on chromosome 19p12 (rs66562280, per-allele OR = 1.55, 95% CI = 1.33-1.81, P = 10^-7.6). A suggestive relatively independent association was detected on chromosome 19p13.3 (rs2656937, OR = 1.53, 95% CI = 1.31-1.78, P = 10^-7.0). Similar associations were seen for these SNPs among the full and part Jews combined. This is the first GWAS conducted for pancreatic cancer in the increased-risk Jewish population. The SNPs rs66562280 and rs2656937 are located in introns of ZNF100-like and ARRDC5, respectively, and are known to alter regulatory motifs of genes that play integral roles in pancreatic carcinogenesis.

Patterns of genomic divergence and introgression between Japanese stickleback species with overlapping breeding habitats

With only a few absolute geographic barriers in marine environments, the factors maintaining reproductive isolation among marine organisms remain elusive. However, spatial structuring in breeding habitat can contribute to reproductive isolation. This is particularly important for marine organisms that migrate to use fresh- or brackish water environments to breed. The Japanese Gasterosteus stickleback species, the Pacific Ocean three-spined stickleback (G. aculeatus) and the Japan Sea stickleback (G. nipponicus) overwinter in the sea, but migrate to rivers for spawning. Although they co-occur at several locations across the Japanese islands, they are reproductively isolated. Our previous studies in Bekanbeushi River showed that the Japan Sea stickleback spawns in the estuary, while the Pacific Ocean stickleback mainly spawns further upstream in freshwater. Overall genomic divergence was very high with many interspersed regions of introgression. Here, we investigated genomic divergence and introgression between the sympatric species in the much shorter Tokotan River, where they share spawning sites. The levels of genome-wide divergence were reduced and introgression was increased, suggesting that habitat isolation substantially contributes to a reduction in gene flow. We also found that genomic regions of introgression were largely shared between the two systems. Furthermore, some regions of introgression were located near loci with a heterozygote advantage for juvenile survival. Taken together, introgression may be partially driven by adaptation in this system. Although, the two species remain clearly genetically differentiated. Regions with low recombination rates showed especially low introgression. Speciation reversal is therefore likely prevented by barriers other than habitat isolation.

Extending long-range phasing and haplotype library imputation algorithms to large and heterogeneous datasets

BACKGROUND

We describe the latest improvements to the long-range phasing (LRP) and haplotype library imputation (HLI) algorithms for successful phasing of both datasets with one million individuals and datasets genotyped using different sets of single nucleotide polymorphisms (SNPs). Previous publicly available implementations of the LRP algorithm implemented in AlphaPhase could not phase large datasets due to the computational cost of defining surrogate parents by exhaustive all-against-all searches. Furthermore, the AlphaPhase implementations of LRP and HLI were not designed to deal with large amounts of missing data that are inherent when using multiple SNP arrays.

METHODS

We developed methods that avoid the need for all-against-all searches by performing LRP on subsets of individuals and then concatenating the results. We also extended LRP and HLI algorithms to enable the use of different sets of markers, including missing values, when determining surrogate parents and identifying haplotypes. We implemented and tested these extensions in an updated version of AlphaPhase, and compared its performance to the software package Eagle2.

RESULTS

A simulated dataset with one million individuals genotyped with the same 6711 SNPs for a single chromosome took less than a day to phase, compared to more than seven days for Eagle2. The percentage of correctly phased alleles at heterozygous loci was 90.2 and 99.9% for AlphaPhase and Eagle2, respectively. A larger dataset with one million individuals genotyped with 49,579 SNPs for a single chromosome took AlphaPhase 23 days to phase, with 89.9% of alleles at heterozygous loci phased correctly. The phasing accuracy was generally lower for datasets with different sets of markers than with one set of markers. For a simulated dataset with three sets of markers, 1.5% of alleles at heterozygous positions were phased incorrectly, compared to 0.4% with one set of markers.

CONCLUSIONS

The improved LRP and HLI algorithms enable AlphaPhase to quickly and accurately phase very large and heterogeneous datasets. AlphaPhase is an order of magnitude faster than the other tested packages, although Eagle2 showed a higher level of phasing accuracy. The speed gain will make phasing achievable for very large genomic datasets in livestock, enabling more powerful breeding and genetics research and application.

Mitochondrial GWAS and association of nuclear - mitochondrial epistasis with BMI in T1DM patients

BACKGROUND

BMI is a strong indicator of complications from type I diabetes, especially under intensive treatment.

METHODS

We have genotyped 435 type 1 diabetics using Illumina Infinium Omni Express Exome-8 v1.4 arrays and performed mitoGWAS on BMI. We identified additive interactions between mitochondrial and nuclear variants in genes associated with mitochondrial functioning MitoCarta2.0 and confirmed and refined the results on external cohorts: the Framingham Heart Study (FHS) and GTEx data. Linear mixed model analysis was performed using the GENESIS package in R/Bioconductor.

RESULTS

We find a borderline significant association between the mitochondrial variant rs28357980, localized to MT-ND2, and BMI (β = - 0.69, p = 0.056). This BMI association was confirmed on 1889 patients from FHS cohort (β = - 0.312, p = 0.047). Next, we searched for additive interactions between mitochondrial and nuclear variants. MT-ND2 variants interacted with variants in the genes SIRT3, ATP5B, CYCS, TFB2M and POLRMT. TFB2M is a mitochondrial transcription factor and together with TFAM creates a transcription promoter complex for the mitochondrial polymerase POLRMT. We have found an interaction between rs3021088 in MT-ND2 and rs6701836 in TFB2M leading to BMI decrease (inter_pval = 0.0241), while interaction of rs3021088 in MT-ND2 and rs41542013 in POLRMT led to BMI increase (inter_pval = 0.0004). The influence of these interactions on BMI was confirmed in external cohorts.

CONCLUSIONS

Here, we have shown that variants in the mitochondrial genome as well as additive interactions between mitochondrial and nuclear SNPs influence BMI in T1DM and general cohorts.

Population-Specific Recombination Maps from Segments of Identity by Descent

Recombination rates vary significantly across the genome, and estimates of recombination rates are needed for downstream analyses such as haplotype phasing and genotype imputation. Existing methods for recombination rate estimation are limited by insufficient amounts of informative genetic data or by high computational cost. We present a method and software, called IBDrecomb, for using segments of identity by descent to infer recombination rates. IBDrecomb can be applied to sequenced population cohorts to obtain high-resolution, population-specific recombination maps. In simulated admixed data, IBDrecomb obtains higher accuracy than admixture-based estimation of recombination rates. When applied to 2,500 simulated individuals, IBDrecomb obtains similar accuracy to a linkage-disequilibrium (LD)-based method applied to 96 individuals (the largest number for which computation is tractable). Compared to LD-based maps, our IBD-based maps have the advantage of estimating recombination rates in the recent past rather than the distant past. We used IBDrecomb to generate new recombination maps for European Americans and for African Americans from TOPMed sequence data from the Framingham Heart Study (1,626 unrelated individuals) and the Jackson Heart Study (2,046 unrelated individuals), and we compare them to LD-based, admixture-based, and family-based maps.

Native American gene flow into Polynesia predating Easter Island settlement

The possibility of voyaging contact between prehistoric Polynesian and Native American populations has long intrigued researchers. Proponents have pointed to the existence of New World crops, such as the sweet potato and bottle gourd, in the Polynesian archaeological record, but nowhere else outside the pre-Columbian Americas1-6, while critics have argued that these botanical dispersals need not have been human mediated7. The Norwegian explorer Thor Heyerdahl controversially suggested that prehistoric South American populations had an important role in the settlement of east Polynesia and particularly of Easter Island (Rapa Nui)2. Several limited molecular genetic studies have reached opposing conclusions, and the possibility continues to be as hotly contested today as it was when first suggested8-12. Here we analyse genome-wide variation in individuals from islands across Polynesia for signs of Native American admixture, analysing 807 individuals from 17 island populations and 15 Pacific coast Native American groups. We find conclusive evidence for prehistoric contact of Polynesian individuals with Native American individuals (around AD 1200) contemporaneous with the settlement of remote Oceania13-15. Our analyses suggest strongly that a single contact event occurred in eastern Polynesia, before the settlement of Rapa Nui, between Polynesian individuals and a Native American group most closely related to the indigenous inhabitants of present-day Colombia.

Oligonucleotide hybridization with magnetic separation assay for multiple SNP phasing

Since humans have two copies of each gene, multiple mutations in different loci may or may not be found on the same strand of DNA (i.e., inherited from one parent). When a person is heterozygous at more than one position, the placement of these mutations, also called the haplotype phase, (i.e., cis for the same strand and trans for different strands) can result in the expression of different amount and type of proteins. In this work, we described an enzyme-free method to phase two single nucleotide polymorphisms (SNPs) using two fluorophore/quencher-labelled probes, where one of which was biotinylated. The fluorescence signal was obtained twice: first, after the addition of the labelled probes and second, after the addition of the magnetic beads. The first signal was shown to be proportional to the total number of SNP A and SNP B present in the target analyte, while the second signal showed a marked decrease of the fluorescence signal from the non-biotinylated probe when the SNPs were in trans, showing that the probe immobilized on the magnetic bead selectively captures targets with SNPs in a cis configuration. We then mimic the nature of the human genome which consists of two haplotype copies of each gene, and showed that 250 nM of the 10 possible pairs of haplotypes could be differentiated using a combination of fluorescence microscopy and fluorescence detection.

Linking protein to phenotype with Mendelian Randomization detects 38 proteins with causal roles in human diseases and traits

To efficiently transform genetic associations into drug targets requires evidence that a particular gene, and its encoded protein, contribute causally to a disease. To achieve this, we employ a three-step proteome-by-phenome Mendelian Randomization (MR) approach. In step one, 154 protein quantitative trait loci (pQTLs) were identified and independently replicated. From these pQTLs, 64 replicated locally-acting variants were used as instrumental variables for proteome-by-phenome MR across 846 traits (step two). When its assumptions are met, proteome-by-phenome MR, is equivalent to simultaneously running many randomized controlled trials. Step 2 yielded 38 proteins that significantly predicted variation in traits and diseases in 509 instances. Step 3 revealed that amongst the 271 instances from GeneAtlas (UK Biobank), 77 showed little evidence of pleiotropy (HEIDI), and 92 evidence of colocalization (eCAVIAR). Results were wide ranging: including, for example, new evidence for a causal role of tyrosine-protein phosphatase non-receptor type substrate 1 (SHPS1; SIRPA) in schizophrenia, and a new finding that intestinal fatty acid binding protein (FABP2) abundance contributes to the pathogenesis of cardiovascular disease. We also demonstrated confirmatory evidence for the causal role of four further proteins (FGF5, IL6R, LPL, LTA) in cardiovascular disease risk.

SNPs in lncRNA Regions and Breast Cancer Risk

Long non-coding RNAs (lncRNAs) play crucial roles in human physiology, and have been found to be associated with various cancers. Transcribed ultraconserved regions (T-UCRs) are a subgroup of lncRNAs conserved in several species, and are often located in cancer-related regions. Breast cancer is the most common cancer in women worldwide and the leading cause of female cancer deaths. We investigated the association of genetic variants in lncRNA and T-UCR regions with breast cancer risk to uncover candidate loci for further analysis. Our focus was on low-penetrance variants that can be discovered in a large dataset. We selected 565 regions of lncRNAs and T-UCRs that are expressed in breast or breast cancer tissue, or show expression correlation to major breast cancer associated genes. We studied the association of single nucleotide polymorphisms (SNPs) in these regions with breast cancer risk in the 122970 case samples and 105974 controls of the Breast Cancer Association Consortium's genome-wide data, and also by in silico functional analyses using Integrated Expression Quantitative trait and in silico prediction of GWAS targets (INQUISIT) and expression quantitative trait loci (eQTL) analysis. The eQTL analysis was carried out using the METABRIC dataset and analyses from GTEx and ncRNA eQTL databases. We found putative breast cancer risk variants (p < 1 × 10-5) targeting the lncRNA GABPB1-AS1 in INQUISIT and eQTL analysis. In addition, putative breast cancer risk associated SNPs (p < 1 × 10-5) in the region of two T-UCRs, uc.184 and uc.313, located in protein coding genes CPEB4 and TIAL1, respectively, targeted these genes in INQUISIT and in eQTL analysis. Other non-coding regions containing SNPs with the defined p-value and highly significant false discovery rate (FDR) for breast cancer risk association were discovered that may warrant further studies. These results suggest candidate lncRNA loci for further research on breast cancer risk and the molecular mechanisms.

Genome-wide association study identifies genetic susceptibility loci and pathways of radiation-induced acute oral mucositis

Background

Radiation-induced oral mucositis (OM) is one of the most common acute complications for head and neck cancer. Severe OM is associated with radiation treatment breaks, which harms successful tumor management. Radiogenomics studies have indicated that genetic variants are associated with adverse effects of radiotherapy.

Methods

A large-scale genome-wide scan was performed in 1467 nasopharyngeal carcinoma patients, including 753 treated with 2D-CRT from Genetic Architecture of the Radiotherapy Toxicity and Prognosis (GARTP) cohort and 714 treated with IMRT (192 from the GARTP and 522 newly recruited). Subgroup analysis by radiotherapy technique was further performed in the top associations. We also performed physical and regulatory mapping of the risk loci and gene set enrichment analysis of the candidate target genes.

Results

We identified 50 associated genomic loci and 64 genes via positional mapping, expression quantitative trait locus (eQTL) mapping, chromatin interaction mapping and gene-based analysis, and 36 of these loci were replicated in subgroup analysis. Interestingly, one of the top loci located in TNKS, a gene relevant to radiation toxicity, was associated with increased OM risk with OR = 3.72 of the lead SNP rs117157809 (95% CI 2.10–6.57; P = 6.33 × 10⁻⁶). Gene set analyses showed that the 64 candidate target genes were enriched in the biological processes of regulating telomere capping and maintenance and telomerase activity (Top P = 7.73 × 10⁻⁷).

Conclusions

These results enhance the biological understanding of radiotherapy toxicity. The association signals enriched in telomere function regulation implicate the potential underlying mechanism and warrant further functional investigation and potential individual radiotherapy applications.

Identification of novel epithelial ovarian cancer loci in women of African ancestry

Women of African ancestry have lower incidence of epithelial ovarian cancer (EOC) yet worse survival compared to women of European ancestry. We conducted a genome-wide association study in African ancestry women with 755 EOC cases, including 537 high-grade serous ovarian carcinomas (HGSOC) and 1,235 controls. We identified four novel loci with suggestive evidence of association with EOC (p < 1 × 10-6 ), including rs4525119 (intronic to AKR1C3), rs7643459 (intronic to LOC101927394), rs4286604 (12 kb 3' of UGT2A2) and rs142091544 (5 kb 5' of WWC1). For HGSOC, we identified six loci with suggestive evidence of association including rs37792 (132 kb 5' of follistatin [FST]), rs57403204 (81 kb 3' of MAGEC1), rs79079890 (LOC105376360 intronic), rs66459581 (5 kb 5' of PRPSAP1), rs116046250 (GABRG3 intronic) and rs192876988 (32 kb 3' of GK2). Among the identified variants, two are near genes known to regulate hormones and diseases of the ovary (AKR1C3 and FST), and two are linked to cancer (AKR1C3 and MAGEC1). In follow-up studies of the 10 identified variants, the GK2 region SNP, rs192876988, showed an inverse association with EOC in European ancestry women (p = 0.002), increased risk of ER positive breast cancer in African ancestry women (p = 0.027) and decreased expression of GK2 in HGSOC tissue from African ancestry women (p = 0.004). A European ancestry-derived polygenic risk score showed positive associations with EOC and HGSOC in women of African ancestry suggesting shared genetic architecture. Our investigation presents evidence of variants for EOC shared among European and African ancestry women and identifies novel EOC risk loci in women of African ancestry.

XAF1 as a modifier of p53 function and cancer susceptibility

Cancer risk is highly variable in carriers of the common TP53-R337H founder allele, possibly due to the influence of modifier genes. Whole-genome sequencing identified a variant in the tumor suppressor XAF1 (E134*/Glu134Ter/rs146752602) in a subset of R337H carriers. Haplotype-defining variants were verified in 203 patients with cancer, 582 relatives, and 42,438 newborns. The compound mutant haplotype was enriched in patients with cancer, conferring risk for sarcoma (P = 0.003) and subsequent malignancies (P = 0.006). Functional analyses demonstrated that wild-type XAF1 enhances transactivation of wild-type and hypomorphic TP53 variants, whereas XAF1-E134* is markedly attenuated in this activity. We propose that cosegregation of XAF1-E134* and TP53-R337H mutations leads to a more aggressive cancer phenotype than TP53-R337H alone, with implications for genetic counseling and clinical management of hypomorphic TP53 mutant carriers.

Distinct subtypes of polycystic ovary syndrome with novel genetic associations: An unsupervised, phenotypic clustering analysis

BACKGROUND

Polycystic ovary syndrome (PCOS) is a common, complex genetic disorder affecting up to 15% of reproductive-age women worldwide, depending on the diagnostic criteria applied. These diagnostic criteria are based on expert opinion and have been the subject of considerable controversy. The phenotypic variation observed in PCOS is suggestive of an underlying genetic heterogeneity, but a recent meta-analysis of European ancestry PCOS cases found that the genetic architecture of PCOS defined by different diagnostic criteria was generally similar, suggesting that the criteria do not identify biologically distinct disease subtypes. We performed this study to test the hypothesis that there are biologically relevant subtypes of PCOS.

METHODS AND FINDINGS

Using biochemical and genotype data from a previously published PCOS genome-wide association study (GWAS), we investigated whether there were reproducible phenotypic subtypes of PCOS with subtype-specific genetic associations. Unsupervised hierarchical cluster analysis was performed on quantitative anthropometric, reproductive, and metabolic traits in a genotyped cohort of 893 PCOS cases (median and interquartile range [IQR]: age = 28 [25-32], body mass index [BMI] = 35.4 [28.2-41.5]). The clusters were replicated in an independent, ungenotyped cohort of 263 PCOS cases (median and IQR: age = 28 [24-33], BMI = 35.7 [28.4-42.3]). The clustering revealed 2 distinct PCOS subtypes: a "reproductive" group (21%-23%), characterized by higher luteinizing hormone (LH) and sex hormone binding globulin (SHBG) levels with relatively low BMI and insulin levels, and a "metabolic" group (37%-39%), characterized by higher BMI, glucose, and insulin levels with lower SHBG and LH levels. We performed a GWAS on the genotyped cohort, limiting the cases to either the reproductive or metabolic subtypes. We identified alleles in 4 loci that were associated with the reproductive subtype at genome-wide significance (PRDM2/KAZN, P = 2.2 × 10-10; IQCA1, P = 2.8 × 10-9; BMPR1B/UNC5C, P = 9.7 × 10-9; CDH10, P = 1.2 × 10-8) and one locus that was significantly associated with the metabolic subtype (KCNH7/FIGN, P = 1.0 × 10-8). We developed a predictive model to classify a separate, family-based cohort of 73 women with PCOS (median and IQR: age = 28 [25-33], BMI = 34.3 [27.8-42.3]) and found that the subtypes tended to cluster in families and that carriers of previously reported rare variants in DENND1A, a gene that regulates androgen biosynthesis, were significantly more likely to have the reproductive subtype of PCOS. Limitations of our study were that only PCOS cases of European ancestry diagnosed by National Institutes of Health (NIH) criteria were included, the sample sizes for the subtype GWAS were small, and the GWAS findings were not replicated.

CONCLUSIONS

In conclusion, we have found reproducible reproductive and metabolic subtypes of PCOS. Furthermore, these subtypes were associated with novel, to our knowledge, susceptibility loci. Our results suggest that these subtypes are biologically relevant because they appear to have distinct genetic architecture. This study demonstrates how phenotypic subtyping can be used to gain additional insights from GWAS data.

Leveraging Genome and Phenome-Wide Association Studies to Investigate Genetic Risk of Acute Lymphoblastic Leukemia

BACKGROUND

Genome-wide association studies (GWAS) of childhood cancers remain limited, highlighting the need for novel analytic strategies. We describe a hybrid GWAS and phenome-wide association study (PheWAS) approach to uncover genotype-phenotype relationships and candidate risk loci, applying it to acute lymphoblastic leukemia (ALL).

METHODS

PheWAS was performed for 12 ALL SNPs identified by prior GWAS and two control SNP-sets using UK Biobank data. PheWAS-traits significantly associated with ALL SNPs compared with control SNPs were assessed for association with ALL risk (959 cases, 2,624 controls) using polygenic score and Mendelian randomization analyses. Trait-associated SNPs were tested for association with ALL risk in single-SNP analyses, with replication in an independent case-control dataset (1,618 cases, 9,409 controls).

RESULTS

Platelet count was the trait most enriched for association with known ALL risk loci. A polygenic score for platelet count (223 SNPs) was not associated with ALL risk (P = 0.82) and Mendelian randomization did not suggest a causal relationship. However, twelve platelet count-associated SNPs were nominally associated with ALL risk in COG data and three were replicated in UK data (rs10058074, rs210142, rs2836441).

CONCLUSIONS

In our hybrid GWAS-PheWAS approach, we identify pleiotropic genetic variation contributing to ALL risk and platelet count. Three SNPs known to influence platelet count were reproducibly associated with ALL risk, implicating genomic regions containing IRF1, proapoptotic protein BAK1, and ERG in platelet production and leukemogenesis.

IMPACT

Incorporating PheWAS data into association studies can leverage genetic pleiotropy to identify cancer risk loci, highlighting the utility of our novel approach.

Family-based genome-wide association study of leprosy in Vietnam

Leprosy is a chronic infectious disease of the skin and peripheral nerves with a strong genetic predisposition. Recent genome-wide approaches have identified numerous common variants associated with leprosy, almost all in the Chinese population. We conducted the first family-based genome-wide association study of leprosy in 622 affected offspring from Vietnam, followed by replication in an independent sample of 1181 leprosy cases and 668 controls of the same ethnic origin. The most significant results were observed within the HLA region, in which six SNPs displayed genome-wide significant associations, all of which were replicated in the independent case/control sample. We investigated the signal in the HLA region in more detail, by conducting a multivariate analysis on the case/control sample of 319 GWAS-suggestive HLA hits for which evidence for replication was obtained. We identified three independently associated SNPs, two located in the HLA class I region (rs1265048: OR = 0.69 [0.58-0.80], combined p-value = 5.53x10-11; and rs114598080: OR = 1.47 [1.46-1.48], combined p-value = 8.77x10-13), and one located in the HLA class II region (rs3187964 (OR = 1.67 [1.55-1.80], combined p-value = 8.35x10-16). We also validated two previously identified risk factors for leprosy: the missense variant rs3764147 in the LACC1 gene (OR = 1.52 [1.41-1.63], combined p-value = 5.06x10-14), and the intergenic variant rs6871626 located close to the IL12B gene (OR = 0.73 [0.61-0.84], combined p-value = 6.44x10-8). These results shed new light on the genetic control of leprosy, by dissecting the influence of HLA SNPs, and validating the independent role of two additional variants in a large Vietnamese sample.

On the origin of European sheep as revealed by the diversity of the Balkan breeds and by optimizing population-genetic analysis tools

Background

In the Neolithic, domestic sheep migrated into Europe and subsequently spread in westerly and northwesterly directions. Reconstruction of these migrations and subsequent genetic events requires a more detailed characterization of the current phylogeographic differentiation.

Results

We collected 50 K single nucleotide polymorphism (SNP) profiles of Balkan sheep that are currently found near the major Neolithic point of entry into Europe, and combined these data with published genotypes from southwest-Asian, Mediterranean, central-European and north-European sheep and from Asian and European mouflons. We detected clines, ancestral components and admixture by using variants of common analysis tools: geography-informative supervised principal component analysis (PCA), breed-specific admixture analysis, across-breed \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_{4}$$\end{document}f4 profiles and phylogenetic analysis of regional pools of breeds. The regional Balkan sheep populations exhibit considerable genetic overlap, but are clearly distinct from the breeds in surrounding regions. The Asian mouflon did not influence the differentiation of the European domestic sheep and is only distantly related to present-day sheep, including those from Iran where the mouflons were sampled. We demonstrate the occurrence, from southeast to northwest Europe, of a continuously increasing ancestral component of up to 20% contributed by the European mouflon, which is assumed to descend from the original Neolithic domesticates. The overall patterns indicate that the Balkan region and Italy served as post-domestication migration hubs, from which wool sheep reached Spain and north Italy with subsequent migrations northwards. The documented dispersal of Tarentine wool sheep during the Roman period may have been part of this process. Our results also reproduce the documented 18th century admixture of Spanish Merino sheep into several central-European breeds.

Conclusions

Our results contribute to a better understanding of the events that have created the present diversity pattern, which is relevant for the management of the genetic resources represented by the European sheep population.

A Rare Mutation in SPLUNC1 Affects Bacterial Adherence and Invasion in Meningococcal Disease

BACKGROUND

Neisseria meningitidis (Nm) is a nasopharyngeal commensal carried by healthy individuals. However, invasive infections occurs in a minority of individuals, with devastating consequences. There is evidence that common polymorphisms are associated with invasive meningococcal disease (IMD), but the contributions of rare variants other than those in the complement system have not been determined.

METHODS

We identified familial cases of IMD in the UK meningococcal disease study and the European Union Life-Threatening Infectious Disease Study. Candidate genetic variants were identified by whole-exome sequencing of 2 patients with familial IMD. Candidate variants were further validated by in vitro assays.

RESULTS

Exomes of 2 siblings with IMD identified a novel heterozygous missense mutation in BPIFA1/SPLUNC1. Sequencing of 186 other nonfamilial cases identified another unrelated IMD patient with the same mutation. SPLUNC1 is an innate immune defense protein expressed in the nasopharyngeal epithelia; however, its role in invasive infections is unknown. In vitro assays demonstrated that recombinant SPLUNC1 protein inhibits biofilm formation by Nm, and impedes Nm adhesion and invasion of human airway cells. The dominant negative mutant recombinant SPLUNC1 (p.G22E) showed reduced antibiofilm activity, increased meningococcal adhesion, and increased invasion of cells, compared with wild-type SPLUNC1.

CONCLUSIONS

A mutation in SPLUNC1 affecting mucosal attachment, biofilm formation, and invasion of mucosal epithelial cells is a new genetic cause of meningococcal disease.

Dissecting human North African gene-flow into its western coastal surroundings

North African history and populations have exerted a pivotal influence on surrounding geographical regions, although scant genetic studies have addressed this issue. Our aim is to understand human historical migrations in the coastal surroundings of North Africa. We built a refined genome-wide dataset of North African populations to unearth the fine-scale genetic structure of the region, using haplotype information. The results suggest that the gene-flow from North Africa into the European Mediterranean coast (Tuscany and the Iberian Peninsula) arrived mainly from the Mediterranean coast of North Africa. In Tuscany, this North African admixture date estimate suggests the movement of peoples during the fall of the Roman Empire around the fourth century. In the Iberian Peninsula, the North African component probably reflects the impact of the Arab expansion since the seventh century and the subsequent expansion of the Christian Kingdoms. By contrast, the North African component in the Canary Islands has a source genetically related to present-day people from the Atlantic North African coast. We also find sub-Saharan gene-flow from the Senegambia region in the Canary Islands. Specifically, we detect a complex signal of admixture involving Atlantic, Senegambian and European sources intermixing around the fifteenth century, soon after the Castilian conquest. Our results highlight the differential genetic influence of North Africa into the surrounding coast and show that specific historical events have not only had a socio-cultural impact but additionally modified the gene pool of the populations.