Submicroscopic deletion in patients with Williams-Beuren syndrome influences expression levels of the nonhemizygous flanking genes

Cyclophosphamide Pharmacogenomic Variation in Cancer Treatment and Its Effect on Bioactivation and Pharmacokinetics

Cyclophosphamide (Cy) is a prodrug that is mainly bioactivated by cytochrome P450 (CYP) 2B6 enzyme. Several other enzymes are also involved in its bioactivation and affect its kinetics. Previous studies have shown the effect of the enzymes' genetic polymorphisms on Cy kinetics and its clinical outcome. These results were controversial primarily because of the involvement of several interacting enzymes in the Cy metabolic pathway, which can also be affected by several clinical factors as well as other drug interactions. In this review article, we present the effect of CYP2B6 polymorphisms on Cy kinetics since it is the main bioactivating enzyme, as well as discussing all previously reported enzymes and clinical factors that can alter Cy efficacy. Additionally, we present explanations for key Cy side effects related to the nature and site of its bioactivation. Finally, we discuss the role of busulphan in conditioning regimens in the Cy metabolic pathway as a clinical example of drug-drug interactions involving several enzymes. By the end of this article, our aim is to have provided a comprehensive summary of Cy pharmacogenomics and the effect on its kinetics. The utility of these findings in the development of new strategies for Cy personalized patient dose adjustment will aid in the future optimization of patient specific Cy dosages and ultimately in improving clinical outcomes. In conclusion, CYP2B6 and several other enzyme polymorphisms can alter Cy kinetics and consequently the clinical outcomes. However, the precise quantification of Cy kinetics in any individual patient is complex as it is clearly under multifactorial genetic control. Additionally, other clinical factors such as the patient's age, diagnosis, concomitant medications, and clinical status should also be considered.

A phenome-wide association study identifies effects of copy-number variation of VNTRs and multicopy genes on multiple human traits

The human genome contains tens of thousands of large tandem repeats and hundreds of genes that show common and highly variable copy-number changes. Due to their large size and repetitive nature, these variable number tandem repeats (VNTRs) and multicopy genes are generally recalcitrant to standard genotyping approaches and, as a result, this class of variation is poorly characterized. However, several recent studies have demonstrated that copy-number variation of VNTRs can modify local gene expression, epigenetics, and human traits, indicating that many have a functional role. Here, using read depth from whole-genome sequencing to profile copy number, we report results of a phenome-wide association study (PheWAS) of VNTRs and multicopy genes in a discovery cohort of ∼35,000 samples, identifying 32 traits associated with copy number of 38 VNTRs and multicopy genes at 1% FDR. We replicated many of these signals in an independent cohort and observed that VNTRs showing trait associations were significantly enriched for expression QTLs with nearby genes, providing strong support for our results. Fine-mapping studies indicated that in the majority (∼90%) of cases, the VNTRs and multicopy genes we identified represent the causal variants underlying the observed associations. Furthermore, several lie in regions where prior SNV-based GWASs have failed to identify any significant associations with these traits. Our study indicates that copy number of VNTRs and multicopy genes contributes to diverse human traits and suggests that complex structural variants potentially explain some of the so-called "missing heritability" of SNV-based GWASs.

Sustained Trichodysplasia Spinulosa Polyomavirus Viremia Illustrating a Primary Disseminated Infection in a Kidney Transplant Recipient

Novel human polyomaviruses (HPyV) have been recently identified in solid organ transplant recipients. Trichodysplasia spinulosa (TS) is a rare disease associated with immunosuppression and induced by a polyomavirus (TSPyV). We report here a case of primary and disseminated TSPyV infection after kidney transplantation with extensive skin lesions, sustained viremia, and high viral loads in urine specimens, anal, nasal and throat swabs, assessed via specific real-time PCR for TSPyV during a follow-up period of 32 months after transplantation. The detection of TSPyV with a high viral load in respiratory and anal swab samples is compatible with viral replication and thus may suggest potential respiratory and oro-fecal routes of transmission.

Integration of genetic, transcriptomic, and clinical data provides insight into 16p11.2 and 22q11.2 CNV genes

BACKGROUND

Deletions and duplications of the multigenic 16p11.2 and 22q11.2 copy number variant (CNV) regions are associated with brain-related disorders including schizophrenia, intellectual disability, obesity, bipolar disorder, and autism spectrum disorder (ASD). The contribution of individual CNV genes to each of these identified phenotypes is unknown, as well as the contribution of these CNV genes to other potentially subtler health implications for carriers. Hypothesizing that DNA copy number exerts most effects via impacts on RNA expression, we attempted a novel in silico fine-mapping approach in non-CNV carriers using both GWAS and biobank data.

METHODS

We first asked whether gene expression level in any individual gene in the CNV region alters risk for a known CNV-associated behavioral phenotype(s). Using transcriptomic imputation, we performed association testing for CNV genes within large genotyped cohorts for schizophrenia, IQ, BMI, bipolar disorder, and ASD. Second, we used a biobank containing electronic health data to compare the medical phenome of CNV carriers to controls within 700,000 individuals in order to investigate the full spectrum of health effects of the CNVs. Third, we used genotypes for over 48,000 individuals within the biobank to perform phenome-wide association studies between imputed expressions of individual 16p11.2 and 22q11.2 genes and over 1500 health traits.

RESULTS

Using large genotyped cohorts, we found individual genes within 16p11.2 associated with schizophrenia (TMEM219, INO80E, YPEL3), BMI (TMEM219, SPN, TAOK2, INO80E), and IQ (SPN), using conditional analysis to identify upregulation of INO80E as the driver of schizophrenia, and downregulation of SPN and INO80E as increasing BMI. We identified both novel and previously observed over-represented traits within the electronic health records of 16p11.2 and 22q11.2 CNV carriers. In the phenome-wide association study, we found seventeen significant gene-trait pairs, including psychosis (NPIPB11, SLX1B) and mood disorders (SCARF2), and overall enrichment of mental traits.

CONCLUSIONS

Our results demonstrate how integration of genetic and clinical data aids in understanding CNV gene function and implicates pleiotropy and multigenicity in CNV biology.

Combinatorial patterns of gene expression changes contribute to variable expressivity of the developmental delay-associated 16p12.1 deletion

BACKGROUND

Recent studies have suggested that individual variants do not sufficiently explain the variable expressivity of phenotypes observed in complex disorders. For example, the 16p12.1 deletion is associated with developmental delay and neuropsychiatric features in affected individuals, but is inherited in > 90% of cases from a mildly-affected parent. While children with the deletion are more likely to carry additional "second-hit" variants than their parents, the mechanisms for how these variants contribute to phenotypic variability are unknown.

METHODS

We performed detailed clinical assessments, whole-genome sequencing, and RNA sequencing of lymphoblastoid cell lines for 32 individuals in five large families with multiple members carrying the 16p12.1 deletion. We identified contributions of the 16p12.1 deletion and "second-hit" variants towards a range of expression changes in deletion carriers and their family members, including differential expression, outlier expression, alternative splicing, allele-specific expression, and expression quantitative trait loci analyses.

RESULTS

We found that the deletion dysregulates multiple autism and brain development genes such as FOXP1, ANK3, and MEF2. Carrier children also showed an average of 5323 gene expression changes compared with one or both parents, which matched with 33/39 observed developmental phenotypes. We identified significant enrichments for 13/25 classes of "second-hit" variants in genes with expression changes, where 4/25 variant classes were only enriched when inherited from the noncarrier parent, including loss-of-function SNVs and large duplications. In 11 instances, including for ZEB2 and SYNJ1, gene expression was synergistically altered by both the deletion and inherited "second-hits" in carrier children. Finally, brain-specific interaction network analysis showed strong connectivity between genes carrying "second-hits" and genes with transcriptome alterations in deletion carriers.

