Integrating common and rare genetic variation in diverse human populations

The female protective effect against autism spectrum disorder

Autism spectrum disorder (ASD) is diagnosed three to four times more frequently in males than in females. Genetic studies of rare variants support a female protective effect (FPE) against ASD. However, sex differences in common inherited genetic risk for ASD are less studied, particularly within families. Leveraging the Danish iPSYCH resource, we found siblings of female ASD cases (n = 1,707) had higher rates of ASD than siblings of male ASD cases (n = 6,270; p < 1.0 × 10-10). In the Simons Simplex and SPARK collections, mothers of ASD cases (n = 7,436) carried more polygenic risk for ASD than fathers of ASD cases (n = 5,926; 0.08 polygenic risk score [PRS] SD; p = 7.0 × 10-7). Further, male unaffected siblings under-inherited polygenic risk (n = 1,519; p = 0.03). Using both epidemiologic and genetic approaches, our findings strongly support an FPE against ASD's common inherited influences.

Multivariate GWAS of psychiatric disorders and their cardinal symptoms reveal two dimensions of cross-cutting genetic liabilities

Understanding which biological pathways are specific versus general across diagnostic categories and levels of symptom severity is critical to improving nosology and treatment of psychopathology. Here, we combine transdiagnostic and dimensional approaches to genetic discovery for the first time, conducting a novel multivariate genome-wide association study of eight psychiatric symptoms and disorders broadly related to mood disturbance and psychosis. We identify two transdiagnostic genetic liabilities that distinguish between common forms of psychopathology versus rarer forms of serious mental illness. Biological annotation revealed divergent genetic architectures that differentially implicated prenatal neurodevelopment and neuronal function and regulation. These findings inform psychiatric nosology and biological models of psychopathology, as they suggest that the severity of mood and psychotic symptoms present in serious mental illness may reflect a difference in kind rather than merely in degree.

Dissection of multiple sclerosis genetics identifies B and CD4+ T cells as driver cell subsets

Background

Multiple sclerosis (MS) is an autoimmune condition of the central nervous system with a well-characterized genetic background. Prior analyses of MS genetics have identified broad enrichments across peripheral immune cells, yet the driver immune subsets are unclear.

Results

We utilize chromatin accessibility data across hematopoietic cells to identify cell type-specific enrichments of MS genetic signals. We find that CD4 T and B cells are independently enriched for MS genetics and further refine the driver subsets to T_h17 and memory B cells, respectively. We replicate our findings in data from untreated and treated MS patients and find that immunomodulatory treatments suppress chromatin accessibility at driver cell types. Integration of statistical fine-mapping and chromatin interactions nominate numerous putative causal genes, illustrating complex interplay between shared and cell-specific genes.

Conclusions

Overall, our study finds that open chromatin regions in CD4 T cells and B cells independently drive MS genetic signals. Our study highlights how careful integration of genetics and epigenetics can provide fine-scale insights into causal cell types and nominate new genes and pathways for disease.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02694-y.

Analysis of MRI-derived spleen iron in the UK Biobank identifies genetic variation linked to iron homeostasis and hemolysis

The spleen plays a key role in iron homeostasis. It is the largest filter of the blood and performs iron reuptake from old or damaged erythrocytes. Despite this role, spleen iron concentration has not been measured in a large, population-based cohort. In this study, we quantify spleen iron in 41,764 participants of the UK Biobank by using magnetic resonance imaging and provide a reference range for spleen iron in an unselected population. Through genome-wide association study, we identify associations between spleen iron and regulatory variation at two hereditary spherocytosis genes, ANK1 and SPTA1. Spherocytosis-causing coding mutations in these genes are associated with lower reticulocyte volume and increased reticulocyte percentage, while these common alleles are associated with increased expression of ANK1 and SPTA1 in blood and with larger reticulocyte volume and reduced reticulocyte percentage. As genetic modifiers, these common alleles may explain mild spherocytosis phenotypes that have been observed clinically. Our genetic study also identifies a signal that co-localizes with a splicing quantitative trait locus for MS4A7, and we show this gene is abundantly expressed in the spleen and in macrophages. The combination of deep learning and efficient image processing enables non-invasive measurement of spleen iron and, in turn, characterization of genetic factors related to the lytic phase of the erythrocyte life cycle and iron reuptake in the spleen.

TOP-LD: A tool to explore linkage disequilibrium with TOPMed whole-genome sequence data

Current publicly available tools that allow rapid exploration of linkage disequilibrium (LD) between markers (e.g., HaploReg and LDlink) are based on whole-genome sequence (WGS) data from 2,504 individuals in the 1000 Genomes Project. Here, we present TOP-LD, an online tool to explore LD inferred with high-coverage (∼30×) WGS data from 15,578 individuals in the NHLBI Trans-Omics for Precision Medicine (TOPMed) program. TOP-LD provides a significant upgrade compared to current LD tools, as the TOPMed WGS data provide a more comprehensive representation of genetic variation than the 1000 Genomes data, particularly for rare variants and in the specific populations that we analyzed. For example, TOP-LD encompasses LD information for 150.3, 62.2, and 36.7 million variants for European, African, and East Asian ancestral samples, respectively, offering 2.6- to 9.1-fold increase in variant coverage compared to HaploReg 4.0 or LDlink. In addition, TOP-LD includes tens of thousands of structural variants (SVs). We demonstrate the value of TOP-LD in fine-mapping at the GGT1 locus associated with gamma glutamyltransferase in the African ancestry participants in UK Biobank. Beyond fine-mapping, TOP-LD can facilitate a wide range of applications that are based on summary statistics and estimates of LD. TOP-LD is freely available online.

Integrated analysis of direct and proxy genome wide association studies highlights polygenicity of Alzheimer's disease outside of the APOE region

Recent meta-analyses combining direct genome-wide association studies (GWAS) with those of family history (GWAX) have indicated very low SNP heritability of Alzheimer's disease (AD). These low estimates may call into question the prospects of continued progress in genetic discovery for AD within the spectrum of common variants. We highlight dramatic downward biases in previous methods, and we validate a novel method for the estimation of SNP heritability via integration of GWAS and GWAX summary data. We apply our method to investigate the genetic architecture of AD using GWAX from UK Biobank and direct case-control GWAS from the International Genomics of Alzheimer's Project (IGAP). We estimate the liability scale common variant SNP heritability of Clinical AD outside of APOE region at ~7-11%, and we project the corresponding estimate for AD pathology to be up to approximately 23%. We estimate that nearly 90% of common variant SNP heritability of Clinical AD exists outside the APOE region. Rare variants not tagged in standard GWAS may account for additional variance. Our results indicate that, while GWAX for AD in UK Biobank may result in greater attenuation of genetic effects beyond that conventionally assumed, it does not introduce appreciable contamination of signal by genetically distinct traits relative to direct case-control GWAS in IGAP. Genetic risk for AD represents a strong effect of APOE superimposed upon a highly polygenic background.

Multi-phenotype analyses of hemostatic traits with cardiovascular events reveal novel genetic associations

BACKGROUND

Multi-phenotype analysis of genetically correlated phenotypes can increase the statistical power to detect loci associated with multiple traits, leading to the discovery of novel loci. This is the first study to date to comprehensively analyze the shared genetic effects within different hemostatic traits, and between these and their associated disease outcomes.

OBJECTIVES

To discover novel genetic associations by combining summary data of correlated hemostatic traits and disease events.

METHODS

Summary statistics from genome wide-association studies (GWAS) from seven hemostatic traits (factor VII [FVII], factor VIII [FVIII], von Willebrand factor [VWF] factor XI [FXI], fibrinogen, tissue plasminogen activator [tPA], plasminogen activator inhibitor 1 [PAI-1]) and three major cardiovascular (CV) events (venous thromboembolism [VTE], coronary artery disease [CAD], ischemic stroke [IS]), were combined in 27 multi-trait combinations using metaUSAT. Genetic correlations between phenotypes were calculated using Linkage Disequilibrium Score Regression (LDSC). Newly associated loci were investigated for colocalization. We considered a significance threshold of 1.85 × 10-9 obtained after applying Bonferroni correction for the number of multi-trait combinations performed (n = 27).

RESULTS

Across the 27 multi-trait analyses, we found 4 novel pleiotropic loci (XXYLT1, KNG1, SUGP1/MAU2, TBL2/MLXIPL) that were not significant in the original individual datasets, were not described in previous GWAS for the individual traits, and that presented a common associated variant between the studied phenotypes.

