The genome-wide landscape of copy number variations in the MUSGEN study provides evidence for a founder effect in the isolated Finnish population

Population history modulates the fitness effects of Copy Number Variation in the Roma

We provide the first whole genome Copy Number Variant (CNV) study addressing Roma, along with reference populations from South Asia, the Middle East and Europe. Using CNV calling software for short-read sequence data, we identified 3171 deletions and 489 duplications. Taking into account the known population history of the Roma, as inferred from whole genome nucleotide variation, we could discern how this history has shaped CNV variation. As expected, patterns of deletion variation, but not duplication, in the Roma followed those obtained from single nucleotide polymorphisms (SNPs). Reduced effective population size resulting in slightly relaxed natural selection may explain our observation of an increase in intronic (but not exonic) deletions within Loss of Function (LoF)-intolerant genes. Over-representation analysis for LoF-intolerant gene sets hosting intronic deletions highlights a substantial accumulation of shared biological processes in Roma, intriguingly related to signaling, nervous system and development features, which may be related to the known profile of private disease in the population. Finally, we show the link between deletions and known trait-related SNPs reported in the genome-wide association study (GWAS) catalog, which exhibited even frequency distributions among the studied populations. This suggests that, in general human populations, the strong association between deletions and SNPs associated to biomedical conditions and traits could be widespread across continental populations, reflecting a common background of potentially disease/trait-related CNVs.

Variation in plant Toll/Interleukin-1 receptor domain protein dependence on ENHANCED DISEASE SUSCEPTIBILITY 1

Toll/Interleukin-1 receptor (TIR) domains are integral to immune systems across all kingdoms. In plants, TIRs are present in nucleotide-binding leucine-rich repeat (NLR) immune receptors, NLR-like, and TIR-only proteins. Although TIR-NLR and TIR signaling in plants require the ENHANCED DISEASE SUSCEPTIBILITY 1 (EDS1) protein family, TIRs persist in species that have no EDS1 members. To assess whether particular TIR groups evolved with EDS1, we searched for TIR-EDS1 co-occurrence patterns. Using a large-scale phylogenetic analysis of TIR domains from 39 algal and land plant species, we identified 4 TIR families that are shared by several plant orders. One group occurred in TIR-NLRs of eudicots and another in TIR-NLRs across eudicots and magnoliids. Two further groups were more widespread. A conserved TIR-only group co-occurred with EDS1 and members of this group elicit EDS1-dependent cell death. In contrast, a maize (Zea mays) representative of TIR proteins with tetratricopeptide repeats was also present in species without EDS1 and induced EDS1-independent cell death. Our data provide a phylogeny-based plant TIR classification and identify TIRs that appear to have evolved with and are dependent on EDS1, while others have EDS1-independent activity.

How far musicality and perfect pitch are derived from genetic factors?

There is an agreement about joint genetic and environmental background of musical reception and performance. Musical abilities tend to cluster in families. The studies done on a random population, twins and families of gifted musicians provided a strong support for genetic contribution. Modern biomolecular techniques exploring linkage analysis, variation of gene copy number, scanning for whole-genome expression helped to identify genes, or chromosome regions associated with musical aptitude. Some studies were focused on rare ability to recognize tone without reference that is known as a perfect pitch where a far ethnic differentiation was established. On the other hand, gene deletion leading to dysfunction in amusical individuals also indicated appropriate loci “by negation.” The strongest support for an association of genes with musicality was provided for genes: AVPR1 (12q14.2), SLC6A4 (17q11.2), GALM (2p22), PCDH7 (4p15.1), GATA2 (3q21.3), and few others as well for 4q22, 4q23, and 8q13–21 chromosome bands.

Convergent evidence for the molecular basis of musical traits

To obtain aggregate evidence for the molecular basis of musical abilities and the effects of music, we integrated gene-level data from 105 published studies across multiple species including humans, songbirds and several other animals and used a convergent evidence method to prioritize the top candidate genes. Several of the identified top candidate genes like EGR1, FOS, ARC, BDNF and DUSP1 are known to be activity-dependent immediate early genes that respond to sensory and motor stimuli in the brain. Several other top candidate genes like MAPK10, SNCA, ARHGAP24, TET2, UBE2D3, FAM13A and NUDT9 are located on chromosome 4q21-q24, on the candidate genomic region for music abilities in humans. Functional annotation analyses showed the enrichment of genes involved in functions like cognition, learning, memory, neuronal excitation and apoptosis, long-term potentiation and CDK5 signaling pathway. Interestingly, all these biological functions are known to be essential processes underlying learning and memory that are also fundamental for musical abilities including recognition and production of sound. In summary, our study prioritized top candidate genes related to musical traits.

Global patterns of large copy number variations in the human genome reveal complexity in chromosome organization