CONCLUSIONS

Our results suggest a potential mechanism for how "second-hit" variants modulate expressivity of complex disorders such as the 16p12.1 deletion through transcriptomic perturbation of gene networks important for early development. Our work further shows that family-based assessments of transcriptome data are highly relevant towards understanding the genetic mechanisms associated with complex disorders.

Core transcriptional networks in Williams syndrome: IGF1-PI3K-AKT-mTOR, MAPK and actin signaling at the synapse echo autism

Gene networks for disorders of social behavior provide the mechanisms critical for identifying therapeutic targets and biomarkers. Large behavioral phenotypic effects of small human deletions make the positive sociality of Williams syndrome (WS) ideal for determining transcriptional networks for social dysfunction currently based on DNA variations for disorders such as autistic spectrum disorder (ASD) and schizophrenia (SCHZ). Consensus on WS networks has been elusive due to the need for larger cohort size, sensitive genome-wide detection and analytic tools. We report a core set of WS network perturbations in a cohort of 58 individuals (34 with typical, 6 atypical deletions and 18 controls). Genome-wide exon-level expression arrays robustly detected changes in differentially expressed gene (DEG) transcripts from WS deleted genes that ranked in the top 11 of 12 122 transcripts, validated by quantitative reverse transcription PCR, RNASeq and western blots. WS DEG's were strictly dosed in the full but not the atypical deletions that revealed a breakpoint position effect on non-deleted CLIP2, a caveat for current phenotypic mapping based on copy number variants. Network analyses tested the top WS DEG's role in the dendritic spine, employing GeneMANIA to harmonize WS DEGs with comparable query gene-sets. The results indicate perturbed actin cytoskeletal signaling analogous to the excitatory dendritic spines. Independent protein-protein interaction analyses of top WS DEGs generated a 100-node graph annotated topologically revealing three interacting pathways, MAPK, IGF1-PI3K-AKT-mTOR/insulin and actin signaling at the synapse. The results indicate striking similarity of WS transcriptional networks to genome-wide association study-based ASD and SCHZ risk suggesting common network dysfunction for these disorders of divergent sociality.

Molecules, Mechanisms, and Disorders of Self-Domestication: Keys for Understanding Emotional and Social Communication from an Evolutionary Perspective

The neural crest hypothesis states that the phenotypic features of the domestication syndrome are due to a reduced number or disruption of neural crest cells (NCCs) migration, as these cells differentiate at their final destinations and proliferate into different tissues whose activity is reduced by domestication. Comparing the phenotypic characteristics of modern and prehistoric man, it is clear that during their recent evolutionary past, humans also went through a process of self-domestication with a simultaneous prolongation of the period of socialization. This has led to the development of social abilities and skills, especially language, as well as neoteny. Disorders of neural crest cell development and migration lead to many different conditions such as Waardenburg syndrome, Hirschsprung disease, fetal alcohol syndrome, DiGeorge and Treacher-Collins syndrome, for which the mechanisms are already relatively well-known. However, for others, such as Williams-Beuren syndrome and schizophrenia that have the characteristics of hyperdomestication, and autism spectrum disorders, and 7dupASD syndrome that have the characteristics of hypodomestication, much less is known. Thus, deciphering the biological determinants of disordered self-domestication has great potential for elucidating the normal and disturbed ontogenesis of humans, as well as for the understanding of evolution of mammals in general.

A de novo 13q31.3 microduplication encompassing the miR-17 ~ 92 cluster results in features mirroring those associated with Feingold syndrome 2

Hemizygosity of the MIR17HG gene encoding the miR-17 ~ 92 cluster is associated with Feingold syndrome 2 characterized by intellectual disability, skeletal abnormalities, short stature, and microcephaly. Here, we report on a female with a de novo 13q31.3 microduplication encompassing MIR17HG but excluding GPC5. She presented developmental delay, skeletal and digital abnormalities, and features such as tall stature and macrocephaly mirroring those of Feingold syndrome 2 patients. The limited extent of the proband's rearrangement to the miR cluster and the corresponding normal expression level of the neighboring GPC5 in her cells, together with previously described data on affected individuals of two families carrying overlapping duplications of the miR-17 ~ 92 cluster that comprise part of GPC5, who likewise presented macrocephaly, developmental delay, as well as skeletal, digital and stature abnormalities, allow to define a new syndrome due to independent microduplication of the miR-17 ~ 92 cluster.

The Prenatal Diagnosis of Seven Fetuses with 7q11.23 Microdeletion or Microduplication

Objective: There is scant information available about fetuses with 7q11.23 copy number variants (CNVs) found during pregnancy. We studied the clinical significance of 7q11.23 CNVs in prenatal diagnosis. Materials and methods: The amniocentesis was performed on pregnant women who underwent ultrasound (US) of fetal abnormalities. After karyotype analysis, CNVs were detected using BACs-on-Beads (BoBs) technique and chromosome microarray analysis (CMA). Results: Of seven fetuses with CNV of 7q11.23, five had microdeletions and two had microduplications. Case 1 had a 7q11.23 microdeletion along with other CNVs. Case 7 was a newborn with a normal phenotype and 7q11.23 microduplication. Conclusion: The CNVs in 7q11.23 results in many clinical manifestations, but the specificity of clinical features is not high. This study demonstrated that BoBs combined with CMA allows prenatal diagnosis of CNVs involving 7q11.23, and provide a clinical basis for prenatal diagnosis and genetic counseling of such CNVs.

The Ligon lintless -2 Short Fiber Mutation Is Located within a Terminal Deletion of Chromosome 18 in Cotton

Extreme elongation distinguishes about one-fourth of cotton (Gossypium sp.) seed epidermal cells as "lint" fibers, useful for the textile industry, from "fuzz" fibers (<5 mm). Ligon lintless-2 (Li ₂ ), a dominant mutation that results in no lint fiber but normal fuzz fiber, offers insight into pathways and mechanisms that differentiate spinnable cotton from its progenitors. A genetic map developed using 1,545 F2 plants showed that marker CISP15 was 0.4 cM from Li ₂ , and "dominant" simple sequence repeat (SSR) markers (i.e. with null alleles in the Li ₂ genotype) SSR7 and SSR18 showed complete linkage with Li ₂ Nonrandom distribution of markers with null alleles suggests that the Li ₂ phenotype results from a 176- to 221-kb deletion of the terminal region of chromosome 18 that may have been masked in prior pooled-sample mapping strategies. The deletion includes 10 genes with putative roles in fiber development. Two Glycosyltransferase Family 1 genes showed striking expression differences during elongation of wild-type versus Li ₂ fiber, and virus-induced silencing of these genes in the wild type induced Li ₂ -like phenotypes. Further, at least 7 of the 10 putative fiber development genes in the deletion region showed higher expression in the wild type than in Li ₂ mutants during fiber development stages, suggesting coordinated regulation of processes in cell wall development and cell elongation, consistent with the hypothesis that some fiber-related quantitative trait loci comprise closely spaced groups of functionally diverse but coordinately regulated genes.

Defining the Critical Region for Intellectual Disability and Brain Malformations in 6q27 Microdeletions

