Germline mutation in POLR2A: a heterogeneous, multi-systemic developmental disorder characterized by transcriptional dysregulation.

De novo germline variation in POLR2A was recently reported to associate with a neurodevelopmental disorder. We report twelve individuals harboring putatively pathogenic de novo or inherited variants in POLR2A, detail their phenotypes, and map all known variants to the domain structure of POLR2A and crystal structure of RNA polymerase II. Affected individuals were ascertained from a local data lake, pediatric genetics clinic, and an online community of families of affected individuals. These include six affected by de novo missense variants (including one previously reported individual), four clinical laboratory samples affected by missense variation with unknown inheritance-with yeast functional assays further supporting altered function-one affected by a de novo in-frame deletion, and one affected by a C-terminal frameshift variant inherited from a largely asymptomatic mother. Recurrently observed phenotypes include ataxia, joint hypermobility, short stature, skin abnormalities, congenital cardiac abnormalities, immune system abnormalities, hip dysplasia, and short Achilles tendons. We report a significantly higher occurrence of epilepsy (8/12, 66.7%) than previously reported (3/15, 20%) (p value = 0.014196; chi-square test) and a lower occurrence of hypotonia (8/12, 66.7%) than previously reported (14/15, 93.3%) (p value = 0.076309). POLR2A-related developmental disorders likely represent a spectrum of related, multi-systemic developmental disorders, driven by distinct mechanisms, converging at a single locus.

Phenotypic expansion in KIF1A-related dominant disorders: A description of novel variants and review of published cases.

KIF1A is a molecular motor for membrane-bound cargo important to the development and survival of sensory neurons. KIF1A dysfunction has been associated with several Mendelian disorders with a spectrum of overlapping phenotypes, ranging from spastic paraplegia to intellectual disability. We present a novel pathogenic in-frame deletion in the KIF1A molecular motor domain inherited by two affected siblings from an unaffected mother with apparent germline mosaicism. We identified eight additional cases with heterozygous, pathogenic KIF1A variants ascertained from a local data lake. Our data provide evidence for the expansion of KIF1A-associated phenotypes to include hip subluxation and dystonia as well as phenotypes observed in only a single case: gelastic cataplexy, coxa valga, and double collecting system. We review the literature and suggest that KIF1A dysfunction is better understood as a single neuromuscular disorder with variable involvement of other organ systems than a set of discrete disorders converging at a single locus.

Gene content evolution in the arthropods.

BACKGROUND: Arthropods comprise the largest and most diverse phylum on Earth and play vital roles in nearly every ecosystem. Their diversity stems in part from variations on a conserved body plan, resulting from and recorded in adaptive changes in the genome. Dissection of the genomic record of sequence change enables broad questions regarding genome evolution to be addressed, even across hyper-diverse taxa within arthropods. RESULTS: Using 76 whole genome sequences representing 21 orders spanning more than 500 million years of arthropod evolution, we document changes in gene and protein domain content and provide temporal and phylogenetic context for interpreting these innovations. We identify many novel gene families that arose early in the evolution of arthropods and during the diversification of insects into modern orders. We reveal unexpected variation in patterns of DNA methylation across arthropods and examples of gene family and protein domain evolution coincident with the appearance of notable phenotypic and physiological adaptations such as flight, metamorphosis, sociality, and chemoperception. CONCLUSIONS: These analyses demonstrate how large-scale comparative genomics can provide broad new insights into the genotype to phenotype map and generate testable hypotheses about the evolution of animal diversity.

Genomic Profiling of Childhood Tumor Patient-Derived Xenograft Models to Enable Rational Clinical Trial Design.

