CNV-association meta-analysis in 191,161 European adults reveals new loci associated with anthropometric traits

Fibin is a crucial mitochondrial regulatory gene in skeletal muscle development

Fin bud initiation factor homolog (Fibin) is a secreted protein that is relatively conserved among species. It is closely related to fin bud development and can regulate a variety of cellular processes. In our previous high-throughput chromosome conformation capture (Hi-C) study of pig embryonic muscle development, it was found that Fibin has high expression and activity during the development of pig primary muscle fibers. Therefore, we speculated Fibin participated in myogenesis severely. Specific deletion of Fibin in mouse skeletal muscle resulted in abnormal primary muscle fiber development during the embryonic period and a substantial decrease in skeletal muscle mass in adulthood. In vitro, knocking out Fibin in C2C12 cells promoted cell proliferation; however, after myogenic induction, cells lacking Fibin had almost no ability to differentiate into myotubes. During myogenic differentiation, loss of Fibin disrupts the normal function of mitochondria and impairs oxidative phosphorylation, resulting in decrease of NADH and FADH in the electron transport chain. Transmission electron microscopy also showed that mitochondrial morphology of Fibin-deficient C2C12 was impaired. In conclusion, our research has unveiled a novel mechanism of myogenesis regulation in mitochondrial function and potential target Fibin, and improved understanding of a broad range of mitochondrial muscle diseases.

Genome-wide copy number variation association study in anorexia nervosa

This study represents the first large-scale investigation of rare (<1% population frequency) copy number variants (CNVs) in anorexia nervosa (AN). Large, rare CNVs are reported to be causally associated with anthropometric traits, neurodevelopmental disorders, and schizophrenia, yet their role in the genetic basis of AN is unclear. Using genome-wide association study (GWAS) array data from the Anorexia Nervosa Genetics Initiative (ANGI), which included 7414 AN case and 5044 controls, we investigated the association of 67 well-established syndromic CNVs and 178 pleiotropic disease-risk dosage-sensitive CNVs with AN. To identify novel CNV regions (CNVRs) that increase the risk of AN, we conducted genome-wide association studies with a focus on rare CNV-breakpoints (CNV-GWAS). We found no net enrichment of rare CNVs, either deletions or duplications, in AN, and none of the well-established syndromic or pleiotropic CNVs had a significant association with AN status. However, the CNV-GWAS found 21 nominally associated CNVRs that contribute to AN risk, covering protein-coding genes implicated in synaptic function, metabolic/mitochondrial factors, and lipid characteristics, like the CD36 (7q21.11) gene, which transports long-chain fatty acids into cells. CNVRs intersecting genes previously related to neurodevelopmental traits include deletions of NRXN1 intron 5 (2p16.3), IMMP2L (7q31.1), and PTPRD (9p23). Overall, given that our study is well powered to detect the CNV burden level reported for schizophrenia, we can conclude that rare CNVs have a limited role in the etiology of AN, as reported for bipolar disorder. Our nominal associations for the 21 discovered CNVRs are consistent with AN being a metabo-psychiatric trait, as demonstrated by the common genetic architecture of AN, and we provide association results to allow for replication in future research.

Disentangling mechanisms behind the pleiotropic effects of proximal 16p11.2 BP4-5 CNVs

Whereas 16p11.2 BP4-5 copy-number variants (CNVs) represent one of the most pleiotropic etiologies of genomic syndromes in both clinical and population cohorts, the mechanisms leading to such pleiotropy remain understudied. Identifying 73 deletion and 89 duplication carrier individuals among unrelated White British UK Biobank participants, we performed a phenome-wide association study (PheWAS) between the region's copy number and 117 complex traits and diseases, mimicking four dosage models. Forty-six phenotypes (39%) were affected by 16p11.2 BP4-5 CNVs, with the deletion-only, mirror, U-shape, and duplication-only models being the best fit for 30, 10, 4, and 2 phenotypes, respectively, aligning with the stronger deleteriousness of the deletion. Upon individually adjusting CNV effects for either body mass index (BMI), height, or educational attainment (EA), we found that sixteen testable deletion-driven associations-primarily with cardiovascular and metabolic traits-were BMI dependent, with EA playing a more subtle role and no association depending on height. Bidirectional Mendelian randomization supported that 13 out of these 16 associations were secondary consequences of the CNV's impact on BMI. For the 23 traits that remained significantly associated upon individual adjustment for mediators, matched-control analyses found that 10 phenotypes, including musculoskeletal traits, liver enzymes, fluid intelligence, platelet count, and pneumonia and acute kidney injury risk, remained associated under strict Bonferroni correction, with 10 additional nominally significant associations. These results paint a complex picture of 16p11.2 BP4-5's pleiotropic pattern that involves direct effects on multiple physiological systems and indirect co-morbidities consequential to the CNV's impact on BMI and EA, acting through trait-specific dosage mechanisms.

The pleiotropic spectrum of proximal 16p11.2 CNVs

Recurrent genomic rearrangements at 16p11.2 BP4-5 represent one of the most common causes of genomic disorders. Originally associated with increased risk for autism spectrum disorder, schizophrenia, and intellectual disability, as well as adiposity and head circumference, these CNVs have since been associated with a plethora of phenotypic alterations, albeit with high variability in expressivity and incomplete penetrance. Here, we comprehensively review the pleiotropy associated with 16p11.2 BP4-5 rearrangements to shine light on its full phenotypic spectrum. Illustrating this phenotypic heterogeneity, we expose many parallels between findings gathered from clinical versus population-based cohorts, which often point to the same physiological systems, and emphasize the role of the CNV beyond neuropsychiatric and anthropometric traits. Revealing the complex and variable clinical manifestations of this CNV is crucial for accurate diagnosis and personalized treatment strategies for carrier individuals. Furthermore, we discuss areas of research that will be key to identifying factors contributing to phenotypic heterogeneity and gaining mechanistic insights into the molecular pathways underlying observed associations, while demonstrating how diversity in affected individuals, cohorts, experimental models, and analytical approaches can catalyze discoveries.

Germline copy number variants and endometrial cancer risk

Known risk loci for endometrial cancer explain approximately one third of familial endometrial cancer. However, the association of germline copy number variants (CNVs) with endometrial cancer risk remains relatively unknown. We conducted a genome-wide analysis of rare CNVs overlapping gene regions in 4115 endometrial cancer cases and 17,818 controls to identify functionally relevant variants associated with disease. We identified a 1.22-fold greater number of CNVs in DNA samples from cases compared to DNA samples from controls (p = 4.4 × 10-63). Under three models of putative CNV impact (deletion, duplication, and loss of function), genome-wide association studies identified 141 candidate gene loci associated (p < 0.01) with endometrial cancer risk. Pathway analysis of the candidate loci revealed an enrichment of genes involved in the 16p11.2 proximal deletion syndrome, driven by a large recurrent deletion (chr16:29,595,483-30,159,693) identified in 0.15% of endometrial cancer cases and 0.02% of control participants. Together, these data provide evidence that rare copy number variants have a role in endometrial cancer susceptibility and that the proximal 16p11.2 BP4-BP5 region contains 25 candidate risk gene(s) that warrant further analysis to better understand their role in human disease.

Whole-genome resequencing uncovers diversity and selective sweep in Kazakh cattle

The Kazakh cattle in the Xinjiang Uygur Autonomous Region of China are highly adaptable and have multiple uses, including milk and meat production, and use as draft animals. They are an excellent original breed that could be enhanced by breeding and hybrid improvement. However, the genomic diversity and signature of selection underlying the germplasm characteristics require further elucidation. Herein, we evaluated 26 Kazakh cattle genomes in comparison with 103 genomes of seven other cattle breeds from regions around the world to assess the Kazakh cattle genetic variability. We revealed that the relatively low linkage disequilibrium at large SNP distances was strongly correlated with the largest effective population size among Kazakh cattle. Using population structural analysis, we next demonstrated a taurine lineage with restricted Bos indicus introgression among Kazakh cattle. Notably, we identified putative selected genes associated with resistance to disease and body size within Kazakh cattle. Together, our findings shed light on the evolutionary history and breeding profile of Kazakh cattle, as well as offering indispensable resources for germplasm resource conservation and crossbreeding program implementation.

Tall stature and gigantism in adult patients with acromegaly

OBJECTIVES

Increased height in patients with acromegaly could be a manifestation of growth hormone (GH) excess before epiphysis closure. The aim of this study was to evaluate the relationship between the height of adult patients with GH excess related to mid-parental height (MPH) and population mean and to find whether taller patients with acromegaly come from tall families.