CONCLUSIONS

The discovery of four novel loci contributes to the understanding of the relationship between hemostasis and CV events and elucidate common genetic factors between these traits.

Using genetic variation to disentangle the complex relationship between food intake and health outcomes

Diet is considered as one of the most important modifiable factors influencing human health, but efforts to identify foods or dietary patterns associated with health outcomes often suffer from biases, confounding, and reverse causation. Applying Mendelian randomization in this context may provide evidence to strengthen causality in nutrition research. To this end, we first identified 283 genetic markers associated with dietary intake in 445,779 UK Biobank participants. We then converted these associations into direct genetic effects on food exposures by adjusting them for effects mediated via other traits. The SNPs which did not show evidence of mediation were then used for MR, assessing the association between genetically predicted food choices and other risk factors, health outcomes. We show that using all associated SNPs without omitting those which show evidence of mediation, leads to biases in downstream analyses (genetic correlations, causal inference), similar to those present in observational studies. However, MR analyses using SNPs which have only a direct effect on the exposure on food exposures provided unequivocal evidence of causal associations between specific eating patterns and obesity, blood lipid status, and several other risk factors and health outcomes.

Assessment of Bidirectional Relationships between Leisure Sedentary Behaviors and Neuropsychiatric Disorders: A Two-Sample Mendelian Randomization Study

(1) Background: Increasing evidence shows that sedentary behaviors are associated with neuropsychiatric disorders (NPDs) and thus may be a modifiable factor to target for the prevention of NPDs. However, the direction and causality for the relationship remain unknown; sedentary behaviors could increase or decrease the risk of NPDs, and/or NPDs may increase or decrease engagement in sedentary behaviors. (2) Methods: This Mendelian randomization (MR) study with two samples included independent genetic variants related to sedentary behaviors (n = 408,815), Alzheimer's disease (AD; n = 63,926), schizophrenia (SCZ; n = 105,318), and major depressive disorder (MDD; n = 500,199), which were extracted from several of the largest non-overlapping genome-wide association studies (GWASs), as instrumental variables. The summarized MR effect sizes from each instrumental variable were combined in an IVW (inverse-variance-weighted) approach, with various approaches (e.g., MR-Egger, weighted median, MR-pleiotropy residual sum and outlier), and sensitivity analyses were performed to identify and remove outliers and assess the horizontal pleiotropy. (3) Results: The MR evidence and linkage disequilibrium score regression revealed a consistent directional association between television watching and MDD (odds ratio (OR), 1.13 for MDD per one standard deviation (SD) increase in mean television watching time; 95% CI, 1.06-1.20; p = 6.80 × 10^-5) and a consistent relationship between computer use and a decrease in the risk of AD (OR, 0.52 for AD per one SD increase in mean computer use time; 95% CI, 0.32-0.84; p = 8.20 × 10^-3). In the reverse direction, MR showed a causal association between a reduced risk of SCZ and an increase in driving time (β, -0.016; 95% CI, -0.027--0.004; p = 8.30 × 10^-3). (4) Conclusions: Using genetic instrumental variables identified from large-scale GWASs, we found robust evidence for a causal relationship between long computer use time and a reduced risk of AD, and for a causal relationship between long television watching time and an increased risk of MDD. In reverse analyses, we found that SCZ was causally associated with reduced driving time. These findings fit in with our observations and prior knowledge as well as emphasizing the importance of distinguishing between different domains of sedentary behaviors in epidemiologic studies of NPDs.

Imaging and genetics in Parkinson's disease: assessment of the GBA1 mutation

INTRODUCTION

Several genetic variants are associated with an increased risk for developing Parkinson's Disease (PD) and limited genotype/phenotype correlation. Specifically, mutations in GBA1, the gene coding for the lysosomal enzyme glucocerebrosidase, are associated with an earlier age of onset and faster disease progression. Given these phenotypic differences associated with GBA1 variants, we explored whether cortical thickness and other biomarkers of neurodegeneration differed in healthy controls and PD patients with and without GBA1 variants.

METHODS

To understand how different GBA1 variants influence PD phenotype early in the disease, we retrieved neuroimaging and biospecimen data from the Parkinson's Progression Markers Initiative database. Using FreeSurfer, we compared T1-weighted MRI images from healthy controls (N = 47) to PD patients with heterozygous N370S (N = 21), heterozygous E326K (N = 18) or heterozygous T369M (N = 8) variants, and GBA1 non-mutation carriers (N = 47).

RESULTS

Cortical thickness in PD patients differed from controls in the parietal cortex, with E365K, T369M variants, and GBA1 non-mutation carriers showing more cortical thinning than N370S variants. Patients with N370S variants had significantly higher serum neurofilament light levels among all groups.

CONCLUSION

Our results demonstrate significant cortical thinning in PD patients independent of genotype in superior parietal and postcentral regions when compared to the controls. They highlight the impact of GBA1 variants on cortical thickness in the parietal cortex. Finally, they suggest that recently diagnosed PD patients with N370S variants have a higher cortical thickness and increased active neurodegeneration when compared to PD patients without GBA1 mutations and PD patients with E326K or T369M variants.

Human brain anatomy reflects separable genetic and environmental components of socioeconomic status

Socioeconomic status (SES) correlates with brain structure, a relation of interest given the long-observed relations of SES to cognitive abilities and health. Yet, major questions remain open, in particular, the pattern of causality that underlies this relation. In an unprecedently large study, here, we assess genetic and environmental contributions to SES differences in neuroanatomy. We first establish robust SES-gray matter relations across a number of brain regions, cortical and subcortical. These regional correlates are parsed into predominantly genetic factors and those potentially due to the environment. We show that genetic effects are stronger in some areas (prefrontal cortex, insula) than others. In areas showing less genetic effect (cerebellum, lateral temporal), environmental factors are likely to be influential. Our results imply a complex interplay of genetic and environmental factors that influence the SES-brain relation and may eventually provide insights relevant to policy.

Large-scale GWAS of food liking reveals genetic determinants and genetic correlations with distinct neurophysiological traits

We present the results of a GWAS of food liking conducted on 161,625 participants from the UK-Biobank. Liking was assessed over 139 specific foods using a 9-point scale. Genetic correlations coupled with structural equation modelling identified a multi-level hierarchical map of food-liking with three main dimensions: "Highly-palatable", "Acquired" and "Low-caloric". The Highly-palatable dimension is genetically uncorrelated from the other two, suggesting that independent processes underlie liking high reward foods. This is confirmed by genetic correlations with MRI brain traits which show with distinct associations. Comparison with the corresponding food consumption traits shows a high genetic correlation, while liking exhibits twice the heritability. GWAS analysis identified 1,401 significant food-liking associations which showed substantial agreement in the direction of effects with 11 independent cohorts. In conclusion, we created a comprehensive map of the genetic determinants and associated neurophysiological factors of food-liking.

Innovative in Silico Approaches for Characterization of Genes and Proteins

Bioinformatics is an amalgamation of biology, mathematics and computer science. It is a science which gathers the information from biology in terms of molecules and applies the informatic techniques to the gathered information for understanding and organizing the data in a useful manner. With the help of bioinformatics, the experimental data generated is stored in several databases available online like nucleotide database, protein databases, GENBANK and others. The data stored in these databases is used as reference for experimental evaluation and validation. Till now several online tools have been developed to analyze the genomic, transcriptomic, proteomics, epigenomics and metabolomics data. Some of them include Human Splicing Finder (HSF), Exonic Splicing Enhancer Mutation taster, and others. A number of SNPs are observed in the non-coding, intronic regions and play a role in the regulation of genes, which may or may not directly impose an effect on the protein expression. Many mutations are thought to influence the splicing mechanism by affecting the existing splice sites or creating a new sites. To predict the effect of mutation (SNP) on splicing mechanism/signal, HSF was developed. Thus, the tool is helpful in predicting the effect of mutations on splicing signals and can provide data even for better understanding of the intronic mutations that can be further validated experimentally. Additionally, rapid advancement in proteomics have steered researchers to organize the study of protein structure, function, relationships, and dynamics in space and time. Thus the effective integration of all of these technological interventions will eventually lead to steering up of next-generation systems biology, which will provide valuable biological insights in the field of research, diagnostic, therapeutic and development of personalized medicine.

SVA retrotransposons and a low copy repeat in humans and great apes: a mobile connection