Terminal microdeletions of the long arm of chromosome 6 are associated with a phenotype that includes multiple brain malformations, intellectual disability, and epilepsy. A 1.7-Mb region has been proposed to contain a gene responsible for the brain anomalies. Here, we present the case of a 12-year-old girl with multiple brain alterations and moderate intellectual disability with a 18-kb deletion in chromosome 6q27, which is smaller than the microdeletions previously described by microarray analysis. We refined the smallest region of overlap possibly associated with the phenotype of brain malformations and intellectual disability to a segment of 325 kb, comprising the DLL1, PSMB1, TBP, and PDCD2 genes since these genes were structurally and/or functionally lost in the smaller deletions described to date. We hypothesize that DLL1 is responsible for brain malformations and possibly interacts with other adjacent genes. The TBP gene encodes a transcription factor which is potentially related to cognitive development. TBP is linked to PSMB1 and PDCD2 in a conserved manner among mammals, suggesting a potential interaction between these genes. In conclusion, the 6q27 microdeletion is a complex syndrome with variable expressivity of brain malformations and intellectual disability phenotypes which are possibly triggered by the 4 genes described and adjacent genes susceptible to gene regulation changes.

An evolving view of copy number variants

Copy number variants (CNVs) are regions of the genome that vary in integer copy number. CNVs, which comprise both amplifications and deletions of DNA sequence, have been identified across all domains of life, from bacteria and archaea to plants and animals. CNVs are an important source of genetic diversity, and can drive rapid adaptive evolution and progression of heritable and somatic human diseases, such as cancer. However, despite their evolutionary importance and clinical relevance, CNVs remain understudied compared to single-nucleotide variants (SNVs). This is a consequence of the inherent difficulties in detecting CNVs at low-to-intermediate frequencies in heterogeneous populations of cells. Here, we discuss molecular methods used to detect CNVs, the limitations associated with using these techniques, and the application of new and emerging technologies that present solutions to these challenges. The goal of this short review and perspective is to highlight aspects of CNV biology that are understudied and define avenues for further research that address specific gaps in our knowledge of these complex alleles. We describe our recently developed method for CNV detection in which a fluorescent gene functions as a single-cell CNV reporter and present key findings from our evolution experiments in Saccharomyces cerevisiae. Using a CNV reporter, we found that CNVs are generated at a high rate and undergo selection with predictable dynamics across independently evolving replicate populations. Many CNVs appear to be generated through DNA replication-based processes that are mediated by the presence of short, interrupted, inverted-repeat sequences. Our results have important implications for the role of CNVs in evolutionary processes and the molecular mechanisms that underlie CNV formation. We discuss the possible extension of our method to other applications, including tracking the dynamics of CNVs in models of human tumors.

Integrative network analysis reveals biological pathways associated with Williams syndrome

BACKGROUND

Williams syndrome (WS) is a neurodevelopmental disorder that has been attributed to heterozygous deletions in chromosome 7q11.23 and exhibits a variety of physical, cognitive, and behavioral features. However, the genetic basis of this phenotypic variability is unclear. In this study, we identified genetic clues underlying these complex phenotypes.

METHODS

Neurobehavioral function was assessed in WS patients and healthy controls. Total RNA was extracted from peripheral blood and subjected to microarray analysis, RNA-sequencing, and qRT-PCR. Weighted gene co-expression network analysis was performed to identify specific alterations related to intermediate disease phenotypes. To functionally interpret each WS-related module, gene ontology and disease-related gene enrichment were examined. We also investigated the micro (mi)RNA expression profiles and miRNA co-expression networks to better explain the regulation of the transcriptome in WS.

RESULTS

Our analysis identified four significant co-expression modules related to intermediate WS phenotypes. Notably, the three upregulated WS-related modules were composed exclusively of genes located outside the 7q11.23 region. They were significantly enriched in genes related to B-cell activation, RNA processing, and RNA transport. BCL11A, which is known for its association with speech disorders and intellectual disabilities, was identified as one of the hub genes in the top WS-related module. Finally, these key upregulated mRNA co-expression modules appear to be inversely correlated with a specific downregulated WS-related miRNA co-expression module.

CONCLUSIONS

Dysregulation of the mRNA/miRNA network involving genes outside of the 7q11.23 region is likely related to the complex phenotypes observed in WS patients.

Whole genome paired-end sequencing elucidates functional and phenotypic consequences of balanced chromosomal rearrangement in patients with developmental disorders

BACKGROUND

Balanced chromosomal rearrangements associated with abnormal phenotype are rare events, but may be challenging for genetic counselling, since molecular characterisation of breakpoints is not performed routinely. We used next-generation sequencing to characterise breakpoints of balanced chromosomal rearrangements at the molecular level in patients with intellectual disability and/or congenital anomalies.

METHODS

Breakpoints were characterised by a paired-end low depth whole genome sequencing (WGS) strategy and validated by Sanger sequencing. Expression study of disrupted and neighbouring genes was performed by RT-qPCR from blood or lymphoblastoid cell line RNA.

RESULTS

Among the 55 patients included (41 reciprocal translocations, 4 inversions, 2 insertions and 8 complex chromosomal rearrangements), we were able to detect 89% of chromosomal rearrangements (49/55). Molecular signatures at the breakpoints suggested that DNA breaks arose randomly and that there was no major influence of repeated elements. Non-homologous end-joining appeared as the main mechanism of repair (55% of rearrangements). A diagnosis could be established in 22/49 patients (44.8%), 15 by gene disruption (KANSL1, FOXP1, SPRED1, TLK2, MBD5, DMD, AUTS2, MEIS2, MEF2C, NRXN1, NFIX, SYNGAP1, GHR, ZMIZ1) and 7 by position effect (DLX5, MEF2C, BCL11B, SATB2, ZMIZ1). In addition, 16 new candidate genes were identified. Systematic gene expression studies further supported these results. We also showed the contribution of topologically associated domain maps to WGS data interpretation.

CONCLUSION

Paired-end WGS is a valid strategy and may be used for structural variation characterisation in a clinical setting.

Downregulation of genes outside the deleted region in individuals with 22q11.2 deletion syndrome

The 22q11.2 deletion syndrome (22q11.2DS) is caused by recurrent hemizygous deletions of chromosome 22q11.2. The phenotype of the syndrome is complex and varies widely among individuals. Little is known about the role of the different genes located in 22q11.2, and we hypothesized that genetic risk factors lying elsewhere in the genome might contribute to the phenotype. Here, we present the whole-genome gene expression data of 11 patients with approximately 3 Mb deletions. Apart from the hemizygous genes mapped to the 22q11.2 region, the TUBA8 and GNAZ genes, neighboring the deleted interval but in normal copy number, showed altered expression. When genes mapped to other chromosomes were considered in the gene expression analysis, a genome-wide dysregulation was observed, with increased or decreased expression levels. The enriched pathways of these genes were related to immune response, a deficiency that is frequently observed in 22q11.2DS patients. We also used the hypothesis-free weighted gene co-expression network analysis (WGCNA), which revealed the co-expression gene network modules with clear connection to mechanisms associated with 22q11.2DS such as immune response and schizophrenia. These findings, combined with the traditional gene expression profile, can be used for the identification of potential pathways and genes not previously considered to be related to the 22q11.2 deletion syndrome.

Copy Number Variations of KLF6 Modulate Gene Transcription and Growth Traits in Chinese Datong Yak (Bos Grunniens)

