Search

The adoption of agriculture triggered a rapid shift towards starch-rich diets in human populations1. Amylase genes facilitate starch digestion, and increased amylase copy number has been observed in some modern human populations with high-starch intake2, although evidence of recent selection is lacking3,4. Here, using 94 long-read haplotype-resolved assemblies and short-read data from approximately 5,600 contemporary and ancient humans, we resolve the diversity and evolutionary history of structural variation at the amylase locus. We find that amylase genes have higher copy numbers in agricultural populations than in fishing, hunting and pastoral populations. We identify 28 distinct amylase structural architectures and demonstrate that nearly identical structures have arisen recurrently on different haplotype backgrounds throughout recent human history. AMY1 and AMY2A genes each underwent multiple duplication/deletion events with mutation rates up to more than 10,000-fold the single-nucleotide polymorphism mutation rate, whereas AMY2B gene duplications share a single origin. Using a pangenome-based approach, we infer structural haplotypes across thousands of humans identifying extensively duplicated haplotypes at higher frequency in modern agricultural populations. Leveraging 533 ancient human genomes, we find that duplication-containing haplotypes (with more gene copies than the ancestral haplotype) have rapidly increased in frequency over the past 12,000 years in West Eurasians, suggestive of positive selection. Together, our study highlights the potential effects of the agricultural revolution on human genomes and the importance of structural variation in human adaptation.

Statin drugs lower blood cholesterol levels for cardiovascular disease prevention. Women are more likely than men to experience adverse statin effects, particularly new-onset diabetes (NOD) and muscle weakness. Here we find that impaired glucose homeostasis and muscle weakness in statin-treated female mice are associated with reduced levels of the omega-3 fatty acid, docosahexaenoic acid (DHA), impaired redox tone, and reduced mitochondrial respiration. Statin adverse effects are prevented in females by administering fish oil as a source of DHA, by reducing dosage of the X chromosome or the Kdm5c gene, which escapes X chromosome inactivation and is normally expressed at higher levels in females than males. As seen in female mice, we find that women experience more severe reductions than men in DHA levels after statin administration, and that DHA levels are inversely correlated with glucose levels. Furthermore, induced pluripotent stem cells from women who developed NOD exhibit impaired mitochondrial function when treated with statin, whereas cells from men do not. These studies identify X chromosome dosage as a genetic risk factor for statin adverse effects and suggest DHA supplementation as a preventive co-therapy.

The 22q11.2 locus contains genes critical for brain development. Reciprocal Copy Number Variations (CNVs) at this locus impact risk for neurodevelopmental and psychiatric disorders. Both 22q11.2 deletions (22qDel) and duplications (22qDup) are associated with autism, but 22qDel uniquely elevates schizophrenia risk. Understanding brain phenotypes associated with these highly penetrant CNVs can provide insights into genetic pathways underlying neuropsychiatric disorders. Human neuroimaging and animal models indicate subcortical brain alterations in 22qDel, yet little is known about developmental differences across specific nuclei between reciprocal 22q11.2 CNV carriers and typically developing (TD) controls. We conducted a longitudinal MRI study in a total of 385 scans from 22qDel (n = 96, scans = 191, 53.1% female), 22qDup (n = 37, scans = 64, 45.9% female), and TD controls (n = 80, scans = 130, 51.2% female), across a wide age range (5.5-49.5 years). Volumes of the thalamus, hippocampus, amygdala, and anatomical subregions were estimated using FreeSurfer, and the linear effects of 22q11.2 gene dosage and non-linear effects of age were characterized with generalized additive mixed models (GAMMs). Positive gene dosage effects (volume increasing with copy number) were observed for total intracranial and whole hippocampus volumes, but not whole thalamus or amygdala volumes. Several amygdala subregions exhibited similar positive effects, with bi-directional effects found across thalamic nuclei. Distinct age-related trajectories were observed across the three groups. Notably, both 22qDel and 22qDup carriers exhibited flattened development of hippocampal CA2/3 subfields relative to TD controls. This study provides novel insights into the impact of 22q11.2 CNVs on subcortical brain structures and their developmental trajectories.

Loss of Cdx2 in vivo leads to stunted development of the allantois, an extraembryonic mesoderm-derived structure critical for nutrient delivery and waste removal in the early embryo. Here, we investigate how CDX2 dose-dependently influences the gene regulatory network underlying extraembryonic mesoderm development. By engineering human induced pluripotent stem cells (hiPSCs) consisting of wild-type (WT), heterozygous (CDX2-Het), and homozygous null CDX2 (CDX2-KO) genotypes, differentiating these cells in a 2D gastruloid model, and subjecting these cells to single-nucleus RNA and ATAC sequencing, we identify several pathways that are dose-dependently regulated by CDX2 including VEGF and non-canonical WNT. snATAC-seq reveals that CDX2-Het cells retain a WT-like chromatin accessibility profile, suggesting accessibility alone is not sufficient to drive this variability in gene expression. Because the loss of CDX2 or TBXT phenocopy one another in vivo, we compared differentially expressed genes in our CDX2-KO to those from TBXT-KO hiPSCs differentiated in an analogous experiment. This comparison identifies several communally misregulated genes that are critical for cytoskeletal integrity and tissue permeability. Together, these results clarify how CDX2 dose-dependently regulates gene expression in the extraembryonic mesoderm and reveal pathways that may underlie the defects in vascular development and allantoic elongation seen in vivo.

IntroductionPeople with Down syndrome (DS) are predisposed to Alzheimer's disease (AD). The amyloid hypothesis informs studies of AD. In AD-DS, but not sporadic AD, increased APP copy number is necessary, defining the APP gene dose hypothesis. Which amyloid precursor protein (APP) products contribute needs to be determined.MethodsBrain levels of full-length protein (fl-hAPP), C-terminal fragments (hCTFs), and amyloid beta (Aβ) peptides were measured in DS, AD-DS, non-demented controls (ND), and sporadic AD cases. The APP gene-dose hypothesis was evaluated in the Dp16 model.ResultsDS and AD-DS differed from ND and AD for all APP products. In AD-DS, Aβ42 and Aβ40 levels exceeded AD. APP products were increased in the Dp16 model; increased APP gene dose was necessary for loss of vulnerable neurons, tau pathology, and activation of astrocytes and microglia.DiscussionIncreases in APP products other than Aβ distinguished AD-DS from AD. Deciphering AD-DS pathogenesis necessitates deciphering which APP products contribute and how.

Sequence variants in gene regulatory regions alter gene expression and contribute to phenotypes of individual cells and the whole organism, including disease susceptibility and progression. Single-nucleotide variants in enhancers or promoters may affect gene transcription by altering transcription factor binding sites. Differential transcription factor binding in heterozygous genomic loci provides a natural source of information on such regulatory variants. We present a novel approach to call the allele-specific transcription factor binding events at single-nucleotide variants in ChIP-Seq data, taking into account the joint contribution of aneuploidy and local copy number variation, that is estimated directly from variant calls. We have conducted a meta-analysis of more than 7 thousand ChIP-Seq experiments and assembled the database of allele-specific binding events listing more than half a million entries at nearly 270 thousand single-nucleotide polymorphisms for several hundred human transcription factors and cell types. These polymorphisms are enriched for associations with phenotypes of medical relevance and often overlap eQTLs, making candidates for causality by linking variants with molecular mechanisms. Specifically, there is a special class of switching sites, where different transcription factors preferably bind alternative alleles, thus revealing allele-specific rewiring of molecular circuitry.