Accelerating cures for children with cancer remains an immediate challenge as a result of extensive oncogenic heterogeneity between and within histologies, distinct molecular mechanisms evolving between diagnosis and relapsed disease, and limited therapeutic options. To systematically prioritize and rationally test novel agents in preclinical murine models, researchers within the Pediatric Preclinical Testing Consortium are continuously developing patient-derived xenografts (PDXs)-many of which are refractory to current standard-of-care treatments-from high-risk childhood cancers. Here, we genomically characterize 261 PDX models from 37 unique pediatric cancers; demonstrate faithful recapitulation of histologies and subtypes; and refine our understanding of relapsed disease. In addition, we use expression signatures to classify tumors for TP53 and NF1 pathway inactivation. We anticipate that these data will serve as a resource for pediatric oncology drug development and will guide rational clinical trial design for children with cancer.

Genetic architecture of laterality defects revealed by whole exome sequencing.

Aberrant left-right patterning in the developing human embryo can lead to a broad spectrum of congenital malformations. The causes of most laterality defects are not known, with variants in established genes accounting for <20% of cases. We sought to characterize the genetic spectrum of these conditions by performing whole-exome sequencing of 323 unrelated laterality cases. We investigated the role of rare, predicted-damaging variation in 1726 putative laterality candidate genes derived from model organisms, pathway analyses, and human phenotypes. We also evaluated the contribution of homo/hemizygous exon deletions and gene-based burden of rare variation. A total of 28 candidate variants (26 rare predicted-damaging variants and 2 hemizygous deletions) were identified, including variants in genes known to cause heterotaxy and primary ciliary dyskinesia (ACVR2B, NODAL, ZIC3, DNAI1, DNAH5, HYDIN, MMP21), and genes without a human phenotype association, but with prior evidence for a role in embryonic laterality or cardiac development. Sanger validation of the latter variants in probands and their parents revealed no de novo variants, but apparent transmitted heterozygous (ROCK2, ISL1, SMAD2), and hemizygous (RAI2, RIPPLY1) variant patterns. Collectively, these variants account for 7.1% of our study subjects. We also observe evidence for an excess burden of rare, predicted loss-of-function variation in PXDNL and BMS1- two genes relevant to the broader laterality phenotype. These findings highlight potential new genes in the development of laterality defects, and suggest extensive locus heterogeneity and complex genetic models in this class of birth defects.

Elucidating the molecular pathogenesis of glioma: integrated germline and somatic profiling of a familial glioma case series.

Background: The genomic characterization of sporadically arising gliomas has delineated molecularly and clinically distinct subclasses of disease. However, less is known about the molecular nature of gliomas that are familial in origin. We performed molecular subtyping of 163 tumor specimens from individuals with a family history of glioma and integrated germline and somatic genomic data to characterize the pathogenesis of 20 tumors in additional detail. Methods: Immunohistochemical analyses were performed on formalin-fixed, paraffin-embedded tumor sections to determine molecular subtypes of glioma. For 20 cases, tumor DNA was exome sequenced on an Illumina HiSeq 2000 platform and copy number profiling was performed on the Illumina HumanOmniExpress BeadChip. Genotypes at glioma risk polymorphisms were determined from germline DNA profiled on the Illumina Infinium OncoArray and deleterious germline mutations were identified from germline sequencing data. Results: All 3 molecular subtypes of sporadic glioma were represented in the overall case series, including molecular glioblastoma (n = 102), oligodendroglioma (n = 21), and astrocytoma (n = 20). Detailed profiling of 20 of these cases showed characteristic subtype-specific alterations at frequencies comparable to sporadic glioma cases. All 20 cases had alterations in genes regulating telomere length. Frequencies of common glioma risk alleles were similar to those among sporadic cases, and correlations between risk alleles and same-gene somatic mutations were not observed. Conclusions: This study illustrates that the molecular characteristics of familial tumors profiled largely recapitulate what is known about sporadic glioma and that both germline and somatic molecular features target common core pathways involved in gliomagenesis. Key Points: 1. Familial and sporadic gliomas display highly comparable molecular landscapes. 2. Germline and somatic molecular events target common core pathways involved in gliomagenesis. 3. Carriage of germline glioma risk variants is not associated with somatic events in the same gene.

Non-malignant respiratory epithelial cells preferentially proliferate from resected non-small cell lung cancer specimens cultured under conditionally reprogrammed conditions.