Copy number variation (CNV) is a significant marker of the genetic and phenotypic diversity among individuals that accounts for complex quantitative traits of phenotype and diseases via modulating gene dosage and disrupting coding regions in the genome. Biochemically, Kruppel-like factor 6 (KLF6) genes plays a significant role in the regulation of cell differentiation and proliferation and muscle development. The aim of this study was to detect the distributions of KLF6 copy number variations (CNVs) in five breeds of domestic yak and to explore their effect on growth traits and gene expression. The data were analyzed by real-time quantitative PCR (qPCR). Our results elucidated that a decreased CNV in the KLF6 gene is more highly associated (p < 0.05) with various growth traits than increased or normal CNVs in six-month-old and five-year-old Datong yak. Nevertheless, negative correlations between the DNA copy number and KLF6 gene expression were observed in the skeletal muscle of adult Datong yak. These results suggest that CNVs of the KLF6 gene could be crucial genomic markers for growth phenotypes of Chinese Datong yak breeds and this finding constitutes the first evidence of the biological role of KLF6 CNVs in Chinese Datong yak breeds.

Activity of Genes with Functions in Human Williams-Beuren Syndrome Is Impacted by Mobile Element Insertions in the Gray Wolf Genome

In canines, transposon dynamics have been associated with a hyper-social behavioral syndrome, although the functional mechanism has yet to be described. We investigate the epigenetic and transcriptional consequences of these behavior-associated mobile element insertions (MEIs) in dogs and Yellowstone gray wolves. We posit that the transposons themselves may not be the causative feature; rather, their transcriptional regulation may exert the functional impact. We survey four outlier transposons associated with hyper-sociability, with the expectation that they are targeted for epigenetic silencing. We predict hyper-methylation of MEIs, suggestive that the epigenetic silencing of and not the MEIs themselves may be driving dysregulation of nearby genes. We found that transposon-derived sequences are significantly hyper-methylated, regardless of their copy number or species. Further, we have assessed transcriptome sequence data and found evidence that MEIs impact the expression levels of six genes (WBSCR17, LIMK1, GTF2I, WBSCR27, BAZ1B, and BCL7B), all of which have known roles in human Williams-Beuren syndrome due to changes in copy number, typically hemizygosity. Although further evidence is needed, our results suggest that a few insertions alter local expression at multiple genes, likely through a cis-regulatory mechanism that excludes proximal methylation.

Segmental duplications: evolution and impact among the current Lepidoptera genomes

Background

Structural variation among genomes is now viewed to be as important as single nucleoid polymorphisms in influencing the phenotype and evolution of a species. Segmental duplication (SD) is defined as segments of DNA with homologous sequence.

Results

Here, we performed a systematic analysis of segmental duplications (SDs) among five lepidopteran reference genomes (Plutella xylostella, Danaus plexippus, Bombyx mori, Manduca sexta and Heliconius melpomene) to understand their potential impact on the evolution of these species. We find that the SDs content differed substantially among species, ranging from 1.2% of the genome in B. mori to 15.2% in H. melpomene. Most SDs formed very high identity (similarity higher than 90%) blocks but had very few large blocks. Comparative analysis showed that most of the SDs arose after the divergence of each linage and we found that P. xylostella and H. melpomene showed more duplications than other species, suggesting they might be able to tolerate extensive levels of variation in their genomes. Conserved ancestral and species specific SD events were assessed, revealing multiple examples of the gain, loss or maintenance of SDs over time. SDs content analysis showed that most of the genes embedded in SDs regions belonged to species-specific SDs (“Unique” SDs). Functional analysis of these genes suggested their potential roles in the lineage-specific evolution. SDs and flanking regions often contained transposable elements (TEs) and this association suggested some involvement in SDs formation. Further studies on comparison of gene expression level between SDs and non-SDs showed that the expression level of genes embedded in SDs was significantly lower, suggesting that structure changes in the genomes are involved in gene expression differences in species.

Conclusions

The results showed that most of the SDs were “unique SDs”, which originated after species formation. Functional analysis suggested that SDs might play different roles in different species. Our results provide a valuable resource beyond the genetic mutation to explore the genome structure for future Lepidoptera research.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-1007-y) contains supplementary material, which is available to authorized users.

Structural variants in genes associated with human Williams-Beuren syndrome underlie stereotypical hypersociability in domestic dogs

Although considerable progress has been made in understanding the genetic basis of morphologic traits (for example, body size and coat color) in dogs and wolves, the genetic basis of their behavioral divergence is poorly understood. An integrative approach using both behavioral and genetic data is required to understand the molecular underpinnings of the various behavioral characteristics associated with domestication. We analyze a 5-Mb genomic region on chromosome 6 previously found to be under positive selection in domestic dog breeds. Deletion of this region in humans is linked to Williams-Beuren syndrome (WBS), a multisystem congenital disorder characterized by hypersocial behavior. We associate quantitative data on behavioral phenotypes symptomatic of WBS in humans with structural changes in the WBS locus in dogs. We find that hypersociability, a central feature of WBS, is also a core element of domestication that distinguishes dogs from wolves. We provide evidence that structural variants in GTF2I and GTF2IRD1, genes previously implicated in the behavioral phenotype of patients with WBS and contained within the WBS locus, contribute to extreme sociability in dogs. This finding suggests that there are commonalities in the genetic architecture of WBS and canine tameness and that directional selection may have targeted a unique set of linked behavioral genes of large phenotypic effect, allowing for rapid behavioral divergence of dogs and wolves, facilitating coexistence with humans.

Diversity and regulatory impact of copy number variation in the primate Macaca fascicularis

BACKGROUND

Copy number variations (CNVs) are a significant source of genetic diversity and commonly found in mammalian genomes. We have generated a genome-wide CNV map for Cynomolgus monkeys (Macaca fascicularis). This crab-eating macaque is the closest animal model to humans that is used in biomedical research.

RESULTS

We show that Cynomolgus monkey CNVs are in general much smaller in size than gene loci and are specific to the population of origin. Genome-wide expression data from five vitally important organs demonstrates that CNVs in close proximity to transcription start sites associate strongly with expression changes. Among these eQTL genes we find an overrepresentation of genes involved in metabolism, receptor activity, and transcription.

CONCLUSION

These results provide evidence that CNVs shape tissue transcriptomes in monkey populations, potentially offering an adaptive advantage. We suggest that this genetic diversity should be taken into account when using Cynomolgus macaques as models.

Modeling Williams syndrome with induced pluripotent stem cells

The development of induced pluripotent stem cells (iPSCs) like never before has opened novel opportunity to study diseases in relevant cell types. In our recent study, Williams syndrome (WS), a rare genetic neurodevelopmental disorder, that is caused by hemizygous deletion of 25-28 genes on chromosome 7, is of interest because of its unique cognitive and social profiles. Little is known about haploinsufficiency effect of those deleted genes on molecular and cellular phenotypes at the neural level due to the lack of relevant human cellular model. Using the cellular reprogramming approach, we reported that WS iPSC-derived neural progenitor cells (NPCs) has increased apoptosis and therefore increased doubling time, which could be rescued by complementation of frizzled 9, one of the genes typically deleted in WS. Moreover, WS iPSC-derived CTIP2-positive pyramidal neurons exhibit morphologic alterations including longer total dendrites and increasing dendritic spine number. In addition, WS iPSC-derived neurons show an increase in calcium transient frequency and synchronized activity likely due to increased number of dendritic spines and synapses. Our work integrated cross-level data from genetics to behavior of WS individuals and revealed altered cellular phenotypes in WS human NPCs and neurons that could be validated in other model systems such as magnetic resonance imaging (MRI) in live subjects and postmortem brain tissues.

Genetic Variants in Isolated Ebstein Anomaly Implicated in Myocardial Development Pathways