Haploinsufficiency of transcriptional regulators causes human congenital heart disease (CHD); however, the underlying CHD gene regulatory network (GRN) imbalances are unknown. Here, we define transcriptional consequences of reduced dosage of the CHD transcription factor, TBX5, in individual cells during cardiomyocyte differentiation from human induced pluripotent stem cells (iPSCs). We discovered highly sensitive dysregulation of TBX5-dependent pathways-including lineage decisions and genes associated with heart development, cardiomyocyte function, and CHD genetics-in discrete subpopulations of cardiomyocytes. Spatial transcriptomic mapping revealed chamber-restricted expression for many TBX5-sensitive transcripts. GRN analysis indicated that cardiac network stability, including vulnerable CHD-linked nodes, is sensitive to TBX5 dosage. A GRN-predicted genetic interaction between Tbx5 and Mef2c, manifesting as ventricular septation defects, was validated in mice. These results demonstrate exquisite and diverse sensitivity to TBX5 dosage in heterogeneous subsets of iPSC-derived cardiomyocytes and predicts candidate GRNs for human CHDs, with implications for quantitative transcriptional regulation in disease.

Tuberculosis (TB) is the largest infectious disease with 10 million new active-TB patients and1.7 million deaths per year. Active-TB is an inflammatory disease and is increasingly viewed as an imbalance of immune responses to M. tb. infection. The mechanisms of a switch from latent infection to active disease is not well worked out but a shift in the immune responses is thought to be responsible. Increasingly, the role of gut microbiota has been described as a major influencer of the immune system. And because the gut is the largest immune organ, we aimed to analyze the gut microbiome in active-TB patients in a TB-endemic country, Pakistan. The study revealed that Ruminococcacea, Enetrobactericeae, Erysipelotrichaceae, Bifidobacterium, etc. were the major genera associated with active-TB, also associated with chronic inflammatory disease. Plasma antibody profiles against several M. tb. antigens, as specific biomarkers for active-TB, correlated closely with the patient gut microbial profiles. Besides, bcoA gene copy number, indicative of the level of butyrate production by the gut microbiome was five-fold lower in TB patients compared to healthy individuals. These findings suggest that gut health in TB patients is compromised, with implications for disease morbidity (e.g., severe weight loss) as well as immune impairment.

Males and females differ in body composition and fat distribution. Using a mouse model that segregates gonadal sex (ovaries and testes) from chromosomal sex (XX and XY), we showed that XX chromosome complement in combination with a high-fat diet led to enhanced weight gain in the presence of male or female gonads. We identified the genomic dosage of Kdm5c, an X chromosome gene that escapes X chromosome inactivation, as a determinant of the X chromosome effect on adiposity. Modulating Kdm5c gene dosage in XX female mice to levels that are normally present in males resulted in reduced body weight, fat content, and food intake to a degree similar to that seen with altering the entire X chromosome dosage. In cultured preadipocytes, the levels of KDM5C histone demethylase influenced chromatin accessibility (ATAC-Seq), gene expression (RNA-Seq), and adipocyte differentiation. Both in vitro and in vivo, Kdm5c dosage influenced gene expression involved in extracellular matrix remodeling, which is critical for adipocyte differentiation and adipose tissue expansion. In humans, adipose tissue KDM5C mRNA levels and KDM5C genetic variants were associated with body mass. These studies demonstrate that the sex-dependent dosage of Kdm5c contributes to male/female differences in adipocyte biology and highlight X-escape genes as a critical component of female physiology.

Non-coding mutations at the far end of a large gene desert surrounding the SOX9 gene result in a human craniofacial disorder called Pierre Robin sequence (PRS). Leveraging a human stem cell differentiation model, we identify two clusters of enhancers within the PRS-associated region that regulate SOX9 expression during a restricted window of facial progenitor development at distances up to 1.45 Mb. Enhancers within the 1.45 Mb cluster exhibit highly synergistic activity that is dependent on the Coordinator motif. Using mouse models, we demonstrate that PRS phenotypic specificity arises from the convergence of two mechanisms: confinement of Sox9 dosage perturbation to developing facial structures through context-specific enhancer activity and heightened sensitivity of the lower jaw to Sox9 expression reduction. Overall, we characterize the longest-range human enhancers involved in congenital malformations, directly demonstrate that PRS is an enhanceropathy, and illustrate how small changes in gene expression can lead to morphological variation.

Marek's disease (MD) is a neoplastic disease of chickens caused by Marek's disease virus (MDV), a member of the subfamily Alphaherpesvirinae Like other alphaherpesviruses, MDV encodes a serine/threonine protein kinase, US3. The functions of US3 have been extensively studied in other alphaherpesviruses; however, the biological functions of MDV US3 and its substrates have not been studied in detail. In this study, we investigated potential cellular pathways that are regulated by MDV US3 and identified chicken CREB (chCREB) as a substrate of MDV US3. We show that wild-type MDV US3, but not kinase-dead US3 (US3-K220A), increases CREB phosphorylation, leading to recruitment of phospho-CREB (pCREB) to the promoter of the CREB-responsive gene and activation of CREB target gene expression. Using US3 deletion and US3 kinase-dead recombinant MDV, we identified US3-responsive MDV genes during infection and found that the majority of US3-responsive genes were located in the MDV repeat regions. Chromatin immunoprecipitation sequencing (ChIP-seq) studies determined that some US3-regulated genes colocalized with Meq (an MDV-encoded oncoprotein) recruitment sites. Chromatin immunoprecipitation-PCR (ChIP-PCR) further confirmed Meq binding to the ICP4/LAT region, which is also regulated by US3. Furthermore, biochemical studies demonstrated that MDV US3 interacts with Meq in transfected cells and MDV-infected chicken embryonic fibroblasts in a phosphorylation-dependent manner. Finally, in vitro kinase studies revealed that Meq is a US3 substrate. MDV US3 thus acts as an upstream kinase of the CREB signaling pathway to regulate the transcription function of the CREB/Meq heterodimer, which targets cellular and viral gene expression.IMPORTANCE MDV is a potent oncogenic herpesvirus that induces T-cell lymphoma in infected chickens. Marek's disease continues to have a significant economic impact on the poultry industry worldwide. US3 encoded by alphaherpesviruses is a multifunctional kinase involved in the regulation of various cellular pathways. Using an MDV genome quantitative reverse transcriptase PCR (qRT-PCR) array and chromatin immunoprecipitation, we elucidated the role of MDV US3 in viral and cellular gene regulation. Our results provide insights into how viral kinase regulates host cell signaling pathways to activate both viral and host gene expression. This is an important step toward understanding host-pathogen interaction through activation of signaling cascades.