Segmental duplications (SDs) constitute a considerable fraction of primate genomes. They contribute to genetic variation and provide raw material for evolution. Groups of SDs are characterized by the presence of shared core duplicons. One of these core duplicons, low copy repeat (lcr)16a, has been shown to be particularly active in the propagation of interspersed SDs in primates. The underlying mechanisms are, however, only partially understood. Alu short interspersed elements (SINEs) are frequently found at breakpoints and have been implicated in the expansion of SDs.

Detailed analysis of lcr16a-containing SDs shows that the hominid-specific SVA (SINE-R-VNTR-Alu) retrotransposon is an integral component of the core duplicon in Asian and African great apes. In orang-utan, it provides breakpoints and contributes to both interchromosomal and intrachromosomal lcr16a mobility by inter-element recombination. Furthermore, the data suggest that in hominines (human, chimpanzee, gorilla) SVA recombination-mediated integration of a circular intermediate is the founding event of a lineage-specific lcr16a expansion. One of the hominine lcr16a copies displays large flanking direct repeats, a structural feature shared by other SDs in the human genome.

Taken together, the results obtained extend the range of SVAs’ contribution to genome evolution from RNA-mediated transduction to DNA-based recombination. In addition, they provide further support for a role of circular intermediates in SD mobilization.

Genetic analyses identify pleiotropy and causality for blood proteins and highlight Wnt/β-catenin signalling in migraine

Migraine is a common complex disorder with a significant polygenic SNP heritability (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${h}_{{SNP}}^{2}$$\end{document}hSNP2). Here we utilise genome-wide association study (GWAS) summary statistics to study pleiotropy between blood proteins and migraine under the polygenic model. We estimate \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${h}_{{SNP}}^{2}$$\end{document}hSNP2 for 4625 blood protein GWASs and identify 325 unique proteins with a significant \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${h}_{{SNP}}^{2}$$\end{document}hSNP2 for use in subsequent genetic analyses. Pleiotropy analyses link 58 blood proteins to migraine risk at genome-wide, gene and/or single-nucleotide polymorphism levels—suggesting shared genetic influences or causal relationships. Notably, the identified proteins are largely distinct from migraine GWAS loci. We show that higher levels of DKK1 and PDGFB, and lower levels of FARS2, GSTA4 and CHIC2 proteins have a significant causal effect on migraine. The risk-increasing effect of DKK1 is particularly interesting—indicating a role for downregulation of β-catenin-dependent Wnt signalling in migraine risk, suggesting Wnt activators that restore Wnt/β-catenin signalling in brain could represent therapeutic tools against migraine.

Understanding of the causes and treatment of migraine is incomplete. Here, the authors detect pleiotropic genetic effects and causal relationships between migraine and 58 proteins that are largely distinct from migraine-associated loci identified by GWAS.

An eigenvalue ratio approach to inferring population structure from whole genome sequencing data

Inference of population structure from genetic data plays an important role in population and medical genetics studies. With the advancement and decreasing cost of sequencing technology, the increasingly available whole genome sequencing data provide much richer information about the underlying population structure. The traditional method (Patterson et al., 2006) originally developed for array-based genotype data for computing and selecting top principal components that capture population structure may not perform well on sequencing data for two reasons. First, the number of genetic variants p is much larger than the sample size n in sequencing data such that the sample-to-marker ratio n/p is nearly zero, violating the assumption of the Tracy-Widom test used in their method. Second, their method might not be able to handle the linkage disequilibrium well in sequencing data. To resolve those two practical issues, we propose a new method called ERStruct to determine the number of top informative principal components based on sequencing data. More specifically, we propose to use the ratio of consecutive eigenvalues as a more robust test statistic, and then we approximate its null distribution using modern random matrix theory. Both simulation studies and applications to two public data sets from the HapMap 3 and the 1000 Genomes Projects demonstrate the empirical performance of our ERStruct method. This article is protected by copyright. All rights reserved.

The sexual brain, genes, and cognition: A machine-predicted brain sex score explains individual differences in cognitive intelligence and genetic influence in young children

Sex impacts the development of the brain and cognition differently across individuals. However, the literature on brain sex dimorphism in humans is mixed. We aim to investigate the biological underpinnings of the individual variability of sexual dimorphism in the brain and its impact on cognitive performance. To this end, we tested whether the individual difference in brain sex would be linked to that in cognitive performance that is influenced by genetic factors in prepubertal children (N = 9,658, ages 9-10 years old; the Adolescent Brain Cognitive Development study). To capture the interindividual variability of the brain, we estimated the probability of being male or female based on the brain morphometry and connectivity features using machine learning (herein called a brain sex score). The models accurately classified the biological sex with a test ROC-AUC of 93.32%. As a result, a greater brain sex score correlated significantly with greater intelligence (p_fdr  < .001, η p 2 $$ {\eta}_p^2 $$  = .011-.034; adjusted for covariates) and higher cognitive genome-wide polygenic scores (GPSs) (p_fdr  < .001, η p 2 $$ {\eta}_p^2 $$  < .005). Structural equation models revealed that the GPS-intelligence association was significantly modulated by the brain sex score, such that a brain with a higher maleness score (or a lower femaleness score) mediated a positive GPS effect on intelligence (indirect effects = .006-.009; p = .002-.022; sex-stratified analysis). The finding of the sex modulatory effect on the gene-brain-cognition relationship presents a likely biological pathway to the individual and sex differences in the brain and cognitive performance in preadolescence.

Computationally repurposing drugs for breast cancer subtypes using a network-based approach

‘De novo’ drug discovery is costly, slow, and with high risk. Repurposing known drugs for treatment of other diseases offers a fast, low-cost/risk and highly-efficient method toward development of efficacious treatments. The emergence of large-scale heterogeneous biomolecular networks, molecular, chemical and bioactivity data, and genomic and phenotypic data of pharmacological compounds is enabling the development of new area of drug repurposing called ‘in silico’ drug repurposing, i.e., computational drug repurposing (CDR). The aim of CDR is to discover new indications for an existing drug (drug-centric) or to identify effective drugs for a disease (disease-centric). Both drug-centric and disease-centric approaches have the common challenge of either assessing the similarity or connections between drugs and diseases. However, traditional CDR is fraught with many challenges due to the underlying complex pharmacology and biology of diseases, genes, and drugs, as well as the complexity of their associations. As such, capturing highly non-linear associations among drugs, genes, diseases by most existing CDR methods has been challenging. We propose a network-based integration approach that can best capture knowledge (and complex relationships) contained within and between drugs, genes and disease data. A network-based machine learning approach is applied thereafter by using the extracted knowledge and relationships in order to identify single and pair of approved or experimental drugs with potential therapeutic effects on different breast cancer subtypes. Indeed, further clinical analysis is needed to confirm the therapeutic effects of identified drugs on each breast cancer subtype.

Common variants contribute to intrinsic human brain functional networks

