Chromosomal variation in lymphoblastoid cell lines

A Burden of Rare Copy Number Variants in Obsessive-Compulsive Disorder

Current genetic research on obsessive-compulsive disorder (OCD) supports contributions to risk specifically from common single nucleotide variants (SNVs), along with rare coding SNVs and small insertion-deletions (indels). The contribution to OCD risk from large, rare copy number variants (CNVs), however, has not been formally assessed at a similar scale. Here we describe an analysis of rare CNVs called from genotype array data in 2,248 deeply phenotyped OCD cases and 3,608 unaffected controls from Sweden and Norway. We found that in general cases carry an elevated burden of large (>30kb, at least 15 probes) CNVs (OR=1.12, P=1.77×10-3). The excess rate of these CNVs in cases versus controls was around 0.07 (95% CI 0.02-0.11, P=2.58×10-3). This signal was largely driven by CNVs overlapping protein-coding regions (OR=1.19, P=3.08×10-4), particularly deletions impacting loss-of-function intolerant genes (pLI>0.995, OR=4.12, P=2.54×10-5). We did not identify any specific locus where CNV burden was associated with OCD case status at genome-wide significance, but we noted non-random recurrence of CNV deletions in cases (permutation P = 2.60×10-3). In cases where sufficient clinical data were available (n=1612) we found that carriers of neurodevelopmental duplications were more likely to have comorbid autism (P<0.001), and that carriers of deletions overlapping neurodevelopmental genes had lower treatment response (P=0.02). The results demonstrate a contribution of large, rare CNVs to OCD risk, and suggest that studies of rare coding variation in OCD would have increased power to identify risk genes if this class of variation were incorporated into formal tests.

A cryptic microdeletion del(12)(p11.21p11.23) within an unbalanced translocation t(7;12)(q21.13;q23.1) implicates new candidate loci for intellectual disability and Kallmann syndrome

In a patient diagnosed with both Kallmann syndrome (KS) and intellectual disability (ID), who carried an apparently balanced translocation t(7;12)(q22;q24)dn, array comparative genomic hybridization (aCGH) disclosed a cryptic heterozygous 4.7 Mb deletion del(12)(p11.21p11.23), unrelated to the translocation breakpoint. This novel discovery prompted us to consider the possibility that the combination of KS and neurological disorder in this patient could be attributed to gene(s) within this specific deletion at 12p11.21-12p11.23, rather than disrupted or dysregulated genes at the translocation breakpoints. To further support this hypothesis, we expanded our study by screening five candidate genes at both breakpoints of the chromosomal translocation in a cohort of 48 KS patients. However, no mutations were found, thus reinforcing our supposition. In order to delve deeper into the characterization of the 12p11.21-12p11.23 region, we enlisted six additional patients with small copy number variations (CNVs) and analyzed eight individuals carrying small CNVs in this region from the DECIPHER database. Our investigation utilized a combination of complementary approaches. Firstly, we conducted a comprehensive phenotypic-genotypic comparison of reported CNV cases. Additionally, we reviewed knockout animal models that exhibit phenotypic similarities to human conditions. Moreover, we analyzed reported variants in candidate genes and explored their association with corresponding phenotypes. Lastly, we examined the interacting genes associated with these phenotypes to gain further insights. As a result, we identified a dozen candidate genes: TSPAN11 as a potential KS candidate gene, TM7SF3, STK38L, ARNTL2, ERGIC2, TMTC1, DENND5B, and ETFBKMT as candidate genes for the neurodevelopmental disorder, and INTS13, REP15, PPFIBP1, and FAR2 as candidate genes for KS with ID. Notably, the high-level expression pattern of these genes in relevant human tissues further supported their candidacy. Based on our findings, we propose that dosage alterations of these candidate genes may contribute to sexual and/or cognitive impairments observed in patients with KS and/or ID. However, the confirmation of their causal roles necessitates further identification of point mutations in these candidate genes through next-generation sequencing.

A microdeletion del(12)(p11.21p11.23) with a cryptic unbalanced translocation t(7;12)(q21.13;q23.1) implicates new candidate loci for intellectual disability and Kallmann syndrome