Ebstein anomaly (EA) is a rare heart defect in which the tricuspid valve is malformed and displaced. The tricuspid valve abnormalities can lead to backflow of blood from the right ventricle to the right atrium, preventing proper circulation of blood to the lungs. Although the etiology of EA is largely unresolved, increased prevalence of EA in those with a family history of congenital heart disease suggests EA has a genetic component. Copy number variants (CNVs) are a major source of genetic variation and have been implicated in a range of congenital heart defect phenotypes. We performed a systematic, genome-wide search for CNVs in 47 isolated EA cases using genotyping microarrays. In addition, we used a custom HaloPlex panel to sequence three known EA genes and 47 candidate EA genes. We identified 35 candidate CNVs in 24 (51%) EA cases. Rare sequence variants in genes associated with cardiomyopathy were identified in 11 (23%) EA cases. Two CNVs near the transcriptional repressor HEY1, a member of the NOTCH signaling pathway, were identified in three unrelated cases. All other candidate CNVs were each identified in a single case. At least 11 of 35 candidate CNVs include genes involved in myocardial development or function, including multiple genes in the BMP signaling pathway. We identified enrichment of gene sets involved in histone modification and cardiomyocyte differentiation, supporting the involvement of the developing myocardium in the etiology of EA. Gene set enrichment analysis also identified ribosomal RNA processing, a potentially novel pathway of altered cardiac development in EA. Our results suggest an altered myocardial program may contribute to abnormal tricuspid valve development in EA. Future studies should investigate abnormal differentiation of cardiomyocytes as a potential etiological factor in EA.

Mechanisms underlying structural variant formation in genomic disorders

With the recent burst of technological developments in genomics, and the clinical implementation of genome-wide assays, our understanding of the molecular basis of genomic disorders, specifically the contribution of structural variation to disease burden, is evolving quickly. Ongoing studies have revealed a ubiquitous role for genome architecture in the formation of structural variants at a given locus, both in DNA recombination-based processes and in replication-based processes. These reports showcase the influence of repeat sequences on genomic stability and structural variant complexity and also highlight the tremendous plasticity and dynamic nature of our genome in evolution, health and disease susceptibility.

Mitochondrial Dysfunction may explain symptom variation in Phelan-McDermid Syndrome

Phelan-McDermid Syndrome (PMS), which is defined by a deletion within 22q13, demonstrates significant phenotypic variation. Given that six mitochondrial genes are located within 22q13, including complex I and IV genes, we hypothesize that mitochondrial complex activity abnormalities may explain phenotypic variation in PMS symptoms. Complex I, II, II + III and IV activity was measured in 51 PMS participants. Caretakers completed questionnaires and provided genetic information through the PMS foundation registry. Complex activity was abnormal in 59% of PMS participants. Abnormalities were found in complex I and IV but not complex II + III and II activity, consistent with disruption of genes within the 22q13 region. However, complex activity abnormalities were not related to specific gene deletions suggesting a "neighboring effect" of regional deletions on adjacent gene expression. A specific combination of symptoms (autism spectrum disorder, developmental regression, failure-to-thrive, exercise intolerance/fatigue) was associated with complex activity abnormalities. 64% of 106 individuals in the PMS foundation registry who did not have complex activity measured also endorsed this pattern of symptoms. These data suggest that mitochondrial abnormalities, specifically abnormalities in complex I and IV activity, may explain some phenotypic variation in PMS individuals. These results point to novel pathophysiology mechanisms and treatment targets for PMS patients.

Human induced pluripotent stem cell derived neurons as a model for Williams-Beuren syndrome

Background

Williams-Beuren Syndrome (WBS) is caused by the microdeletion of approximately 25 genes on chromosome 7q11.23, and is characterized by a spectrum of cognitive and behavioural features.

Results

We generated cortical neurons from a WBS individual and unaffected (WT) control by directed differentiation of induced pluripotent stem cells (iPSCs). Single cell mRNA analyses and immunostaining demonstrated very efficient production of differentiated cells expressing markers of mature neurons of mixed subtypes and from multiple cortical layers. We found that there was a profound alteration in action potentials, with significantly prolonged WBS repolarization times and a WBS deficit in voltage-activated K⁺ currents. Miniature excitatory synaptic currents were normal, indicating that unitary excitatory synaptic transmission was not altered. Gene expression profiling identified 136 negatively enriched gene sets in WBS compared to WT neurons including gene sets involved in neurotransmitter receptor activity, synaptic assembly, and potassium channel complexes.

Conclusions

Our findings provide insight into gene dysregulation and electrophysiological defects in WBS patient neurons.

Electronic supplementary material

The online version of this article (doi:10.1186/s13041-015-0168-0) contains supplementary material, which is available to authorized users.

Quantitative analysis of chromatin interaction changes upon a 4.3 Mb deletion at mouse 4E2

Background

Circular chromosome conformation capture (4C) has provided important insights into three dimensional (3D) genome organization and its critical impact on the regulation of gene expression. We developed a new quantitative framework based on polymer physics for the analysis of paired-end sequencing 4C (PE-4Cseq) data. We applied this strategy to the study of chromatin interaction changes upon a 4.3 Mb DNA deletion in mouse region 4E2.

Results

A significant number of differentially interacting regions (DIRs) and chromatin compaction changes were detected in the deletion chromosome compared to a wild-type (WT) control. Selected DIRs were validated by 3D DNA FISH experiments, demonstrating the robustness of our pipeline. Interestingly, significant overlaps of DIRs with CTCF/Smc1 binding sites and differentially expressed genes were observed.

Conclusions

Altogether, our PE-4Cseq analysis pipeline provides a comprehensive characterization of DNA deletion effects on chromatin structure and function.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2137-5) contains supplementary material, which is available to authorized users.

Conserved piRNA Expression from a Distinct Set of piRNA Cluster Loci in Eutherian Mammals

The Piwi pathway is deeply conserved amongst animals because one of its essential functions is to repress transposons. However, many Piwi-interacting RNAs (piRNAs) do not base-pair to transposons and remain mysterious in their targeting function. The sheer number of piRNA cluster (piC) loci in animal genomes and infrequent piRNA sequence conservation also present challenges in determining which piC loci are most important for development. To address this question, we determined the piRNA expression patterns of piC loci across a wide phylogenetic spectrum of animals, and reveal that most genic and intergenic piC loci evolve rapidly in their capacity to generate piRNAs, regardless of known transposon silencing function. Surprisingly, we also uncovered a distinct set of piC loci with piRNA expression conserved deeply in Eutherian mammals. We name these loci Eutherian-Conserved piRNA cluster (ECpiC) loci. Supporting the hypothesis that conservation of piRNA expression across ~100 million years of Eutherian evolution implies function, we determined that one ECpiC locus generates abundant piRNAs antisense to the STOX1 transcript, a gene clinically associated with preeclampsia. Furthermore, we confirmed reduced piRNAs in existing mouse mutations at ECpiC-Asb1 and -Cbl, which also display spermatogenic defects. The Asb1 mutant testes with strongly reduced Asb1 piRNAs also exhibit up-regulated gene expression profiles. These data indicate ECpiC loci may be specially adapted to support Eutherian reproduction.

Animal genomes from flies to humans contain many hundreds of non-coding elements called Piwi-interacting RNAs (piRNAs) cluster loci (piC loci). Some of these elements generate piRNAs that direct the silencing of transposable elements, which are pervasive genetic parasites. However, we lack an understanding of the targeting function for the remaining bulk of piRNAs because their loci are not complementarity to transposable elements. In addition, the field does not know if all piC loci are quickly evolving, or if some piC loci might be deeply conserved in piRNA expression, an indication of its potentially functional importance. Our study confirms the highly rapid evolution in piRNA expression capacity for the majority of piC loci in flies and mammals, with many clade- and species-specific piC loci expression patterns. In spite of this, we also discover a cohort of piC loci that are deeply conserved in piRNA expression from the human to the dog, a significantly broad phylogenetic spectrum of eutherian mammals. However, this conservation in piRNA expression ends at non-eutherian mammals like marsupials and monotremes. Existing mutations in two of these Eutherian-Conserved piC (ECpiC) loci impair mouse reproduction and abrogate piRNA production. Therefore, we suggest these ECpiC loci are conserved for piRNA expression due to their important function in eutherian reproduction and stand out as prime candidates for future genetic studies.