The human brain forms functional networks of correlated activity, which have been linked with both cognitive and clinical outcomes. However, the genetic variants affecting brain function are largely unknown. Here, we used resting-state functional magnetic resonance images from 47,276 individuals to discover and validate common genetic variants influencing intrinsic brain activity. We identified 45 new genetic regions associated with brain functional signatures (P < 2.8 × 10^-11), including associations to the central executive, default mode, and salience networks involved in the triple-network model of psychopathology. A number of brain activity-associated loci colocalized with brain disorders (e.g., the APOE ε4 locus with Alzheimer's disease). Variation in brain function was genetically correlated with brain disorders, such as major depressive disorder and schizophrenia. Together, our study provides a step forward in understanding the genetic architecture of brain functional networks and their genetic links to brain-related complex traits and disorders.

Genetic variants associated with longitudinal changes in brain structure across the lifespan

Human brain structure changes throughout the lifespan. Altered brain growth or rates of decline are implicated in a vast range of psychiatric, developmental and neurodegenerative diseases. In this study, we identified common genetic variants that affect rates of brain growth or atrophy in what is, to our knowledge, the first genome-wide association meta-analysis of changes in brain morphology across the lifespan. Longitudinal magnetic resonance imaging data from 15,640 individuals were used to compute rates of change for 15 brain structures. The most robustly identified genes GPR139, DACH1 and APOE are associated with metabolic processes. We demonstrate global genetic overlap with depression, schizophrenia, cognitive functioning, insomnia, height, body mass index and smoking. Gene set findings implicate both early brain development and neurodegenerative processes in the rates of brain changes. Identifying variants involved in structural brain changes may help to determine biological pathways underlying optimal and dysfunctional brain development and aging.

Early Tumor-Immune Microenvironmental Remodeling and Response to First-Line Fluoropyrimidine and Platinum Chemotherapy in Advanced Gastric Cancer

Chemotherapy is ubiquitous in first-line treatment of advanced gastric cancer, yet responses are heterogeneous, and little is known about mediators of chemotherapy response. To move forward, an understanding of the effects of standard chemotherapy on the tumor-immune microenvironment (TME) is needed. Coupling whole-exome sequencing, bulk RNA and single-cell transcriptomics from paired pretreatment and on-treatment samples in treatment-naïve patients with HER2-positive and HER2-negative gastric cancer, we define features associated with response to platinum-based chemotherapy. Response was associated with on-treatment TME remodeling including natural killer (NK) cell recruitment, decreased tumor-associated macrophages, M1-macrophage repolarization, and increased effector T-cell infiltration. Among chemotherapy nonresponders, we observed low/absent PD-L1 expression or modulation, on-treatment increases in Wnt signaling, B-cell infiltration, and LAG3-expressing T cells coupled to an exodus of dendritic cells. We did not observe significant genomic changes in early on-treatment sampling. We provide a map of on-treatment TME modulation with standard chemotherapy and nominate candidate future approaches.

SIGNIFICANCE

Using paired pretreatment and on-treatment samples during standard first-line chemotherapy, we identify chemotherapy-induced NK-cell infiltration, macrophage repolarization, and increased antigen presentation among responders. Increased LAG3 expression and decreased dendritic cell abundance were seen in nonresponders, emphasizing remodeling of the TME during chemotherapy response and resistance. This article is highlighted in the In This Issue feature, p. 873.

Genome-wide identification of the shared genetic basis of cannabis and cigarette smoking and schizophrenia implicates NCAM1 and neuronal abnormality

OBJECTIVES

Confirming the existence and composition of the shared genetic basis of Schizophrenia and cannabis and cigarette smoking has critical values for the clinical prevention and intervention of psychosis.

METHODS

To achieve this goal, we leveraged Genome-Wide summary statistics of Schizophrenia (n = 99,934), cigarette smoking (n = 518,633) and cannabis usage (n = 162,082). We applied Causal Analysis Using Summary Effect Estimates (CAUSE) and genomic structural equation modeling (GenomicSEM) to quantify the contribution of a common genetic factor of cannabis and cigarette smoking and schizophrenia (referred to as SCZ_SMO), then identified genome-wide loci that made up SCZ_SMO.

RESULTS

We estimated that SCZ_SMO explained 8.6% of Schizophrenia heritability (Z score <-2.5 in CAUSE, p<10^-20 in Genomic SEM). There were 20 independent loci showing association with SCZ_SMO at the genome-wide threshold of p<5 × 10^-8. At the top locus on chromosome 11, fine-mapping identified rs7945073 (posterior inclusion probability =0.12, p = 2.24 × 10^-32) as the top risk variants. Gene-level association and fine-mapping highlighted NCAM1, PHC2, and SEMA6D as risk genes of SCZ_SMO. Other risk genes were enriched in cortex, neuron, and dendritic spines (adjusted p<0.05). SCZ_SMO showed significant positive correlation (p<10^-6) with the genetic risk of attention deficit hyperactivity disorder (r = 0.50), lifestyle problems (r = 0.83), social deprivation (r = 0.58) and all-cause pregnant loss (r = 0.60).

CONCLUSION

Our result provided new evidence on the shared genetic basis model for the association between Schizophrenia and smoking and provided genetic and biological insights into their shared mechanism.

Epigenome-Wide Analysis Reveals DNA Methylation Alteration in ZFP57 and Its Target RASGFR2 in a Mexican Population Cohort with Autism

Autism Spectrum Disorders (ASD) comprise a group of heterogeneous and complex neurodevelopmental disorders. Genetic and environmental factors contribute to ASD etiology. DNA methylation is particularly relevant for ASD due to its mediating role in the complex interaction between genotype and environment and has been implicated in ASD pathophysiology. The lack of diversity in DNA methylation studies in ASD individuals is remarkable. Since genetic and environmental factors are likely to vary across populations, the study of underrepresented populations is necessary to understand the molecular alterations involved in ASD and the risk factors underlying these changes. This study explored genome-wide differences in DNA methylation patterns in buccal epithelium cells between Mexican ASD patients (n = 27) and age-matched typically developing (TD: n = 15) children. DNA methylation profiles were evaluated with the Illumina 450k array. We evaluated the interaction between sex and ASD and found a differentially methylated region (DMR) over the 5′UTR region of ZFP57 and one of its targets, RASGRF2. These results match previous findings in brain tissue, which may indicate that ZFP57 could be used as a proxy for DNA methylation in different tissues. This is the first study performed in a Mexican, and subsequently, Latin American, population that evaluates DNA methylation in ASD patients.

The relationships between women's reproductive factors: a Mendelian randomisation analysis

BACKGROUND

Women's reproductive factors include their age at menarche and menopause, the age at which they start and stop having children and the number of children they have. Studies that have linked these factors with disease risk have largely investigated individual reproductive factors and have not considered the genetic correlation and total interplay that may occur between them. This study aimed to investigate the nature of the relationships between eight female reproductive factors.

METHODS

We used data from the UK Biobank and genetic consortia with data available for the following reproductive factors: age at menarche, age at menopause, age at first birth, age at last birth, number of births, being parous, age first had sexual intercourse and lifetime number of sexual partners. Linkage disequilibrium score regression (LDSC) was performed to investigate the genetic correlation between reproductive factors. We then applied Mendelian randomisation (MR) methods to estimate the causal relationships between these factors. Sensitivity analyses were used to investigate directionality of the effects, test for evidence of pleiotropy and account for sample overlap.

RESULTS

LDSC indicated that most reproductive factors are genetically correlated (rg range: |0.06-0.94|), though there was little evidence for genetic correlations between lifetime number of sexual partners and age at last birth, number of births and ever being parous (rg < 0.01). MR revealed potential causal relationships between many reproductive factors, including later age at menarche (1 SD increase) leading to a later age at first sexual intercourse (beta (B) = 0.09 SD, 95% confidence intervals (CI) = 0.06,0.11), age at first birth (B = 0.07 SD, CI = 0.04,0.10), age at last birth (B = 0.06 SD, CI = 0.04,0.09) and age at menopause (B = 0.06 SD, CI = 0.03,0.10). Later age at first birth was found to lead to a later age at menopause (B = 0.21 SD, CI = 0.13,0.29), age at last birth (B = 0.72 SD, CI = 0.67, 0.77) and a lower number of births (B = -0.38 SD, CI = -0.44, -0.32).

CONCLUSION

This study presents evidence that women's reproductive factors are genetically correlated and causally related. Future studies examining the health sequelae of reproductive factors should consider a woman's entire reproductive history, including the causal interplay between reproductive factors.

Genetic mosaicism in the human brain: from lineage tracing to neuropsychiatric disorders

Genetic mosaicism is the result of the accumulation of somatic mutations in the human genome starting from the first postzygotic cell generation and continuing throughout the whole life of an individual. The rapid development of next-generation and single-cell sequencing technologies is now allowing the study of genetic mosaicism in normal tissues, revealing unprecedented insights into their clonal architecture and physiology. The somatic variant repertoire of an adult human neuron is the result of somatic mutations that accumulate in the brain by different mechanisms and at different rates during development and ageing. Non-pathogenic developmental mutations function as natural barcodes that once identified in deep bulk or single-cell sequencing can be used to retrospectively reconstruct human lineages. This approach has revealed novel insights into the clonal structure of the human brain, which is a mosaic of clones traceable to the early embryo that contribute differentially to the brain and distinct areas of the cortex. Some of the mutations happening during development, however, have a pathogenic effect and can contribute to some epileptic malformations of cortical development and autism spectrum disorder. In this Review, we discuss recent findings in the context of genetic mosaicism and their implications for brain development and disease.

Towards the detection of copy number variation from single sperm sequencing in cattle

BACKGROUND

Copy number variation (CNV) has been routinely studied using bulk-cell sequencing. However, CNV is not well studied on the single-cell level except for humans and a few model organisms.

RESULTS

We sequenced 143 single sperms of two Holstein bulls, from which we predicted CNV events using 14 single sperms with deep sequencing. We then compared the CNV results derived from single sperms with the bulk-cell sequencing of one bull's family trio of diploid genomes. As a known CNV hotspot, segmental duplications were also predicted using the bovine ARS-UCD1.2 genome. Although the trio CNVs validated only some single sperm CNVs, they still showed a distal chromosomal distribution pattern and significant associations with segmental duplications and satellite repeats.

CONCLUSION

Our preliminary results pointed out future research directions and highlighted the importance of uniform whole genome amplification, deep sequence coverage, and dedicated software pipelines for CNV detection using single cell sequencing data.

Positive selection acts on regulatory genetic variants in populations of European ancestry that affect ALDH2 gene expression

ALDH2 is a key enzyme in alcohol metabolism that protects cells from acetaldehyde toxicity. Using iHS, iSAFE and FST statistics, we identified regulatory acting variants affecting ALDH2 gene expression under positive selection in populations of European ancestry. Several SNPs (rs3184504, rs4766578, rs10774625, rs597808, rs653178, rs847892, rs2013002) that function as eQTLs for ALDH2 in various tissues showed evidence of strong positive selection. Very large pairwise FST values indicated high genetic differentiation at these loci between populations of European ancestry and populations of other global ancestries. Estimating the timing of positive selection on the beneficial alleles suggests that these variants were recently adapted approximately 3000-3700 years ago. The derived beneficial alleles are in complete linkage disequilibrium with the derived ALDH2 promoter variant rs886205, which is associated with higher transcriptional activity. The SNPs rs4766578 and rs847892 are located in binding sequences for the transcription factor HNF4A, which is an important regulatory element of ALDH2 gene expression. In contrast to the missense variant ALDH2 rs671 (ALDH2*2), which is common only in East Asian populations and is associated with greatly reduced enzyme activity and alcohol intolerance, the beneficial alleles of the regulatory variants identified in this study are associated with increased expression of ALDH2. This suggests adaptation of Europeans to higher alcohol consumption.

The role of single-nucleotide polymorphisms of the 8q24 chromosome region in patients with concomitant bladder and prostate cancer

OBJECTIVE

To assess whether single-nucleotide polymorphisms (SNPs) of the 8q24 chromosome region are associated with recurrence-free survival (RFS) after radical cystoprostatectomy (RC) in patients with concomitant bladder (BC) and prostate cancer (PC).

MATERIALS AND METHODS

A cohort of thirty-six patients treated with RC and pelvic lymph node dissection and histologically exhibited invasive BC and incidental PC. Using Sanger sequencing, a total of seven SNPs in the androgen-responsive element of the promoter region of the following genes were assessed in tumor-free lymph nodes and correlated with oncological outcomes: PSCA (rs2294008, rs2978974, rs1045531, rs3736001), MYC (rs6983267), FXBO32 (rs7830622), and MIR151A (rs14974929). The median follow-up was 26 months (range: 4-68).

RESULTS

In a dominant model, patients exhibiting rs2978974 as a minor allelic variant of the PSCA gene had worse RFS (32 vs. 75%, p = 0.015). No associations were found for the other SNPs.

CONCLUSIONS

These data suggest that the rs2978974 of the PSCA gene correlates with inferior BC-specific RFS after RC and should be further evaluated in larger studies.

Fine-Scale Genetic Structure in the United Arab Emirates Reflects Endogamous and Consanguineous Culture, Population History, and Geography

The indigenous population of the United Arab Emirates (UAE) has a unique demographic and cultural history. Its tradition of endogamy and consanguinity is expected to produce genetic homogeneity and partitioning of gene pools while population movements and intercontinental trade are likely to have contributed to genetic diversity. Emiratis and neighboring populations of the Middle East have been underrepresented in the population genetics literature with few studies covering the broader genetic history of the Arabian Peninsula. Here, we genotyped 1,198 individuals from the seven Emirates using 1.7 million markers and by employing haplotype-based algorithms and admixture analyses, we reveal the fine-scale genetic structure of the Emirati population. Shared ancestry and gene flow with neighboring populations display their unique geographic position while increased intra- versus inter-Emirati kinship and sharing of uniparental haplogroups, reflect the endogamous and consanguineous cultural traditions of the Emirates and their tribes.

Genetic Risk for Smoking: Disentangling Interplay Between Genes and Socioeconomic Status

This study aims to disentangle the contribution of genetic liability, educational attainment (EA), and their overlap and interaction in lifetime smoking. We conducted genome-wide association studies (GWASs) in UK Biobank (N = 394,718) to (i) capture variants for lifetime smoking, (ii) variants for EA, and (iii) variants that contribute to lifetime smoking independently from EA ('smoking-without-EA'). Based on the GWASs, three polygenic scores (PGSs) were created for individuals from the Netherlands Twin Register (NTR, N = 17,805) and the Netherlands Mental Health Survey and Incidence Study-2 (NEMESIS-2, N = 3090). We tested gene-environment (G × E) interactions between each PGS, neighborhood socioeconomic status (SES) and EA on lifetime smoking. To assess if the PGS effects were specific to smoking or had broader implications, we repeated the analyses with measures of mental health. After subtracting EA effects from the smoking GWAS, the SNP-based heritability decreased from 9.2 to 7.2%. The genetic correlation between smoking and SES characteristics was reduced, whereas overlap with smoking traits was less affected by subtracting EA. The PGSs for smoking, EA, and smoking-without-EA all predicted smoking. For mental health, only the PGS for EA was a reliable predictor. There were suggestions for G × E for some relationships, but there were no clear patterns per PGS type. This study showed that the genetic architecture of smoking has an EA component in addition to other, possibly more direct components. PGSs based on EA and smoking-without-EA had distinct predictive profiles. This study shows how disentangling different models of genetic liability and interplay can contribute to our understanding of the etiology of smoking.

Interaction Testing and Polygenic Risk Scoring to Estimate the Association of Common Genetic Variants With Treatment Resistance in Schizophrenia

IMPORTANCE

About 20% to 30% of people with schizophrenia have psychotic symptoms that do not respond adequately to first-line antipsychotic treatment. This clinical presentation, chronic and highly disabling, is known as treatment-resistant schizophrenia (TRS). The causes of treatment resistance and their relationships with causes underlying schizophrenia are largely unknown. Adequately powered genetic studies of TRS are scarce because of the difficulty in collecting data from well-characterized TRS cohorts.

OBJECTIVE

To examine the genetic architecture of TRS through the reassessment of genetic data from schizophrenia studies and its validation in carefully ascertained clinical samples.

DESIGN, SETTING, AND PARTICIPANTS

Two case-control genome-wide association studies (GWASs) of schizophrenia were performed in which the case samples were defined as individuals with TRS (n = 10 501) and individuals with non-TRS (n = 20 325). The differences in effect sizes for allelic associations were then determined between both studies, the reasoning being such differences reflect treatment resistance instead of schizophrenia. Genotype data were retrieved from the CLOZUK and Psychiatric Genomics Consortium (PGC) schizophrenia studies. The output was validated using polygenic risk score (PRS) profiling of 2 independent schizophrenia cohorts with TRS and non-TRS: a prevalence sample with 817 individuals (Cardiff Cognition in Schizophrenia [CardiffCOGS]) and an incidence sample with 563 individuals (Genetics Workstream of the Schizophrenia Treatment Resistance and Therapeutic Advances [STRATA-G]).

MAIN OUTCOMES AND MEASURES

GWAS of treatment resistance in schizophrenia. The results of the GWAS were compared with complex polygenic traits through a genetic correlation approach and were used for PRS analysis on the independent validation cohorts using the same TRS definition.

RESULTS

The study included a total of 85 490 participants (48 635 [56.9%] male) in its GWAS stage and 1380 participants (859 [62.2%] male) in its PRS validation stage. Treatment resistance in schizophrenia emerged as a polygenic trait with detectable heritability (1% to 4%), and several traits related to intelligence and cognition were found to be genetically correlated with it (genetic correlation, 0.41-0.69). PRS analysis in the CardiffCOGS prevalence sample showed a positive association between TRS and a history of taking clozapine (r2 = 2.03%; P = .001), which was replicated in the STRATA-G incidence sample (r2 = 1.09%; P = .04).

CONCLUSIONS AND RELEVANCE

In this GWAS, common genetic variants were differentially associated with TRS, and these associations may have been obscured through the amalgamation of large GWAS samples in previous studies of broadly defined schizophrenia. Findings of this study suggest the validity of meta-analytic approaches for studies on patient outcomes, including treatment resistance.

Assessing the contribution of rare variants to complex trait heritability from whole-genome sequence data

Analyses of data from genome-wide association studies on unrelated individuals have shown that, for human traits and diseases, approximately one-third to two-thirds of heritability is captured by common SNPs. However, it is not known whether the remaining heritability is due to the imperfect tagging of causal variants by common SNPs, in particular whether the causal variants are rare, or whether it is overestimated due to bias in inference from pedigree data. Here we estimated heritability for height and body mass index (BMI) from whole-genome sequence data on 25,465 unrelated individuals of European ancestry. The estimated heritability was 0.68 (standard error 0.10) for height and 0.30 (standard error 0.10) for body mass index. Low minor allele frequency variants in low linkage disequilibrium (LD) with neighboring variants were enriched for heritability, to a greater extent for protein-altering variants, consistent with negative selection. Our results imply that rare variants, in particular those in regions of low linkage disequilibrium, are a major source of the still missing heritability of complex traits and disease.

Association with HLA-DRβ1 position 37 distinguishes juvenile Dermatomyositis from adult-onset myositis

Objectives

Juvenile dermatomyositis (JDM) is a rare, severe autoimmune disease and the most common idiopathic inflammatory myopathy (IIM) of children. JDM and adult-onset dermatomyositis (DM) have similar clinical, biological and serological features, although these features differ in prevalence between childhood-onset and adult-onset disease, suggesting age of disease onset may influence pathogenesis. Therefore, a JDM-focused genetic analysis was performed using the largest collection of JDM samples to date.

Methods

Caucasian JDM samples (n = 952) obtained via international collaboration were genotyped using the Illumina HumanCoreExome chip. Additional non-assayed HLA loci and genome-wide SNPs were imputed.

Results

HLA-DRB1*03:01 was confirmed as the classical HLA allele most strongly associated with JDM (OR 1.66; 95% CI 1.46, 1.89; P = 1.4 × 10−14), with an independent association at HLA-C*02:02 (OR = 1.74; 95% CI 1.42, 2.13, P = 7.13 × 10−8). Analyses of amino acid positions within HLA-DRB1 indicated the strongest association was at position 37 (omnibus P = 3.3 × 10−19), with suggestive evidence this association was independent of position 74 (omnibus P = 5.1 × 10−5), the position most strongly associated with adult-onset DM. Conditional analyses also suggested the association at position 37 of HLA-DRB1 was independent of some alleles of the Caucasian HLA 8.1 ancestral haplotype (AH8.1) such as HLA-DQB1*02:01 (OR = 1.62; 95% CI 1.36, 1.93; P = 8.70 × 10−8), but not HLA-DRB1*03:01 (OR = 1.49; 95% CR 1.24, 1.80; P = 2.24 × 10−5). No associations outside the HLA region were identified.

Conclusions

Our findings confirm previous associations with AH8.1 and HLA-DRB1*03:01, HLA-C*02:02 and identify a novel association with amino acid position 37 within HLA-DRB1 which may distinguish JDM from adult DM.

Investigating the effect of sexual behaviour on oropharyngeal cancer risk: a methodological assessment of Mendelian randomization

BACKGROUND

Human papilloma virus infection is known to influence oropharyngeal cancer (OPC) risk, likely via sexual transmission. However, sexual behaviour has been correlated with other risk factors including smoking and alcohol, meaning independent effects are difficult to establish. We aimed to evaluate the causal effect of sexual behaviour on the risk of OPC using Mendelian randomization (MR).

METHODS

Genetic variants robustly associated with age at first sex (AFS) and the number of sexual partners (NSP) were used to perform both univariable and multivariable MR analyses with summary data on 2641 OPC cases and 6585 controls, obtained from the largest available genome-wide association studies (GWAS). Given the potential for genetic pleiotropy, we performed a number of sensitivity analyses: (i) MR methods to account for horizontal pleiotropy, (ii) MR of sexual behaviours on positive (cervical cancer and seropositivity for Chlamydia trachomatis) and negative control outcomes (lung and oral cancer), (iii) Causal Analysis Using Summary Effect estimates (CAUSE), to account for correlated and uncorrelated horizontal pleiotropic effects, (iv) multivariable MR analysis to account for the effects of smoking, alcohol, risk tolerance and educational attainment.

RESULTS

In univariable MR, we found evidence supportive of an effect of both later AFS (IVW OR = 0.4, 95%CI (0.3, 0.7), per standard deviation (SD), p = < 0.001) and increasing NSP (IVW OR = 2.2, 95%CI (1.3, 3.8) per SD, p = < 0.001) on OPC risk. These effects were largely robust to sensitivity analyses accounting for horizontal pleiotropy. However, negative control analysis suggested potential violation of the core MR assumptions and subsequent CAUSE analysis implicated pleiotropy of the genetic instruments used to proxy sexual behaviours. Finally, there was some attenuation of the univariable MR results in the multivariable models (AFS IVW OR = 0.7, 95%CI (0.4, 1.2), p = 0.21; NSP IVW OR = 0.9, 95%CI (0.5 1.7), p = 0.76).

CONCLUSIONS

Despite using genetic variants strongly related sexual behaviour traits in large-scale GWAS, we found evidence for correlated pleiotropy. This emphasizes a need for multivariable approaches and the triangulation of evidence when performing MR of complex behavioural traits.

Comprehensive Statistical and Bioinformatics Analysis in the Deciphering of Putative Mechanisms by Which Lipid-Associated GWAS Loci Contribute to Coronary Artery Disease

The study was designed to evaluate putative mechanisms by which lipid-associated loci identified by genome-wide association studies (GWAS) are involved in the molecular pathogenesis of coronary artery disease (CAD) using a comprehensive statistical and bioinformatics analysis. A total of 1700 unrelated individuals of Slavic origin from the Central Russia, including 991 CAD patients and 709 healthy controls were examined. Sixteen lipid-associated GWAS loci were selected from European studies and genotyped using the MassArray-4 system. The polymorphisms were associated with plasma lipids such as total cholesterol (rs12328675, rs4846914, rs55730499, and rs838880), LDL-cholesterol (rs3764261, rs55730499, rs1689800, and rs838880), HDL-cholesterol (rs3764261) as well as carotid intima-media thickness/CIMT (rs12328675, rs11220463, and rs1689800). Polymorphisms such as rs4420638 of APOC1 (p = 0.009), rs55730499 of LPA (p = 0.0007), rs3136441 of F2 (p < 0.0001), and rs6065906 of PLTP (p = 0.002) showed significant associations with the risk of CAD, regardless of sex, age, and body mass index. A majority of the observed associations were successfully replicated in large independent cohorts. Bioinformatics analysis allowed establishing (1) phenotype-specific and shared epistatic gene–gene and gene–smoking interactions contributing to all studied cardiovascular phenotypes; (2) lipid-associated GWAS loci might be allele-specific binding sites for transcription factors from gene regulatory networks controlling multifaceted molecular mechanisms of atherosclerosis.

Regenerating zebrafish scales express a subset of evolutionary conserved genes involved in human skeletal disease

Background

Scales are mineralised exoskeletal structures that are part of the dermal skeleton. Scales have been mostly lost during evolution of terrestrial vertebrates whilst bony fish have retained a mineralised dermal skeleton in the form of fin rays and scales. Each scale is a mineralised collagen plate that is decorated with both matrix-building and resorbing cells. When removed, an ontogenetic scale is quickly replaced following differentiation of the scale pocket-lining cells that regenerate a scale. Processes promoting de novo matrix formation and mineralisation initiated during scale regeneration are poorly understood. Therefore, we performed transcriptomic analysis to determine gene networks and their pathways involved in dermal scale regeneration.

Results

We defined the transcriptomic profiles of ontogenetic and regenerating scales of zebrafish and identified 604 differentially expressed genes (DEGs). These were enriched for extracellular matrix, ossification, and cell adhesion pathways, but not in enamel or dentin formation processes indicating that scales are reminiscent to bone. Hypergeometric tests involving monogenetic skeletal disorders showed that DEGs were strongly enriched for human orthologues that are mutated in low bone mass and abnormal bone mineralisation diseases (P< 2× 10⁻³). The DEGs were also enriched for human orthologues associated with polygenetic skeletal traits, including height (P< 6× 10⁻⁴), and estimated bone mineral density (eBMD, P< 2× 10⁻⁵). Zebrafish mutants of two human orthologues that were robustly associated with height (COL11A2, P=6× 10⁻²⁴) or eBMD (SPP1, P=6× 10⁻²⁰) showed both exo- and endo- skeletal abnormalities as predicted by our genetic association analyses; col11a2^Y228X/Y228X mutants showed exoskeletal and endoskeletal features consistent with abnormal growth, whereas spp1^P160X/P160X mutants predominantly showed mineralisation defects.

Conclusion

We show that scales have a strong osteogenic expression profile comparable to other elements of the dermal skeleton, enriched in genes that favour collagen matrix growth. Despite the many differences between scale and endoskeletal developmental processes, we also show that zebrafish scales express an evolutionarily conserved sub-population of genes that are relevant to human skeletal disease.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01209-8.

Detecting selection using extended haplotype homozygosity (EHH)-based statistics in unphased or unpolarized data

Analysis of population genetic data often includes a search for genomic regions with signs of recent positive selection. One of such approaches involves the concept of extended haplotype homozygosity (EHH) and its associated statistics. These statistics typically require phased haplotypes, and some of them necessitate polarized variants. Here, we unify and extend previously proposed modifications to loosen these requirements. We compare the modified versions with the original ones by measuring the false discovery rate in simulated whole-genome scans and by quantifying the overlap of inferred candidate regions in empirical data. We find that phasing information is indispensable for accurate estimation of within-population statistics (for all but very large samples) and of cross-population statistics for small samples. Ancestry information, in contrast, is of lesser importance for both types of statistic. Our publicly available R package rehh incorporates the modified statistics presented here.

Chromosomal microarray analysis of 410 Han Chinese patients with autism spectrum disorder or unexplained intellectual disability and developmental delay

Copy number variants (CNVs) are recognized as a crucial genetic cause of neurodevelopmental disorders (NDDs). Chromosomal microarray analysis (CMA), the first-tier diagnostic test for individuals with NDDs, has been utilized to detect CNVs in clinical practice, but most reports are still from populations of European ancestry. To contribute more worldwide clinical genomics data, we investigated the genetic etiology of 410 Han Chinese patients with NDDs (151 with autism and 259 with unexplained intellectual disability (ID) and developmental delay (DD)) using CMA (Affymetrix) after G-banding karyotyping. Among all the NDD patients, 109 (26.6%) carried clinically relevant CNVs or uniparental disomies (UPDs), and 8 (2.0%) had aneuploidies (6 with trisomy 21 syndrome, 1 with 47,XXY, 1 with 47,XYY). In total, we found 129 clinically relevant CNVs and UPDs, including 32 CNVs in 30 ASD patients, and 92 CNVs and 5 UPDs in 79 ID/DD cases. When excluding the eight patients with aneuploidies, the diagnostic yield of pathogenic and likely pathogenic CNVs and UPDs was 20.9% for all NDDs (84/402), 3.3% in ASD (5/151), and 31.5% in ID/DD (79/251). When aneuploidies were included, the diagnostic yield increased to 22.4% for all NDDs (92/410), and 33.6% for ID/DD (87/259). We identified a de novo CNV in 14.9% (60/402) of subjects with NDDs. Interestingly, a higher diagnostic yield was observed in females (31.3%, 40/128) compared to males (16.1%, 44/274) for all NDDs (P = 4.8 × 10^-4), suggesting that a female protective mechanism exists for deleterious CNVs and UPDs.

Ancestry- and sex-specific effects underlying inguinal hernia susceptibility identified in a multiethnic genome-wide association study meta-analysis

Inguinal hernias are some of the most frequently diagnosed conditions in clinical practice and inguinal hernia repair is the most common procedure performed by general surgeons. Studies of inguinal hernias in non-European populations are lacking, though it is expected that such studies could identify novel loci. Further, the cumulative lifetime incidence of inguinal hernia is nine times greater in men than women, however, it is not clear why this difference exists. We conducted a genome-wide association meta-analysis of inguinal hernia risk across 513 120 individuals (35 774 cases and 477 346 controls) of Hispanic/Latino, African, Asian and European descent, with replication in 728 418 participants (33 491 cases and 694 927 controls) from the 23andMe, Inc dataset. We identified 63 genome-wide significant loci (P < 5 × 10⁻⁸), including 41 novel. Ancestry-specific analyses identified two loci (LYPLAL1-AS1/SLC30A10 and STXBP6-NOVA1) in African ancestry individuals. Sex-stratified analyses identified two loci (MYO1D and ZBTB7C) that are specific to women, and four (EBF2, EMX2/RAB11FIP2, VCL and FAM9A/FAM9B) that are specific to men. Functional experiments demonstrated that several of the associated regions (EFEMP1 and LYPLAL1-SLC30A10) function as enhancers and show differential activity between risk and reference alleles. Our study highlights the importance of large-scale genomic studies in ancestrally diverse populations for identifying ancestry-specific inguinal hernia susceptibility loci and provides novel biological insights into inguinal hernia etiology.

PopHumanVar: an interactive application for the functional characterization and prioritization of adaptive genomic variants in humans

Adaptive challenges that humans faced as they expanded across the globe left specific molecular footprints that can be decoded in our today's genomes. Different sets of metrics are used to identify genomic regions that have undergone selection. However, there are fewer methods capable of pinpointing the allele ultimately responsible for this selection. Here, we present PopHumanVar, an interactive online application that is designed to facilitate the exploration and thorough analysis of candidate genomic regions by integrating both functional and population genomics data currently available. PopHumanVar generates useful summary reports of prioritized variants that are putatively causal of recent selective sweeps. It compiles data and graphically represents different layers of information, including natural selection statistics, as well as functional annotations and genealogical estimations of variant age, for biallelic single nucleotide variants (SNVs) of the 1000 Genomes Project phase 3. Specifically, PopHumanVar amasses SNV-based information from GEVA, SnpEFF, GWAS Catalog, ClinVar, RegulomeDB and DisGeNET databases, as well as accurate estimations of iHS, nSL and iSAFE statistics. Notably, PopHumanVar can successfully identify known causal variants of frequently reported candidate selection regions, including EDAR in East-Asians, ACKR1 (DARC) in Africans and LCT/MCM6 in Europeans. PopHumanVar is open and freely available at https://pophumanvar.uab.cat.

Mendelian randomization study reveals a population-specific putative causal effect of type 2 diabetes in risk of cataract

BACKGROUND

The epidemiological association between type 2 diabetes and cataract has been well established. However, it remains unclear whether the two diseases share a genetic basis, and if so, whether this reflects a putative causal relationship.

METHODS

We used East Asian population-based genome-wide association studies (GWAS) summary statistics of type 2 diabetes (Ncase = 36 614, Ncontrol = 155 150) and cataract (Ncase = 24 622, Ncontrol = 187 831) to comprehensively investigate the shared genetics between the two diseases. We performed: (i) linkage disequilibrium score regression (LDSC) and heritability estimation from summary statistics (ρ-HESS) to estimate the genetic correlation and local genetic correlation pattern between type 2 diabetes and cataract; (ii) multiple Mendelian randomization (MR) analyses to infer the putative causality between type 2 diabetes and cataract; and (iii) summary-data-based Mendelian randomization (SMR) to identify candidate risk genes underling the putative causality. Moreover, to investigate the extent of the population-specific genetic effect size underlying the shared genetics between type 2 diabetes and cataract, we applied the same analytical pipeline to perform a comparative analysis on European population-based GWAS of type 2 diabetes (Ncase = 62 892, Ncontrol = 596 424) and cataract (Ncase = 5045, Ncontrol = 356 096).

RESULTS

Using East Asian population-based GWAS summary data, we observed a strong genetic correlation [rg = 0.58, 95% confidence interval (CI) = 0.33, 0.83), P-value = 5.60 × 10-6] between type 2 diabetes and cataract. Both ρ-HESS and multiple MR methods consistently showed a putative causal effect of type 2 diabetes on cataract, with estimated liability-scale MR odds ratios (ORs) at around 1.10 (95% CI = 1.06, 1.17). In contrast, no evidence supports a causal effect of cataract on type 2 diabetes. SMR analysis identified two novel genes MIR4453HG (βSMR = -0.34, 95% CI = -0.46, -0.22, P-value = 6.41 × 10-8) and KCNK17 (βSMR = -0.07, 95% CI = -0.09, -0.05, P-value = 2.49 × 10-10), whose expression levels were likely involved in the putative causality of type 2 diabetes on cataract. On the contrary, our comparative analysis on European population provided universally weak evidence on the genetic correlation and causal relationship between the two diseases.