METHODS

This is a single-centre, observational study involving 135 consecutive patients with acromegaly diagnosed as adults and no family history of GH excess. We established three categories for height for patients with acromegaly: normal stature, tall stature (TS, height above the 97th percentile (1.88 standard deviations (SD)) to <3 SD for gender- and country-specific data or as a height which was greater than 1.5 SD but less than 2 SD above the MPH) and gigantism (height which was greater than 3 SD) above the gender- and country-specific mean or greater than 2 SD above MPH).

RESULTS

Thirteen percent (17/135) of patients (53% females) met the criteria for gigantism, 10% (14/135) fulfilled the criteria for TS (57% females). Parents and adult siblings were not taller than the population mean.

CONCLUSION

In a group of 135 consecutive adult patients with acromegaly, 23% had increased height based on country-specific and MPH data: 13% presented with gigantism while 10% had TS. The frequency of gigantism and TS in patients diagnosed with GH excess as adults is not higher in males than in females. Patients with acromegaly come from normal-stature families.

The association between DNA methylation and human height and a prospective model of DNA methylation-based height prediction

As a vital anthropometric characteristic, human height information not only helps to understand overall developmental status and genetic risk factors, but is also important for forensic DNA phenotyping. We utilized linear regression analysis to test the association between each CpG probe and the height phenotype. Next, we designed a methylation sequencing panel targeting 959 CpGs and subsequent height inference models were constructed for the Chinese population. A total of 11,730 height-associated sites were identified. By employing KPCA and deep neural networks, a prediction model was developed, of which the cross-validation RMSE, MAE and R2 were 5.62 cm, 4.45 cm and 0.64, respectively. Genetic factors could explain 39.4% of the methylation level variance of sites used in the height inference models. Collectively, we demonstrated an association between height and DNA methylation status through an EWAS analysis. Targeted methylation sequencing of only 959 CpGs combined with deep learning techniques could provide a model to estimate human height with higher accuracy than SNP-based prediction models.

Rare copy-number variants as modulators of common disease susceptibility

BACKGROUND

Copy-number variations (CNVs) have been associated with rare and debilitating genomic disorders (GDs) but their impact on health later in life in the general population remains poorly described.

METHODS

Assessing four modes of CNV action, we performed genome-wide association scans (GWASs) between the copy-number of CNV-proxy probes and 60 curated ICD-10 based clinical diagnoses in 331,522 unrelated white British UK Biobank (UKBB) participants with replication in the Estonian Biobank.

RESULTS

We identified 73 signals involving 40 diseases, all of which indicating that CNVs increased disease risk and caused earlier onset. We estimated that 16% of these associations are indirect, acting by increasing body mass index (BMI). Signals mapped to 45 unique, non-overlapping regions, nine of which being linked to known GDs. Number and identity of genes affected by CNVs modulated their pathogenicity, with many associations being supported by colocalization with both common and rare single-nucleotide variant association signals. Dissection of association signals provided insights into the epidemiology of known gene-disease pairs (e.g., deletions in BRCA1 and LDLR increased risk for ovarian cancer and ischemic heart disease, respectively), clarified dosage mechanisms of action (e.g., both increased and decreased dosage of 17q12 impacted renal health), and identified putative causal genes (e.g., ABCC6 for kidney stones). Characterization of the pleiotropic pathological consequences of recurrent CNVs at 15q13, 16p13.11, 16p12.2, and 22q11.2 in adulthood indicated variable expressivity of these regions and the involvement of multiple genes. Finally, we show that while the total burden of rare CNVs-and especially deletions-strongly associated with disease risk, it only accounted for ~ 0.02% of the UKBB disease burden. These associations are mainly driven by CNVs at known GD CNV regions, whose pleiotropic effect on common diseases was broader than anticipated by our CNV-GWAS.

CONCLUSIONS

Our results shed light on the prominent role of rare CNVs in determining common disease susceptibility within the general population and provide actionable insights for anticipating later-onset comorbidities in carriers of recurrent CNVs.

Chromosomal deletions on 16p11.2 encompassing SH2B1 are associated with accelerated metabolic disease

New approaches are needed to treat people whose obesity and type 2 diabetes (T2D) are driven by specific mechanisms. We investigate a deletion on chromosome 16p11.2 (breakpoint 2-3 [BP2-3]) encompassing SH2B1, a mediator of leptin and insulin signaling. Phenome-wide association scans in the UK (N = 502,399) and Estonian (N = 208,360) biobanks show that deletion carriers have increased body mass index (BMI; p = 1.3 × 10-10) and increased rates of T2D. Compared with BMI-matched controls, deletion carriers have an earlier onset of T2D, with poorer glycemic control despite higher medication usage. Cystatin C, a biomarker of kidney function, is significantly elevated in deletion carriers, suggesting increased risk of renal impairment. In a Mendelian randomization study, decreased SH2B1 expression increases T2D risk (p = 8.1 × 10-6). We conclude that people with 16p11.2 BP2-3 deletions have early, complex obesity and T2D and may benefit from therapies that enhance leptin and insulin signaling.

A Review of Anthropometric Measurements for Saudi Adults and Elderly, Directions for Future Work and Recommendations to Establish Saudi Guidelines in Line with the Saudi 2030 Vision

Body weight is a significant risk factor for the disease burden of noncommunicable diseases (NCDs). Anthropometric measurements are the first step in determining NCDs risk, and clinicians must have access to valid cutoffs. This study aims to review the literature of Saudi national guidelines and studies previously conducted in Saudi Arabia (SA) and to provide insights and recommendations to establish national guidelines in anthropometric measurements for Saudi adults/elderly in line with the Saudi 2030 Vision. In total, 163 studies were included, and 12 of them contributed to the development of specific anthropometric cutoffs. Cutoffs for metabolic syndrome, waist circumference, and body mass index were established in Saudi adults. However, limited studies were conducted in the elderly. This review warrants establishing standard cutoffs of Saudi adult anthropometrics to avoid over/underreporting of malnutrition and adiposity. This review will help policymakers and the Ministry of Health to establish national guidelines and standard cutoffs to be used in SA for anthropometric measurements that may assist in detecting malnutrition and NCDs.

Functional characterization of human genomic variation linked to polygenic diseases

The burden of human disease lies predominantly in polygenic diseases. Since the early 2000s, genome-wide association studies (GWAS) have identified genetic variants and loci associated with complex traits. These have ranged from variants in coding sequences to mutations in regulatory regions, such as promoters and enhancers, as well as mutations affecting mediators of mRNA stability and other downstream regulators, such as 5' and 3'-untranslated regions (UTRs), long noncoding RNA (lncRNA), and miRNA. Recent research advances in genetics have utilized a combination of computational techniques, high-throughput in vitro and in vivo screening modalities, and precise genome editing to impute the function of diverse classes of genetic variants identified through GWAS. In this review, we highlight the vastness of genomic variants associated with polygenic disease risk and address recent advances in how genetic tools can be used to functionally characterize them.

Keeping the balance: Trade-offs between human brain evolution, autism, and schizophrenia

The unique qualities of the human brain are a product of a complex evolutionary process. Evolution, famously described by François Jacob as a “tinkerer,” builds upon existing genetic elements by modifying and repurposing them for new functions. Genetic changes in DNA may lead to the emergence of new genes or cause altered gene expression patterns. Both gene and regulatory element mutations may lead to new functions. Yet, this process may lead to side-effects. An evolutionary trade-off occurs when an otherwise beneficial change, which is important for evolutionary success and is under strong positive selection, concurrently results in a detrimental change in another trait. Pleiotropy occurs when a gene affects multiple traits. Antagonistic pleiotropy is a phenomenon whereby a genetic variant leads to an increase in fitness at one life-stage or in a specific environment, but simultaneously decreases fitness in another respect. Therefore, it is conceivable that the molecular underpinnings of evolution of highly complex traits, including brain size or cognitive ability, under certain conditions could result in deleterious effects, which would increase the susceptibility to psychiatric or neurodevelopmental diseases. Here, we discuss possible trade-offs and antagonistic pleiotropies between evolutionary change in a gene sequence, dosage or activity and the susceptibility of individuals to autism spectrum disorders and schizophrenia. We present current knowledge about genes and alterations in gene regulatory landscapes, which have likely played a role in establishing human-specific traits and have been implicated in those diseases.