In an apparently balanced translocation t(7;12)(q22;q24)dn exhibiting both Kallmann syndrome (KS) and intellectual disability (ID), we detected a cryptic heterozygous 4.7 Mb del(12)(p11.21p11.23) unrelated to the translocation breakpoint. This new finding raised the possibility that KS combined with neurological disorder in this patient could be caused by gene(s) within this deletion at 12p11.21-12p11.23 instead of disrupted or dysregulated genes at the genomic breakpoints. Screening of five candidate genes at both breakpoints in 48 KS patients we recruited found no mutation, corroborating our supposition. To substantiate this hypothesis further, we recruited six additional subjects with small CNVs and analyzed eight individuals carrying small CNVs in this region from DECIPHER to dissect 12p11.21-12p11.23. We used multiple complementary approaches including a phenotypic-genotypic comparison of reported cases, a review of knockout animal models recapitulating the human phenotypes, and analyses of reported variants in the interacting genes with corresponding phenotypes. The results identified one potential KS candidate gene ( TSPAN11 ), seven candidate genes for the neurodevelopmental disorder ( TM7SF3 , STK38L , ARNTL2 , ERGIC2 , TMTC1 , DENND5B , and ETFBKMT ), and four candidate genes for KS with ID ( INTS13 , REP15 , PPFIBP1 , and FAR2 ). The high-level expression pattern in the relevant human tissues further suggested the candidacy of these genes. We propose that the dosage alterations of the candidate genes may contribute to sexual and/or cognitive impairment in patients with KS and/or ID. Further identification of point mutations through next generation sequencing will be necessary to confirm their causal roles.

Functional analysis of structural variants in single cells using Strand-seq

Somatic structural variants (SVs) are widespread in cancer, but their impact on disease evolution is understudied due to a lack of methods to directly characterize their functional consequences. We present a computational method, scNOVA, which uses Strand-seq to perform haplotype-aware integration of SV discovery and molecular phenotyping in single cells by using nucleosome occupancy to infer gene expression as a readout. Application to leukemias and cell lines identifies local effects of copy-balanced rearrangements on gene deregulation, and consequences of SVs on aberrant signaling pathways in subclones. We discovered distinct SV subclones with dysregulated Wnt signaling in a chronic lymphocytic leukemia patient. We further uncovered the consequences of subclonal chromothripsis in T cell acute lymphoblastic leukemia, which revealed c-Myb activation, enrichment of a primitive cell state and informed successful targeting of the subclone in cell culture, using a Notch inhibitor. By directly linking SVs to their functional effects, scNOVA enables systematic single-cell multiomic studies of structural variation in heterogeneous cell populations.

A high-throughput real-time PCR tissue-of-origin test to distinguish blood from lymphoblastoid cell line DNA for (epi)genomic studies

Lymphoblastoid cell lines (LCLs) derive from blood infected in vitro by Epstein–Barr virus and were used in several genetic, transcriptomic and epigenomic studies. Although few changes were shown between LCL and blood genotypes (SNPs) validating their use in genetics, more were highlighted for other genomic features and/or in their transcriptome and epigenome. This could render them less appropriate for these studies, notably when blood DNA could still be available. Here we developed a simple, high-throughput and cost-effective real-time PCR approach allowing to distinguish blood from LCL DNA samples based on the presence of EBV relative load and rearranged T-cell receptors γ and β. Our approach was able to achieve 98.5% sensitivity and 100% specificity on DNA of known origin (458 blood and 316 LCL DNA). It was further applied to 1957 DNA samples from the CEPH Aging cohort comprising DNA of uncertain origin, identifying 784 blood and 1016 LCL DNA. A subset of these DNA was further analyzed with an epigenetic clock indicating that DNA extracted from blood should be preferred to LCL for DNA methylation-based age prediction analysis. Our approach could thereby be a powerful tool to ascertain the origin of DNA in old collections prior to (epi)genomic studies.

seGMM: A New Tool for Gender Determination From Massively Parallel Sequencing Data

In clinical genetic testing, checking the concordance between self-reported gender and genotype-inferred gender from genomic data is a significant quality control measure because mismatched gender due to sex chromosomal abnormalities or misregistration of clinical information can significantly affect molecular diagnosis and treatment decisions. Targeted gene sequencing (TGS) is widely recommended as a first-tier diagnostic step in clinical genetic testing. However, the existing gender-inference tools are optimized for whole genome and whole exome data and are not adequate and accurate for analyzing TGS data. In this study, we validated a new gender-inference tool, seGMM, which uses unsupervised clustering (Gaussian mixture model) to determine the gender of a sample. The seGMM tool can also identify sex chromosomal abnormalities in samples by aligning the sequencing reads from the genotype data. The seGMM tool consistently demonstrated >99% gender-inference accuracy in a publicly available 1,000-gene panel dataset from the 1,000 Genomes project, an in-house 785 hearing loss gene panel dataset of 16,387 samples, and a 187 autism risk gene panel dataset from the Autism Clinical and Genetic Resources in China (ACGC) database. The performance and accuracy of seGMM was significantly higher for the targeted gene sequencing (TGS), whole exome sequencing (WES), and whole genome sequencing (WGS) datasets compared to the other existing gender-inference tools such as PLINK, seXY, and XYalign. The results of seGMM were confirmed by the short tandem repeat analysis of the sex chromosome marker gene, amelogenin. Furthermore, our data showed that seGMM accurately identified sex chromosomal abnormalities in the samples. In conclusion, the seGMM tool shows great potential in clinical genetics by determining the sex chromosomal karyotypes of samples from massively parallel sequencing data with high accuracy.