CONCLUSIONS

Our results provided robust evidence supporting a putative causal effect of type 2 diabetes on the risk of cataract in East Asians, and revealed potential genetic heterogeneity in the shared genetics underlying type 2 diabetes and cataract between East Asians and Europeans. These findings posed new paths on guiding the prevention and early-stage diagnosis of cataract in type 2 diabetes patients.

Limb development genes underlie variation in human fingerprint patterns

Fingerprints are of long-standing practical and cultural interest, but little is known about the mechanisms that underlie their variation. Using genome-wide scans in Han Chinese cohorts, we identified 18 loci associated with fingerprint type across the digits, including a genetic basis for the long-recognized “pattern-block” correlations among the middle three digits. In particular, we identified a variant near EVI1 that alters regulatory activity and established a role for EVI1 in dermatoglyph patterning in mice. Dynamic EVI1 expression during human development supports its role in shaping the limbs and digits, rather than influencing skin patterning directly. Trans-ethnic meta-analysis identified 43 fingerprint-associated loci, with nearby genes being strongly enriched for general limb development pathways. We also found that fingerprint patterns were genetically correlated with hand proportions. Taken together, these findings support the key role of limb development genes in influencing the outcome of fingerprint patterning.

•

GWAS identifies variants associated with fingerprint type across all digits