The "conditionally reprogrammed cells" (CRC) method, using a Rho kinase inhibitor and irradiated mouse fibroblast cells has been described for the efficient growth of cells from malignant and non-malignant samples from primary tumor and non-malignant sites. Using the CRC method, four institutions independently cultured tumor tissues from 48 non-small cell lung cancers (NSCLC, mostly from primary resected tumors) and 22 non-malignant lungs. We found that epithelial cells could be cultured from tumor and non-malignant lung. However, epithelial cells cultured from tumors had features of non-malignant respiratory epithelial cells which include: 1) among 22 mutations found in the original tumors only two mutations were found in the CRC cultures with reduced frequency (31% to 13% and 92% to 15% from original tumor and CRC culture respectively); 2) copy number variation was analyzed in 9 tumor and their CRC cultures and only diploid patterns were found in CRC cultures; 3) mRNA expression profiles were similar to those of normal respiratory epithelial cells; and 4) co-culture of tumor and non-malignant lung epithelial cells resulted in mostly non-malignant cells. We conclude that CRC method is a highly selective and useful method for the growth of non-malignant respiratory epithelial cells from tumor specimens and only occasionally do such CRC cultures contain a small subpopulation of cancer cells marked by oncogenic mutations. While our findings are restricted to resected primary NSCLC, they indicated the necessity to fully characterize all CRC cultures and the need to develop culture technology that facilitates the growth of primary lung cancers.

Evidence for stabilizing selection on codon usage in chromosomal rearrangements of Drosophila pseudoobscura.

There has been a renewed interest in investigating the role of stabilizing selection acting on genome-wide traits such as codon usage bias. Codon bias, when synonymous codons are used at unequal frequencies, occurs in a wide variety of taxa. Standard evolutionary models explain the maintenance of codon bias through a balance of genetic drift, mutation and weak purifying selection. The efficacy of selection is expected to be reduced in regions of suppressed recombination. Contrary to observations in Drosophila melanogaster, some recent studies have failed to detect a relationship between the recombination rate, intensity of selection acting at synonymous sites, and the magnitude of codon bias as predicted under these standard models. Here, we examined codon bias in 2798 protein coding loci on the third chromosome of D. pseudoobscura using whole-genome sequences of 47 individuals, representing five common third chromosome gene arrangements. Fine-scale recombination maps were constructed using more than 1 million segregating sites. As expected, recombination was demonstrated to be significantly suppressed between chromosome arrangements, allowing for a direct examination of the relationship between recombination, selection, and codon bias. As with other Drosophila species, we observe a strong mutational bias away from the most frequently used codons. We find the rate of synonymous and nonsynonymous polymorphism is variable between different amino acids. However, we do not observe a reduction in codon bias or the strength of selection in regions of suppressed recombination as expected. Instead, we find that the interaction between weak stabilizing selection and mutational bias likely plays a role in shaping the composition of synonymous codons across the third chromosome in D. pseudoobscura.

Comparative validation of the D. melanogaster modENCODE transcriptome annotation.

Accurate gene model annotation of reference genomes is critical for making them useful. The modENCODE project has improved the D. melanogaster genome annotation by using deep and diverse high-throughput data. Since transcriptional activity that has been evolutionarily conserved is likely to have an advantageous function, we have performed large-scale interspecific comparisons to increase confidence in predicted annotations. To support comparative genomics, we filled in divergence gaps in the Drosophila phylogeny by generating draft genomes for eight new species. For comparative transcriptome analysis, we generated mRNA expression profiles on 81 samples from multiple tissues and developmental stages of 15 Drosophila species, and we performed cap analysis of gene expression in D. melanogaster and D. pseudoobscura. We also describe conservation of four distinct core promoter structures composed of combinations of elements at three positions. Overall, each type of genomic feature shows a characteristic divergence rate relative to neutral models, highlighting the value of multispecies alignment in annotating a target genome that should prove useful in the annotation of other high priority genomes, especially human and other mammalian genomes that are rich in noncoding sequences. We report that the vast majority of elements in the annotation are evolutionarily conserved, indicating that the annotation will be an important springboard for functional genetic testing by the Drosophila community.