Adult Height, 22q11.2 Deletion Extent, and Short Stature in 22q11.2 Deletion Syndrome

The 22q11.2 deletion syndrome (22q11.2DS) manifests as a wide range of medical conditions across a number of systems. Pediatric growth deficiency with some catch-up growth is reported, but there are few studies of final adult height. We aimed to investigate how final adult height in 22q11.2DS compared with general population norms, and to examine predictors of short stature in in a cohort of 397 adults with 22q11.2DS (aged 17.6–76.3 years) with confirmed typical 22q11.2 microdeletion (overlapping the LCR22A to LCR22B region). We defined short stature as <3rd percentile using population norms. For the subset (n = 314, 79.1%) with 22q11.2 deletion extent, we used a binomial logistic regression model to predict short stature in 22q11.2DS, accounting for effects of sex, age, ancestry, major congenital heart disease (CHD), moderate-to-severe intellectual disability (ID), and 22q11.2 deletion extent. Adult height in 22q11.2DS showed a normal distribution but with a shift to the left, compared with population norms. Those with short stature represented 22.7% of the 22q11.2DS sample, 7.6-fold greater than population expectations (p < 0.0001). In the regression model, moderate-to-severe ID, major CHD, and the common LCR22A-LCR22D (A-D) deletion were significant independent risk factors for short stature while accounting for other factors (model p = 0.0004). The results suggest that the 22q11.2 microdeletion has a significant effect on final adult height distribution, and on short stature with effects appearing to arise from reduced gene dosage involving both the proximal and distal sub-regions of the A-D region. Future studies involving larger sample sizes with proximal nested 22q11.2 deletions, longitudinal lifetime data, parental heights, and genotype data will be valuable.

Influences of rare copy-number variation on human complex traits

The human genome contains hundreds of thousands of regions harboring copy-number variants (CNV). However, the phenotypic effects of most such polymorphisms are unknown because only larger CNVs have been ascertainable from SNP-array data generated by large biobanks. We developed a computational approach leveraging haplotype sharing in biobank cohorts to more sensitively detect CNVs. Applied to UK Biobank, this approach accounted for approximately half of all rare gene inactivation events produced by genomic structural variation. This CNV call set enabled a detailed analysis of associations between CNVs and 56 quantitative traits, identifying 269 independent associations (p < 5 × 10-8) likely to be causally driven by CNVs. Putative target genes were identifiable for nearly half of the loci, enabling insights into dosage sensitivity of these genes and uncovering several gene-trait relationships. These results demonstrate the ability of haplotype-informed analysis to provide insights into the genetic basis of human complex traits.

Omics-informed CNV calls reduce false-positive rates and improve power for CNV-trait associations

Copy-number variations (CNV) are believed to play an important role in a wide range of complex traits, but discovering such associations remains challenging. While whole-genome sequencing (WGS) is the gold-standard approach for CNV detection, there are several orders of magnitude more samples with available genotyping microarray data. Such array data can be exploited for CNV detection using dedicated software (e.g., PennCNV); however, these calls suffer from elevated false-positive and -negative rates. In this study, we developed a CNV quality score that weights PennCNV calls (pCNVs) based on their likelihood of being true positive. First, we established a measure of pCNV reliability by leveraging evidence from multiple omics data (WGS, transcriptomics, and methylomics) obtained from the same samples. Next, we built a predictor of omics-confirmed pCNVs, termed omics-informed quality score (OQS), using only PennCNV software output parameters. Promisingly, OQS assigned to pCNVs detected in close family members was up to 35% higher than the OQS of pCNVs not carried by other relatives (p < 3.0 × 10-90), outperforming other scores. Finally, in an association study of four anthropometric traits in 89,516 Estonian Biobank samples, the use of OQS led to a relative increase in the trait variance explained by CNVs of up to 56% compared with published quality filtering methods or scores. Overall, we put forward a flexible framework to improve any CNV detection method leveraging multi-omics evidence, applied it to improve PennCNV calls, and demonstrated its utility by improving the statistical power for downstream association analyses.

CNest: A novel copy number association discovery method uncovers 862 new associations from 200,629 whole-exome sequence datasets in the UK Biobank

Copy number variation (CNV) is known to influence human traits, having a rich history of research into common and rare genetic disease, and although CNV is accepted as an important class of genomic variation, progress on copy-number-based genome-wide association studies (GWASs) from next-generation sequencing (NGS) data has been limited. Here we present a novel method for large-scale copy number analysis from NGS data generating robust copy number estimates and allowing copy number GWASs (CN-GWASs) to be performed genome-wide in discovery mode. We provide a detailed analysis in the UK Biobank resource and a specifically designed software package. We use these methods to perform CN-GWAS analysis across 78 human traits, discovering over 800 genetic associations that are likely to contribute strongly to trait distributions. Finally, we compare CNV and SNP association signals across the same traits and samples, defining specific CNV association classes.

A cross-disorder dosage sensitivity map of the human genome

Rare copy-number variants (rCNVs) include deletions and duplications that occur infrequently in the global human population and can confer substantial risk for disease. In this study, we aimed to quantify the properties of haploinsufficiency (i.e., deletion intolerance) and triplosensitivity (i.e., duplication intolerance) throughout the human genome. We harmonized and meta-analyzed rCNVs from nearly one million individuals to construct a genome-wide catalog of dosage sensitivity across 54 disorders, which defined 163 dosage sensitive segments associated with at least one disorder. These segments were typically gene dense and often harbored dominant dosage sensitive driver genes, which we were able to prioritize using statistical fine-mapping. Finally, we designed an ensemble machine-learning model to predict probabilities of dosage sensitivity (pHaplo & pTriplo) for all autosomal genes, which identified 2,987 haploinsufficient and 1,559 triplosensitive genes, including 648 that were uniquely triplosensitive. This dosage sensitivity resource will provide broad utility for human disease research and clinical genetics.

The individual and global impact of copy-number variants on complex human traits

The impact of copy-number variations (CNVs) on complex human traits remains understudied. We called CNVs in 331,522 UK Biobank participants and performed genome-wide association studies (GWASs) between the copy number of CNV-proxy probes and 57 continuous traits, revealing 131 signals spanning 47 phenotypes. Our analysis recapitulated well-known associations (e.g., 1q21 and height), revealed the pleiotropy of recurrent CNVs (e.g., 26 and 16 traits for 16p11.2-BP4-BP5 and 22q11.21, respectively), and suggested gene functionalities (e.g., MARF1 in female reproduction). Forty-eight CNV signals (38%) overlapped with single-nucleotide polymorphism (SNP)-GWASs signals for the same trait. For instance, deletion of PDZK1, which encodes a urate transporter scaffold protein, decreased serum urate levels, while deletion of RHD, which encodes the Rhesus blood group D antigen, associated with hematological traits. Other signals overlapped Mendelian disorder regions, suggesting variable expressivity and broad impact of these loci, as illustrated by signals mapping to Rotor syndrome (SLCO1B1/3), renal cysts and diabetes syndrome (HNF1B), or Charcot-Marie-Tooth (PMP22) loci. Total CNV burden negatively impacted 35 traits, leading to increased adiposity, liver/kidney damage, and decreased intelligence and physical capacity. Thirty traits remained burden associated after correcting for CNV-GWAS signals, pointing to a polygenic CNV architecture. The burden negatively correlated with socio-economic indicators, parental lifespan, and age (survivorship proxy), suggesting a contribution to decreased longevity. Together, our results showcase how studying CNVs can expand biological insights, emphasizing the critical role of this mutational class in shaping human traits and arguing in favor of a continuum between Mendelian and complex diseases.

SavvyCNV: Genome-wide CNV calling from off-target reads

Identifying copy number variants (CNVs) can provide diagnoses to patients and provide important biological insights into human health and disease. Current exome and targeted sequencing approaches cannot detect clinically and biologically-relevant CNVs outside their target area. We present SavvyCNV, a tool which uses off-target read data from exome and targeted sequencing data to call germline CNVs genome-wide. Up to 70% of sequencing reads from exome and targeted sequencing fall outside the targeted regions. We have developed a new tool, SavvyCNV, to exploit this 'free data' to call CNVs across the genome. We benchmarked SavvyCNV against five state-of-the-art CNV callers using truth sets generated from genome sequencing data and Multiplex Ligation-dependent Probe Amplification assays. SavvyCNV called CNVs with high precision and recall, outperforming the five other tools at calling CNVs genome-wide, using off-target or on-target reads from targeted panel and exome sequencing. We then applied SavvyCNV to clinical samples sequenced using a targeted panel and were able to call previously undetected clinically-relevant CNVs, highlighting the utility of this tool within the diagnostic setting. SavvyCNV outperforms existing tools for calling CNVs from off-target reads. It can call CNVs genome-wide from targeted panel and exome data, increasing the utility and diagnostic yield of these tests. SavvyCNV is freely available at https://github.com/rdemolgen/SavvySuite.