•

Fingerprint-associated genes are strongly enriched for limb development functions

•

Evi1 alters dermatoglyphs in mice by modulating limb rather than skin development

•

Fingerprint patterns are genetically correlated with hand and finger proportions

Are Alzheimer's and coronary artery diseases genetically related to longevity?

INTRODUCTION

In the last decade researchers have attempted to investigate the shared genetic architecture of longevity and age-related diseases and assess whether the increased longevity in certain people is due to protective alleles in the risk genes for a particular condition or whether there are specific "longevity" genes increasing the lifespan independently of age-related conditions' risk genes. The aim of this study was to investigate the shared genetic component between longevity and two age-related conditions.

METHODS

We performed a cross-trait meta-analysis of publicly available genome-wide data for Alzheimer's disease, coronary artery disease and longevity using a subset-based approach provided by the R package ASSET.

RESULTS

Despite the lack of strong genetic correlation between longevity and the two diseases, we identified 38 genome-wide significant lead SNPs across 22 independent genomic loci. Of them 6 were found to be potentially shared among the three traits mapping to genes including DAB2IP, DNM2, FCHO1, CLPTM1, and SNRPD2. We also identified 19 novel genome-wide associations for the individual traits in this study. Functional annotations and biological pathway enrichment analyses suggested that pleiotropic variants are involved in clathrin-mediated endocytosis and plasma lipoprotein and neurotransmitter clearance processes.