Genomic instability contributes to tumorigenesis through the amplification and deletion of cancer driver genes. DNA copy number (CN) profiling of ensembles of tumors allows a thermodynamic analysis of the profile for each tumor. The free energy of the distribution of CNs is found to be a monotonically increasing function of the average chromosomal ploidy. The dependence is universal across several cancer types. Surprisal analysis distinguishes two main known subgroups: tumors with cells that have or have not undergone whole-genome duplication (WGD). The analysis uncovers that CN states having a narrower distribution are energetically more favorable toward the WGD transition. Surprisal analysis also determines the deviations from a fully stable-state distribution. These deviations reflect constraints imposed by tumor fitness selection pressures. The results point to CN changes that are more common in high-ploidy tumors and thus support altered selection pressures upon WGD.

A key goal of whole-genome sequencing for studies of human genetics is to interrogate all forms of variation, including single-nucleotide variants, small insertion or deletion (indel) variants and structural variants. However, tools and resources for the study of structural variants have lagged behind those for smaller variants. Here we used a scalable pipeline1 to map and characterize structural variants in 17,795 deeply sequenced human genomes. We publicly release site-frequency data to create the largest, to our knowledge, whole-genome-sequencing-based structural variant resource so far. On average, individuals carry 2.9 rare structural variants that alter coding regions; these variants affect the dosage or structure of 4.2 genes and account for 4.0-11.2% of rare high-impact coding alleles. Using a computational model, we estimate that structural variants account for 17.2% of rare alleles genome-wide, with predicted deleterious effects that are equivalent to loss-of-function coding alleles; approximately 90% of such structural variants are noncoding deletions (mean 19.1 per genome). We report 158,991 ultra-rare structural variants and show that 2% of individuals carry ultra-rare megabase-scale structural variants, nearly half of which are balanced or complex rearrangements. Finally, we infer the dosage sensitivity of genes and noncoding elements, and reveal trends that relate to element class and conservation. This work will help to guide the analysis and interpretation of structural variants in the era of whole-genome sequencing.

Humanin is a member of a new family of peptides that are encoded by short open reading frames within the mitochondrial genome. It is conserved in animals and is both neuroprotective and cytoprotective. Here we report that in C. elegans the overexpression of humanin is sufficient to increase lifespan, dependent on daf-16/Foxo. Humanin transgenic mice have many phenotypes that overlap with the worm phenotypes and, similar to exogenous humanin treatment, have increased protection against toxic insults. Treating middle-aged mice twice weekly with the potent humanin analogue HNG, humanin improves metabolic healthspan parameters and reduces inflammatory markers. In multiple species, humanin levels generally decline with age, but here we show that levels are surprisingly stable in the naked mole-rat, a model of negligible senescence. Furthermore, in children of centenarians, who are more likely to become centenarians themselves, circulating humanin levels are much greater than age-matched control subjects. Further linking humanin to healthspan, we observe that humanin levels are decreased in human diseases such as Alzheimer's disease and MELAS (Mitochondrial Encephalopathy, Lactic Acidosis, and Stroke-like episodes). Together, these studies are the first to demonstrate that humanin is linked to improved healthspan and increased lifespan.

The invasive behavior of glioblastoma is essential to its aggressive potential. Here, we show that pleckstrin homology domain interacting protein (PHIP), acting through effects on the force transduction layer of the focal adhesion complex, drives glioblastoma motility and invasion. Immunofluorescence analysis localized PHIP to the leading edge of glioblastoma cells, together with several focal adhesion proteins: vinculin (VCL), talin 1 (TLN1), integrin beta 1 (ITGB1), as well as phosphorylated forms of paxillin (pPXN) and focal adhesion kinase (pFAK). Confocal microscopy specifically localized PHIP to the force transduction layer, together with TLN1 and VCL. Immunoprecipitation revealed a physical interaction between PHIP and VCL. Targeted suppression of PHIP resulted in significant down-regulation of these focal adhesion proteins, along with zyxin (ZYX), and produced profoundly disorganized stress fibers. Live-cell imaging of glioblastoma cells overexpressing a ZYX-GFP construct demonstrated a role for PHIP in regulating focal adhesion dynamics. PHIP silencing significantly suppressed the migratory and invasive capacity of glioblastoma cells, partially restored following TLN1 or ZYX cDNA overexpression. PHIP knockdown produced substantial suppression of tumor growth upon intracranial implantation, as well as significantly reduced microvessel density and secreted VEGF levels. PHIP copy number was elevated in the classical glioblastoma subtype and correlated with elevated EGFR levels. These results demonstrate PHIP's role in regulating the actin cytoskeleton, focal adhesion dynamics, and tumor cell motility, and identify PHIP as a key driver of glioblastoma migration and invasion.

Sex chromosome dosage (SCD) variation increases risk for neuropsychiatric impairment, which may reflect direct SCD effects on brain organization. Here, we 1) map cumulative X- and Y-chromosome dosage effects on regional cortical thickness (CT) and investigate potential functional implications of these effects using Neurosynth, 2) test if this map is organized by patterns of CT covariance that are evident in health, and 3) characterize SCD effects on CT covariance itself. We modeled SCD effects on CT and CT covariance for 308 equally sized regions of the cortical sheet using structural neuroimaging data from 301 individuals with varying numbers of sex chromosomes (169 euploid, 132 aneuploid). Mounting SCD increased CT in the rostral frontal cortex and decreased CT in the lateral temporal cortex, bilaterally. Regions targeted by SCD were associated with social functioning, language processing, and comprehension. Cortical regions with a similar degree of SCD-sensitivity showed heightened CT covariance in health. Finally, greater SCD also increased covariance among regions similarly affected by SCD. Our study both 1) develops novel methods for comparing typical and disease-related structural covariance networks in the brain and 2) uses these techniques to resolve and identify organizing principles for SCD effects on regional cortical anatomy and anatomical covariance.

Cancer develops through a process of somatic evolution1,2. Sequencing data from a single biopsy represent a snapshot of this process that can reveal the timing of specific genomic aberrations and the changing influence of mutational processes3. Here, by whole-genome sequencing analysis of 2,658 cancers as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA)4, we reconstruct the life history and evolution of mutational processes and driver mutation sequences of 38 types of cancer. Early oncogenesis is characterized by mutations in a constrained set of driver genes, and specific copy number gains, such as trisomy 7 in glioblastoma and isochromosome 17q in medulloblastoma. The mutational spectrum changes significantly throughout tumour evolution in 40% of samples. A nearly fourfold diversification of driver genes and increased genomic instability are features of later stages. Copy number alterations often occur in mitotic crises, and lead to simultaneous gains of chromosomal segments. Timing analyses suggest that driver mutations often precede diagnosis by many years, if not decades. Together, these results determine the evolutionary trajectories of cancer, and highlight opportunities for early cancer detection.