Natural variation in genome architecture among 205 Drosophila melanogaster Genetic Reference Panel lines.

The Drosophila melanogaster Genetic Reference Panel (DGRP) is a community resource of 205 sequenced inbred lines, derived to improve our understanding of the effects of naturally occurring genetic variation on molecular and organismal phenotypes. We used an integrated genotyping strategy to identify 4,853,802 single nucleotide polymorphisms (SNPs) and 1,296,080 non-SNP variants. Our molecular population genomic analyses show higher deletion than insertion mutation rates and stronger purifying selection on deletions. Weaker selection on insertions than deletions is consistent with our observed distribution of genome size determined by flow cytometry, which is skewed toward larger genomes. Insertion/deletion and single nucleotide polymorphisms are positively correlated with each other and with local recombination, suggesting that their nonrandom distributions are due to hitchhiking and background selection. Our cytogenetic analysis identified 16 polymorphic inversions in the DGRP. Common inverted and standard karyotypes are genetically divergent and account for most of the variation in relatedness among the DGRP lines. Intriguingly, variation in genome size and many quantitative traits are significantly associated with inversions. Approximately 50% of the DGRP lines are infected with Wolbachia, and four lines have germline insertions of Wolbachia sequences, but effects of Wolbachia infection on quantitative traits are rarely significant. The DGRP complements ongoing efforts to functionally annotate the Drosophila genome. Indeed, 15% of all D. melanogaster genes segregate for potentially damaged proteins in the DGRP, and genome-wide analyses of quantitative traits identify novel candidate genes. The DGRP lines, sequence data, genotypes, quality scores, phenotypes, and analysis and visualization tools are publicly available.

Epistasis dominates the genetic architecture of Drosophila quantitative traits.

Epistasis-nonlinear genetic interactions between polymorphic loci-is the genetic basis of canalization and speciation, and epistatic interactions can be used to infer genetic networks affecting quantitative traits. However, the role that epistasis plays in the genetic architecture of quantitative traits is controversial. Here, we compared the genetic architecture of three Drosophila life history traits in the sequenced inbred lines of the Drosophila melanogaster Genetic Reference Panel (DGRP) and a large outbred, advanced intercross population derived from 40 DGRP lines (Flyland). We assessed allele frequency changes between pools of individuals at the extremes of the distribution for each trait in the Flyland population by deep DNA sequencing. The genetic architecture of all traits was highly polygenic in both analyses. Surprisingly, none of the SNPs associated with the traits in Flyland replicated in the DGRP and vice versa. However, the majority of these SNPs participated in at least one epistatic interaction in the DGRP. Despite apparent additive effects at largely distinct loci in the two populations, the epistatic interactions perturbed common, biologically plausible, and highly connected genetic networks. Our analysis underscores the importance of epistasis as a principal factor that determines variation for quantitative traits and provides a means to uncover genetic networks affecting these traits. Knowledge of epistatic networks will contribute to our understanding of the genetic basis of evolutionarily and clinically important traits and enhance predictive ability at an individualized level in medicine and agriculture.

The Drosophila melanogaster Genetic Reference Panel.

A major challenge of biology is understanding the relationship between molecular genetic variation and variation in quantitative traits, including fitness. This relationship determines our ability to predict phenotypes from genotypes and to understand how evolutionary forces shape variation within and between species. Previous efforts to dissect the genotype-phenotype map were based on incomplete genotypic information. Here, we describe the Drosophila melanogaster Genetic Reference Panel (DGRP), a community resource for analysis of population genomics and quantitative traits. The DGRP consists of fully sequenced inbred lines derived from a natural population. Population genomic analyses reveal reduced polymorphism in centromeric autosomal regions and the X chromosome, evidence for positive and negative selection, and rapid evolution of the X chromosome. Many variants in novel genes, most at low frequency, are associated with quantitative traits and explain a large fraction of the phenotypic variance. The DGRP facilitates genotype-phenotype mapping using the power of Drosophila genetics.