DISCUSSION

In summary, we have been able to advance in the knowledge of the genetic overlap existing among longevity and the two most common age-related disorders.

Working memory and reaction time variability mediate the relationship between polygenic risk and ADHD traits in a general population sample

Endophenotypes are heritable and quantifiable traits indexing genetic liability for a disorder. Here, we examined three potential endophenotypes, working memory function, response inhibition, and reaction time variability, for attention-deficit hyperactivity disorder (ADHD) measured as a dimensional latent trait in a large general population sample derived from the Adolescent Brain Cognitive Development^SM Study. The genetic risk for ADHD was estimated using polygenic risk scores (PRS) whereas ADHD traits were quantified as a dimensional continuum using Bartlett factor score estimates, derived from Attention Problems items from the Child Behaviour Checklist and Effortful Control items from the Early Adolescent Temperament Questionnaire-Revised. The three candidate cognitive endophenotypes were quantified using task-based performance measures. Higher ADHD PRSs were associated with higher ADHD traits, as well as poorer working memory performance and increased reaction time variability. Lower working memory performance, poorer response inhibition, and increased reaction time variability were associated with more pronounced ADHD traits. Working memory and reaction time variability partially statistically mediated the relationship between ADHD PRS and ADHD traits, explaining 14% and 16% of the association, respectively. The mediation effect was specific to the genetic risk for ADHD and did not generalise to genetic risk for four other major psychiatric disorders. Together, these findings provide robust evidence from a large general population sample that working memory and reaction time variability can be considered endophenotypes for ADHD that mediate the relationship between ADHD PRS and ADHD traits.