Global patterns of copy number variations (CNVs) in chromosomes are required to understand the dynamics of genome organization and complexity. For this study, analysis was performed using the Affymetrix Genome-Wide Human SNP Array 6.0 chip and CytoScan High-Density arrays. We identified a total of 44 109 CNVs from 1715 genomes with a mean of 25 CNVs in an individual, which established the first drafts of population-specific CNV maps providing a rationale for prioritizing chromosomal regions. About 19 905 ancient CNVs were identified across all chromosomes and populations at varying frequencies. CNV count, and sometimes CNV size, contributed to the bulk CNV size of the chromosome. Population specific lengthening and shortening of chromosomal length was observed. Sex bias for CNV presence was largely dependent on ethnicity. Lower CNV inheritance rate was observed for India, compared to YRI and CEU. A total of 33 candidate CNV hotspots from 5382 copy number (CN) variable region (CNVR) clusters were identified. Population specific CNV distribution patterns in p and q arms disturbed the assumption that CNV counts in the p arm are less common compared to long arms, and the CNV occurrence and distribution in chromosomes is length independent. This study unraveled the force of independent evolutionary dynamics on genome organization and complexity across chromosomes and populations.

The landscape of copy number variations in Finnish families with autism spectrum disorders

Rare de novo and inherited copy number variations (CNVs) have been implicated in autism spectrum disorder (ASD) risk. However, the genetic underpinnings of ASD remain unknown in more than 80% of cases. Therefore, identification of novel candidate genes and corroboration of known candidate genes may broaden the horizons of determining genetic risk alleles, and subsequent development of diagnostic testing. Here, using genotyping arrays, we characterized the genetic architecture of rare CNVs (<1% frequency) in a Finnish case-control dataset. Unsurprisingly, ASD cases harbored a significant excess of rare, large (>1 Mb) CNVs and rare, exonic CNVs. The exonic rare de novo CNV rate (∼22.5%) seemed higher compared to previous reports. We identified several CNVs in well-known ASD regions including GSTM1-5, DISC1, FHIT, RBFOX1, CHRNA7, 15q11.2, 15q13.2-q13.3, 17q12, and 22q11.21. Additionally, several novel candidate genes (BDKRB1, BDKRB2, AP2M1, SPTA1, PTH1R, CYP2E1, PLCD3, F2RL1, UQCRC2, LILRB3, RPS9, and COL11A2) were identified through gene prioritization. The majority of these genes belong to neuroactive ligand-receptor interaction pathways, and calcium signaling pathways, thus suggesting that a subset of these novel candidate genes may contribute to ASD risk. Furthermore, several metabolic pathways like caffeine metabolism, drug metabolism, retinol metabolism, and calcium-signaling pathway were found to be affected by the rare exonic ASD CNVs. Additionally, biological processes such as bradykinin receptor activity, endoderm formation and development, and oxidoreductase activity were enriched among the rare exonic ASD CNVs. Overall, our findings may add data about new genes and pathways that contribute to the genetic architecture of ASD.

Global spectrum of copy number variations reveals genome organizational plasticity and proposes new migration routes

Global spectrum of CNVs is required to catalog variations to provide a high-resolution on the dynamics of genome-organization and human migration. In this study, we performed genome-wide genotyping using high-resolution arrays and identified 44,109 CNVs from 1,715 genomes across 12 populations. The study unraveled the force of independent evolutionary dynamics on genome-organizational plasticity across populations. We demonstrated the use of CNV tool to study human migration and identified a second major settlement establishing new migration routes in addition to existing ones.

Impact of copy number variations burden on coding genome in humans using integrated high resolution arrays

Copy number variations (CNVs) alter the transcriptional and translational levels of genes by disrupting the coding structure and this burden of CNVs seems to be a significant contributor to phenotypic variations. Therefore it was necessary to assess the complexities of CNV burden on the coding genome. A total of 1715 individuals from 12 populations were used for CNV analysis in the present investigation. Analysis was performed using Affymetrix Genome-Wide Human SNP Array 6·0 chip and CytoScan High-Density arrays. CNVs were more frequently observed in the coding region than in the non-coding region. CNVs were observed vastly more frequently in the coding region than the non-coding region. CNVs were found to be enriched in the regions containing functional genes (83-96%) compared with the regions containing pseudogenes (4-17%). CNVs across the genome of an individual showed multiple hits across many genes, whose proteins interact physically and function under the same pathway. We identified varying numbers of proteins and degrees of interactions within protein complexes of single individual genomes. This study represents the first draft of a population-specific CNV genes map as well as a cross-populational map. The complex relationship of CNVs on genes and their physically interacting partners unravels many complexities involved in phenotype expression. This study identifies four mechanisms contributing to the complexities caused by the presence of multiple CNVs across many genes in the coding part of the genome.

Clinical interpretation of CNVs with cross-species phenotype data

BACKGROUND