Microsatellite Markers in Biobanking: A New Multiplexed Assay

Microsatellites, or MSATs, offer a fast and cost-effective way for biobanks to establish a biospecimen genetic profile. Importantly, this genetic profile can be used to authenticate multiple submissions derived from the same individual as well as biospecimens derived from the same original sample submission over time. While the Certificate of Confidentiality provided by the National Institutes of Health offers some meaningful protection to prevent the disclosure of potentially identifiable information to entities within the United States, we consider, in this study, the potential to offer additional protection to participants who choose to donate to biobanks by minimizing the use of forensic Combined DNA Index System (CODIS) MSAT markers in biobanking. To this end, we report the design and validation of a new multiplexed MSAT assay that does not include CODIS markers for use in biobanking operations and quality control management.

The Master Regulator Protein BAZ2B Can Reprogram Human Hematopoietic Lineage-Committed Progenitors into a Multipotent State

Bi-species, fusion-mediated, somatic cell reprogramming allows precise, organism-specific tracking of unknown lineage drivers. The fusion of Tcf7l1^-/- murine embryonic stem cells with EBV-transformed human B cell lymphocytes, leads to the generation of bi-species heterokaryons. Human mRNA transcript profiling at multiple time points permits the tracking of the reprogramming of B cell nuclei to a multipotent state. Interrogation of a human B cell regulatory network with gene expression signatures identifies 8 candidate master regulator proteins. Of these 8 candidates, ectopic expression of BAZ2B, from the bromodomain family, efficiently reprograms hematopoietic committed progenitors into a multipotent state and significantly enhances their long-term clonogenicity, stemness, and engraftment in immunocompromised mice. Unbiased systems biology approaches let us identify the early driving events of human B cell reprogramming.

Integrated Genomic, Epigenomic, and Expression Analyses of Ovarian Cancer Cell Lines

To improve our understanding of ovarian cancer, we performed genome-wide analyses of 45 ovarian cancer cell lines. Given the challenges of genomic analyses of tumors without matched normal samples, we developed approaches for detection of somatic sequence and structural changes and integrated these with epigenetic and expression alterations. Alterations not previously implicated in ovarian cancer included amplification or overexpression of ASXL1 and H3F3B, deletion or underexpression of CDC73 and TGF-beta receptor pathway members, and rearrangements of YAP1-MAML2 and IKZF2-ERBB4. Dose-response analyses to targeted therapies revealed unique molecular dependencies, including increased sensitivity of tumors with PIK3CA and PPP2R1A alterations to PI3K inhibitor GNE-493, MYC amplifications to PARP inhibitor BMN673, and SMAD3/4 alterations to MEK inhibitor MEK162. Genome-wide rearrangements provided an improved measure of sensitivity to PARP inhibition. This study provides a comprehensive and broadly accessible resource of molecular information for the development of therapeutic avenues in ovarian cancer.

The overall survival of patients with late-stage ovarian cancer is dismal. To identify therapeutic opportunities, Papp et al. integrate genomic, epigenomic, and expression analyses to provide a resource of molecular abnormalities in ovarian cancer cell lines and use these to identify tumors sensitive to PARP, MEK, and PI3K inhibitors.

DNA copy number variations in children with vesicoureteral reflux and urinary tract infections

Vesicoureteral reflux (VUR) is a complex, heritable disorder. Genome-wide linkage analyses of families affected by VUR have revealed multiple genomic loci linked to VUR. These loci normally harbor a number of genes whose biologically functional variant is yet to be identified. DNA copy number variations (CNVs) have not been extensively studied at high resolution in VUR patients. In this study, we performed array comparative genomic hybridization (aCGH) on a cohort of patients with a history of both VUR and urinary tract infection (UTI) with the objective of identifying genetic variations responsible for VUR and/or UTI susceptibility. UTI/VUR-associated CNVs were identified by aCGH results from the 192 Randomized Intervention for Children With Vesicoureteral Reflux (RIVUR) patients compared to 683 controls. Rare, large CNVs that are likely pathogenic and lead to VUR development were identified using stringent analysis criteria. Because UTI is a common affliction with multiple risk factors, we utilized standard analysis to identify potential disease-modifying CNVs that can contribute to UTI risk. Gene ontology analysis identified that CNVs in innate immunity and development genes were enriched in RIVUR patients. CNVs affecting innate immune genes may contribute to UTI susceptibility in VUR patients and may provide the first step in assisting clinical medicine in determining adverse outcome risk in children with VUR.