Evaluating brain structure traits as endophenotypes using polygenicity and discoverability

Human brain structure traits have been hypothesized to be broad endophenotypes for neuropsychiatric disorders, implying that brain structure traits are comparatively "closer to the underlying biology." Genome-wide association studies from large sample sizes allow for the comparison of common variant genetic architectures between traits to test the evidence supporting this claim. Endophenotypes, compared to neuropsychiatric disorders, are hypothesized to have less polygenicity, with greater effect size of each susceptible SNP, requiring smaller sample sizes to discover them. Here, we compare polygenicity and discoverability of brain structure traits, neuropsychiatric disorders, and other traits (91 in total) to directly test this hypothesis. We found reduced polygenicity (FDR = 0.01) and increased discoverability (FDR = 3.68 × 10-9 ) of cortical brain structure traits, as compared to aggregated estimates of multiple neuropsychiatric disorders. We predict that ~8 M individuals will be required to explain the full heritability of cortical surface area by genome-wide significant SNPs, whereas sample sizes over 20 M will be required to explain the full heritability of depression. In conclusion, our findings are consistent with brain structure satisfying the higher power criterion of endophenotypes.

Donor Chimerism Study by Single Nucleotide Polymorphism using SYBR green based Real Time PCR

Objective:

To optimize and evaluate a real time PCR of Single Nucleotide Polymorphism by SYBR Green method for detection of donor chimerism after haematopoietic stem cell transplantation.

Methods:

This descriptive study was conducted at Genetic Resource Centre (GRC) Lab Rawalpindi from Oct 2017 - Dec 2019. A total of twenty patients of post haematopoietic stem cell transplant with various haematological disorders were studied to see the status of donor chimerism by using SNP real time PCR using SYBR Green method and short tandem repeat PCR. These patients had undergone allogeneic HSCT from HLA-matched sibling donors at Pakistan Institute of Medical Science and Armed Forces Bone Marrow Transplant Centre.

Results:

Real time PCR using SYBR Green was able to detect significant amount of chimerism in all 20 patients having undergone HSCT. Regarding precision of the real time PCR assay the mean value of donor chimerism was 94.1% (SD 3.96) and by STR PCR it was 95.1% (SD 1.41). The assay was found to be sensitive with a detection limit of <1%.

Conclusion:

Our results demonstrate that SNP analysis by SYBR Green real time PCR may be used for the evaluation of chimerism status in patients having undergone HSCT with a sensitivity of <1%. Hence donor chimerism by this sensitive method can be used in monitoring of chimerism in post-transplant patients with various haematological disorders.

Polygenic architecture and cardiovascular risk of familial combined hyperlipidemia

BACKGROUND AND AIMS

Familial combined hyperlipidemia (FCHL) is one of the most common inherited lipid phenotypes, characterized by elevated plasma concentrations of apolipoprotein B-100 and triglycerides. The genetic inheritance of FCHL remains poorly understood. The goals of this study were to investigate the polygenetic architecture and cardiovascular risk associated with FCHL.

METHODS AND RESULTS

We identified individuals with an FCHL phenotype among 349,222 unrelated participants of European ancestry in the UK Biobank using modified versions of 5 different diagnostic criteria. The prevalence of the FCHL phenotype was 11.44% (n = 39,961), 5.01% (n = 17,485), 1.48% (n = 5,153), 1.10% (n = 3,838), and 0.48% (n = 1,688) according to modified versions of the Consensus Conference, Dutch, Mexico, Brunzell, and Goldstein criteria, respectively. We performed discovery, case-control genome-wide association studies for these different FCHL criteria and identified 175 independent loci associated with FCHL at genome-wide significance. We investigated the association of genetic and clinical risk with FCHL and found that polygenic susceptibility to hypercholesterolemia or hypertriglyceridemia and features of metabolic syndrome were associated with greater prevalence of FCHL. Participants with an FCHL phenotype had a similar risk of incident coronary artery disease compared to participants with monogenic familial hypercholesterolemia (adjusted hazard ratio vs controls [95% confidence interval]: 2.72 [2.31-3.21] and 1.90 [1.30-2.78]).

CONCLUSIONS

These results suggest that, rather than being a single genetic entity, the FCHL phenotype represents a polygenic susceptibility to dyslipidemia in combination with metabolic abnormalities. The cardiovascular risk associated with an FCHL phenotype is similar to that of monogenic familial hypercholesterolemia, despite being ∼5x more common.