Cytochrome P450 Oxidoreductase Influences CYP2B6 Activity in Cyclophosphamide Bioactivation

Introduction

Cyclophosphamide is commonly used as an important component in conditioning prior to hematopoietic stem cell transplantation, a curative treatment for several hematological diseases. Cyclophosphamide is a prodrug activated mainly by cytochrome P450 2B6 (CYP2B6) in the liver. A high degree of inter- and intra-individual variation in cyclophosphamide kinetics has been reported in several studies.

Materials and Methods

Hydroxylation of cyclophosphamide was investigated in vitro using three microsomal batches of CYP2B6*1 with different ratios of POR/CYP expression levels. Twenty patients undergoing hematopoietic stem cell transplantation were also included in the study. All patients received an i.v. infusion of cyclophosphamide (60 mg/kg/day, for two days) as a part of their conditioning. Blood samples were collected from each patient before cyclophosphamide infusion, 6 h after the first dose and before and 6 h after the second dose. POR gene expression was measured by mRNA analysis and the pharmacokinetics of cyclophosphamide and its active metabolite were determined.

Results

A strong correlation between the in vitro intrinsic clearance of cyclophosphamide and the POR/CYP ratio was found. The apparent K_m for CYP2B6.1 was almost constant (3-4 mM), while the CL_int values were proportional to the POR/CYP ratio (3-34 μL/min/nmol CYP). In patients, the average expression of the POR gene in blood was significantly (P <0.001) up-regulated after cyclophosphamide infusion, with high inter-individual variations and significant correlation with the concentration ratio of the active metabolite 4-hydroxy-cyclophosphamide/cyclophosphamide. Nine patients were carriers for POR*28; four patients had relatively high POR expression.

Conclusions

This investigation shows for the first time that POR besides CYP2B6 can influence cyclophosphamide metabolism. Our results indicate that not only CYPs are important, but also POR expression and/or activity may influence cyclophosphamide bioactivation, affecting therapeutic efficacy and treatment related toxicity and hence on clinical outcome. Thus, both POR and CYP genotype and expression levels may have to be taken into account when personalizing treatment schedules to achieve optimal therapeutic drug plasma concentrations of cyclophosphamide.

Mechanisms of Origin, Phenotypic Effects and Diagnostic Implications of Complex Chromosome Rearrangements

Complex chromosome rearrangements (CCRs) are currently defined as structural genome variations that involve more than 2 chromosome breaks and result in exchanges of chromosomal segments. They are thought to be extremely rare, but their detection rate is rising because of improvements in molecular cytogenetic technology. Their population frequency is also underestimated, since many CCRs may not elicit a phenotypic effect. CCRs may be the result of fork stalling and template switching, microhomology-mediated break-induced repair, breakage-fusion-bridge cycles, or chromothripsis. Patients with chromosomal instability syndromes show elevated rates of CCRs due to impaired DNA double-strand break responses during meiosis. Therefore, the putative functions of the proteins encoded by ATM, BLM, WRN, ATR, MRE11, NBS1, and RAD51 in preventing CCRs are discussed. CCRs may exert a pathogenic effect by either (1) gene dosage-dependent mechanisms, e.g. haploinsufficiency, (2) mechanisms based on disruption of the genomic architecture, such that genes, parts of genes or regulatory elements are truncated, fused or relocated and thus their interactions disturbed - these mechanisms will predominantly affect gene expression - or (3) mixed mutation mechanisms in which a CCR on one chromosome is combined with a different type of mutation on the other chromosome. Such inferred mechanisms of pathogenicity need corroboration by mRNA sequencing. Also, future studies with in vitro models, such as inducible pluripotent stem cells from patients with CCRs, and transgenic model organisms should substantiate current inferences regarding putative pathogenic effects of CCRs. The ramifications of the growing body of information on CCRs for clinical and experimental genetics and future treatment modalities are briefly illustrated with 2 cases, one of which suggests KDM4C (JMJD2C) as a novel candidate gene for mental retardation.

A Potential Contributory Role for Ciliary Dysfunction in the 16p11.2 600 kb BP4-BP5 Pathology

The 16p11.2 600 kb copy-number variants (CNVs) are associated with mirror phenotypes on BMI, head circumference, and brain volume and represent frequent genetic lesions in autism spectrum disorders (ASDs) and schizophrenia. Here we interrogated the transcriptome of individuals carrying reciprocal 16p11.2 CNVs. Transcript perturbations correlated with clinical endophenotypes and were enriched for genes associated with ASDs, abnormalities of head size, and ciliopathies. Ciliary gene expression was also perturbed in orthologous mouse models, raising the possibility that ciliary dysfunction contributes to 16p11.2 pathologies. In support of this hypothesis, we found structural ciliary defects in the CA1 hippocampal region of 16p11.2 duplication mice. Moreover, by using an established zebrafish model, we show genetic interaction between KCTD13, a key driver of the mirrored neuroanatomical phenotypes of the 16p11.2 CNV, and ciliopathy-associated genes. Overexpression of BBS7 rescues head size and neuroanatomical defects of kctd13 morphants, whereas suppression or overexpression of CEP290 rescues phenotypes induced by KCTD13 under- or overexpression, respectively. Our data suggest that dysregulation of ciliopathy genes contributes to the clinical phenotypes of these CNVs.

The impact of human copy number variation on gene expression

Recent years have witnessed a flurry of important technological and methodological developments in the discovery and analysis of copy number variations (CNVs), which are increasingly enabling the systematic evaluation of their impact on a broad range of phenotypes from molecular-level (intermediate) traits to higher-order clinical phenotypes. Like single nucleotide variants in the human genome, CNVs have been linked to complex traits in humans, including disease and drug response. These recent developments underscore the importance of incorporating complex forms of genetic variation into disease mapping studies and promise to transform our understanding of genome function and the genetic basis of disease. Here we review some of the findings that have emerged from transcriptome studies of CNVs facilitated by the rapid advances in -omics technologies and corresponding methodologies.

Deciphering the pathogenic consequences of chromosomal aberrations in human genetic disease

Chromosomal aberrations include translocations, deletions, duplications, inversions, aneuploidies and complex rearrangements. They underlie genetic disease in roughly 15% of patients with multiple congenital abnormalities and/or mental retardation (MCA/MR). In genetic diagnostics, the pathogenicity of chromosomal aberrations in these patients is typically assessed based on criteria such as phenotypic similarity to other patients with the same or overlapping aberration, absence in healthy individuals, de novo occurrence, and protein coding gene content. However, a thorough understanding of the molecular mechanisms that lead to MCA/MR as a result of chromosome aberrations is often lacking. Chromosome aberrations can affect one or more genes in a complex manner, such as by changing the regulation of gene expression, by disrupting exons, and by creating fusion genes. The precise delineation of breakpoints by whole-genome sequencing enables the construction of local genomic architecture and facilitates the prediction of the molecular determinants of the patient's phenotype. Here, we review current methods for breakpoint identification and their impact on the interpretation of chromosome aberrations in patients with MCA/MR. In addition, we discuss opportunities to dissect disease mechanisms based on large-scale genomic technologies and studies in model organisms.

The effect of homozygous deletion of the BBOX1 and Fibin genes on carnitine level and acyl carnitine profile

Background