Congenital eye anomalies: More mosaic than thought?

The eye is a sensory organ that primarily captures light and provides the sense of sight, as well as delivering non-visual light information involving biological rhythms and neurophysiological activities to the brain. Since the early 1990s, rapid advances in molecular biology have enabled the identification of developmental genes, genes responsible for human congenital diseases, and relevant genes of mutant animals with various anomalies. In this review, we first look at the development of the eye, and we highlight seminal reports regarding archetypal gene defects underlying three developmental ocular disorders in humans: (1) holoprosencephaly (HPE), with cyclopia being exhibited in the most severe cases; (2) microphthalmia, anophthalmia, and coloboma (MAC) phenotypes; and (3) anterior segment dysgenesis (ASDG), known as Peters anomaly and its related disorders. The recently developed methods, such as next-generation sequencing and genome editing techniques, have aided the discovery of gene mutations in congenital eye diseases and gene functions in normal eye development. Finally, we discuss Pax6-genome edited mosaic eyes and propose that somatic mosaicism in developmental gene mutations should be considered a causal factor for variable phenotypes, sporadic cases, and de novo mutations in human developmental disorders.

Extensive epigenetic and transcriptomic variability between genetically identical human B-lymphoblastoid cells with implications in pharmacogenomics research

Genotyped human B-lymphoblastoid cell lines (LCLs) are widely used models in mapping quantitative trait loci for chromatin features, gene expression, and drug response. The extent of genotype-independent functional genomic variability of the LCL model, although largely overlooked, may inform association study design. In this study, we use flow cytometry, chromatin immunoprecipitation sequencing and mRNA sequencing to study surface marker patterns, quantify genome-wide chromatin changes (H3K27ac) and transcriptome variability, respectively, among five isogenic LCLs derived from the same individual. Most of the studied LCLs were non-monoclonal and had mature B cell phenotypes. Strikingly, nearly one-fourth of active gene regulatory regions showed significantly variable H3K27ac levels, especially enhancers, among which several were classified as clustered enhancers. Large, contiguous genomic regions showed signs of coordinated activity change. Regulatory differences were mirrored by mRNA expression changes, preferentially affecting hundreds of genes involved in specialized cellular processes including immune and drug response pathways. Differential expression of DPYD, an enzyme involved in 5-fluorouracil (5-FU) catabolism, was associated with variable LCL growth inhibition mediated by 5-FU. The extent of genotype-independent functional genomic variability might highlight the need to revisit study design strategies for LCLs in pharmacogenomics.

The Malaria-Protective Human Glycophorin Structural Variant DUP4 Shows Somatic Mosaicism and Association with Hemoglobin Levels

Glycophorin A and glycophorin B are red blood cell surface proteins and are both receptors for the parasite Plasmodium falciparum, which is the principal cause of malaria in sub-Saharan Africa. DUP4 is a complex structural genomic variant that carries extra copies of a glycophorin A-glycophorin B fusion gene and has a dramatic effect on malaria risk by reducing the risk of severe malaria by up to 40%. Using fiber-FISH and Illumina sequencing, we validate the structural arrangement of the glycophorin locus in the DUP4 variant and reveal somatic variation in copy number of the glycophorin B-glycophorin A fusion gene. By developing a simple, specific, PCR-based assay for DUP4, we show that the DUP4 variant reaches a frequency of 13% in the population of a malaria-endemic village in south-eastern Tanzania. We genotype a substantial proportion of that village and demonstrate an association of DUP4 genotype with hemoglobin levels, a phenotype related to malaria, using a family-based association test. Taken together, we show that DUP4 is a complex structural variant that may be susceptible to somatic variation and show that DUP4 is associated with a malarial-related phenotype in a longitudinally followed population.

Genetic and genomic stability across lymphoblastoid cell line expansions

Objective

Lymphoblastoid cell lines are widely used in genetic and genomic studies. Previous work has characterized variant stability in transformed culture and across culture passages. Our objective was to extend this work to evaluate single nucleotide polymorphism and structural variation across cell line expansions, which are commonly used in biorepository distribution. Our study used DNA and cell lines sampled from six research participants. We assayed genome-wide genetic variants and inferred structural variants for DNA extracted from blood, from transformed cell cultures, and from three generations of expansions.