Effects of copy number variations on brain structure and risk for psychiatric illness: Large-scale studies from the ENIGMA working groups on CNVs

The Enhancing NeuroImaging Genetics through Meta-Analysis copy number variant (ENIGMA-CNV) and 22q11.2 Deletion Syndrome Working Groups (22q-ENIGMA WGs) were created to gain insight into the involvement of genetic factors in human brain development and related cognitive, psychiatric and behavioral manifestations. To that end, the ENIGMA-CNV WG has collated CNV and magnetic resonance imaging (MRI) data from ~49,000 individuals across 38 global research sites, yielding one of the largest studies to date on the effects of CNVs on brain structures in the general population. The 22q-ENIGMA WG includes 12 international research centers that assessed over 533 individuals with a confirmed 22q11.2 deletion syndrome, 40 with 22q11.2 duplications, and 333 typically developing controls, creating the largest-ever 22q11.2 CNV neuroimaging data set. In this review, we outline the ENIGMA infrastructure and procedures for multi-site analysis of CNVs and MRI data. So far, ENIGMA has identified effects of the 22q11.2, 16p11.2 distal, 15q11.2, and 1q21.1 distal CNVs on subcortical and cortical brain structures. Each CNV is associated with differences in cognitive, neurodevelopmental and neuropsychiatric traits, with characteristic patterns of brain structural abnormalities. Evidence of gene-dosage effects on distinct brain regions also emerged, providing further insight into genotype-phenotype relationships. Taken together, these results offer a more comprehensive picture of molecular mechanisms involved in typical and atypical brain development. This "genotype-first" approach also contributes to our understanding of the etiopathogenesis of brain disorders. Finally, we outline future directions to better understand effects of CNVs on brain structure and behavior.

Medical manifestations and health care utilization among adult MyCode participants with neurodevelopmental psychiatric copy number variants

PURPOSE

Recurrent pathogenic copy number variants (pCNVs) have large-effect impacts on brain function and represent important etiologies of neurodevelopmental psychiatric disorders (NPDs), including autism and schizophrenia. Patterns of health care utilization in adults with pCNVs have gone largely unstudied and are likely to differ in significant ways from those of children.

METHODS

We compared the prevalence of NPDs and electronic health record-based medical conditions in 928 adults with 26 pCNVs to a demographically-matched cohort of pCNV-negative controls from >135,000 patient-participants in Geisinger's MyCode Community Health Initiative. We also evaluated 3 quantitative health care utilization measures (outpatient, inpatient, and emergency department visits) in both groups.

RESULTS

Adults with pCNVs (24.9%) were more likely than controls (16.0%) to have a documented NPD. They had significantly higher rates of several chronic diseases, including diabetes (29.3% in participants with pCNVs vs 20.4% in participants without pCNVs) and dementia (2.2% in participants with pCNVs vs 1.0% participants without pCNVs), and twice as many annual emergency department visits.

CONCLUSION

These findings highlight the potential for genetic information-specifically, pCNVs-to inform the study of health care outcomes and utilization in adults. If, as our findings suggest, adults with pCNVs have poorer health and require disproportionate health care resources, early genetic diagnosis paired with patient-centered interventions may help to anticipate problems, improve outcomes, and reduce the associated economic burden.

Comprehensive characterization of copy number variation (CNV) called from array, long- and short-read data

BACKGROUND

SNP arrays, short- and long-read genome sequencing are genome-wide high-throughput technologies that may be used to assay copy number variants (CNVs) in a personal genome. Each of these technologies comes with its own limitations and biases, many of which are well-known, but not all of them are thoroughly quantified.

RESULTS

We assembled an ensemble of public datasets of published CNV calls and raw data for the well-studied Genome in a Bottle individual NA12878. This assembly represents a variety of methods and pipelines used for CNV calling from array, short- and long-read technologies. We then performed cross-technology comparisons regarding their ability to call CNVs. Different from other studies, we refrained from using the golden standard. Instead, we attempted to validate the CNV calls by the raw data of each technology.

CONCLUSIONS

Our study confirms that long-read platforms enable recalling CNVs in genomic regions inaccessible to arrays or short reads. We also found that the reproducibility of a CNV by different pipelines within each technology is strongly linked to other CNV evidence measures. Importantly, the three technologies show distinct public database frequency profiles, which differ depending on what technology the database was built on.

Lessons Learned From Neuroimaging Studies of Copy Number Variants: A Systematic Review

Pathogenic copy number variants (CNVs) and aneuploidies alter gene dosage and are associated with neurodevelopmental psychiatric disorders such as autism spectrum disorder and schizophrenia. Brain mechanisms mediating genetic risk for neurodevelopmental psychiatric disorders remain largely unknown, but there is a rapid increase in morphometry studies of CNVs using T1-weighted structural magnetic resonance imaging. Studies have been conducted one mutation at a time, leaving the field with a complex catalog of brain alterations linked to different genomic loci. Our aim was to provide a systematic review of neuroimaging phenotypes across CNVs associated with developmental psychiatric disorders including autism and schizophrenia. We included 76 structural magnetic resonance imaging studies on 20 CNVs at the 15q11.2, 22q11.2, 1q21.1 distal, 16p11.2 distal and proximal, 7q11.23, 15q11-q13, and 22q13.33 (SHANK3) genomic loci as well as aneuploidies of chromosomes X, Y, and 21. Moderate to large effect sizes on global and regional brain morphometry are observed across all genomic loci, which is in line with levels of symptom severity reported for these variants. This is in stark contrast with the much milder neuroimaging effects observed in idiopathic psychiatric disorders. Data also suggest that CNVs have independent effects on global versus regional measures as well as on cortical surface versus thickness. Findings highlight a broad diversity of regional morphometry patterns across genomic loci. This heterogeneity of brain patterns provides insight into the weak effects reported in magnetic resonance imaging studies of cognitive dimension and psychiatric conditions. Neuroimaging studies across many more variants will be required to understand links between gene function and brain morphometry.

Clinical evaluation of patients with a neuropsychiatric risk copy number variant

Several copy number variants (CNVs) have been identified to confer high risk for a range of neuropsychiatric conditions. Because of advances in genetic testing within clinical settings, patients are increasingly receiving diagnoses of copy number variant genomic disorders. However, clinical guidelines surrounding assessment and management are limited. This review synthesises recent research and makes preliminary recommendations regarding the clinical evaluation of patients with neuropsychiatric risk CNVs. We recommend multi-system assessment beyond the initial referral reason, recognition of the potential need for co-ordinated multidisciplinary care, and that interventions take account of relevant multimorbidity. The frequently complex needs of patients with CNVs across the life-course pose challenges for many health care systems and may be best provided for by the establishment of specialist clinics.

The copy number variation and stroke (CaNVAS) risk and outcome study

BACKGROUND AND PURPOSE

The role of copy number variation (CNV) variation in stroke susceptibility and outcome has yet to be explored. The Copy Number Variation and Stroke (CaNVAS) Risk and Outcome study addresses this knowledge gap.

METHODS

Over 24,500 well-phenotyped IS cases, including IS subtypes, and over 43,500 controls have been identified, all with readily available genotyping on GWAS and exome arrays, with case measures of stroke outcome. To evaluate CNV-associated stroke risk and stroke outcome it is planned to: 1) perform Risk Discovery using several analytic approaches to identify CNVs that are associated with the risk of IS and its subtypes, across the age-, sex- and ethnicity-spectrums; 2) perform Risk Replication and Extension to determine whether the identified stroke-associated CNVs replicate in other ethnically diverse datasets and use biomarker data (e.g. methylation, proteomic, RNA, miRNA, etc.) to evaluate how the identified CNVs exert their effects on stroke risk, and lastly; 3) perform outcome-based Replication and Extension analyses of recent findings demonstrating an inverse relationship between CNV burden and stroke outcome at 3 months (mRS), and then determine the key CNV drivers responsible for these associations using existing biomarker data.

RESULTS