To date, over 1000 melanocytic neoplasms, spanning all stages of tumorigenesis, have been sequenced, offering detailed views into their -omic landscapes. This has coincided with advances in genetic engineering technologies that allow molecular biologists to edit the human genome with extreme precision and new mouse models to simulate disease progression. In this review, we describe how these technologies are being harnessed to provide insights into the evolution of melanoma at an unprecedented resolution, revealing that prior models of melanoma evolution, in which pathways are turned 'on' or 'off' in a binary fashion during the run-up to melanoma, are oversimplified.

Humans and orcas are among the very rare species that have a prolonged post-reproductive lifespan (PRLS), during which the aging process continues. Reactive oxygen species (ROS) derived from mitochondria and from the NADPH oxidase (NOX) enzymes of innate immune cells are known to contribute to aging, with the former thought to be dominant. CD33-related-Siglecs are immune receptors that recognize self-associated-molecular-patterns and modulate NOX-derived-ROS. We herewith demonstrate a strong correlation of lifespan with CD33rSIGLEC gene number in 26 species, independent of body weight or phylogeny. The correlation is stronger when considering total CD33rSIGLEC gene number rather than those encoding inhibitory and activating subsets, suggesting that lifetime balancing of ROS is important. Combining independent lines of evidence including the short half-life and spontaneous activation of neutrophils, we calculate that even without inter-current inflammation, a major source of lifetime ROS exposure may actually be neutrophil NOX-derived. However, genomes of human supercentenarians (>110 years) do not harbor a significantly higher number of functional CD33rSIGLEC genes. Instead, lifespan correlation with CD33rSIGLEC gene number was markedly strengthened by excluding the post-reproductive lifespan of humans and orcas (R2  = 0.83; P

Multi-omics datasets represent distinct aspects of the central dogma of molecular biology. Such high-dimensional molecular profiles pose challenges to data interpretation and hypothesis generation. ActivePathways is an integrative method that discovers significantly enriched pathways across multiple datasets using statistical data fusion, rationalizes contributing evidence and highlights associated genes. As part of the ICGC/TCGA Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium, which aggregated whole genome sequencing data from 2658 cancers across 38 tumor types, we integrated genes with coding and non-coding mutations and revealed frequently mutated pathways and additional cancer genes with infrequent mutations. We also analyzed prognostic molecular pathways by integrating genomic and transcriptomic features of 1780 breast cancers and highlighted associations with immune response and anti-apoptotic signaling. Integration of ChIP-seq and RNA-seq data for master regulators of the Hippo pathway across normal human tissues identified processes of tissue regeneration and stem cell regulation. ActivePathways is a versatile method that improves systems-level understanding of cellular organization in health and disease through integration of multiple molecular datasets and pathway annotations.

Mitochondrial DNA copy number (mtDNA-CN), a measure of the number of mitochondrial genomes per cell, is a minimally invasive proxy measure for mitochondrial function and has been associated with several aging-related diseases. Although quantitative real-time PCR (qPCR) is the current gold standard method for measuring mtDNA-CN, mtDNA-CN can also be measured from genotyping microarray probe intensities and DNA sequencing read counts. To conduct a comprehensive examination on the performance of these methods, we use known mtDNA-CN correlates (age, sex, white blood cell count, Duffy locus genotype, incident cardiovascular disease) to evaluate mtDNA-CN calculated from qPCR, two microarray platforms, as well as whole genome (WGS) and whole exome sequence (WES) data across 1,085 participants from the Atherosclerosis Risk in Communities (ARIC) study and 3,489 participants from the Multi-Ethnic Study of Atherosclerosis (MESA). We observe mtDNA-CN derived from WGS data is significantly more associated with known correlates compared to all other methods (p < 0.001). Additionally, mtDNA-CN measured from WGS is on average more significantly associated with traits by 5.6 orders of magnitude and has effect size estimates 5.8 times more extreme than the current gold standard of qPCR. We further investigated the role of DNA extraction method on mtDNA-CN estimate reproducibility and found mtDNA-CN estimated from cell lysate is significantly less variable than traditional phenol-chloroform-isoamyl alcohol (p = 5.44x10-4) and silica-based column selection (p = 2.82x10-7). In conclusion, we recommend the field moves towards more accurate methods for mtDNA-CN, as well as re-analyze trait associations as more WGS data becomes available from larger initiatives such as TOPMed.

Tumor DNA sequencing data can be interpreted by computational methods that analyze genomic heterogeneity to infer evolutionary dynamics. A growing number of studies have used these approaches to link cancer evolution with clinical progression and response to therapy. Although the inference of tumor phylogenies is rapidly becoming standard practice in cancer genome analyses, standards for evaluating them are lacking. To address this need, we systematically assess methods for reconstructing tumor subclonality. First, we elucidate the main algorithmic problems in subclonal reconstruction and develop quantitative metrics for evaluating them. Then we simulate realistic tumor genomes that harbor all known clonal and subclonal mutation types and processes. Finally, we benchmark 580 tumor reconstructions, varying tumor read depth, tumor type and somatic variant detection. Our analysis provides a baseline for the establishment of gold-standard methods to analyze tumor heterogeneity.

Sequences encoding Olduvai protein domains (formerly DUF1220) show the greatest human lineage-specific increase in copy number of any coding region in the genome and have been associated, in a dosage-dependent manner, with brain size, cognitive aptitude, autism, and schizophrenia. Tandem intragenic duplications of a three-domain block, termed the Olduvai triplet, in four NBPF genes in the chromosomal 1q21.1-0.2 region, are primarily responsible for the striking human-specific copy number increase. Interestingly, most of the Olduvai triplets are adjacent to, and transcriptionally coregulated with, three human-specific NOTCH2NL genes that have been shown to promote cortical neurogenesis. Until now, the underlying genomic events that drove the Olduvai hyperamplification in humans have remained unexplained. Here, we show that the presence or absence of an alternative first exon of the Olduvai triplet perfectly discriminates between amplified (58/58) and unamplified (0/12) triplets. We provide sequence and breakpoint analyses that suggest the alternative exon was produced by an nonallelic homologous recombination-based mechanism involving the duplicative transposition of an existing Olduvai exon found in the CON3 domain, which typically occurs at the C-terminal end of NBPF genes. We also provide suggestive in vitro evidence that the alternative exon may promote instability through a putative G-quadraplex (pG4)-based mechanism. Lastly, we use single-molecule optical mapping to characterize the intragenic structural variation observed in NBPF genes in 154 unrelated individuals and 52 related individuals from 16 families and show that the presence of pG4-containing Olduvai triplets is strongly correlated with high levels of Olduvai copy number variation. These results suggest that the same driver of genomic instability that allowed the evolutionarily recent, rapid, and extreme human-specific Olduvai expansion remains highly active in the human genome.

Multiparametric magnetic resonance imaging (mpMRI) has transformed the management of localized prostate cancer by improving identification of clinically significant disease at diagnosis. Approximately 20% of primary prostate tumors are invisible to mpMRI, and we hypothesize that this invisibility reflects fundamental molecular properties of the tumor. We therefore profiled the genomes and transcriptomes of 40 International Society of Urological Pathology grade 2 tumors: 20 mpMRI-invisible (Prostate Imaging-Reporting and Data System [PI-RADS] v2

Large structural variants (SVs) in the human genome are difficult to detect and study by conventional sequencing technologies. With long-range genome analysis platforms, such as optical mapping, one can identify large SVs (>2 kb) across the genome in one experiment. Analyzing optical genome maps of 154 individuals from the 26 populations sequenced in the 1000 Genomes Project, we find that phylogenetic population patterns of large SVs are similar to those of single nucleotide variations in 86% of the human genome, while ~2% of the genome has high structural complexity. We are able to characterize SVs in many intractable regions of the genome, including segmental duplications and subtelomeric, pericentromeric, and acrocentric areas. In addition, we discover ~60 Mb of non-redundant genome content missing in the reference genome sequence assembly. Our results highlight the need for a comprehensive set of alternate haplotypes from different populations to represent SV patterns in the genome.

Natural killer (NK) cells have increasingly become a target of interest for immunotherapies. NK cells express killer immunoglobulin-like receptors (KIRs), which play a vital role in immune response to tumors by detecting cellular abnormalities. The genomic region encoding the 16 KIR genes displays high polymorphic variability in human populations, making it difficult to resolve individual genotypes based on next generation sequencing data. As a result, the impact of polymorphic KIR variation on cancer phenotypes has been understudied. Currently, labor-intensive, experimental techniques are used to determine an individual's KIR gene copy number profile. Here, we develop an algorithm to determine the germline copy number of KIR genes from whole exome sequencing data and apply it to a cohort of nearly 5000 cancer patients. We use a k-mer based approach to capture sequences unique to specific genes, count their occurrences in the set of reads derived from an individual and compare the individual's k-mer distribution to that of the population. Copy number results demonstrate high concordance with population copy number expectations. Our method reveals that the burden of inhibitory KIR genes is associated with survival in two tumor types, highlighting the potential importance of KIR variation in understanding tumor development and response to immunotherapy.

Mitochondrial-derived peptides (MDPs) are rapidly emerging therapeutic targets to combat development of neurodegenerative diseases. SHLP2 (small humanin-like peptide 2) is a newly discovered MDP that is coded from the MT-RNR2 (Mitochondrially encoded 16S rRNA) gene in mitochondrial DNA (mtDNA). In the current study, we examined the biological consequences of treatment with exogenously-added SHLP2 in an in vitro human transmitochondrial age-related macular degeneration (AMD) ARPE-19 cell model. In AMD cells, we observed significant down-regulation of the MDP-coding MT-RNR2 gene, and remarkably reduced levels of all five oxidative phosphorylation (OXPHOS) complex I-V protein subunits that are involved in the electron transport chain; these results suggested mitochondrial toxicity and abnormal OXPHOS complex protein subunits' levels in AMD cells. However, treatment of AMD cells with SHLP2: (1) restored the normal levels of OXPHOS complex protein subunits, (2) prevented loss of viable cells and mitochondria, (3) increased the number of mtDNA copies, (4) induced anti-apoptotic effects, and (5) attenuated amyloid-β-induced cellular and mitochondrial toxicity. Cumulatively, our findings established the protective role of SHLP2 in AMD cells in vitro. In conclusion, this novel study supports the merit of SHLP2 in the treatment of AMD, a primary retinal disease that is a leading cause of blindness among the elderly population in the United States as well as worldwide.

A fundamental question in the biology of sex differences has eluded direct study in humans: How does sex-chromosome dosage (SCD) shape genome function? To address this, we developed a systematic map of SCD effects on gene function by analyzing genome-wide expression data in humans with diverse sex-chromosome aneuploidies (XO, XXX, XXY, XYY, and XXYY). For sex chromosomes, we demonstrate a pattern of obligate dosage sensitivity among evolutionarily preserved X-Y homologs and update prevailing theoretical models for SCD compensation by detecting X-linked genes that increase expression with decreasing X- and/or Y-chromosome dosage. We further show that SCD-sensitive sex-chromosome genes regulate specific coexpression networks of SCD-sensitive autosomal genes with critical cellular functions and a demonstrable potential to mediate previously documented SCD effects on disease. These gene coexpression results converge with analysis of transcription factor binding site enrichment and measures of gene expression in murine knockout models to spotlight the dosage-sensitive X-linked transcription factor ZFX as a key mediator of SCD effects on wider genome expression. Our findings characterize the effects of SCD broadly across the genome, with potential implications for human phenotypic variation.

We elucidated genomic and transcriptomic changes that accompany the evolution of melanoma from pre-malignant lesions by sequencing DNA and RNA from primary melanomas and their adjacent precursors, as well as matched primary tumors and regional metastases. In total, we analyzed 230 histopathologically distinct areas of melanocytic neoplasia from 82 patients. Somatic alterations sequentially induced mitogen-activated protein kinase (MAPK) pathway activation, upregulation of telomerase, modulation of the chromatin landscape, G1/S checkpoint override, ramp-up of MAPK signaling, disruption of the p53 pathway, and activation of the PI3K pathway; no mutations were specifically associated with metastatic progression, as these pathways were perturbed during the evolution of primary melanomas. UV radiation-induced point mutations steadily increased until melanoma invasion, at which point copy-number alterations also became prevalent.

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.

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.

Genome-wide association studies (GWAS) have identified over 100 risk loci for schizophrenia, but the causal mechanisms remain largely unknown. We performed a transcriptome-wide association study (TWAS) integrating a schizophrenia GWAS of 79,845 individuals from the Psychiatric Genomics Consortium with expression data from brain, blood, and adipose tissues across 3,693 primarily control individuals. We identified 157 TWAS-significant genes, of which 35 did not overlap a known GWAS locus. Of these 157 genes, 42 were associated with specific chromatin features measured in independent samples, thus highlighting potential regulatory targets for follow-up. Suppression of one identified susceptibility gene, mapk3, in zebrafish showed a significant effect on neurodevelopmental phenotypes. Expression and splicing from the brain captured most of the TWAS effect across all genes. This large-scale connection of associations to target genes, tissues, and regulatory features is an essential step in moving toward a mechanistic understanding of GWAS.

Variants of TREM2 are associated with Alzheimer's disease (AD). To study whether increasing TREM2 gene dosage could modify the disease pathogenesis, we developed BAC transgenic mice expressing human TREM2 (BAC-TREM2) in microglia. We found that elevated TREM2 expression reduced amyloid burden in the 5xFAD mouse model. Transcriptomic profiling demonstrated that increasing TREM2 levels conferred a rescuing effect, which includes dampening the expression of multiple disease-associated microglial genes and augmenting downregulated neuronal genes. Interestingly, 5xFAD/BAC-TREM2 mice showed further upregulation of several reactive microglial genes linked to phagocytosis and negative regulation of immune cell activation. Moreover, these mice showed enhanced process ramification and phagocytic marker expression in plaque-associated microglia and reduced neuritic dystrophy. Finally, elevated TREM2 gene dosage led to improved memory performance in AD models. In summary, our study shows that a genomic transgene-driven increase in TREM2 expression reprograms microglia responsivity and ameliorates neuropathological and behavioral deficits in AD mouse models.

Our growing awareness of the microbial world's importance and diversity contrasts starkly with our limited understanding of its fundamental structure. Despite recent advances in DNA sequencing, a lack of standardized protocols and common analytical frameworks impedes comparisons among studies, hindering the development of global inferences about microbial life on Earth. Here we present a meta-analysis of microbial community samples collected by hundreds of researchers for the Earth Microbiome Project. Coordinated protocols and new analytical methods, particularly the use of exact sequences instead of clustered operational taxonomic units, enable bacterial and archaeal ribosomal RNA gene sequences to be followed across multiple studies and allow us to explore patterns of diversity at an unprecedented scale. The result is both a reference database giving global context to DNA sequence data and a framework for incorporating data from future studies, fostering increasingly complete characterization of Earth's microbial diversity.

Several common alleles in the oxytocin receptor gene (OXTR) are associated with altered brain function in reward circuitry in neurotypical adults and may increase risk for autism spectrum disorders (ASD). Yet, it is currently unknown how variation in the OXTR relates to brain functioning in individuals with ASD, and, critically, whether neural endophenotypes vary as a function of aggregate genetic risk. Here, for we believe the first time, we use a multi-locus approach to examine how genetic variation across several OXTR single-nucleotide polymorphisms (SNPs) affect functional connectivity of the brain's reward network. Using data from 41 children with ASD and 41 neurotypical children, we examined functional connectivity of the nucleus accumbens (NAcc) - a hub of the reward network - focusing on how connectivity varies with OXTR risk-allele dosage. Youth with ASD showed reduced NAcc connectivity with other areas in the reward circuit as a function of increased OXTR risk-allele dosage, as well as a positive association between risk-allele dosage and symptom severity, whereas neurotypical youth showed increased NAcc connectivity with frontal brain regions involved in mentalizing. In addition, we found that increased NAcc-frontal cortex connectivity in typically developing youth was related to better scores on a standardized measure of social functioning. Our results indicate that cumulative genetic variation on the OXTR impacts reward system connectivity in both youth with ASD and neurotypical controls. By showing differential genetic effects on neuroendophenotypes, these pathways elucidate mechanisms of vulnerability versus resilience in carriers of disease-associated risk alleles.

Reciprocal chromosomal rearrangements at the 22q11.2 locus are associated with elevated risk of neurodevelopmental disorders. The 22q11.2 deletion confers the highest known genetic risk for schizophrenia, but a duplication in the same region is strongly associated with autism and is less common in schizophrenia cases than in the general population. Here we conducted the first study of 22q11.2 gene dosage effects on brain structure in a sample of 143 human subjects: 66 with 22q11.2 deletions (22q-del; 32 males), 21 with 22q11.2 duplications (22q-dup; 14 males), and 56 age- and sex-matched controls (31 males). 22q11.2 gene dosage varied positively with intracranial volume, gray and white matter volume, and cortical surface area (deletion < control < duplication). In contrast, gene dosage varied negatively with mean cortical thickness (deletion > control > duplication). Widespread differences were observed for cortical surface area with more localized effects on cortical thickness. These diametric patterns extended into subcortical regions: 22q-dup carriers had a significantly larger right hippocampus, on average, but lower right caudate and corpus callosum volume, relative to 22q-del carriers. Novel subcortical shape analysis revealed greater radial distance (thickness) of the right amygdala and left thalamus, and localized increases and decreases in subregions of the caudate, putamen, and hippocampus in 22q-dup relative to 22q-del carriers. This study provides the first evidence that 22q11.2 is a genomic region associated with gene-dose-dependent brain phenotypes. Pervasive effects on cortical surface area imply that this copy number variant affects brain structure early in the course of development.SIGNIFICANCE STATEMENT Probing naturally occurring reciprocal copy number variation in the genome may help us understand mechanisms underlying deviations from typical brain and cognitive development. The 22q11.2 genomic region is particularly susceptible to chromosomal rearrangements and contains many genes crucial for neuronal development and migration. Not surprisingly, reciprocal genomic imbalances at this locus confer some of the highest known genetic risks for developmental neuropsychiatric disorders. Here we provide the first evidence that brain morphology differs meaningfully as a function of reciprocal genomic variation at the 22q11.2 locus. Cortical thickness and surface area were affected in opposite directions with more widespread effects of gene dosage on cortical surface area.

Sequencing key cancer-driver genes using formalin-fixed, paraffin-embedded (FFPE) cancer tissues is becoming the standard for identifying the best treatment regimen. However, about 25% of all samples are rejected for genetic analyses for reasons that include too little tissue to extract enough high quality DNA. One way to overcome this is to do whole-genome amplification (WGA) in clinical samples, but only limited studies have tested different WGA methods in FFPE cancer specimens using targeted next-generation sequencing (NGS). We therefore tested the two most commonly used WGA methods, multiple displacement amplification (MDA-Qiagen REPLI-g kit) and the hybrid or modified PCR-based method (Sigma/Rubicon Genomics Inc. GenomePlex kit) in FFPE normal and tumor tissue specimens. For the normalized copy number analysis, the FFPE process caused none or very minimal bias. Variations in copy number were minimal in samples amplified using the GenomePlex kit, but they were statistically significantly higher in samples amplified using the REPLI-g kit. The pattern was similar for variant allele frequencies across the samples, which was minimal for the GenomePlex kit but highly variable for the REPLI-g kit. These findings suggest that each WGA method should be tested thoroughly before using it for clinical cancer samples.

Recent studies demonstrate that the invasion and metastasis of gastric cancer (GC) is closely associated with a multi-subunit vacuolar H+-ATPase (V-ATPase). Here we investigated the expression and role of the human ATP6V1A gene that encodes the catalytic subunit A of V-ATPase in GC. We found that ATP6V1A expression level is significantly elevated in GCs compared to normals, but GC patients with higher expression levels of ATP6V1A have a better prognosis. Genomic analysis revealed that APT6V1A copy number is gained in a small fraction of GC patients and lost in a minimum number. Moreover, the ATP6V1A copy number was positively correlated with its mRNA level. To explore additional mechanisms by which ATP6V1A overexpressed in GCs, we investigated the relationship between transcription factor YY1 and ATP6V1A, and found that mRNA expression of YY1 had significant correlation with that of ATP6V1A. To validate that YY1 transcriptionally regulates ATP6V1A, we discovered that the ATP6V1A core promoter region contains three YY1 binding sites. Moreover, RNAi-mediated knockdown of YY1 in GC cells significantly decreased ATP6V1A mRNA and protein expression, while YY1 overexpression increased ATP6V1A expression level. In conclusion, YY1 may play an important regulatory role in ATP6V1A expression with potential mechanistic and clinical implications in GC.

The cohesin complex is an evolutionarily conserved multi-subunit protein complex which regulates sister chromatid cohesion during mitosis and meiosis. Additionally, the cohesin complex regulates DNA replication, DNA repair, and transcription. The core of the complex consists of four subunits: SMC1A, SMC3, RAD21, and STAG1/2. Loss-of-function mutations in many of these proteins have been implicated in human developmental disorders collectively termed cohesinopathies. Through clinical exome sequencing (CES) of an 8-year-old girl with a clinical history of global developmental delay, microcephaly, microtia with hearing loss, language delay, ADHD, and dysmorphic features, we describe a heterozygous de novo variant (c.205C>T; p.(Arg69*)) in the integral cohesin structural protein, STAG2. This variant is associated with decreased STAG2 protein expression. The analyses of metaphase spreads did not exhibit premature sister chromatid separation; however, delayed sister chromatid cohesion was observed. To further support the pathogenicity of STAG2 variants, we identified two additional female cases from the DECIPHER research database with mutations in STAG2 and phenotypes similar to our patient. Interestingly, the clinical features of these three cases are remarkably similar to those observed in other well-established cohesinopathies. Herein, we suggest that STAG2 is a dosage-sensitive gene and that heterozygous loss-of-function variants lead to a cohesinopathy.

Genomic aberrations have been identified in many human pluripotent stem cell (hPSC) cultures. Commonly observed duplications in portions of chromosomes 12p and 17q have been associated with increases in genetic instability and resistance to apoptosis, respectively. However, the phenotypic consequences related to sporadic mutations have not been evaluated to date. Here, we report on the effects of a single-copy deletion of the chr17p13.1 region, a sporadic mutation that spontaneously arose independently in several subclones of a human embryonic stem cell culture. Compared to cells with two normal copies of chr17p13.1 ("wild-type"), the cells with a single-copy deletion of this region ("mutant") displayed a selective advantage when exposed to stressful conditions, and retained a higher percentage of cells expressing the pluripotency marker POU5F1/OCT4 after 2 weeks of in vitro differentiation. Knockdown of TP53, which is a gene encompassed by the deleted region, in wild-type cells mimicked the chr17p13.1 deletion phenotype. Thus, sporadic mutations in hPSCs can have phenotypic effects that may impact their utility for clinical applications. Stem Cells 2017;35:872-885.

Many bacterial genomes are highly variable but nonetheless are typically published as a single assembled genome. Experiments tracking bacterial genome evolution have not looked at the variation present at a given point in time. Here, we analyzed the mouse-passaged Helicobacter pylori strain SS1 and its parent PMSS1 to assess intra- and intergenomic variability. Using high sequence coverage depth and experimental validation, we detected extensive genome plasticity within these H. pylori isolates, including movement of the transposable element IS607, large and small inversions, multiple single nucleotide polymorphisms, and variation in cagA copy number. The cagA gene was found as 1 to 4 tandem copies located off the cag island in both SS1 and PMSS1; this copy number variation correlated with protein expression. To gain insight into the changes that occurred during mouse adaptation, we also compared SS1 and PMSS1 and observed 46 differences that were distinct from the within-genome variation. The most substantial was an insertion in cagY, which encodes a protein required for a type IV secretion system function. We detected modifications in genes coding for two proteins known to affect mouse colonization, the HpaA neuraminyllactose-binding protein and the FutB α-1,3 lipopolysaccharide (LPS) fucosyltransferase, as well as genes predicted to modulate diverse properties. In sum, our work suggests that data from consensus genome assemblies from single colonies may be misleading by failing to represent the variability present. Furthermore, we show that high-depth genomic sequencing data of a population can be analyzed to gain insight into the normal variation within bacterial strains.IMPORTANCE Although it is well known that many bacterial genomes are highly variable, it is nonetheless traditional to refer to, analyze, and publish "the genome" of a bacterial strain. Variability is usually reduced ("only sequence from a single colony"), ignored ("just publish the consensus"), or placed in the "too-hard" basket ("analysis of raw read data is more robust"). Now that whole-genome sequences are regularly used to assess virulence and track outbreaks, a better understanding of the baseline genomic variation present within single strains is needed. Here, we describe the variability seen in typical working stocks and colonies of pathogen Helicobacter pylori model strains SS1 and PMSS1 as revealed by use of high-coverage mate pair next-generation sequencing (NGS) and confirmed by traditional laboratory techniques. This work demonstrates that reliance on a consensus assembly as "the genome" of a bacterial strain may be misleading.

The healthy human brain is a mosaic of varied genomes. Long interspersed element-1 (LINE-1 or L1) retrotransposition is known to create mosaicism by inserting L1 sequences into new locations of somatic cell genomes. Using a machine learning-based, single-cell sequencing approach, we discovered that somatic L1-associated variants (SLAVs) are composed of two classes: L1 retrotransposition insertions and retrotransposition-independent L1-associated variants. We demonstrate that a subset of SLAVs comprises somatic deletions generated by L1 endonuclease cutting activity. Retrotransposition-independent rearrangements in inherited L1s resulted in the deletion of proximal genomic regions. These rearrangements were resolved by microhomology-mediated repair, which suggests that L1-associated genomic regions are hotspots for somatic copy number variants in the brain and therefore a heritable genetic contributor to somatic mosaicism. We demonstrate that SLAVs are present in crucial neural genes, such as DLG2 (also called PSD93), and affect 44-63% of cells of the cells in the healthy brain.

Multiple myeloma is incurable by standard approaches because of inevitable relapse and development of treatment resistance in all patients. In our prior work, we identified a panel of macropinocytosing human monoclonal antibodies against CD46, a negative regulator of the innate immune system, and constructed antibody-drug conjugates (ADCs). In this report, we show that an anti-CD46 ADC (CD46-ADC) potently inhibited proliferation in myeloma cell lines with little effect on normal cells. CD46-ADC also potently eliminated myeloma growth in orthometastatic xenograft models. In primary myeloma cells derived from bone marrow aspirates, CD46-ADC induced apoptosis and cell death, but did not affect the viability of nontumor mononuclear cells. It is of clinical interest that the CD46 gene resides on chromosome 1q, which undergoes genomic amplification in the majority of relapsed myeloma patients. We found that the cell surface expression level of CD46 was markedly higher in patient myeloma cells with 1q gain than in those with normal 1q copy number. Thus, genomic amplification of CD46 may serve as a surrogate for target amplification that could allow patient stratification for tailored CD46-targeted therapy. Overall, these findings indicate that CD46 is a promising target for antibody-based treatment of multiple myeloma, especially in patients with gain of chromosome 1q.

The cellular proteome is a complex microcosm of structural and regulatory networks that requires continuous surveillance and modification to meet the dynamic needs of the cell. It is therefore crucial that the protein flux of the cell remains in balance to ensure proper cell function. Genetic alterations that range from chromosome imbalance to oncogene activation can affect the speed, fidelity and capacity of protein biogenesis and degradation systems, which often results in proteome imbalance. An improved understanding of the causes and consequences of proteome imbalance is helping to reveal how these systems can be targeted to treat diseases such as cancer.

ObjectiveChromosomal aberrations are frequently associated with birth defects and pregnancy losses. Trisomy 13, Trisomy 18 and Trisomy 21 are the most common, clinically relevant fetal aneusomies. This study used a transcriptomics approach to identify the molecular signatures at the maternal-fetal interface in each aneuploidy.MethodsWe profiled placental gene expression (13-22 weeks) in T13 (n = 4), T18 (n = 4) and T21 (n = 8), and in euploid pregnancies (n = 4).ResultsWe found differentially expressed transcripts (≥2-fold) in T21 (n = 160), T18 (n = 80) and T13 (n = 125). The majority were upregulated and most of the misexpressed genes were not located on the relevant trisomic chromosome, suggesting genome-wide dysregulation. A smaller number of the differentially expressed transcripts were encoded on the trisomic chromosome, suggesting gene dosage. In T21,

Copy number variants at the 16p11.2 chromosomal locus are associated with several neuropsychiatric disorders, including autism, schizophrenia, bipolar disorder, attention-deficit hyperactivity disorder, and speech and language disorders. A gene dosage dependence has been suggested, with 16p11.2 deletion carriers demonstrating higher body mass index and head circumference, and 16p11.2 duplication carriers demonstrating lower body mass index and head circumference. Here, we use diffusion tensor imaging to elucidate this reciprocal relationship in white matter organization, showing widespread increases of fractional anisotropy throughout the supratentorial white matter in pediatric deletion carriers and, in contrast, extensive decreases of white matter fractional anisotropy in pediatric and adult duplication carriers. We find associations of these white matter alterations with cognitive and behavioral impairments. We further demonstrate the value of imaging metrics for characterizing the copy number variant phenotype by employing linear discriminant analysis to predict the gene dosage status of the study subjects. These results show an effect of 16p11.2 gene dosage on white matter microstructure, and further suggest that opposite changes in diffusion tensor imaging metrics can lead to similar cognitive and behavioral deficits. Given the large effect sizes found in this study, our results support the view that specific genetic variations are more strongly associated with specific brain alterations than are shared neuropsychiatric diagnoses. Hum Brain Mapp 37:2833-2848, 2016. © 2016 Wiley Periodicals, Inc.

PTEN loss is a promising prognostic and predictive biomarker in prostate cancer. Because it occurs most commonly via PTEN gene deletion, we developed a clinical-grade, automated, and inexpensive immunohistochemical assay to detect PTEN loss. We studied the sensitivity and specificity of PTEN immunohistochemistry relative to four-color fluorescence in situ hybridization (FISH) for detection of PTEN gene deletion in a multi-institutional cohort of 731 primary prostate tumors. Intact PTEN immunostaining was 91% specific for the absence of PTEN gene deletion (549/602 tumors with two copies of the PTEN gene by FISH showed intact expression of PTEN by immunohistochemistry) and 97% sensitive for the presence of homozygous PTEN gene deletion (absent PTEN protein expression by immunohistochemistry in 65/67 tumors with homozygous deletion). PTEN immunohistochemistry was 65% sensitive for the presence of hemizygous PTEN gene deletion, with protein loss in 40/62 hemizygous tumors. We reviewed the 53 cases where immunohistochemistry showed PTEN protein loss and FISH showed two intact copies of the PTEN gene. On re-review, there was ambiguous immunohistochemistry loss in 6% (3/53) and failure to analyze the same tumor area by both methods in 34% (18/53). Of the remaining discordant cases, 41% (13/32) revealed hemizygous (n=8) or homozygous (n=5) PTEN gene deletion that was focal in most cases (11/13). The remaining 19 cases had two copies of the PTEN gene detected by FISH, representing truly discordant cases. Our automated PTEN immunohistochemistry assay is a sensitive method for detection of homozygous PTEN gene deletions. Immunohistochemistry screening is particularly useful to identify cases with heterogeneous PTEN gene deletion in a subset of tumor glands. Mutations, small insertions, or deletions and/or epigenetic or microRNA-mediated mechanisms may lead to PTEN protein loss in tumors with normal or hemizygous PTEN gene copy number.

The physiological functions of natural killer (NK) cells in human immunity and reproduction depend upon diverse interactions between killer cell immunoglobulin-like receptors (KIRs) and their HLA class I ligands: HLA-A, HLA-B, and HLA-C. The genomic regions containing the KIR and HLA class I genes are unlinked, structurally complex, and highly polymorphic. They are also strongly associated with a wide spectrum of diseases, including infections, autoimmune disorders, cancers, and pregnancy disorders, as well as the efficacy of transplantation and other immunotherapies. To facilitate study of these extraordinary genes, we developed a method that captures, sequences, and analyzes the 13 KIR genes and HLA-A, HLA-B, and HLA-C from genomic DNA. We also devised a bioinformatics pipeline that attributes sequencing reads to specific KIR genes, determines copy number by read depth, and calls high-resolution genotypes for each KIR gene. We validated this method by using DNA from well-characterized cell lines, comparing it to established methods of HLA and KIR genotyping, and determining KIR genotypes from 1000 Genomes sequence data. This identified 116 previously uncharacterized KIR alleles, which were all demonstrated to be authentic by sequencing from source DNA via standard methods. Analysis of just two KIR genes showed that 22% of the 1000 Genomes individuals have a previously uncharacterized allele or a structural variant. The method we describe is suited to the large-scale analyses that are needed for characterizing human populations and defining the precise HLA and KIR factors associated with disease. The methods are applicable to other highly polymorphic genes.

The three-dimensional organization of the genome plays important roles in regulating the functional output of the genome and even in the maintenance of epigenetic inheritance and genome stability. Here, we review and compare a number of newly developed methods-especially those that utilize the CRISPR (clustered regularly interspaced short palindromic repeat)-Cas9 (CRISPR-associated protein 9) system-that enable the direct visualization of specific, endogenous DNA sequences in living cells. We also discuss the practical considerations in implementing the CRISPR imaging technique to achieve sufficient signal-to-background levels, high specificity, and high labeling efficiency. These DNA labeling methods enable tracking of the copy number, localization, and movement of genomic elements, and we discuss the potential applications of these methods in understanding the searching and targeting mechanism of the Cas9-sgRNA complex, investigating chromosome organization, and visualizing genome instability and rearrangement.

To provide a detailed analysis of the molecular components and underlying mechanisms associated with ovarian cancer, we performed a comprehensive mass-spectrometry-based proteomic characterization of 174 ovarian tumors previously analyzed by The Cancer Genome Atlas (TCGA), of which 169 were high-grade serous carcinomas (HGSCs). Integrating our proteomic measurements with the genomic data yielded a number of insights into disease, such as how different copy-number alternations influence the proteome, the proteins associated with chromosomal instability, the sets of signaling pathways that diverse genome rearrangements converge on, and the ones most associated with short overall survival. Specific protein acetylations associated with homologous recombination deficiency suggest a potential means for stratifying patients for therapy. In addition to providing a valuable resource, these findings provide a view of how the somatic genome drives the cancer proteome and associations between protein and post-translational modification levels and clinical outcomes in HGSC. VIDEO ABSTRACT.