Results

Single nucleotide variation was stable between DNA and expanded cell lines (ranging from 99.90 to 99.98% concordance). Structural variation was less consistent across expansions (median 33% concordance) with a noticeable decrease in later expansions. In summary, we demonstrate consistency between SNPs assayed from whole blood DNA and LCL DNA; however, more caution should be taken in using LCL DNA to study structural variation.

Electronic supplementary material

The online version of this article (10.1186/s13104-018-3664-3) contains supplementary material, which is available to authorized users.

Discovery biology of neuropsychiatric syndromes (DBNS): a center for integrating clinical medicine and basic science

BACKGROUND

There is emerging evidence that there are shared genetic, environmental and developmental risk factors in psychiatry, that cut across traditional diagnostic boundaries. With this background, the Discovery biology of neuropsychiatric syndromes (DBNS) proposes to recruit patients from five different syndromes (schizophrenia, bipolar disorder, obsessive compulsive disorder, Alzheimer's dementia and substance use disorders), identify those with multiple affected relatives, and invite these families to participate in this study. The families will be assessed: 1) To compare neuro-endophenotype measures between patients, first degree relatives (FDR) and healthy controls., 2) To identify cellular phenotypes which differentiate the groups., 3) To examine the longitudinal course of neuro-endophenotype measures., 4) To identify measures which correlate with outcome, and 5) To create a unified digital database and biorepository.

METHODS

The identification of the index participants will occur at well-established specialty clinics. The selected individuals will have a strong family history (with at least another affected FDR) of mental illness. We will also recruit healthy controls without family history of such illness. All recruited individuals (N = 4500) will undergo brief clinical assessments and a blood sample will be drawn for isolation of DNA and peripheral blood mononuclear cells (PBMCs). From among this set, a subset of 1500 individuals (300 families and 300 controls) will be assessed on several additional assessments [detailed clinical assessments, endophenotype measures (neuroimaging- structural and functional, neuropsychology, psychophysics-electroencephalography, functional near infrared spectroscopy, eye movement tracking)], with the intention of conducting repeated measurements every alternate year. PBMCs from this set will be used to generate lymphoblastoid cell lines, and a subset of these would be converted to induced pluripotent stem cell lines and also undergo whole exome sequencing.

DISCUSSION

We hope to identify unique and overlapping brain endophenotypes for major psychiatric syndromes. In a proportion of subjects, we expect these neuro-endophenotypes to progress over time and to predict treatment outcome. Similarly, cellular assays could differentiate cell lines derived from such groups. The repository of biomaterials as well as digital datasets of clinical parameters, will serve as a valuable resource for the broader scientific community who wish to address research questions in the area.

Free-living human cells reconfigure their chromosomes in the evolution back to uni-cellularity

Cells of multi-cellular organisms evolve toward uni-cellularity in the form of cancer and, if humans intervene, continue to evolve in cell culture. During this process, gene dosage relationships may evolve in novel ways to cope with the new environment and may regress back to the ancestral uni-cellular state. In this context, the evolution of sex chromosomes vis-a-vis autosomes is of particular interest. Here, we report the chromosomal evolution in ~ 600 cancer cell lines. Many of them jettisoned either Y or the inactive X; thus, free-living male and female cells converge by becoming ‘de-sexualized’. Surprisingly, the active X often doubled, accompanied by the addition of one haploid complement of autosomes, leading to an X:A ratio of 2:3 from the extant ratio of 1:2. Theoretical modeling of the frequency distribution of X:A karyotypes suggests that the 2:3 ratio confers a higher fitness and may reflect aspects of sex chromosome evolution.

Multicellular life relies on a group of cells working together for a common interest. To study these cells, researchers take them out of the organism and grow them in the laboratory. Instead of growing as part of organs and tissues, the cells normally have a free-living lifestyle. Because multicellular life evolved from single-celled organisms, laboratory-grown cells can be considered as life forms that are evolving backward from a multicellular to a single-celled existence.

Normally, the cells that make up most of the tissues in the human body have 22 pairs of chromosomes known as autosomes and a pair of sex chromosomes. The cells of women have two X sex chromosomes, one of which is inactive, while those of men have one X and one Y chromosome. However, free-living single cells do not need to distinguish between male and female cells.

Xu, Peng, Chen et al. have now studied the chromosomes of cancer cells taken from over 600 people and grown in the laboratory. As the cells evolved in response to their free-living lifestyle, they became ‘de-sexualized’; male cells lost their Y chromosome, while female cells abandoned their inactive X chromosome. The cells then evolved toward a new state in which they possessed two active X chromosomes and three sets of autosomes. This new configuration suggests that the current X chromosome to autosome ratio may not be optimal for fitness and hence sheds some light on how mammalian sex chromosomes evolved.

It is currently thought that as cancerous tumors grow, their cells evolve to favor their own interests over the common interests of the rest of the organism. In this way, they develop characteristics more like those of single cells. Further research is therefore needed to investigate whether changes occur to the chromosomes of cancer cells growing within the body, and whether this gives them an advantage over normal cells.

Identification of balanced chromosomal rearrangements previously unknown among participants in the 1000 Genomes Project: implications for interpretation of structural variation in genomes and the future of clinical cytogenetics

Purpose

Recent studies demonstrate that whole-genome sequencing (WGS) enables detection of cryptic rearrangements in apparently balanced chromosomal rearrangements (also known as balanced chromosomal abnormalities, BCAs) previously identified by conventional cytogenetic methods. We aimed to assess our analytical tool for detecting BCAs in The 1000 Genomes Project without knowing affected bands.

Methods

The 1000 Genomes Project provides an unprecedented integrated map of structural variants in phenotypically normal subjects, but there is no information on potential inclusion of subjects with apparently BCAs akin to those traditionally detected in diagnostic cytogenetics laboratories. We applied our analytical tool to 1,166 genomes from the 1000 Genomes Project with sufficient physical coverage (8.25-fold).

Results

Our approach detected four reciprocal balanced translocations and four inversions ranging in size from 57.9 kb to 13.3 Mb, all of which were confirmed by cytogenetic methods and PCR studies. One of DNAs has a subtle translocation that is not readily identified by chromosome analysis due to similar banding patterns and size of exchanged segments, and another results in disruption of all transcripts of an OMIM gene.

Conclusions

Our study demonstrates the extension of utilizing low-coverage WGS for unbiased detection of BCAs including translocations and inversions previously unknown in the 1000 Genomes Project.

Integrated Genomic Medicine: A Paradigm for Rare Diseases and Beyond

Individualized medicine, or the tailoring of therapeutic interventions to a patient's unique genetic, biochemical, physiological, exposure and behavioral profile, has been enhanced, if not enabled, by modern biomedical technologies such as high-throughput DNA sequencing platforms, induced pluripotent stem cell assays, biomarker discovery protocols, imaging modalities, and wireless monitoring devices. Despite successes in the isolated use of these technologies, however, it is arguable that their combined and integrated use in focused studies of individual patients is the best way to not only tailor interventions for those patients, but also shed light on treatment strategies for patients with similar conditions. This is particularly true for individuals with rare diseases since, by definition, they will require study without recourse to other individuals, or at least without recourse to many other individuals. Such integration and focus will require new biomedical scientific paradigms and infrastructure, including the creation of databases harboring study results, the formation of dedicated multidisciplinary research teams and new training programs. We consider the motivation and potential for such integration, point out areas in need of improvement, and argue for greater emphasis on improving patient health via technological innovations, not merely improving the technologies themselves. We also argue that the paradigm described can, in theory, be extended to the study of individuals with more common diseases.

Rare Copy Number Variants in NRXN1 and CNTN6 Increase Risk for Tourette Syndrome

Tourette syndrome (TS) is a model neuropsychiatric disorder thought to arise from abnormal development and/or maintenance of cortico-striato-thalamo-cortical circuits. TS is highly heritable, but its underlying genetic causes are still elusive, and no genome-wide significant loci have been discovered to date. We analyzed a European ancestry sample of 2,434 TS cases and 4,093 ancestry-matched controls for rare (< 1% frequency) copy-number variants (CNVs) using SNP microarray data. We observed an enrichment of global CNV burden that was prominent for large (> 1 Mb), singleton events (OR = 2.28, 95% CI [1.39-3.79], p = 1.2 × 10-3) and known, pathogenic CNVs (OR = 3.03 [1.85-5.07], p = 1.5 × 10-5). We also identified two individual, genome-wide significant loci, each conferring a substantial increase in TS risk (NRXN1 deletions, OR = 20.3, 95% CI [2.6-156.2]; CNTN6 duplications, OR = 10.1, 95% CI [2.3-45.4]). Approximately 1% of TS cases carry one of these CNVs, indicating that rare structural variation contributes significantly to the genetic architecture of TS.

Differences between the genomes of lymphoblastoid cell lines and blood-derived samples

Lymphoblastoid cell lines (LCLs) represent a convenient research tool for expanding the amount of biologic material available from an individual. LCLs are commonly used as reference materials, most notably from the Genome in a Bottle Consortium. However, the question remains how faithfully LCL-derived genome assemblies represent the germline genome of the donor individual as compared to the genome assemblies derived from peripheral blood mononuclear cells. We present an in-depth comparison of a large collection of LCL- and peripheral blood mononuclear cell-derived genomes in terms of distributions of coverage and copy number alterations. We found significant differences in the depth of coverage and copy number calls, which may be driven by differential replication timing. Importantly, these copy number changes preferentially affect regions closer to genes and with higher GC content. This suggests that genomic studies based on LCLs may display locus-specific biases, and that conclusions based on analysis of depth of coverage and copy number variation may require further scrutiny.

Unrevealed mosaicism in the next-generation sequencing era

Mosaicism refers to the presence in an individual of normal and abnormal cells that are genotypically distinct and are derived from a single zygote. The incidence of mosaicism events in the human body is underestimated as the genotypes in the mosaic ratio, especially in the low-grade mosaicism, stay unrevealed. This review summarizes various research outcomes and diagnostic questions in relation to different types of mosaicism. The impact of both tested biological material and applied method on the mosaicism detection rate is especially highlighted. As next-generation sequencing technologies constitute a promising methodological solution in mosaicism detection in the coming years, revisions in current diagnostic protocols are necessary to increase the detection rate of the unrevealed mosaicism events. Since mosaicism identification is a complex process, numerous examples of multistep mosaicism investigations are presented and discussed.

Transcriptome outlier analysis implicates schizophrenia susceptibility genes and enriches putatively functional rare genetic variants

We searched a gene expression dataset comprised of 634 schizophrenia (SZ) cases and 713 controls for expression outliers (i.e., extreme tails of the distribution of transcript expression values) with SZ cases overrepresented compared with controls. These outlier genes were enriched for brain expression and for genes known to be associated with neurodevelopmental disorders. SZ cases showed higher outlier burden (i.e., total outlier events per subject) than controls for genes within copy number variants (CNVs) associated with SZ or neurodevelopmental disorders. Outlier genes were enriched for CNVs and for rare putative regulatory variants, but this only explained a small proportion of the outlier subjects, highlighting the underlying presence of additional genetic and potentially, epigenetic mechanisms.

Copy number variants encompassing Mendelian disease genes in a large multigenerational family segregating bipolar disorder

BACKGROUND

Bipolar affective disorder (BP) is a common, highly heritable psychiatric disorder characterized by periods of depression and mania. Using dense SNP genotype data, we characterized CNVs in 388 members of an Old Order Amish Pedigree with bipolar disorder. We identified CNV regions arising from common ancestral mutations by utilizing the pedigree information. By combining this analysis with whole genome sequence data in the same individuals, we also explored the role of compound heterozygosity.

RESULTS

Here we describe 541 inherited CNV regions, of which 268 are rare in a control population of European origin but present in a large number of Amish individuals. In addition, we highlight a set of CNVs found at higher frequencies in BP individuals, and within genes known to play a role in human development and disease. As in prior reports, we find no evidence for an increased burden of CNVs in BP individuals, but we report a trend towards a higher burden of CNVs in known Mendelian disease loci in bipolar individuals (BPI and BPII, p = 0.06).

CONCLUSIONS

We conclude that CNVs may be contributing factors in the phenotypic presentation of mood disorders and co-morbid medical conditions in this family. These results reinforce the hypothesis of a complex genetic architecture underlying BP disorder, and suggest that the role of CNVs should continue to be investigated in BP data sets.

Upregulation of TFAM and mitochondria copy number in human lymphoblastoid cells

Mitochondria are central to several physiological and pathological conditions in humans. In the present study, we performed copy number analysis of nuclear encoded mitochondrial genes, in peripheral blood mononuclear cells (PBMCs) and its representative lymphoblastoid cells (LCLs). We have observed hyper diploid copies of mitochondrial transcription factor A (TFAM) gene in the LCLs along with increased mtDNA copy number, mitochondrial mass, intracellular ROS and mitochondrial membrane potential, suggesting elevated mitochondrial biogenesis in LCLs. Gene expression analysis confirmed TFAM over-expression in LCLs when compared to PBMC. Based on our observation, we suggest that increased copy number of TFAM gene upregulates its expression, increases mtDNA copy numbers and protects it from oxidative stress induced damage in the transformed LCLs.

Who are the Okinawans? Ancestry, genome diversity, and implications for the genetic study of human longevity from a geographically isolated population

Isolated populations have advantages for genetic studies of longevity from decreased haplotype diversity and long-range linkage disequilibrium. This permits smaller sample sizes without loss of power, among other utilities. Little is known about the genome of the Okinawans, a potential population isolate, recognized for longevity. Therefore, we assessed genetic diversity, structure, and admixture in Okinawans, and compared this with Caucasians, Chinese, Japanese, and Africans from HapMap II, genotyped on the same Affymetrix GeneChip Human Mapping 500K array. Principal component analysis, haplotype coverage, and linkage disequilibrium decay revealed a distinct Okinawan genome-more homogeneity, less haplotype diversity, and longer range linkage disequilibrium. Population structure and admixture analyses utilizing 52 global reference populations from the Human Genome Diversity Cell Line Panel demonstrated that Okinawans clustered almost exclusively with East Asians. Sibling relative risk (λs) analysis revealed that siblings of Okinawan centenarians have 3.11 times (females) and 3.77 times (males) more likelihood of centenarianism. These findings suggest that Okinawans are genetically distinct and share several characteristics of a population isolate, which are prone to develop extreme phenotypes (eg, longevity) from genetic drift, natural selection, and population bottlenecks. These data support further exploration of genetic influence on longevity in the Okinawans.

Vitamin D receptor ChIP-seq in primary CD4+ cells: relationship to serum 25-hydroxyvitamin D levels and autoimmune disease

BACKGROUND

Vitamin D insufficiency has been implicated in autoimmunity. ChIP-seq experiments using immune cell lines have shown that vitamin D receptor (VDR) binding sites are enriched near regions of the genome associated with autoimmune diseases. We aimed to investigate VDR binding in primary CD4+ cells from healthy volunteers.

METHODS

We extracted CD4+ cells from nine healthy volunteers. Each sample underwent VDR ChIP-seq. Our results were analyzed in relation to published ChIP-seq and RNA-seq data in the Genomic HyperBrowser. We used MEMEChIP for de novo motif discovery. 25-Hydroxyvitamin D levels were measured using liquid chromatography-tandem mass spectrometry and samples were divided into vitamin D sufficient (25(OH)D ≥75 nmol/L) and insufficient/deficient (25(OH)D <75 nmol/L) groups.

RESULTS

We found that the amount of VDR binding is correlated with the serum level of 25-hydroxyvitamin D (r = 0.92, P= 0.0005). In vivo VDR binding sites are enriched for autoimmune disease associated loci, especially when 25-hydroxyvitamin D levels (25(OH)D) were sufficient (25(OH)D ≥75: 3.13-fold, P<0.0001; 25(OH)D <75: 2.76-fold, P<0.0001; 25(OH)D ≥75 enrichment versus 25(OH)D <75 enrichment: P= 0.0002). VDR binding was also enriched near genes associated specifically with T-regulatory and T-helper cells in the 25(OH)D ≥75 group. MEME ChIP did not identify any VDR-like motifs underlying our VDR ChIP-seq peaks.

CONCLUSION

Our results show a direct correlation between in vivo 25-hydroxyvitamin D levels and the number of VDR binding sites, although our sample size is relatively small. Our study further implicates VDR binding as important in gene-environment interactions underlying the development of autoimmunity and provides a biological rationale for 25-hydroxyvitamin D sufficiency being based at 75 nmol/L. Our results also suggest that VDR binding in response to physiological levels of vitamin D occurs predominantly in a VDR motif-independent manner.

Sensitive and specific detection of mosaic chromosomal abnormalities using the Parent-of-Origin-based Detection (POD) method

Background

Mosaic somatic alterations are present in all multi-cellular organisms, but the physiological effects of low-level mosaicism are largely unknown. Most mosaic alterations remain undetectable with current analytical approaches, although the presence of such alterations is increasingly implicated as causative for disease.

Results

Here, we present the Parent-of-Origin-based Detection (POD) method for chromosomal abnormality detection in trio-based SNP microarray data. Our software implementation, triPOD, was benchmarked using a simulated dataset, outperformed comparable software for sensitivity of abnormality detection, and displayed substantial improvement in the detection of low-level mosaicism while maintaining comparable specificity. Examples of low-level mosaic abnormalities from a large autism dataset demonstrate the benefits of the increased sensitivity provided by triPOD. The triPOD analyses showed robustness across multiple types of Illumina microarray chips. Two large, clinically-relevant datasets were characterized and compared.

Conclusions

Our method and software provide a significant advancement in the ability to detect low-level mosaic abnormalities, thereby opening new avenues for research into the implications of mosaicism in pathogenic and non-pathogenic processes.