The results of an initial CNV evaluation of 50 samples from each participating dataset are presented demonstrating that the existing GWAS and exome chip data are excellent for the planned CNV analyses. Further, some samples will require additional considerations for analysis, however such samples can readily be identified, as demonstrated by a sample demonstrating clonal mosaicism.

CONCLUSION

The CaNVAS study will cost-effectively leverage the numerous advantages of using existing case-control data sets, exploring the relationships between CNV and IS and its subtypes, and outcome at 3 months, in both men and women, in those of African and European-Caucasian descent, this, across the entire adult-age spectrum.

Association of structural variation with cardiometabolic traits in Finns

The contribution of genome structural variation (SV) to quantitative traits associated with cardiometabolic diseases remains largely unknown. Here, we present the results of a study examining genetic association between SVs and cardiometabolic traits in the Finnish population. We used sensitive methods to identify and genotype 129,166 high-confidence SVs from deep whole-genome sequencing (WGS) data of 4,848 individuals. We tested the 64,572 common and low-frequency SVs for association with 116 quantitative traits and tested candidate associations using exome sequencing and array genotype data from an additional 15,205 individuals. We discovered 31 genome-wide significant associations at 15 loci, including 2 loci at which SVs have strong phenotypic effects: (1) a deletion of the ALB promoter that is greatly enriched in the Finnish population and causes decreased serum albumin level in carriers (p = 1.47 × 10-54) and is also associated with increased levels of total cholesterol (p = 1.22 × 10-28) and 14 additional cholesterol-related traits, and (2) a multi-allelic copy number variant (CNV) at PDPR that is strongly associated with pyruvate (p = 4.81 × 10-21) and alanine (p = 6.14 × 10-12) levels and resides within a structurally complex genomic region that has accumulated many rearrangements over evolutionary time. We also confirmed six previously reported associations, including five led by stronger signals in single nucleotide variants (SNVs) and one linking recurrent HP gene deletion and cholesterol levels (p = 6.24 × 10-10), which was also found to be strongly associated with increased glycoprotein level (p = 3.53 × 10-35). Our study confirms that integrating SVs in trait-mapping studies will expand our knowledge of genetic factors underlying disease risk.

1q21.1 distal copy number variants are associated with cerebral and cognitive alterations in humans

Low-frequency 1q21.1 distal deletion and duplication copy number variant (CNV) carriers are predisposed to multiple neurodevelopmental disorders, including schizophrenia, autism and intellectual disability. Human carriers display a high prevalence of micro- and macrocephaly in deletion and duplication carriers, respectively. The underlying brain structural diversity remains largely unknown. We systematically called CNVs in 38 cohorts from the large-scale ENIGMA-CNV collaboration and the UK Biobank and identified 28 1q21.1 distal deletion and 22 duplication carriers and 37,088 non-carriers (48% male) derived from 15 distinct magnetic resonance imaging scanner sites. With standardized methods, we compared subcortical and cortical brain measures (all) and cognitive performance (UK Biobank only) between carrier groups also testing for mediation of brain structure on cognition. We identified positive dosage effects of copy number on intracranial volume (ICV) and total cortical surface area, with the largest effects in frontal and cingulate cortices, and negative dosage effects on caudate and hippocampal volumes. The carriers displayed distinct cognitive deficit profiles in cognitive tasks from the UK Biobank with intermediate decreases in duplication carriers and somewhat larger in deletion carriers-the latter potentially mediated by ICV or cortical surface area. These results shed light on pathobiological mechanisms of neurodevelopmental disorders, by demonstrating gene dose effect on specific brain structures and effect on cognitive function.

Polygenic burden has broader impact on health, cognition, and socioeconomic outcomes than most rare and high-risk copy number variants