Carnitine is a key molecule in energy metabolism that helps transport activated fatty acids into the mitochondria. Its homeostasis is achieved through oral intake, renal reabsorption and de novo biosynthesis. Unlike dietary intake and renal reabsorption, the importance of de novo biosynthesis pathway in carnitine homeostasis remains unclear, due to lack of animal models and description of a single patient defective in this pathway.

Case presentation

We identified by array comparative genomic hybridization a 42 months-old girl homozygote for a 221 Kb interstitial deletions at 11p14.2, that overlaps the genes encoding Fibin and butyrobetaine-gamma 2-oxoglutarate dioxygenase 1 (BBOX1), an enzyme essential for the biosynthesis of carnitine de novo. She presented microcephaly, speech delay, growth retardation and minor facial anomalies. The levels of almost all evaluated metabolites were normal. Her serum level of free carnitine was at the lower limit of the reference range, while her acylcarnitine to free carnitine ratio was normal.

Conclusions

We present an individual with a completely defective carnitine de novo biosynthesis. This condition results in mildly decreased free carnitine level, but not in clinical manifestations characteristic of carnitine deficiency disorders, suggesting that dietary carnitine intake and renal reabsorption are sufficient to carnitine homeostasis. Our results also demonstrate that haploinsufficiency of BBOX1 and/or Fibin is not associated with Primrose syndrome as previously suggested.

Novel H3K4me3 marks are enriched at human- and chimpanzee-specific cytogenetic structures

Human and chimpanzee genomes are 98.8% identical within comparable sequences. However, they differ structurally in nine pericentric inversions, one fusion that originated human chromosome 2, and content and localization of heterochromatin and lineage-specific segmental duplications. The possible functional consequences of these cytogenetic and structural differences are not fully understood and their possible involvement in speciation remains unclear. We show that subtelomeric regions—regions that have a species-specific organization, are more divergent in sequence, and are enriched in genes and recombination hotspots—are significantly enriched for species-specific histone modifications that decorate transcription start sites in different tissues in both human and chimpanzee. The human lineage-specific chromosome 2 fusion point and ancestral centromere locus as well as chromosome 1 and 18 pericentric inversion breakpoints showed enrichment of human-specific H3K4me3 peaks in the prefrontal cortex. Our results reveal an association between plastic regions and potential novel regulatory elements.

The impact of chromosomal rearrangements on regulation of gene expression

The effects that coding region single-nucleotide polymorphisms or mutations have on gene expression have been well documented, predominantly owing to their association with disease. The effects of structural chromosomal rearrangements are also receiving increasing attention with the development of new techniques that allow accurate, high-resolution data, whether genomic interaction or transcriptome data, to be generated right down to the single-cell level. Over the past 18 months, these advances in experimental techniques have been used to further confirm and delineate the substantial effects that chromosome rearrangements can have on the regulation of gene expression and provide evidence of direct links between the two.

Delineation of candidate genes responsible for structural brain abnormalities in patients with terminal deletions of chromosome 6q27

Patients with terminal deletions of chromosome 6q present with structural brain abnormalities including agenesis of corpus callosum, hydrocephalus, periventricular nodular heterotopia, and cerebellar malformations. The 6q27 region harbors genes that are important for the normal development of brain and delineation of a critical deletion region for structural brain abnormalities may lead to a better genotype-phenotype correlation. We conducted a detailed clinical and molecular characterization of seven unrelated patients with deletions involving chromosome 6q27. All patients had structural brain abnormalities. Using array comparative genomic hybridization, we mapped the size, extent, and genomic content of these deletions. The smallest region of overlap spans 1.7 Mb and contains DLL1, THBS2, PHF10, and C6orf70 (ERMARD) that are plausible candidates for the causation of structural brain abnormalities. Our study reiterates the importance of 6q27 region in normal development of brain and helps identify putative genes in causation of structural brain anomalies.

NADPH-cytochrome P450 reductase is regulated by all-trans retinoic acid and by 1,25-dihydroxyvitamin D3 in human acute myeloid leukemia cells

Acute myeloid leukemia (AML) cell lines can be driven to differentiate to monocyte-like cells by 1,25- dihydroxyvitamin D₃ (1,25D) and to granulocyte-like cells by all-trans-retinoic acid (ATRA). Both compounds activate their specific intracellular receptors, vitamin D receptor (VDR) and retinoic acid receptors (RARs) respectively. Inside the cells 1,25D is degraded to calcitrioic acid by a mitochondrial enzyme CYP24A1, while ATRA is degraded to several polar metabolites by CYP26. NADPH-cytochrome P450 oxidoreductase (POR) is a membrane-bound enzyme required for electron transfer to cytochrome P450 (CYP), vital in the processes of the metabolism of drugs and steroid production in humans. In this paper we report that POR in AML cells, from both cell lines and patients, is upregulated by ATRA and by 1,25D at the level of mRNA and protein. Partial silencing of POR in HL60 cells resulted in augmented differentiation response to 1,25D.

Associations of MYH3 gene copy number variations with transcriptional expression and growth traits in Chinese cattle

Copy number variations (CNVs) have been recently recognized as another important genetic variability complementary to single nucleotide polymorphisms (SNPs). Compelling evidence has indicated that CNVs are responsible for phenotypic traits by changing the copy numbers of functional genes. Myosin heavy chain 3 (MYH3) gene is a critical regulatory factor in skeletal muscle development, and has been detected in the CNVs region by comparative genomic hybridization (CGH) array. This study was conducted to validate and detect the distribution of MYH3 copy numbers (relative to Angus cattle) in four Chinese cattle breeds (NY, QC, LX, and CY), and further to investigate the associations of the copy number changes with its transcriptional expression and cattle growth traits. Substantial genetic differences of MYH3 copy numbers were identified between NY and the other three breeds (P<0.01). The copy numbers of MYH3 gene presented the positive correlations with the transcript level of MYH3 gene in both fetal and adult skeletal muscles (P<0.05). Statistical analysis revealed that CNVs of MYH3 gene were significantly associated with growth traits of NY cattle, and the individuals with copy number gain showed better phenotypes than the loss and/or median groups (P<0.05). This study firstly attempted to establish the correlations between CNVs of candidate genes and growth traits, and our results suggested that the CNVs of MYH3 gene may be utilized as the potential markers for economic traits in selection breeding programs of Chinese cattle.

Structural variation-associated expression changes are paralleled by chromatin architecture modifications

Copy number variants (CNVs) influence the expression of genes that map not only within the rearrangement, but also to its flanks. To assess the possible mechanism(s) underlying this “neighboring effect”, we compared intrachromosomal interactions and histone modifications in cell lines of patients affected by genomic disorders and control individuals. Using chromosome conformation capture (4C-seq), we observed that a set of genes flanking the Williams-Beuren Syndrome critical region (WBSCR) were often looping together. The newly identified interacting genes include AUTS2, mutations of which are associated with autism and intellectual disabilities. Deletion of the WBSCR disrupts the expression of this group of flanking genes, as well as long-range interactions between them and the rearranged interval. We also pinpointed concomitant changes in histone modifications between samples.

We conclude that large genomic rearrangements can lead to chromatin conformation changes that extend far away from the structural variant, thereby possibly modulating expression globally and modifying the phenotype.

GEO Series accession number: GSE33784, GSE33867.

Direct visualization of the highly polymorphic RNU2 locus in proximity to the BRCA1 gene