Clinical evaluation of CNVs identified via techniques such as array comparative genome hybridisation (aCGH) involves the inspection of lists of known and unknown duplications and deletions with the goal of distinguishing pathogenic from benign CNVs. A key step in this process is the comparison of the individual's phenotypic abnormalities with those associated with Mendelian disorders of the genes affected by the CNV. However, because often there is not much known about these human genes, an additional source of data that could be used is model organism phenotype data. Currently, almost 6000 genes in mouse and zebrafish are, when knocked out, associated with a phenotype in the model organism, but no disease is known to be caused by mutations in the human ortholog. Yet, searching model organism databases and comparing model organism phenotypes with patient phenotypes for identifying novel disease genes and medical evaluation of CNVs is hindered by the difficulty in integrating phenotype information across species and the lack of appropriate software tools.

METHODS

Here, we present an integrated ranking scheme based on phenotypic matching, degree of overlap with known benign or pathogenic CNVs and the haploinsufficiency score for the prioritisation of CNVs responsible for a patient's clinical findings.

RESULTS

We show that this scheme leads to significant improvements compared with rankings that do not exploit phenotypic information. We provide a software tool called PhenogramViz, which supports phenotype-driven interpretation of aCGH findings based on multiple data sources, including the integrated cross-species phenotype ontology Uberpheno, in order to visualise gene-to-phenotype relations.

CONCLUSIONS

Integrating and visualising cross-species phenotype information on the affected genes may help in routine diagnostics of CNVs.

Compilation of copy number variants identified in phenotypically normal and parous Japanese women

With increasing public concern about infertility and the frequent involvement of chromosomal anomalies in miscarriage, analyses of copy number variations (CNVs) have been used to identify the genomic regions responsible for each process of childbearing. Although associations between CNVs and diseases have been reported, many CNVs have also been identified in healthy individuals. Like other types of mutations, phenotypically indefinite CNVs may have been retained and accumulated during anthropogenesis. Therefore to distinguish causative variants from other variants is a formidable task. Furthermore, because previous studies have predominantly focused on European and African populations, comprehensive detection of common Asian CNVs is eagerly awaited. Here, using a high-resolution genotyping array and samples from 411 Japanese women with normal parity without significant complications, we have compiled 1043 copy number variable regions. In total, the collected regions cover 164 Mb, or up to 0.5% of the genome. The copy number differences in these regions may be irrelevant not only to infertility but also to a wide range of diseases. The utility of this resource in reducing the candidate pathogenetic variants, especially in Japanese subjects, is also demonstrated.

Copy number variants in a highly inbred Iberian porcine strain

We carried out a comprehensive genomic analysis of porcine copy number variants (CNVs) based on whole-genome SNP genotyping data and provided new measures of genomic diversity (number, length and distribution of CNV events) for a highly inbred strain (the Guadyerbas strain). This strain represents one of the most ancient surviving populations of the Iberian breed, and it is currently in serious danger of extinction. CNV detection was conducted on the complete Guadyerbas population, adjusted for genomic waves, and used strict quality criteria, pedigree information and the latest porcine genome annotation. The analysis led to the detection of 65 CNV regions (CNVRs). These regions cover 0.33% of the autosomal genome of this particular strain. Twenty-nine of these CNVRs were identified here for the first time. The relatively low number of detected CNVRs is in line with the low variability and high inbreeding estimated previously for this Iberian strain using pedigree, microsatellite or SNP data. A comparison across different porcine studies has revealed that more than half of these regions overlap with previously identified CNVRs or multicopy regions. Also, a preliminary analysis of CNV detection using whole-genome sequence data for four Guadyerbas pigs showed overlapping for 16 of the CNVRs, supporting their reliability. Some of the identified CNVRs contain relevant functional genes (e.g., the SCD and USP15 genes), which are worth being further investigated because of their importance in determining the quality of Iberian pig products. The CNVR data generated could be useful for improving the porcine genome annotation.

Influence of music on steroid hormones and the relationship between receptor polymorphisms and musical ability: a pilot study

Studies have shown that music confers plasticity to the brain. In a preliminary pilot study, we examined the effect of music listening on steroid hormones and the relationship between steroid hormone receptor polymorphisms and musical ability. Twenty-one subjects (10 males and 11 females) were recruited and divided into musically talented and control groups. The subjects selected (1) music they preferred (chill-inducing music) and (2) music they did not like. Before and after the experiments, saliva was collected to measure the levels of steroid hormones such as testosterone, estradiol, and cortisol. DNA was also isolated from the saliva samples to determine the androgen receptor (AR) and arginine vasopressin receptor 1A genotypes. Advanced Measures of Music Audiation (AMMA) was used to determine the musical ability of the subjects. With both types of music, the cortisol levels decreased significantly in both sexes. The testosterone (T) levels declined in males when they listened to both types of music. In females, the T levels increased in those listening to chill-inducing music but declined when they listened to music they disliked. However, these differences were not significant. The 17-beta estradiol levels increased in males with both types of music, whereas the levels increased with chill-inducing music but declined with disliked music in females. The AMMA scores were higher for the short repeat length-type AR than for the long repeat length-type. Comparisons of AR polymorphisms and T levels before the experiments showed that the T levels were within the low range in the short repeat length-type group and there was a positive relationship with the repeat length, although it was not significant. This is the first study conducted in humans to analyze the relationships between the AR gene, T levels, and musical ability.