Copy number variants (CNVs) are associated with syndromic and severe neurological and psychiatric disorders (SNPDs), such as intellectual disability, epilepsy, schizophrenia, and bipolar disorder. Although considered high-impact, CNVs are also observed in the general population. This presents a diagnostic challenge in evaluating their clinical significance. To estimate the phenotypic differences between CNV carriers and non-carriers regarding general health and well-being, we compared the impact of SNPD-associated CNVs on health, cognition, and socioeconomic phenotypes to the impact of three genome-wide polygenic risk score (PRS) in two Finnish cohorts (FINRISK, n = 23,053 and NFBC1966, n = 4895). The focus was on CNV carriers and PRS extremes who do not have an SNPD diagnosis. We identified high-risk CNVs (DECIPHER CNVs, risk gene deletions, or large [>1 Mb] CNVs) in 744 study participants (2.66%), 36 (4.8%) of whom had a diagnosed SNPD. In the remaining 708 unaffected carriers, we observed lower educational attainment (EA; OR = 0.77 [95% CI 0.66-0.89]) and lower household income (OR = 0.77 [0.66-0.89]). Income-associated CNVs also lowered household income (OR = 0.50 [0.38-0.66]), and CNVs with medical consequences lowered subjective health (OR = 0.48 [0.32-0.72]). The impact of PRSs was broader. At the lowest extreme of PRS for EA, we observed lower EA (OR = 0.31 [0.26-0.37]), lower-income (OR = 0.66 [0.57-0.77]), lower subjective health (OR = 0.72 [0.61-0.83]), and increased mortality (Cox's HR = 1.55 [1.21-1.98]). PRS for intelligence had a similar impact, whereas PRS for schizophrenia did not affect these traits. We conclude that the majority of working-age individuals carrying high-risk CNVs without SNPD diagnosis have a modest impact on morbidity and mortality, as well as the limited impact on income and educational attainment, compared to individuals at the extreme end of common genetic variation. Our findings highlight that the contribution of traditional high-risk variants such as CNVs should be analyzed in a broader genetic context, rather than evaluated in isolation.

Impact of schizophrenia genetic liability on the association between schizophrenia and physical illness: data-linkage study

BACKGROUND

Individuals with schizophrenia are at higher risk of physical illnesses, which are a major contributor to their 20-year reduced life expectancy. It is currently unknown what causes the increased risk of physical illness in schizophrenia.

AIMS

To link genetic data from a clinically ascertained sample of individuals with schizophrenia to anonymised National Health Service (NHS) records. To assess (a) rates of physical illness in those with schizophrenia, and (b) whether physical illness in schizophrenia is associated with genetic liability.

METHOD

We linked genetic data from a clinically ascertained sample of individuals with schizophrenia (Cardiff Cognition in Schizophrenia participants, n = 896) to anonymised NHS records held in the Secure Anonymised Information Linkage (SAIL) databank. Physical illnesses were defined from the General Practice Database and Patient Episode Database for Wales. Genetic liability for schizophrenia was indexed by (a) rare copy number variants (CNVs), and (b) polygenic risk scores.

RESULTS

Individuals with schizophrenia in SAIL had increased rates of epilepsy (standardised rate ratio (SRR) = 5.34), intellectual disability (SRR = 3.11), type 2 diabetes (SRR = 2.45), congenital disorders (SRR = 1.77), ischaemic heart disease (SRR = 1.57) and smoking (SRR = 1.44) in comparison with the general SAIL population. In those with schizophrenia, carrier status for schizophrenia-associated CNVs and neurodevelopmental disorder-associated CNVs was associated with height (P = 0.015-0.017), with carriers being 7.5-7.7 cm shorter than non-carriers. We did not find evidence that the increased rates of poor physical health outcomes in schizophrenia were associated with genetic liability for the disorder.

CONCLUSIONS

This study demonstrates the value of and potential for linking genetic data from clinically ascertained research studies to anonymised health records. The increased risk for physical illness in schizophrenia is not caused by genetic liability for the disorder.

Large Copy-Number Variants in UK Biobank Caused by Clonal Hematopoiesis May Confound Penetrance Estimates

Large copy-number variants (CNVs) are strongly associated with both developmental delay and cancer, but the type of disease depends strongly on when and where the mutation occurred, i.e., germline versus somatic. We used microarray data from UK Biobank to investigate the prevalence and penetrance of large autosomal CNVs and chromosomal aneuploidies using a standard CNV detection algorithm not designed for detecting mosaic variants. We found 160 individuals that carry >10 Mb copy number changes, including 56 with whole chromosome aneuploidies. Nineteen (12%) individuals had a diagnosis of Down syndrome or other developmental disorder, while 84 (52.5%) individuals had a diagnosis of hematological malignancies or chronic myeloproliferative disorders. Notably, there was no evidence of mosaicism in the blood for many of these large CNVs, so they could easily be mistaken for germline alleles even when caused by somatic mutations. We therefore suggest that somatic mutations associated with blood cancers may result in false estimates of rare variant penetrance from population biobanks.

Genome-Wide Detection of CNVs and Association With Body Weight in Sheep Based on 600K SNP Arrays

Copy number variations (CNVs) are important genomic structural variations and can give rise to significant phenotypic diversity. Herein, we used high-density 600K SNP arrays to detect CNVs in two synthetic lines of sheep (DS and SHH) and in Hu sheep (a local Chinese breed). A total of 919 CNV regions (CNVRs) were detected with a total length of 48.17 Mb, accounting for 1.96% of the sheep genome. These CNVRs consisted of 730 gains, 102 losses, and 87 complex CNVRs. These CNVRs were significantly enriched in the segmental duplication (SD) region. A CNVR-based cluster analysis of the three breeds revealed that the DS and SHH breeds share a close genetic relationship. Functional analysis revealed that some genes in these CNVRs were also significantly enriched in the olfactory transduction pathway (oas04740), including members of the OR gene family such as OR6C76, OR4Q2, and OR4K14. Using association analyses and previous gene annotations, we determined that a subset of identified genes was likely to be associated with body weight, including FOXF2, MAPK12, MAP3K11, STRBP, and C14orf132. Together, these results offer valuable information that will guide future efforts to explore the genetic basis for body weight in sheep.

Association test using Copy Number Profile Curves (CONCUR) enhances power in rare copy number variant analysis

Copy number variants (CNVs) are the gain or loss of DNA segments in the genome that can vary in dosage and length. CNVs comprise a large proportion of variation in human genomes and impact health conditions. To detect rare CNV associations, kernel-based methods have been shown to be a powerful tool due to their flexibility in modeling the aggregate CNV effects, their ability to capture effects from different CNV features, and their accommodation of effect heterogeneity. To perform a kernel association test, a CNV locus needs to be defined so that locus-specific effects can be retained during aggregation. However, CNV loci are arbitrarily defined and different locus definitions can lead to different performance depending on the underlying effect patterns. In this work, we develop a new kernel-based test called CONCUR (i.e., copy number profile curve-based association test) that is free from a definition of locus and evaluates CNV-phenotype associations by comparing individuals' copy number profiles across the genomic regions. CONCUR is built on the proposed concepts of "copy number profile curves" to describe the CNV profile of an individual, and the "common area under the curve (cAUC) kernel" to model the multi-feature CNV effects. The proposed method captures the effects of CNV dosage and length, accounts for the numerical nature of copy numbers, and accommodates between- and within-locus etiological heterogeneity without the need to define artificial CNV loci as required in current kernel methods. In a variety of simulation settings, CONCUR shows comparable or improved power over existing approaches. Real data analyses suggest that CONCUR is well powered to detect CNV effects in the Swedish Schizophrenia Study and the Taiwan Biobank.

Concordance rate between copy number variants detected using either high- or medium-density single nucleotide polymorphism genotype panels and the potential of imputing copy number variants from flanking high density single nucleotide polymorphism haplotypes in cattle

Background

The trading of individual animal genotype information often involves only the exchange of the called genotypes and not necessarily the additional information required to effectively call structural variants. The main aim here was to determine if it is possible to impute copy number variants (CNVs) using the flanking single nucleotide polymorphism (SNP) haplotype structure in cattle. While this objective was achieved using high-density genotype panels (i.e., 713,162 SNPs), a secondary objective investigated the concordance of CNVs called with this high-density genotype panel compared to CNVs called from a medium-density panel (i.e., 45,677 SNPs in the present study). This is the first study to compare CNVs called from high-density and medium-density SNP genotypes from the same animals. High (and medium-density) genotypes were available on 991 Holstein-Friesian, 1015 Charolais, and 1394 Limousin bulls. The concordance between CNVs called from the medium-density and high-density genotypes were calculated separately for each animal. A subset of CNVs which were called from the high-density genotypes was selected for imputation. Imputation was carried out separately for each breed using a set of high-density SNPs flanking the midpoint of each CNV. A CNV was deemed to be imputed correctly when the called copy number matched the imputed copy number.

Results

For 97.0% of CNVs called from the high-density genotypes, the corresponding genomic position on the medium-density of the animal did not contain a called CNV. The average accuracy of imputation for CNV deletions was 0.281, with a standard deviation of 0.286. The average accuracy of imputation of the CNV normal state, i.e. the absence of a CNV, was 0.982 with a standard deviation of 0.022. Two CNV duplications were imputed in the Charolais, a single CNV duplication in the Limousins, and a single CNV duplication in the Holstein-Friesians; in all cases the CNV duplications were incorrectly imputed.

Conclusion

The vast majority of CNVs called from the high-density genotypes were not detected using the medium-density genotypes. Furthermore, CNVs cannot be accurately predicted from flanking SNP haplotypes, at least based on the imputation algorithms routinely used in cattle, and using the SNPs currently available on the high-density genotype panel.

Dose response of the 16p11.2 distal copy number variant on intracranial volume and basal ganglia

Carriers of large recurrent copy number variants (CNVs) have a higher risk of developing neurodevelopmental disorders. The 16p11.2 distal CNV predisposes carriers to e.g., autism spectrum disorder and schizophrenia. We compared subcortical brain volumes of 12 16p11.2 distal deletion and 12 duplication carriers to 6882 non-carriers from the large-scale brain Magnetic Resonance Imaging collaboration, ENIGMA-CNV. After stringent CNV calling procedures, and standardized FreeSurfer image analysis, we found negative dose-response associations with copy number on intracranial volume and on regional caudate, pallidum and putamen volumes (β = -0.71 to -1.37; P < 0.0005). In an independent sample, consistent results were obtained, with significant effects in the pallidum (β = -0.95, P = 0.0042). The two data sets combined showed significant negative dose-response for the accumbens, caudate, pallidum, putamen and ICV (P = 0.0032, 8.9 × 10-6, 1.7 × 10-9, 3.5 × 10-12 and 1.0 × 10-4, respectively). Full scale IQ was lower in both deletion and duplication carriers compared to non-carriers. This is the first brain MRI study of the impact of the 16p11.2 distal CNV, and we demonstrate a specific effect on subcortical brain structures, suggesting a neuropathological pattern underlying the neurodevelopmental syndromes.

Functional and population genetic features of copy number variations in two dairy cattle populations

BACKGROUND

Copy Number Variations (CNVs) are gain or loss of DNA segments that are known to play a role in shaping a wide range of phenotypes. In this study, we used two dairy cattle populations, Holstein Friesian and Jersey, to discover CNVs using the Illumina BovineHD Genotyping BeadChip aligned to the ARS-UCD1.2 assembly. The discovered CNVs were investigated for their functional impact and their population genetics features.

RESULTS

We discovered 14,272 autosomal CNVs, which were aggregated into 1755 CNV regions (CNVR) from 451 animals. These CNVRs together cover 2.8% of the bovine autosomes. The assessment of the functional impact of CNVRs showed that rare CNVRs (MAF < 0.01) are more likely to overlap with genes, than common CNVRs (MAF ≥ 0.05). The Population differentiation index (Fst) based on CNVRs revealed multiple highly diverged CNVRs between the two breeds. Some of these CNVRs overlapped with candidate genes such as MGAM and ADAMTS17 genes, which are related to starch digestion and body size, respectively. Lastly, linkage disequilibrium (LD) between CNVRs and BovineHD BeadChip SNPs was generally low, close to 0, although common deletions (MAF ≥ 0.05) showed slightly higher LD (r2 = ~ 0.1 at 10 kb distance) than the rest. Nevertheless, this LD is still lower than SNP-SNP LD (r2 = ~ 0.5 at 10 kb distance).

CONCLUSIONS

Our analyses showed that CNVRs detected using BovineHD BeadChip arrays are likely to be functional. This finding indicates that CNVs can potentially disrupt the function of genes and thus might alter phenotypes. Also, the population differentiation index revealed two candidate genes, MGAM and ADAMTS17, which hint at adaptive evolution between the two populations. Lastly, low CNVR-SNP LD implies that genetic variation from CNVs might not be fully captured in routine animal genetic evaluation, which relies solely on SNP markers.

TASUKE+: a web-based platform for exploring genome-wide association studies results and large-scale resequencing data

Recent revolutionary advancements in sequencing technologies have made it possible to obtain mass quantities of genome-scale sequence data in a cost-effective manner and have drastically altered molecular biological studies. To utilize these sequence data, genome-wide association studies (GWASs) have become increasingly important. Hence, there is an urgent need to develop a visualization tool that enables efficient data retrieval, integration of GWAS results with diverse information and rapid public release of such large-scale genotypic and phenotypic data. We developed a web-based genome browser TASUKE+ (https://tasuke.dna.affrc.go.jp/), which is equipped with the following functions: (i) interactive GWAS results visualization with genome resequencing data and annotation information, (ii) PCR primer design, (iii) phylogenetic tree reconstruction and (iv) data sharing via the web. GWAS results can be displayed in parallel with polymorphism data, read depths and annotation information in an interactive and scalable manner. Users can design PCR primers for polymorphic sites of interest. In addition, a molecular phylogenetic tree of any region can be reconstructed so that the overall relationship among the examined genomes can be understood intuitively at a glance. All functions are implemented through user-friendly web-based interfaces so that researchers can easily share data with collaborators in remote places without extensive bioinformatics knowledge.

Genetic Imbalance Is Associated With Functional Outcome After Ischemic Stroke

Background and Purpose- We sought to explore the effect of genetic imbalance on functional outcome after ischemic stroke (IS). Methods- Copy number variation was identified in high-density single-nucleotide polymorphism microarray data of IS patients from the CADISP (Cervical Artery Dissection and Ischemic Stroke Patients) and SiGN (Stroke Genetics Network)/GISCOME (Genetics of Ischaemic Stroke Functional Outcome) networks. Genetic imbalance, defined as total number of protein-coding genes affected by copy number variations in an individual, was compared between patients with favorable (modified Rankin Scale score of 0-2) and unfavorable (modified Rankin Scale score of ≥3) outcome after 3 months. Subgroup analyses were confined to patients with imbalance affecting ohnologs-a class of dose-sensitive genes, or to those with imbalance not affecting ohnologs. The association of imbalance with outcome was analyzed by logistic regression analysis, adjusted for age, sex, stroke subtype, stroke severity, and ancestry. Results- The study sample comprised 816 CADISP patients (age 44.2±10.3 years) and 2498 SiGN/GISCOME patients (age 67.7±14.2 years). Outcome was unfavorable in 122 CADISP and 889 SiGN/GISCOME patients. Multivariate logistic regression analysis revealed that increased genetic imbalance was associated with less favorable outcome in both samples (CADISP: P=0.0007; odds ratio=0.89; 95% CI, 0.82-0.95 and SiGN/GISCOME: P=0.0036; odds ratio=0.94; 95% CI, 0.91-0.98). The association was independent of age, sex, stroke severity on admission, stroke subtype, and ancestry. On subgroup analysis, imbalance affecting ohnologs was associated with outcome (CADISP: odds ratio=0.88; 95% CI, 0.80-0.95 and SiGN/GISCOME: odds ratio=0.93; 95% CI, 0.89-0.98) whereas imbalance without ohnologs lacked such an association. Conclusions- Increased genetic imbalance was associated with poorer functional outcome after IS in both study populations. Subgroup analysis revealed that this association was driven by presence of ohnologs in the respective copy number variations, suggesting a causal role of the deleterious effects of genetic imbalance.

Quantifying the polygenic contribution to variable expressivity in eleven rare genetic disorders

Rare genetic disorders (RGDs) often exhibit significant clinical variability among affected individuals, a disease characteristic termed variable expressivity. Recently, the aggregate effect of common variation, quantified as polygenic scores (PGSs), has emerged as an effective tool for predictions of disease risk and trait variation in the general population. Here, we measure the effect of PGSs on 11 RGDs including four sex-chromosome aneuploidies (47,XXX; 47,XXY; 47,XYY; 45,X) that affect height; two copy-number variant (CNV) disorders (16p11.2 deletions and duplications) and a Mendelian disease (melanocortin 4 receptor deficiency (MC4R)) that affect BMI; and two Mendelian diseases affecting cholesterol: familial hypercholesterolemia (FH; LDLR and APOB) and familial hypobetalipoproteinemia (FHBL; PCSK9 and APOB). Our results demonstrate that common, polygenic factors of relevant complex traits frequently contribute to variable expressivity of RGDs and that PGSs may be a useful metric for predicting clinical severity in affected individuals and for risk stratification.

Germline 16p11.2 Microdeletion Predisposes to Neuroblastoma

Neuroblastoma is a cancer of the developing sympathetic nervous system. It is diagnosed in 600-700 children per year in the United States and accounts for 12% of pediatric cancer deaths. Despite recent advances in our understanding of this malignancy's complex genetic architecture, the contribution of rare germline variants remains undefined. Here, we conducted a genome-wide analysis of large (>500 kb), rare (<1%) germline copy number variants (CNVs) in two independent, multi-ethnic cohorts totaling 5,585 children with neuroblastoma and 23,505 cancer-free control children. We identified a 550-kb deletion on chromosome 16p11.2 significantly enriched in neuroblastoma cases (0.39% of cases and 0.03% of controls; p = 3.34 × 10^-9). Notably, this CNV corresponds to a known microdeletion syndrome that affects approximately one in 3,000 children and confers risk for diverse developmental phenotypes including autism spectrum disorder and other neurodevelopmental disorders. The CNV had a substantial impact on neuroblastoma risk, with an odds ratio of 13.9 (95% confidence interval = 5.8-33.4). The association remained significant when we restricted our analysis to individuals of European ancestry in order to mitigate potential confounding by population stratification (0.42% of cases and 0.03% of controls; p = 4.10 × 10^-8). We used whole-genome sequencing (WGS) to validate the deletion in paired germline and tumor DNA from 12 cases. Finally, WGS of four parent-child trios revealed that the deletion primarily arose de novo without maternal or paternal bias. This finding expands the clinical phenotypes associated with 16p11.2 microdeletion syndrome to include cancer, and it suggests that disruption of the 16p11.2 region may dysregulate neurodevelopmental pathways that influence both neurological phenotypes and neuroblastoma.

Phenome-wide Burden of Copy-Number Variation in the UK Biobank

Copy-number variations (CNVs) represent a significant proportion of the genetic differences between individuals and many CNVs associate causally with syndromic disease and clinical outcomes. Here, we characterize the landscape of copy-number variation and their phenome-wide effects in a sample of 472,228 array-genotyped individuals from the UK Biobank. In addition to population-level selection effects against genic loci conferring high mortality, we describe genetic burden from potentially pathogenic and previously uncharacterized CNV loci across more than 3,000 quantitative and dichotomous traits, with separate analyses for common and rare classes of variation. Specifically, we highlight the effects of CNVs at two well-known syndromic loci 16p11.2 and 22q11.2, previously uncharacterized variation at 9p23, and several genic associations in the context of acute coronary artery disease and high body mass index. Our data constitute a deeply contextualized portrait of population-wide burden of copy-number variation, as well as a series of dosage-mediated genic associations across the medical phenome.

Benefits of the Mediterranean diet: Epidemiological and molecular aspects

More than 50 years after the Seven Countries Study, a large number of epidemiological studies have explored the relationship between the Mediterranean diet (MD) and health, through observational, case-control, some longitudinal and a few experimental studies. The overall results show strong evidence suggesting a protective effect of the MD mainly on the risk of cardiovascular disease (CVD) and certain types of cancer. The beneficial effects have been attributed to the types of food consumed, total dietary pattern, components in the food, cooking techniques, eating behaviors and lifestyle behaviors, among others. The aim of this article is to review and summarize the knowledge derived from the literature focusing on the benefits of the MD on health, including those that have been extensively investigated (CVD, cancer) along with more recent issues such as mental health, immunity, quality of life, etc. The review begins with a brief description of the MD and its components. Then we present a review of studies evaluating metabolic biomarkers and genotypes in relation to the MD. Other sections are dedicated to observation and intervention studies for various pathologies. Finally, some insights into the relationship between the MD and sustainability are explored. In conclusion, the research undertaken on metabolomics approaches has identified potential markers for certain MD components and patterns, but more investigation is needed to obtain valid measures. Further evaluation of gene-MD interactions are also required to better understand the mechanisms by which the MD diet exerts its beneficial effects on health. Observation and intervention studies, particularly PREDIMED, have provided invaluable data on the benefits of the MD for a wide range of chronic diseases. However further research is needed to explore the effects of other lifestyle components associated with Mediterranean populations, its environmental impact, as well as the MD extrapolation to non-Mediterranean contexts.

Genome-wide Analysis of Common Copy Number Variation and Epithelial Ovarian Cancer Risk

BACKGROUND

Germline DNA copy number variation (CNV) is a ubiquitous source of genetic variation and remains largely unexplored in association with epithelial ovarian cancer (EOC) risk.

METHODS

CNV was quantified in the DNA of approximately 3,500 cases and controls genotyped with the Illumina 610k and HumanOmni2.5M arrays. We performed a genome-wide association study of common (>1%) CNV regions (CNVRs) with EOC and high-grade serous (HGSOC) risk and, using The Cancer Genome Atlas (TCGA), performed in silico analyses of tumor-gene expression.

RESULTS

Three CNVRs were associated (P < 0.01) with EOC risk: two large (∼100 kb) regions within the 610k set and one small (<5 kb) region with the higher resolution 2.5M data. Large CNVRs included a duplication at LILRA6 (OR = 2.57; P = 0.001) and a deletion at CYP2A7 (OR = 1.90; P = 0.007) that were strongly associated with HGSOC risk (OR = 3.02; P = 8.98 × 10^-5). Somatic CYP2A7 alterations correlated with EGLN2 expression in tumors (P = 2.94 × 10^-47). An intronic ERBB4/HER4 deletion was associated with reduced EOC risk (OR = 0.33; P = 9.5 × 10^-2), and somatic deletions correlated with ERBB4 downregulation (P = 7.05 × 10^-5). Five CNVRs were associated with HGSOC, including two reduced-risk deletions: one at 1p36.33 (OR = 0.28; P = 0.001) that correlated with lower CDKIIA expression in TCGA tumors (P = 2.7 × 10^-7), and another at 8p21.2 (OR = 0.52; P = 0.002) that was present somatically where it correlated with lower GNRH1 expression (P = 5.9 × 10^-5).

CONCLUSIONS

Though CNV appears to not contribute largely to EOC susceptibility, a number of low-to-common frequency variants may influence the risk of EOC and tumor-gene expression.

IMPACT

Further research on CNV and EOC susceptibility is warranted, particularly with CNVs estimated from high-density arrays.

Established and emerging strategies to crack the genetic code of obesity

Tremendous progress has been made in the genetic elucidation of obesity over the past two decades, driven largely by technological, methodological and organizational innovations. Current strategies for identifying obesity-predisposing loci/genes, including cytogenetics, linkage analysis, homozygosity mapping, admixture mapping, candidate gene studies, genome-wide association studies, custom genotyping arrays, whole-exome sequencing and targeted exome sequencing, have achieved differing levels of success, and the identified loci in aggregate explain only a modest fraction of the estimated heritability of obesity. This review outlines the successes and limitations of these approaches and proposes novel strategies, including the use of exceptionally large sample sizes, the study of diverse ethnic groups and deep phenotypes and the application of innovative methods and study designs, to identify the remaining obesity-predisposing genes. The use of both established and emerging strategies has the potential to crack the genetic code of obesity in the not-too-distant future. The resulting knowledge is likely to yield improvements in obesity prediction, prevention and care.

Probe-based association analysis identifies several deletions associated with average daily gain in beef cattle

BACKGROUND

Average daily gain (ADG) is an important trait that contributes to the production efficiency and economic benefits in the beef cattle industry. The molecular mechanisms of ADG have not yet been fully explored because most recent association studies for ADG are based on SNPs or haplotypes. We reported a systematic CNV discovery and association analysis for ADG in Chinese Simmental beef cattle.

RESULTS

Our study identified 4912 nonredundant CNVRs with a total length of ~ 248.7 Mb, corresponding to ~ 8.9% of the cattle genome. Using probe-based CNV association, we identified 24 and 12 significant SNP probes within five deletions and two duplications for ADG, respectively. Among them, we found one common deletion with 89 kb imbedded in LHFPL Tetraspan Subfamily Member 6 (LHFPL6) at 22.9 Mb on BTA12, which has high frequency (12.9%) dispersing across population. CNV selection test using V_ST statistic suggested this common deletion may be under positive selection in Chinese Simmental cattle. Moreover, this deletion was not overlapped with any candidate SNP for ADG compared with previous SNPs-based association studies, suggesting its important role for ADG. In addition, we identified one rare deletion near gene Growth Factor Receptor-bound Protein 10 (GRB10) at 5.1 Mb on BTA4 for ADG using both probe-based association and region-based approaches.

CONCLUSIONS

Our results provided some valuable insights to elucidate the genetic basis of ADG in beef cattle, and these findings offer an alternative perspective to understand the genetic mechanism of complex traits in terms of copy number variations in farm animals.

Insulin translates unfavourable lifestyle into obesity

Lifestyle factors conferring increased diabetes risk are associated with elevated basal insulin levels (hyperinsulinaemia). The latter predicts later obesity in children and adolescents.A causal role of hyperinsulinaemia for adipose tissue growth is probable because pharmacological reduction of insulin secretion lowers body weight in people who are obese. Genetic inactivation of insulin gene alleles in mice also lowers their systemic insulin levels and prevents or ameliorates high-fat diet-induced obesity. Hyperinsulinaemia causes weight gain because of a physiological property of insulin. Insulin levels that are on the high side of normal, or which are slightly elevated, are sufficient to suppress lipolysis and promote lipogenesis in adipocytes. The effect of insulin on glucose transport or hepatic glucose production requires six or two times higher hormone levels, respectively.It seems justified to suggest a lifestyle that avoids high insulin levels in order to limit anabolic fat tissue activity.

Effects of pathogenic CNVs on physical traits in participants of the UK Biobank

Background

Copy number variants (CNVs) have been shown to increase risk for physical anomalies, developmental, psychiatric and medical disorders. Some of them have been associated with changes in weight, height, and other physical traits. As most studies have been performed on children and young people, these effects of CNVs in middle-aged and older people are not well established. The UK Biobank recruited half a million adults who provided a variety of physical measurements. We called all CNVs from the Affymetrix microarrays and selected a set of 54 CNVs implicated as pathogenic (including their reciprocal deletions/duplications) and that were found in five or more persons. Linear regression analysis was used to establish their association with 16 physical traits relevant to human health.

Results

396,725 participants of white British or Irish descent (excluding first-degree relatives) passed our quality control filters. Out of the 864 CNV/trait associations, 214 were significant at a false discovery rate of 0.1, most of them novel. Many of these traits increase risk for adverse health outcomes: e.g. increases in weight, waist-to-hip ratio, pulse rate and body fat composition. Deletions at 16p11.2, 16p12.1, NRXN1 and duplications at 16p13.11 and 22q11.2 produced the highest numbers of significant associations. Five CNVs produced average changes of over one standard deviation for the 16 traits, compared to controls: deletions at 16p11.2 and 22q11.2, and duplications at 3q29, the Williams-Beuren and Potocki-Lupski regions. CNVs at 1q21.1, 2q13, 16p11.2 and 16p11.2 distal, 16p12.1, 17p12 and 17q12 demonstrated one or more mirror image effects of deletions versus duplications.

Conclusions

Carriers of many CNVs should be monitored for physical traits that increase morbidity and mortality. Genes within these CNVs can give insights into biological processes and therapeutic interventions.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5292-7) contains supplementary material, which is available to authorized users.

The causes and consequences of pituitary gigantism

In the general population, height is determined by a complex interplay between genetic and environmental factors. Pituitary gigantism is a rare but very important subgroup of patients with excessive height, as it has an identifiable and clinically treatable cause. The disease is caused by chronic growth hormone and insulin-like growth factor 1 secretion from a pituitary somatotrope adenoma that forms before the closure of the epiphyses. If not controlled effectively, this hormonal hypersecretion could lead to extremely elevated final adult height. The past 10 years have seen marked advances in the understanding of pituitary gigantism, including the identification of genetic causes in ~50% of cases, such as mutations in the AIP gene or chromosome Xq26.3 duplications in X-linked acrogigantism syndrome. Pituitary gigantism has a male preponderance, and patients usually have large pituitary adenomas. The large tumour size, together with the young age of patients and frequent resistance to medical therapy, makes the management of pituitary gigantism complex. Early diagnosis and rapid referral for effective therapy appear to improve outcomes in patients with pituitary gigantism; therefore, a high level of clinical suspicion and efficient use of diagnostic resources is key to controlling overgrowth and preventing patients from reaching very elevated final adult heights.