Although the breast cancer susceptibility gene BRCA1 is one of the most extensively characterized genetic loci, much less is known about its upstream variable number tandem repeat element, the RNU2 locus. RNU2 encodes the U2 small nuclear RNA, an essential splicing element, but this locus is missing from the human genome assembly due to the inherent difficulty in the assembly of repetitive sequences. To fill the gap between RNU2 and BRCA1, we have reconstructed the physical map of this region by re-examining genomic clone sequences of public databases, which allowed us to precisely localize the RNU2 array 124 kb telomeric to BRCA1. We measured by performing FISH analyses on combed DNA for the first time the exact number of repeats carried by each of the two alleles in 41 individuals and found a range of 6-82 copies and a level of heterozygosity of 98%. The precise localisation of the RNU2 locus in the genome reference assembly and the implementation of a new technical tool to study it will make the detailed exploration of this locus possible. This recently neglected macrosatellite could be valuable for evaluating the potential role of structural variations in disease due to its location next to a major cancer susceptibility gene.

The human WBSCR22 protein is involved in the biogenesis of the 40S ribosomal subunits in mammalian cells

The human WBSCR22 protein was previously shown to be up-regulated in invasive breast cancer and its ectopic expression enhances tumor cell survival in the vasculature. In the current study, we show that the WBSCR22 protein is important for cell growth. Knock-down of WBSCR22 with siRNA results in slower growth of WBSCR22-depleted cells. Treatment with siWBSCR22 causes defects in the processing of pre-rRNAs and reduces the level of free 40S ribosomal subunit, suggesting that WBSCR22 is involved in ribosome small subunit biosynthesis. The human WBSCR22 partially complements the growth of WBSCR22 yeast homologue, bud23 deletion mutant suggesting that the human WBSCR22 is a functional homologue of yeast Bud23. WBSCR22 is localized throughout the cell nucleus and is not stably associated with ribosomal subunits within the cell nucleus. We also show that the WBSCR22 protein level is decreased in lymphoblastoid cell lines derived from William-Beuren Syndrome (WBS) patients compared to healthy controls. Our data suggest that the WBSCR22 protein is a ribosome biogenesis factor involved in the biosynthesis of 40S ribosomal particles in mammalian cells.

Presence of autism, hyperserotonemia, and severe expressive language impairment in Williams-Beuren syndrome

BACKGROUND

Deletion of the Williams-Beuren syndrome (WBS) critical region (WBSCR), at 7q11.23, causes a developmental disorder commonly characterized by hypersociability and excessive talkativeness and often considered the opposite behavioral phenotype to autism. Duplication of the WBSCR leads to severe delay in expressive language. Gene-dosage effects on language development at 7q11.23 have been hypothesized.

METHODS

Molecular characterization of the WBSCR was performed by fluorescence in situ hybridization and high-resolution single-nucleotide polymorphism array in two individuals with severe autism enrolled in a genetic study of autism who showed typical WBS facial dysmorphism on systematic clinical genetic examination. The serotonin transporter promoter polymorphism (5-HTTLPR, locus SLC6A4) was genotyped. Platelet serotonin levels and urinary 6-sulfatoxymelatonin excretion were measured. Behavioral and cognitive phenotypes were examined.

RESULTS

The two patients had common WBSCR deletions between proximal and medial low copy repeat clusters, met diagnostic criteria for autism and displayed severe impairment in communication, including a total absence of expressive speech. Both patients carried the 5-HTTLPR ss genotype and exhibited platelet hyperserotonemia and low melatonin production.

CONCLUSIONS

Our observations indicate that behaviors and neurochemical phenotypes typically associated with autism can occur in patients with common WBSCR deletions. The results raise intriguing questions about phenotypic heterogeneity in WBS and regarding genetic and/or environmental factors interacting with specific genes at 7q11.23 sensitive to dosage alterations that can influence the development of social communication skills. Thus, the influence of WBSCR genes on social communication expression might be dramatically modified by other genes, such as 5-HTTLPR, known to influence the severity of social communication impairments in autism, or by environmental factors, such as hyperserotonemia, given that hyperserotonemia is found in WBS associated with autism but not in WBS without autism. In this regard, WBS provides a potentially fruitful model with which to develop integrated genetic, cognitive, behavioral and neurochemical approaches to study genotype-phenotype correlations, possible gene-environment interactions and genetic background effects. The results underscore the importance of considering careful clinical and molecular genetic examination of individuals diagnosed with autism.

Validation of microarray data in human lymphoblasts shows a role of the ubiquitin-proteasome system and NF-kB in the pathogenesis of Down syndrome

BACKGROUND

Down syndrome (DS) is a complex disorder caused by the trisomy of either the entire, or a critical region of chromosome 21 (21q22.1-22.3). Despite representing the most common cause of mental retardation, the molecular bases of the syndrome are still largely unknown.

METHODS

To better understand the pathogenesis of DS, we analyzed the genome-wide transcription profiles of lymphoblastoid cell lines (LCLs) from six DS and six euploid individuals and investigated differential gene expression and pathway deregulation associated with trisomy 21. Connectivity map and PASS-assisted exploration were used to identify compounds whose molecular signatures counteracted those of DS lymphoblasts and to predict their therapeutic potential. An experimental validation in DS LCLs and fetal fibroblasts was performed for the most deregulated GO categories, i.e. the ubiquitin mediated proteolysis and the NF-kB cascade.

RESULTS

We show, for the first time, that the level of protein ubiquitination is reduced in human DS cell lines and that proteasome activity is increased in both basal conditions and oxidative microenvironment. We also provide the first evidence that NF-kB transcription levels, a paradigm of gene expression control by ubiquitin-mediated degradation, is impaired in DS due to reduced IkB-alfa ubiquitination, increased NF-kB inhibitor (IkB-alfa) and reduced p65 nuclear fraction. Finally, the DSCR1/DYRK1A/NFAT genes were analysed. In human DS LCLs, we confirmed the presence of increased protein levels of DSCR1 and DYRK1A, and showed that the levels of the transcription factor NFATc2 were decreased in DS along with a reduction of its nuclear translocation upon induction of calcium fluxes.

CONCLUSIONS

The present work offers new perspectives to better understand the pathogenesis of DS and suggests a rationale for innovative approaches to treat some pathological conditions associated to DS.

Phenotypic impact of genomic structural variation: insights from and for human disease

Genomic structural variants have long been implicated in phenotypic diversity and human disease, but dissecting the mechanisms by which they exert their functional impact has proven elusive. Recently however, developments in high-throughput DNA sequencing and chromosomal engineering technology have facilitated the analysis of structural variants in human populations and model systems in unprecedented detail. In this Review, we describe how structural variants can affect molecular and cellular processes, leading to complex organismal phenotypes, including human disease. We further present advances in delineating disease-causing elements that are affected by structural variants, and we discuss future directions for research on the functional consequences of structural variants.

A 1.3-mb 7q11.23 atypical deletion identified in a cohort of patients with williams-beuren syndrome

Williams-Beuren syndrome is a rare multisystem neurodevelopmental disorder caused by a 1.55-1.84-Mb hemizygous deletion on chromosome 7q11.23. The classical phenotype consists of characteristic facial features, supravalvular aortic stenosis, intellectual disability, overfriendliness, and visuospatial impairment. So far, 26-28 genes have been shown to contribute to the multisystem phenotype associated with Williams-Beuren syndrome. Among them, haploinsufficiency of the ELN gene has been shown to cause the cardiovascular anomalies. Identification of patients with atypical deletions has provided valuable information for genotype-phenotype correlation, in which other genes such as LIMK1,CLIP2, GTF2IRD1, or GTF2I have been correlated with specific cognitive profiles or craniofacial features. Here, we report the clinical and molecular characteristics of a patient with an atypical deletion that does not include the GTF2I gene and only partially includes the GTF2IRD1 gene.