The genome of the colonial hydroid Hydractinia reveals that their stem cells use a toolkit of evolutionarily shared genes with all animals

Hydractinia is a colonial marine hydroid that shows remarkable biological properties, including the capacity to regenerate its entire body throughout its lifetime, a process made possible by its adult migratory stem cells, known as i-cells. Here, we provide an in-depth characterization of the genomic structure and gene content of two Hydractinia species, Hydractinia symbiolongicarpus and Hydractinia echinata, placing them in a comparative evolutionary framework with other cnidarian genomes. We also generated and annotated a single-cell transcriptomic atlas for adult male H. symbiolongicarpus and identified cell-type markers for all major cell types, including key i-cell markers. Orthology analyses based on the markers revealed that Hydractinia's i-cells are highly enriched in genes that are widely shared amongst animals, a striking finding given that Hydractinia has a higher proportion of phylum-specific genes than any of the other 41 animals in our orthology analysis. These results indicate that Hydractinia's stem cells and early progenitor cells may use a toolkit shared with all animals, making it a promising model organism for future exploration of stem cell biology and regenerative medicine. The genomic and transcriptomic resources for Hydractinia presented here will enable further studies of their regenerative capacity, colonial morphology, and ability to distinguish self from nonself.

Random forest classifiers trained on simulated data enable accurate short read-based genotyping of structural variants in the alpha globin region at Chr16p13.3

In regions where reads don't align well to a reference, it is generally difficult to characterize structural variation using short read sequencing. Here, we utilize machine learning classifiers and short sequence reads to genotype structural variants in the alpha globin locus on chromosome 16, a medically-relevant region that is challenging to genotype in individuals. Using models trained only with simulated data, we accurately genotype two hard-to-distinguish deletions in two separate human cohorts. Furthermore, population allele frequencies produced by our methods across a wide set of ancestries agree more closely with previously-determined frequencies than those obtained using currently available genotyping software.

The genome of the colonial hydroid Hydractinia reveals their stem cells utilize a toolkit of evolutionarily shared genes with all animals

Hydractinia is a colonial marine hydroid that exhibits remarkable biological properties, including the capacity to regenerate its entire body throughout its lifetime, a process made possible by its adult migratory stem cells, known as i-cells. Here, we provide an in-depth characterization of the genomic structure and gene content of two Hydractinia species, H. symbiolongicarpus and H. echinata, placing them in a comparative evolutionary framework with other cnidarian genomes. We also generated and annotated a single-cell transcriptomic atlas for adult male H. symbiolongicarpus and identified cell type markers for all major cell types, including key i-cell markers. Orthology analyses based on the markers revealed that Hydractinia's i-cells are highly enriched in genes that are widely shared amongst animals, a striking finding given that Hydractinia has a higher proportion of phylum-specific genes than any of the other 41 animals in our orthology analysis. These results indicate that Hydractinia's stem cells and early progenitor cells may use a toolkit shared with all animals, making it a promising model organism for future exploration of stem cell biology and regenerative medicine. The genomic and transcriptomic resources for Hydractinia presented here will enable further studies of their regenerative capacity, colonial morphology, and ability to distinguish self from non-self.

Exome-wide assessment of isolated biliary atresia: A report from the National Birth Defects Prevention Study using child-parent trios and a case-control design to identify novel rare variants

The etiology of biliary atresia (BA) is unknown, but recent studies suggest a role for rare protein-altering variants (PAVs). Exome sequencing data from the National Birth Defects Prevention Study on 54 child-parent trios, one child-mother duo, and 1513 parents of children with other birth defects were analyzed. Most (91%) cases were isolated BA. We performed (1) a trio-based analysis to identify rare de novo, homozygous, and compound heterozygous PAVs and (2) a case-control analysis using a sequence kernel-based association test to identify genes enriched with rare PAVs. While we replicated previous findings on PKD1L1, our results do not suggest that recurrent de novo PAVs play important roles in BA susceptibility. In fact, our finding in NOTCH2, a disease gene associated with Alagille syndrome, highlights the difficulty in BA diagnosis. Notably, IFRD2 has been implicated in other gastrointestinal conditions and warrants additional study. Overall, our findings strengthen the hypothesis that the etiology of BA is complex.

Exome sequencing identifies genetic variants in anophthalmia and microphthalmia

Anophthalmia and microphthalmia (A/M) are rare birth defects affecting up to 2 per 10,000 live births. These conditions are manifested by the absence of an eye or reduced eye volumes within the orbit leading to vision loss. Although clinical case series suggest a strong genetic component in A/M, few systematic investigations have been conducted on potential genetic contributions owing to low population prevalence. To overcome this challenge, we utilized DNA samples and data collected as part of the National Birth Defects Prevention Study (NBDPS). The NBDPS employed multi-center ascertainment of infants affected by A/M. We performed exome sequencing on 67 family trios and identified numerous genes affected by rare deleterious nonsense and missense variants in this cohort, including de novo variants. We identified 9 nonsense changes and 86 missense variants that are absent from the reference human population (Genome Aggregation Database), and we suggest that these are high priority candidate genes for A/M. We also performed literature curation, single cell transcriptome comparisons, and molecular pathway analysis on the candidate genes and performed protein structure modeling to determine the potential pathogenic variant consequences on PAX6 in this disease.

Exome sequencing identifies variants in infants with sacral agenesis

BACKGROUND

Sacral agenesis (SA) consists of partial or complete absence of the caudal end of the spine and often presents with additional birth defects. Several studies have examined gene variants for syndromic forms of SA, but only one has examined exomes of children with non-syndromic SA.

METHODS

Using buccal cell specimens from families of children with non-syndromic SA, exomes of 28 child-parent trios (eight with and 20 without a maternal diagnosis of pregestational diabetes) and two child-father duos (neither with diagnosis of maternal pregestational diabetes) were exome sequenced.

RESULTS

Three children had heterozygous missense variants in ID1 (Inhibitor of DNA Binding 1), with CADD scores >20 (top 1% of deleterious variants in the genome); two children inherited the variant from their fathers and one from the child's mother. Rare missense variants were also detected in PDZD2 (PDZ Domain Containing 2; N = 1) and SPTBN5 (Spectrin Beta, Non-erythrocytic 5; N = 2), two genes previously suggested to be associated with SA etiology. Examination of variants with autosomal recessive and X-linked recessive inheritance identified five and two missense variants, respectively. Compound heterozygous variants were identified in several genes. In addition, 12 de novo variants were identified, all in different genes in different children.

CONCLUSIONS

To our knowledge, this is the first study reporting a possible association between ID1 and non-syndromic SA. Although maternal pregestational diabetes has been strongly associated with SA, the missense variants in ID1 identified in two of three children were paternally inherited. These findings add to the knowledge of gene variants associated with non-syndromic SA and provide data for future studies.

The complete sequence of a human genome

Since its initial release in 2000, the human reference genome has covered only the euchromatic fraction of the genome, leaving important heterochromatic regions unfinished. Addressing the remaining 8% of the genome, the Telomere-to-Telomere (T2T) Consortium presents a complete 3.055 billion-base pair sequence of a human genome, T2T-CHM13, that includes gapless assemblies for all chromosomes except Y, corrects errors in the prior references, and introduces nearly 200 million base pairs of sequence containing 1956 gene predictions, 99 of which are predicted to be protein coding. The completed regions include all centromeric satellite arrays, recent segmental duplications, and the short arms of all five acrocentric chromosomes, unlocking these complex regions of the genome to variational and functional studies.

KLF3 and PAX6 are candidate driver genes in late-stage, MSI-hypermutated endometrioid endometrial carcinomas

Endometrioid endometrial carcinomas (EECs) are the most common histological subtype of uterine cancer. Late-stage disease is an adverse prognosticator for EEC. The purpose of this study was to analyze EEC exome mutation data to identify late-stage-specific statistically significantly mutated genes (SMGs), which represent candidate driver genes potentially associated with disease progression. We exome sequenced 15 late-stage (stage III or IV) non-ultramutated EECs and paired non-tumor DNAs; somatic variants were called using Strelka, Shimmer, SomaticSniper and MuTect. Additionally, somatic mutation calls were extracted from The Cancer Genome Atlas (TCGA) data for 66 late-stage and 270 early-stage (stage I or II) non-ultramutated EECs. MutSigCV (v1.4) was used to annotate SMGs in the two late-stage cohorts and to derive p-values for all mutated genes in the early-stage cohort. To test whether late-stage SMGs are statistically significantly mutated in early-stage tumors, q-values for late-stage SMGs were re-calculated from the MutSigCV (v1.4) early-stage p-values, adjusting for the number of late-stage SMGs tested. We identified 14 SMGs in the combined late-stage EEC cohorts. When the 14 late-stage SMGs were examined in the TCGA early-stage data, only Krüppel-like factor 3 (KLF3) and Paired box 6 (PAX6) failed to reach significance as early-stage SMGs, despite the inclusion of enough early-stage cases to ensure adequate statistical power. Within TCGA, nonsynonymous mutations in KLF3 and PAX6 were, respectively, exclusive or nearly exclusive to the microsatellite instability (MSI)-hypermutated molecular subgroup and were dominated by insertions-deletions at homopolymer tracts. In conclusion, our findings are hypothesis-generating and suggest that KLF3 and PAX6, which encode transcription factors, are MSI target genes and late-stage-specific SMGs in EEC.

Low mutation burden and frequent loss of CDKN2A/B and SMARCA2, but not PRC2, define premalignant neurofibromatosis type 1-associated atypical neurofibromas

BACKGROUND

Neurofibromatosis type 1 (NF1) is a tumor-predisposition disorder caused by germline mutations in NF1. NF1 patients have an 8-16% lifetime risk of developing a malignant peripheral nerve sheath tumor (MPNST), a highly aggressive soft-tissue sarcoma, often arising from preexisting benign plexiform neurofibromas (PNs) and atypical neurofibromas (ANFs). ANFs are distinct from both PN and MPNST, representing an intermediate step in malignant transformation.

METHODS

In the first comprehensive genomic analysis of ANF originating from multiple patients, we performed tumor/normal whole-exome sequencing (WES) of 16 ANFs. In addition, we conducted WES of 3 MPNSTs, copy-number meta-analysis of 26 ANFs and 28 MPNSTs, and whole transcriptome sequencing analysis of 5 ANFs and 5 MPNSTs.

RESULTS

We identified a low number of mutations (median 1, range 0-5) in the exomes of ANFs (only NF1 somatic mutations were recurrent), and frequent deletions of CDKN2A/B (69%) and SMARCA2 (42%). We determined that polycomb repressor complex 2 (PRC2) genes EED and SUZ12 were frequently mutated, deleted, or downregulated in MPNSTs but not in ANFs. Our pilot gene expression study revealed upregulated NRAS, MDM2, CCND1/2/3, and CDK4/6 in ANFs and MPNSTs, and overexpression of EZH2 in MPNSTs only.

CONCLUSIONS

The PN-ANF transition is primarily driven by the deletion of CDKN2A/B. Further progression from ANF to MPNST likely involves broad chromosomal rearrangements and frequent inactivation of the PRC2 genes, loss of the DNA repair genes, and copy-number increase of signal transduction and cell-cycle and pluripotency self-renewal genes.

Clinical and genomic analysis of a large Chinese family with familial cortical myoclonic tremor with epilepsy and SAMD12 intronic repeat expansion

Objective

Our goal was to perform detailed clinical and genomic analysis of a large multigenerational Chinese family with 21 individuals showing symptoms of Familial Cortical Myoclonic Tremor with Epilepsy (FCMTE) that we have followed for over 20 years.

Methods

Patients were subjected to clinical evaluation, routine EEG, and structural magnetic resonance imaging. Whole exome sequencing, repeat-primed PCR, long-range PCR, and PacBio sequencing were performed to characterize the disease-causing mutation in this family.

Results

All evaluated patients manifested adult-onset seizures and presented with progressive myoclonic postural tremors starting after the third or fourth decade of life. Seizures typically diminished markedly in frequency with implementation of antiseizure medications but did not completely cease. The electroencephalogram of affected individuals showed generalized or multifocal spikes and slow wave complexes. An expansion of TTTTA motifs with addition of TTTCA motifs in intron 4 of SAMD12 was identified to segregate with the disease phenotype in this family. Furthermore, we found that the mutant allele is unstable and can undergo both contraction and expansion by changes in the number of repeat motifs each time it is passed to the next generation. The size of mutant allele varied from 5 to 5.5 kb with 549-603 copies of TTTTA and 287-343 copies of TTTCA repeat motifs in this family.

Significance

Our study provides a detailed description of clinical progression of FCMTE symptoms and its management with antiseizure medications. Our method of repeat analysis by PacBio sequencing of long-range PCR products does not require high-quality DNA and hence can be easily applied to other families to elucidate any correlation between the repeat size and phenotypic variables, such as, age of onset, and severity of symptoms.

Telomere-to-telomere assembly of a complete human X chromosome

After two decades of improvements, the current human reference genome (GRCh38) is the most accurate and complete vertebrate genome ever produced. However, no single chromosome has been finished end to end, and hundreds of unresolved gaps persist1,2. Here we present a human genome assembly that surpasses the continuity of GRCh382, along with a gapless, telomere-to-telomere assembly of a human chromosome. This was enabled by high-coverage, ultra-long-read nanopore sequencing of the complete hydatidiform mole CHM13 genome, combined with complementary technologies for quality improvement and validation. Focusing our efforts on the human X chromosome3, we reconstructed the centromeric satellite DNA array (approximately 3.1 Mb) and closed the 29 remaining gaps in the current reference, including new sequences from the human pseudoautosomal regions and from cancer-testis ampliconic gene families (CT-X and GAGE). These sequences will be integrated into future human reference genome releases. In addition, the complete chromosome X, combined with the ultra-long nanopore data, allowed us to map methylation patterns across complex tandem repeats and satellite arrays. Our results demonstrate that finishing the entire human genome is now within reach, and the data presented here will facilitate ongoing efforts to complete the other human chromosomes.

Comparative clinical and genomic analysis of neurofibromatosis type 2-associated cranial and spinal meningiomas

Neurofibromatosis type 2 (NF2) is an autosomal dominant Mendelian tumor predisposition disorder caused by germline pathogenic variants in the tumor suppressor NF2. Meningiomas are the second most common neoplasm in NF2, often occurring in multiple intracranial and spinal locations within the same patient. In this prospective longitudinal study, we assessed volumes and growth rates of ten spinal and ten cranial benign meningiomas in seven NF2 patients that concluded with surgical resection and performed whole-exome sequencing and copy-number variant (CNV) analysis of the tumors. Our comparison of the volume and the growth rate of NF2-associated spinal and cranial meningiomas point to the differences in timing of tumor initiation and/or to the differences in tumor progression (e.g., non-linear, saltatory growth) at these two anatomical locations. Genomic investigation of these tumors revealed that somatic inactivation of NF2 is the principal and perhaps the only driver of tumor initiation; and that tumor progression likely occurs via accumulation of CNVs, rather than point mutations. Results of this study contribute to a better understanding of NF2-associated meningiomas clinical behavior and their genetic underpinnings.

Antibody Lineages with Vaccine-Induced Antigen-Binding Hotspots Develop Broad HIV Neutralization

The vaccine-mediated elicitation of antibodies (Abs) capable of neutralizing diverse HIV-1 strains has been a long-standing goal. To understand how broadly neutralizing antibodies (bNAbs) can be elicited, we identified, characterized, and tracked five neutralizing Ab lineages targeting the HIV-1-fusion peptide (FP) in vaccinated macaques over time. Genetic and structural analyses revealed two of these lineages to belong to a reproducible class capable of neutralizing up to 59% of 208 diverse viral strains. B cell analysis indicated each of the five lineages to have been initiated and expanded by FP-carrier priming, with envelope (Env)-trimer boosts inducing cross-reactive neutralization. These Abs had binding-energy hotspots focused on FP, whereas several FP-directed Abs induced by immunization with Env trimer-only were less FP-focused and less broadly neutralizing. Priming with a conserved subregion, such as FP, can thus induce Abs with binding-energy hotspots coincident with the target subregion and capable of broad neutralization.

Admixture mapping identifies genetic regions associated with blood pressure phenotypes in African Americans

Hypertension occurs at a higher rate in African Americans than in European Americans. Based on the assumption that causal variants are more frequently found on DNA segments inherited from the ancestral population with higher disease risk, we employed admixture mapping to identify genetic loci with excess local African ancestry associated with blood pressure. Chromosomal regions 1q21.2–21.3, 4p15.1, 19q12 and 20p13 were significantly associated with diastolic blood pressure (β = 5.28, -7.94, -6.82 and 5.89, P-value = 6.39E-04, 2.07E-04, 6.56E-05 and 5.04E-04, respectively); 1q21.2–21.3 and 19q12 were also significantly associated with mean arterial pressure (β = 5.86 and -6.40, P-value = 5.32E-04 and 6.37E-04, respectively). We further selected SNPs that had large allele frequency differences within these regions and tested their association with blood pressure. SNP rs4815428 was significantly associated with diastolic blood pressure after Bonferroni correction (β = -2.42, P-value = 9.57E-04), and it partially explained the admixture mapping signal at 20p13. SNPs rs771205 (β = -1.99, P-value = 3.37E-03), rs3126067, rs2184953 and rs58001094 (the latter three exhibit strong linkage disequilibrium, β = -2.3, P-value = 1.4E-03) were identified to be significantly associated with mean arterial pressure, and together they fully explained the admixture signal at 1q21.2–21.3. Although no SNP at 4p15.1 showed large ancestral allele frequency differences in our dataset, we detected association at low-frequency African-specific variants that mapped predominantly to the gene PCDH7, which is most highly expressed in aorta. Our results suggest that these regions may harbor genetic variants that contribute to the different prevalence of hypertension.

Mutations that prevent caspase cleavage of RIPK1 cause autoinflammatory disease

Receptor Interacting Protein Kinase 1 (RIPK1) is a key regulator of innate immune signalling pathways. To ensure an optimal inflammatory response, RIPK1 is post-translationally regulated by well characterised ubiquitylation and phosphorylation events, as well as caspase-8 mediated cleavage^1–7. The physiological relevance of this cleavage remains unclear, though it is believed to inhibit activation of RIPK3 and necroptosis⁸. Here we show that heterozygous missense mutations p.D324N, p.D324H and p.D324Y prevent caspase cleavage of RIPK1 in humans and result in early-onset periodic fever episodes and severe intermittent lymphadenopathy, a condition we designate ‘Cleavage-resistant RIPK1-Induced Autoinflammatory’ (CRIA) syndrome. To define the mechanism for this disease we generated a cleavage-resistant Ripk1^D325A mutant mouse strain. While Ripk1^-/- mice die postnatally from systemic inflammation, Ripk1^D325A/D325A mice died during embryogenesis. Embryonic lethality was completely prevented by combined loss of Casp8 and Ripk3 but not by loss of Ripk3 or Mlkl alone. Loss of RIPK1 kinase activity also prevented Ripk1^D325A/D325A embryonic lethality, however the mice died before weaning from multi organ inflammation in a RIPK3 dependent manner. Consistently, Ripk1^D325A/D325A and Ripk1^D325A/+ cells were hypersensitive to RIPK3 dependent TNF-induced apoptosis and necroptosis. Heterozygous Ripk1^D325A/+ mice were viable and grossly normal, but were hyper-responsive to inflammatory stimuli in vivo. Our results demonstrate the importance of caspase-mediated RIPK1 cleavage during embryonic development and show that caspase cleavage of RIPK1 not only inhibits necroptosis but maintains inflammatory homeostasis throughout life.

Exome sequencing of family trios from the National Birth Defects Prevention Study: Tapping into a rich resource of genetic and environmental data

BACKGROUND

The National Birth Defects Prevention Study (NBDPS) is a multisite, population-based, case-control study of genetic and nongenetic risk factors for major structural birth defects. Eligible women had a pregnancy affected by a birth defect or a liveborn child without a birth defect between 1997 and 2011. They were invited to complete a telephone interview to collect pregnancy exposure data and were mailed buccal cell collection kits to collect specimens from themselves, their child (if living), and their child's father. Over 23,000 families representing more than 30 major structural birth defects provided DNA specimens.

METHODS

To evaluate their utility for exome sequencing (ES), specimens from 20 children with colonic atresia were studied. Evaluations were conducted on specimens collected using cytobrushes stored and transported in open versus closed packaging, on native genomic DNA (gDNA) versus whole genome amplified (WGA) products and on a library preparation protocol adapted to low amounts of DNA.

RESULTS

The DNA extracted from brushes in open packaging yielded higher quality sequence data than DNA from brushes in closed packaging. Quality metrics of sequenced gDNA were consistently higher than metrics from corresponding WGA products and were consistently high when using a low input protocol.

CONCLUSIONS

This proof-of-principle study established conditions under which ES can be applied to NBDPS specimens. Successful sequencing of exomes from well-characterized NBDPS families indicated that this unique collection can be used to investigate the roles of genetic variation and gene-environment interaction effects in birth defect etiologies, providing a valuable resource for birth defect researchers.

De novo assembly of the goldfish (Carassius auratus) genome and the evolution of genes after whole-genome duplication

For over a thousand years, the common goldfish (Carassius auratus) was raised throughout Asia for food and as an ornamental pet. As a very close relative of the common carp (Cyprinus carpio), goldfish share the recent genome duplication that occurred approximately 14 million years ago in their common ancestor. The combination of centuries of breeding and a wide array of interesting body morphologies provides an exciting opportunity to link genotype to phenotype and to understand the dynamics of genome evolution and speciation. We generated a high-quality draft sequence and gene annotations of a "Wakin" goldfish using 71X PacBio long reads. The two subgenomes in goldfish retained extensive synteny and collinearity between goldfish and zebrafish. However, genes were lost quickly after the carp whole-genome duplication, and the expression of 30% of the retained duplicated gene diverged substantially across seven tissues sampled. Loss of sequence identity and/or exons determined the divergence of the expression levels across all tissues, while loss of conserved noncoding elements determined expression variance between different tissues. This assembly provides an important resource for comparative genomics and understanding the causes of goldfish variants.

Single Cell Sequencing of the Pineal Gland: The Next Chapter

The analysis of pineal cell biology has undergone remarkable development as techniques have become available which allow for sequencing of entire transcriptomes and, most recently, the sequencing of the transcriptome of individual cells. Identification of at least nine distinct cell types in the rat pineal gland has been made possible, allowing identification of the precise cells of origin and expression of transcripts for the first time. Here the history and current state of knowledge generated by these transcriptomic efforts is reviewed, with emphasis on the insights suggested by the findings.

Brief Report: Whole-Exome Sequencing to Identify Rare Variants and Gene Networks That Increase Susceptibility to Scleroderma in African Americans

OBJECTIVE

Whole-exome sequencing (WES) studies in systemic sclerosis (SSc) patients of European American (EA) ancestry have identified variants in the ATP8B4 gene and enrichment of variants in genes in the extracellular matrix (ECM)-related pathway that increase SSc susceptibility. This study was undertaken to evaluate the association of the ATP8B4 gene and the ECM-related pathway with SSc in a cohort of African American (AA) patients.

METHODS

SSc patients of AA ancestry were enrolled from 23 academic centers across the US under the Genome Research in African American Scleroderma Patients consortium. Unrelated AA individuals without serologic evidence of autoimmunity who were enrolled in the Howard University Family Study were used as unaffected controls. Functional variants in genes reported in the 2 WES studies in EA patients with SSc were selected for gene association testing using the optimized sequence kernel association test (SKAT-O) and pathway analysis by Ingenuity Pathway Analysis in 379 patients and 411 controls.

RESULTS

Principal components analysis demonstrated that the patients and controls had similar ancestral backgrounds, with roughly equal proportions of mean European admixture. Using SKAT-O, we examined the association of individual genes that were previously reported in EA patients and none remained significant, including ATP8B4 (P = 0.98). However, we confirmed the previously reported association of the ECM-related pathway with enrichment of variants within the COL13A1, COL18A1, COL22A1, COL4A3, COL4A4, COL5A2, PROK1, and SERPINE1 genes (corrected P = 1.95 × 10^-4 ).

CONCLUSION

In the largest genetic study in AA patients with SSc to date, our findings corroborate the role of functional variants that aggregate in a fibrotic pathway and increase SSc susceptibility.

Common genetic causes of holoprosencephaly are limited to a small set of evolutionarily conserved driver genes of midline development coordinated by TGF-β, hedgehog, and FGF signaling

Here, we applied targeted capture to examine 153 genes representative of all the major vertebrate developmental pathways among 333 probands to rank their relative significance as causes for holoprosencephaly (HPE). We now show that comparisons of variant transmission versus nontransmission among 136 HPE Trios indicates some reported genes now lack confirmation, while novel genes are implicated. Furthermore, we demonstrate that variation of modest intrinsic effect can synergize with these driver mutations as gene modifiers.

Joubert Syndrome: Ophthalmological Findings in Correlation with Genotype and Hepatorenal Disease in 99 Patients Prospectively Evaluated at a Single Center

PURPOSE

Joubert syndrome (JS) is caused by mutations in >34 genes that encode proteins involved with primary (nonmotile) cilia and the cilium basal body. This study describes the varying ocular phenotypes in JS patients, with correlation to systemic findings and genotype.

DESIGN

Patients were systematically and prospectively examined at the National Institutes of Health (NIH) Clinical Center in the setting of a dedicated natural history clinical trial.

PARTICIPANTS

Ninety-nine patients with JS examined at a single center.

METHODS

All patients underwent genotyping for JS, followed by complete age-appropriate ophthalmic examinations at the NIH Clinical Center, including visual acuity (VA), fixation behavior, lid position, motility assessment, slit-lamp biomicroscopy, dilated fundus examination with an indirect ophthalmoscope, and retinoscopy. Color and fundus autofluorescence imaging, Optos wide-field photography (Dunfermline, Scotland, UK), and electroretinography (ERG) were performed when possible.

MAIN OUTCOME MEASURES

The VA (with longitudinal follow-up where possible), ptosis, extraocular muscle function, retinal and optic nerve status, and retinal function as measured by ERG.

RESULTS

Among patients with JS with quantifiable VA (68/99), values ranged from 0 logarithm of the minimum angle of resolution (logMAR) (Snellen 20/20) to 1.5 logMAR (Snellen 20/632). Strabismus (71/98), nystagmus (66/99), oculomotor apraxia (60/77), ptosis (30/98), coloboma (28/99), retinal degeneration (20/83), and optic nerve atrophy (8/86) were identified.

CONCLUSIONS

We recommend regular monitoring for ophthalmological manifestations of JS beginning soon after birth or diagnosis. We demonstrate delayed visual development and note that the amblyogenic time frame may last significantly longer in JS than is typical. In general, patients with coloboma were less likely to display retinal degeneration, and those with retinal degeneration did not have coloboma. Severe retinal degeneration that is early and aggressive is seen in disease caused by specific genes, such as CEP290- and AHI1-associated JS. Retinal degeneration in INPP5E-, MKS1-, and NPHP1-associated JS was generally milder. Finally, ptosis surgery can be helpful in a subset of patients with JS; decisions as to timing and benefit/risk ratio need to be made on an individual basis according to expert consultation.

Author Correction: Applications and efficiencies of the first cat 63K DNA array

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has been fixed in the paper.

Applications and efficiencies of the first cat 63K DNA array

The development of high throughput SNP genotyping technologies has improved the genetic dissection of simple and complex traits in many species including cats. The properties of feline 62,897 SNPs Illumina Infinium iSelect DNA array are described using a dataset of over 2,000 feline samples, the most extensive to date, representing 41 cat breeds, a random bred population, and four wild felid species. Accuracy and efficiency of the array’s genotypes and its utility in performing population-based analyses were evaluated. Average marker distance across the array was 37,741 Kb, and across the dataset, only 1% (625) of the markers exhibited poor genotyping and only 0.35% (221) showed Mendelian errors. Marker polymorphism varied across cat breeds and the average minor allele frequency (MAF) of all markers across domestic cats was 0.21. Population structure analysis confirmed a Western to Eastern structural continuum of cat breeds. Genome-wide linkage disequilibrium ranged from 50–1,500 Kb for domestic cats and 750 Kb for European wildcats (Felis silvestris silvestris). Array use in trait association mapping was investigated under different modes of inheritance, selection and population sizes. The efficient array design and cat genotype dataset continues to advance the understanding of cat breeds and will support monogenic health studies across feline breeds and populations.

A Neutralizing Antibody Recognizing Primarily N-Linked Glycan Targets the Silent Face of the HIV Envelope

Virtually the entire surface of the HIV-1-envelope trimer is recognized by neutralizing antibodies, except for a highly glycosylated region at the center of the "silent face" on the gp120 subunit. From an HIV-1-infected donor, #74, we identified antibody VRC-PG05, which neutralized 27% of HIV-1 strains. The crystal structure of the antigen-binding fragment of VRC-PG05 in complex with gp120 revealed an epitope comprised primarily of N-linked glycans from N262, N295, and N448 at the silent face center. Somatic hypermutation occurred preferentially at antibody residues that interacted with these glycans, suggesting somatic development of glycan recognition. Resistance to VRC-PG05 in donor #74 involved shifting of glycan-N448 to N446 or mutation of glycan-proximal residue E293. HIV-1 neutralization can thus be achieved at the silent face center by glycan-recognizing antibody; along with other known epitopes, the VRC-PG05 epitope completes coverage by neutralizing antibody of all major exposed regions of the prefusion closed trimer.

A comprehensive approach to identification of pathogenic FANCA variants in Fanconi anemia patients and their families

Fanconi anemia (FA) is a rare recessive DNA repair deficiency resulting from mutations in one of at least 22 genes. Two-thirds of FA families harbor mutations in FANCA. To genotype patients in the International Fanconi Anemia Registry (IFAR) we employed multiple methodologies, screening 216 families for FANCA mutations. We describe identification of 57 large deletions and 261 sequence variants, in 159 families. All but seven families harbored distinct combinations of two mutations demonstrating high heterogeneity. Pathogenicity of the 18 novel missense variants was analyzed functionally by determining the ability of the mutant cDNA to improve the survival of a FANCA-null cell line when treated with MMC. Overexpressed pathogenic missense variants were found to reside in the cytoplasm, and nonpathogenic in the nucleus. RNA analysis demonstrated that two variants (c.522G > C and c.1565A > G), predicted to encode missense variants, which were determined to be nonpathogenic by a functional assay, caused skipping of exons 5 and 16, respectively, and are most likely pathogenic. We report 48 novel FANCA sequence variants. Defining both variants in a large patient cohort is a major step toward cataloging all FANCA variants, and permitting studies of genotype-phenotype correlations.

Prospective Evaluation of Kidney Disease in Joubert Syndrome

BACKGROUND AND OBJECTIVES

Joubert syndrome is a genetically heterogeneous ciliopathy associated with >30 genes. The characteristics of kidney disease and genotype-phenotype correlations have not been evaluated in a large cohort at a single center.

DESIGN, SETTING, PARTICIPANTS, & MEASUREMENTS

We evaluated 97 individuals with Joubert syndrome at the National Institutes of Health Clinical Center using abdominal ultrasonography, blood and urine chemistries, and DNA sequencing.

RESULTS

Patients were ages 0.6-36 years old (mean of 9.0±7.6 years old); 41 were female. Mutations were identified in 19 genes in 92 patients; two thirds of the mutations resided in six genes: TMEM67, C5orf42, CC2D2A, CEP290, AHI1, and KIAA0586. Kidney disease was detected in 30%, most commonly in association with the following genes: CEP290 (six of six), TMEM67 (11 of 22), and AHI1 (three of six). No kidney disease was identified in patients with mutations in C5orf42 (zero of 15) or KIAA0586 (zero of six). Prenatal ultrasonography of kidneys was normal in 72% of patients with kidney disease. Specific types of kidney disease included nephronophthisis (31%), an overlap phenotype of autosomal recessive polycystic kidney disease/nephronophthisis (35%), unilateral multicystic dysplastic kidney (10%), and indeterminate-type cystic kidney disease (24%). Early-onset hypertension occurred in 24% of patients with kidney disease. Age at ESRD (n=13) ranged from 6 to 24 years old (mean of 11.3±4.8 years old).

CONCLUSIONS

Kidney disease occurs in up to one third of patients with Joubert syndrome, most commonly in those with mutations in CEP290, TMEM67, and AHI1. Patients with mutations in C5orf42 or KIAA0586 are less likely to develop kidney disease. Prenatal ultrasonography is a poor predictor of kidney involvement in Joubert syndrome. Unilateral multicystic dysplastic kidney and autosomal recessive polycystic kidney disease-like enlarged kidneys with early-onset hypertension can be part of the Joubert syndrome kidney phenotype.

The FOXA2 transcription factor is frequently somatically mutated in uterine carcinosarcomas and carcinomas

BACKGROUND

Uterine carcinosarcomas (UCSs) are a rare but clinically aggressive form of cancer. They are biphasic tumors consisting of both epithelial and sarcomatous components. The majority of uterine carcinosarcomas are clonal, with the carcinomatous cells undergoing metaplasia to give rise to the sarcomatous component. The objective of the current study was to identify novel somatically mutated genes in UCSs.

METHODS

We whole exome sequenced paired tumor and nontumor DNAs from 14 UCSs and orthogonally validated 464 somatic variants using Sanger sequencing. Fifteen genes that were somatically mutated in at least 2 tumor exomes were Sanger sequenced in another 39 primary UCSs.

RESULTS

Overall, among 53 UCSs in the current study, the most frequently mutated of these 15 genes were tumor protein p53 (TP53) (75.5%), phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit alpha (PIK3CA) (34.0%), protein phosphatase 2, regulatory subunit A, alpha (PPP2R1A) (18.9%), F-box and WD repeat domain containing 7 (FBXW7) (18.9%), chromodomain helicase DNA binding protein 4 (CHD4) (17.0%), and forkhead box A2 (FOXA2) (15.1%). FOXA2 has not previously been implicated in UCSs and was predominated by frameshift and nonsense mutations. One UCS with a FOXA2 frameshift mutation expressed truncated FOXA2 protein by immunoblotting. Sequencing of FOXA2 in 160 primary endometrial carcinomas revealed somatic mutations in 5.7% of serous, 22.7% of clear cell, 9% of endometrioid, and 11.1% of mixed endometrial carcinomas, the majority of which were frameshift mutations.

CONCLUSIONS

Collectively, the findings of the current study provide compelling genetic evidence that FOXA2 is a pathogenic driver gene in the etiology of primary uterine cancers, including UCSs. Cancer 2018;124:65-73. © 2017 American Cancer Society.

Assessing the spectrum of germline variation in Fanconi anemia genes among patients with head and neck carcinoma before age 50

BACKGROUND

Patients with Fanconi anemia (FA) have an increased risk for head and neck squamous cell carcinoma (HNSCC). The authors sought to determine the prevalence of undiagnosed FA and FA carriers among patients with HNSCC as well as an age cutoff for FA genetic screening.

METHODS

Germline DNA samples from 417 patients with HNSCC aged <50 years were screened for sequence variants by targeted next-generation sequencing of the entire length of 16 FA genes.

RESULTS

The sequence revealed 194 FA gene variants in 185 patients (44%). The variant spectrum was comprised of 183 nonsynonymous point mutations, 9 indels, 1 large deletion, and 1 synonymous variant that was predicted to effect splicing. One hundred eight patients (26%) had at least 1 rare variant that was predicted to be damaging, and 57 (14%) had at least 1 rare variant that was predicted to be damaging and had been previously reported. Fifteen patients carried 2 rare variants or an X-linked variant in an FA gene. Overall, an age cutoff for FA screening was not identified among young patients with HNSCC, because there were no significant differences in mutation rates when patients were stratified by age, tumor site, ethnicity, smoking status, or human papillomavirus status. However, an increased burden, or mutation load, of FA gene variants was observed in carriers of the genes FA complementation group D2 (FANCD2), FANCE, and FANCL in the HNSCC patient cohort relative to the 1000 Genomes population.

CONCLUSIONS

FA germline functional variants offer a novel area of study in HNSCC tumorigenesis. FANCE and FANCL, which are components of the core complex, are known to be responsible for the recruitment and ubiquitination, respectively, of FANCD2, a critical step in the FA DNA repair pathway. In the current cohort, the increased mutation load of FANCD2, FANCE, and FANCL variants among younger patients with HNSCC indicates the importance of the FA pathway in HNSCC. Cancer 2017;123:3943-54. © 2017 American Cancer Society.

Gene-Specific Substitution Profiles Describe the Types and Frequencies of Amino Acid Changes during Antibody Somatic Hypermutation

Somatic hypermutation (SHM) plays a critical role in the maturation of antibodies, optimizing recognition initiated by recombination of V(D)J genes. Previous studies have shown that the propensity to mutate is modulated by the context of surrounding nucleotides and that SHM machinery generates biased substitutions. To investigate the intrinsic mutation frequency and substitution bias of SHMs at the amino acid level, we analyzed functional human antibody repertoires and developed mGSSP (method for gene-specific substitution profile), a method to construct amino acid substitution profiles from next-generation sequencing-determined B cell transcripts. We demonstrated that these gene-specific substitution profiles (GSSPs) are unique to each V gene and highly consistent between donors. We also showed that the GSSPs constructed from functional antibody repertoires are highly similar to those constructed from antibody sequences amplified from non-productively rearranged passenger alleles, which do not undergo functional selection. This suggests the types and frequencies, or mutational space, of a majority of amino acid changes sampled by the SHM machinery to be well captured by GSSPs. We further observed the rates of mutational exchange between some amino acids to be both asymmetric and context dependent and to correlate weakly with their biochemical properties. GSSPs provide an improved, position-dependent alternative to standard substitution matrices, and can be utilized to developing software for accurately modeling the SHM process. GSSPs can also be used for predicting the amino acid mutational space available for antigen-driven selection and for understanding factors modulating the maturation pathways of antibody lineages in a gene-specific context. The mGSSP method can be used to build, compare, and plot GSSPs¹; we report the GSSPs constructed for 69 common human V genes (DOI: 10.6084/m9.figshare.3511083) and provide high-resolution logo plots for each (DOI: 10.6084/m9.figshare.3511085).

Somatic mutation profiles of clear cell endometrial tumors revealed by whole exome and targeted gene sequencing

BACKGROUND

The molecular pathogenesis of clear cell endometrial cancer (CCEC), a tumor type with a relatively unfavorable prognosis, is not well defined. We searched exome-wide for novel somatically mutated genes in CCEC and assessed the mutational spectrum of known and candidate driver genes in a large cohort of cases.

METHODS

We conducted whole exome sequencing of paired tumor-normal DNAs from 16 cases of CCEC (12 CCECs and the CCEC components of 4 mixed histology tumors). Twenty-two genes-of-interest were Sanger-sequenced from another 47 cases of CCEC. Microsatellite instability (MSI) and microsatellite stability (MSS) were determined by genotyping 5 mononucleotide repeats.

RESULTS

Two tumor exomes had relatively high mutational loads and MSI. The other 14 tumor exomes were MSS and had 236 validated nonsynonymous or splice junction somatic mutations among 222 protein-encoding genes. Among the 63 cases of CCEC in this study, we identified frequent somatic mutations in TP53 (39.7%), PIK3CA (23.8%), PIK3R1 (15.9%), ARID1A (15.9%), PPP2R1A (15.9%), SPOP (14.3%), and TAF1 (9.5%), as well as MSI (11.3%). Five of 8 mutations in TAF1, a gene with no known role in CCEC, localized to the putative histone acetyltransferase domain and included 2 recurrently mutated residues. Based on patterns of MSI and mutations in 7 genes, CCEC subsets molecularly resembled serous endometrial cancer (SEC) or endometrioid endometrial cancer (EEC).

CONCLUSION

Our findings demonstrate molecular similarities between CCEC and SEC and EEC and implicate TAF1 as a novel candidate CCEC driver gene. Cancer 2017;123:3261-8. © 2017 American Cancer Society.

Clinical and molecular phenotyping of a child with Hermansky-Pudlak syndrome-7, an uncommon genetic type of HPS

PURPOSE

Hermansky-Pudlak syndrome (HPS) is a rare inherited disorder with ten reported genetic types; each type has defects in subunits of either Adaptor Protein-3 complex or Biogenesis of Lysosome-related Organelles Complex (BLOC)-1, -2, or -3. Very few patients with BLOC-1 deficiency (HPS-7, -8, and -9 types) have been diagnosed. We report results of comprehensive clinical testing and molecular analyses of primary fibroblasts from a new case of HPS-7.

RESULTS

A 6-year old Paraguayan male presented with hypopigmentation, ocular albinism, nystagmus, reduced visual acuity, and easy bruising. He also experienced delayed motor and language development as a very young child; head and chest trauma resulted in intracranial hemorrhage with subsequent right hemiparesis and lung scarring. There was no clinical evidence of immunodeficiency or colitis. Whole mount transmission electron microscopy revealed absent platelet delta granules; platelet aggregation testing was abnormal. Exome sequencing revealed a homozygous nonsense mutation in the Dystrobrevin binding protein 1 (DTNBP1) gene [NM_032122.4: c.307C>T; p.Gln103*], previously reported in a Portuguese adult. The gene encodes the dysbindin subunit of BLOC-1. Dysbindin protein expression was negligible in our patient's dermal fibroblasts, while his DTNBP1 mRNA level was similar to that of a normal control.

CONCLUSIONS

Comprehensive clinical evaluation of the first pediatric case reported with HPS-7 reveals oculocutaneous albinism and platelet storage pool deficiency; his phenotype is consistent with findings in other patients with BLOC-1 disorders. This patient's markedly reduced Dysbindin protein expression in HPS-7 resulted from a mechanism other than nonsense mediated decay.

Mutations in KIAA0753 cause Joubert syndrome associated with growth hormone deficiency

Joubert syndrome and related disorders (JSRD) are a heterogeneous group of ciliopathies defined based on the mid-hindbrain abnormalities that result in the characteristic "molar tooth sign" on brain imaging. The core clinical findings of JSRD are hypotonia, developmental delay, abnormal eye movements and breathing abnormalities. To date, more than 30 JSRD genes that encode proteins important for structure and/or function of cilia have been identified. Here, we present 2 siblings with Joubert syndrome associated with growth hormone deficiency. Whole exome sequencing of the family identified compound heterozygous mutations in KIAA0753, i.e., a missense mutation (p.Arg257Gly) and an intronic mutation (c.2359-1G>C). The intronic mutation alters normal splicing by activating a cryptic acceptor splice site in exon 16. The novel acceptor site skips nine nucleotides, deleting three amino acids from the protein coding frame. KIAA0753 (OFIP) is a centrosome and pericentriolar satellite protein, previously not known to cause Joubert syndrome. We present comprehensive clinical descriptions of the Joubert syndrome patients as well as the cellular phenotype of defective ciliogenesis in the patients' fibroblasts.

First insight into the somatic mutation burden of neurofibromatosis type 2-associated grade I and grade II meningiomas: a case report comprehensive genomic study of two cranial meningiomas with vastly different clinical presentation

Background

Neurofibromatosis type 2 (NF2) is a rare autosomal dominant nervous system tumor predisposition disorder caused by constitutive inactivation of one of the two copies of NF2. Meningiomas affect about one half of NF2 patients, and are associated with a higher disease burden. Currently, the somatic mutation landscape in NF2-associated meningiomas remains largely unexamined.

Case presentation

Here, we present an in-depth genomic study of benign and atypical meningiomas, both from a single NF2 patient. While the grade I tumor was asymptomatic, the grade II tumor exhibited an unusually high growth rate: expanding to 335 times its initial volume within one year. The genomes of both tumors were examined by whole-exome sequencing (WES) complemented with spectral karyotyping (SKY) and SNP-array copy-number analyses. To better understand the clonal composition of the atypical meningioma, the tumor was divided in four sections and each section was investigated independently. Both tumors had second copy inactivation of NF2, confirming the central role of the gene in meningioma formation. The genome of the benign tumor closely resembled that of a normal diploid cell and had only one other deleterious mutation (EPHB3). In contrast, the chromosomal architecture of the grade II tumor was highly re-arranged, yet uniform among all analyzed fragments, implying that this large and fast growing tumor was composed of relatively few clones. Besides multiple gains and losses, the grade II meningioma harbored numerous chromosomal translocations. WES analysis of the atypical tumor identified deleterious mutations in two genes: ADAMTSL3 and CAPN5 in all fragments, indicating that the mutations were present in the cell undergoing fast clonal expansion

Conclusions

This is the first WES study of NF2-associated meningiomas. Besides second NF2 copy inactivation, we found low somatic burden in both tumors and high level of genomic instability in the atypical meningioma. Genomic instability resulting in altered gene dosage and compromised structural integrity of multiple genes may be the primary reason of the high growth rate for the grade II tumor. Further study of ADAMTSL3 and CAPN5 may lead to elucidation of their molecular implications in meningioma pathogenesis.

Electronic supplementary material

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

Biallelic Mutations in DNAJC12 Cause Hyperphenylalaninemia, Dystonia, and Intellectual Disability

Phenylketonuria (PKU, phenylalanine hydroxylase deficiency), an inborn error of metabolism, can be detected through newborn screening for hyperphenylalaninemia (HPA). Most individuals with HPA harbor mutations in the gene encoding phenylalanine hydroxylase (PAH), and a small proportion (2%) exhibit tetrahydrobiopterin (BH4) deficiency with additional neurotransmitter (dopamine and serotonin) deficiency. Here we report six individuals from four unrelated families with HPA who exhibited progressive neurodevelopmental delay, dystonia, and a unique profile of neurotransmitter deficiencies without mutations in PAH or BH4 metabolism disorder-related genes. In these six affected individuals, whole-exome sequencing (WES) identified biallelic mutations in DNAJC12, which encodes a heat shock co-chaperone family member that interacts with phenylalanine, tyrosine, and tryptophan hydroxylases catalyzing the BH4-activated conversion of phenylalanine into tyrosine, tyrosine into L-dopa (the precursor of dopamine), and tryptophan into 5-hydroxytryptophan (the precursor of serotonin), respectively. DNAJC12 was undetectable in fibroblasts from the individuals with null mutations. PAH enzyme activity was reduced in the presence of DNAJC12 mutations. Early treatment with BH4 and/or neurotransmitter precursors had dramatic beneficial effects and resulted in the prevention of neurodevelopmental delay in the one individual treated before symptom onset. Thus, DNAJC12 deficiency is a preventable and treatable cause of intellectual disability that should be considered in the early differential diagnosis when screening results are positive for HPA. Sequencing of DNAJC12 may resolve any uncertainty and should be considered in all children with unresolved HPA.

CELSR2, encoding a planar cell polarity protein, is a putative gene in Joubert syndrome with cortical heterotopia, microophthalmia, and growth hormone deficiency

Joubert syndrome is a ciliopathy characterized by a specific constellation of central nervous system malformations that result in the pathognomonic "molar tooth sign" on imaging. More than 27 genes are associated with Joubert syndrome, but some patients do not have mutations in any of these genes. Celsr1, Celsr2, and Celsr3 are the mammalian orthologues of the drosophila planar cell polarity protein, flamingo; they play important roles in neural development, including axon guidance, neuronal migration, and cilium polarity. Here, we report bi-allelic mutations in CELSR2 in a Joubert patient with cortical heterotopia, microophthalmia, and growth hormone deficiency. © 2017 Wiley Periodicals, Inc.

Phospholipase A2-activating protein is associated with a novel form of leukoencephalopathy

Leukoencephalopathies are a group of white matter disorders related to abnormal formation, maintenance, and turnover of myelin in the central nervous system. These disorders of the brain are categorized according to neuroradiological and pathophysiological criteria. Herein, we have identified a unique form of leukoencephalopathy in seven patients presenting at ages 2 to 4 months with progressive microcephaly, spastic quadriparesis, and global developmental delay. Clinical, metabolic, and imaging characterization of seven patients followed by homozygosity mapping and linkage analysis were performed. Next generation sequencing, bioinformatics, and segregation analyses followed, to determine a loss of function sequence variation in the phospholipase A₂-activating protein encoding gene (PLAA). Expression and functional studies of the encoded protein were performed and included measurement of prostaglandin E₂ and cytosolic phospholipase A₂ activity in membrane fractions of fibroblasts derived from patients and healthy controls. Plaa-null mice were generated and prostaglandin E₂ levels were measured in different tissues. The novel phenotype of our patients segregated with a homozygous loss-of-function sequence variant, causing the substitution of leucine at position 752 to phenylalanine, in PLAA, which causes disruption of the protein's ability to induce prostaglandin E₂ and cytosolic phospholipase A₂ synthesis in patients' fibroblasts. Plaa-null mice were perinatal lethal with reduced brain levels of prostaglandin E₂ The non-functional phospholipase A₂-activating protein and the associated neurological phenotype, reported herein for the first time, join other complex phospholipid defects that cause leukoencephalopathies in humans, emphasizing the importance of this axis in white matter development and maintenance.

Biallelic SCN10A mutations in neuromuscular disease and epileptic encephalopathy

Objectives

Two consanguineous families, one of Sudanese ethnicity presenting progressive neuromuscular disease, severe cognitive impairment, muscle weakness, upper motor neuron lesion, anhydrosis, facial dysmorphism, and recurrent seizures and the other of Egyptian ethnicity presenting with neonatal hypotonia, bradycardia, and recurrent seizures, were evaluated for the causative gene mutation.

Methods and Results

Homozygosity mapping and whole exome sequencing (WES) identified damaging homozygous variants in SCN10A, namely c.4514C>T; p.Thr1505Met in the first family and c.4735C>T; p.Arg1579* in the second family. A third family, of Western European descent, included a child with febrile infection‐related epilepsy syndrome (FIRES) who also had compound heterozygous missense mutations in SCN10A, namely, c.3482T>C; p.Met1161Thr and c.4709C>A; p.Thr1570Lys. A search for SCN10A variants in three consortia datasets (EuroEPINOMICS, Epi4K/EPGP, Autism/dbGaP) identified an additional five individuals with compound heterozygous variants. A Hispanic male with infantile spasms [c.2842G>C; p.Val948Leu and c.1453C>T; p.Arg485Cys], and a Caucasian female with Lennox–Gastaut syndrome [c.1529C>T; p.Pro510Leu and c.4984G>A; p.Gly1662Ser] in the epilepsy databases and three in the autism databases with [c.4009T>A; p.Ser1337Thr and c.1141A>G; p.Ile381Val], [c.2972C>T; p.Pro991Leu and c.2470C>T; p.His824Tyr], and [c.4009T>A; p.Ser1337Thr and c.2052G>A; p.Met684Ile].

Interpretation

SCN10A is a member of the voltage‐gated sodium channel (VGSC) gene family. Sodium channels are responsible for the instigation and proliferation of action potentials in central and peripheral nervous systems. Heterozygous mutations in VGSC genes cause a wide range of epileptic and peripheral nervous system disorders. This report presents autosomal recessive mutations in SCN10A that may be linked to epilepsy‐related phenotypes, Lennox–Gastaut syndrome, infantile spasms, and Autism Spectrum Disorder.

Mapping Complex Traits in a Diversity Outbred F1 Mouse Population Identifies Germline Modifiers of Metastasis in Human Prostate Cancer

It is unclear how standing genetic variation affects the prognosis of prostate cancer patients. To provide one controlled answer to this problem, we crossed a dominant, penetrant mouse model of prostate cancer to Diversity Outbred mice, a collection of animals that carries over 40 million SNPs. Integration of disease phenotype and SNP variation data in 493 F1 males identified a metastasis modifier locus on Chromosome 8 (LOD = 8.42); further analysis identified the genes Rwdd4, Cenpu, and Casp3 as functional effectors of this locus. Accordingly, analysis of over 5,300 prostate cancer patient samples revealed correlations between the presence of genetic variants at these loci, their expression levels, cancer aggressiveness, and patient survival. We also observed that ectopic overexpression of RWDD4 and CENPU increased the aggressiveness of two human prostate cancer cell lines. In aggregate, our approach demonstrates how well-characterized genetic variation in mice can be harnessed in conjunction with systems genetics approaches to identify and characterize germline modifiers of human disease processes.

Stringent comparative sequence analysis reveals SOX10 as a putative inhibitor of glial cell differentiation

Background

The transcription factor SOX10 is essential for all stages of Schwann cell development including myelination. SOX10 cooperates with other transcription factors to activate the expression of key myelin genes in Schwann cells and is therefore a context-dependent, pro-myelination transcription factor. As such, the identification of genes regulated by SOX10 will provide insight into Schwann cell biology and related diseases. While genome-wide studies have successfully revealed SOX10 target genes, these efforts mainly focused on myelinating stages of Schwann cell development. We propose that less-biased approaches will reveal novel functions of SOX10 outside of myelination.

Results

We developed a stringent, computational-based screen for genome-wide identification of SOX10 response elements. Experimental validation of a pilot set of predicted binding sites in multiple systems revealed that SOX10 directly regulates a previously unreported alternative promoter at SOX6, which encodes a transcription factor that inhibits glial cell differentiation. We further explored the utility of our computational approach by combining it with DNase-seq analysis in cultured Schwann cells and previously published SOX10 ChIP-seq data from rat sciatic nerve. Remarkably, this analysis enriched for genomic segments that map to loci involved in the negative regulation of gliogenesis including SOX5, SOX6, NOTCH1, HMGA2, HES1, MYCN, ID4, and ID2. Functional studies in Schwann cells revealed that: (1) all eight loci are expressed prior to myelination and down-regulated subsequent to myelination; (2) seven of the eight loci harbor validated SOX10 binding sites; and (3) seven of the eight loci are down-regulated upon repressing SOX10 function.

Conclusions

Our computational strategy revealed a putative novel function for SOX10 in Schwann cells, which suggests a model where SOX10 activates the expression of genes that inhibit myelination during non-myelinating stages of Schwann cell development. Importantly, the computational and functional datasets we present here will be valuable for the study of transcriptional regulation, SOX protein function, and glial cell biology.

Electronic supplementary material

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

Identification of a novel mutation in HPS6 in a patient with hemophilia B and oculocutaneous albinism

PURPOSE

Hemophilia B, an X-linked disease, manifests with recurrent soft tissue bleeding episodes. Hermansky-Pudlak syndrome, a rare autosomal recessive disorder, is characterized by oculocutaneous albinism and an increased tendency to bleed due to a platelet storage pool defect. We report a novel mutation in HPS6 in a Caucasian man with hemophilia B and oculocutaneous albinism.

RESULTS

The patient was diagnosed with hemophilia B at age 4months due to recurrent soft tissue bleeding episodes, and he was also diagnosed with Hermansky-Pudlak syndrome at 32years of age due to unexplained oculocutaneous albinism. His factor IX level was markedly reduced at 13%; whole exome and Sanger sequencing showed the Durham mutation in F9 (NM_000133.3). The diagnosis of Hermansky-Pudlak syndrome subtype 6 was established by demonstrating absence of platelet delta granules on whole mount electron microscopy, an abnormal secondary wave in platelet aggregation studies, and a novel homozygous c.1114 C>T (p.Arg372*) mutation in HPS6 (NM_024747.5) on exome analysis and Sanger sequencing. Clinical phenotyping revealed no evidence of recurrent or unusual infections, interstitial lung disease or pulmonary fibrosis, or neurological disorders. The patient was treated with fresh frozen plasma, recombinant factor IX, and aminocaproic acid. Treatment with desmopressin was added to his regimen after he was diagnosed with Hermansky-Pudlak syndrome. Treatment of bleeding episodes results in effective hemostasis, and the patient has not required platelet or blood product transfusions.

CONCLUSIONS

This report highlights the need to consider Hermansky-Pudlak syndrome as an etiology of oculocutaneous albinism even in patients with known hematologic disorders associated with bleeding. Identification of a novel mutation in HPS6 in an individual with hemophilia B shows that, although quite rare, patients may be diagnosed with two independent inherited bleeding disorders. No evidence of lung disease was found in this adult patient with Hermansky-Pudlak syndrome subtype 6.

Alternative Isoform Analysis of Ttc8 Expression in the Rat Pineal Gland Using a Multi-Platform Sequencing Approach Reveals Neural Regulation

Alternative isoform regulation (AIR) vastly increases transcriptome diversity and plays an important role in numerous biological processes and pathologies. However, the detection and analysis of isoform-level differential regulation is difficult, particularly in the face of complex and incompletely-annotated transcriptomes. Here we have used Illumina short-read/high-throughput RNA-Seq to identify 55 genes that exhibit neurally-regulated AIR in the pineal gland, and then used two other complementary experimental platforms to further study and characterize the Ttc8 gene, which is involved in Bardet-Biedl syndrome and non-syndromic retinitis pigmentosa. Use of the JunctionSeq analysis tool led to the detection of several novel exons and splice junctions in this gene, including two novel alternative transcription start sites which were found to display disproportionately strong neurally-regulated differential expression in several independent experiments. These high-throughput sequencing results were validated and augmented via targeted qPCR and long-read Pacific Biosciences SMRT sequencing. We confirmed the existence of numerous novel splice junctions and the selective upregulation of the two novel start sites. In addition, we identified more than 20 novel isoforms of the Ttc8 gene that are co-expressed in this tissue. By using information from multiple independent platforms we not only greatly reduce the risk of errors, biases, and artifacts influencing our results, we also are able to characterize the regulation and splicing of the Ttc8 gene more deeply and more precisely than would be possible via any single platform. The hybrid method outlined here represents a powerful strategy in the study of the transcriptome.

The cnidarian Hydractinia echinata employs canonical and highly adapted histones to pack its DNA

Background

Cnidarians are a group of early branching animals including corals, jellyfish and hydroids that are renowned for their high regenerative ability, growth plasticity and longevity. Because cnidarian genomes are conventional in terms of protein-coding genes, their remarkable features are likely a consequence of epigenetic regulation. To facilitate epigenetics research in cnidarians, we analysed the histone complement of the cnidarian model organism Hydractinia echinata using phylogenomics, proteomics, transcriptomics and mRNA in situ hybridisations.

Results

We find that the Hydractinia genome encodes 19 histones and analyse their spatial expression patterns, genomic loci and replication-dependency. Alongside core and other replication-independent histone variants, we find several histone replication-dependent variants, including a rare replication-dependent H3.3, a female germ cell-specific H2A.X and an unusual set of five H2B variants, four of which are male germ cell-specific. We further confirm the absence of protamines in Hydractinia.

Conclusions

Since no protamines are found in hydroids, we suggest that the novel H2B variants are pivotal for sperm DNA packaging in this class of Cnidaria. This study adds to the limited number of full histone gene complements available in animals and sets a comprehensive framework for future studies on the role of histones and their post-translational modifications in cnidarian epigenetics. Finally, it provides insight into the evolution of spermatogenesis.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-016-0085-1) contains supplementary material, which is available to authorized users.

TMEM231 Gene Conversion Associated with Joubert and Meckel-Gruber Syndromes in the Same Family

Joubert and Meckel-Gruber syndromes (JS and MGS) are ciliopathies with overlapping features. JS patients manifest the "molar tooth sign" on brain imaging and variable eye, kidney, and liver disease. MGS presents with polycystic kidneys, occipital encephalocele, and polydactyly; it is typically perinatally fatal. Both syndromes are genetically heterogeneous; some genes cause either syndrome. Here, we report two brothers married to unrelated women. The first brother had three daughters with JS and a son with polycystic kidneys who died at birth. The second brother's wife had a fetal demise due to MGS. Whole exome sequencing identified TMEM231 NM_001077416.2: c.784G>A; p.(Asp262Asn) in all children and the wife of the first brother; the second brother's wife had a c.406T>G;p.(Trp136Gly) change. In-depth analysis uncovered a rare gene conversion event in TMEM231, leading to loss of exon 4, in all the affected children of first brother. We believe that the combination of this gene conversion with different missense mutations led to a spectrum of phenotypes that span JS and MGS.

Assessing the reproducibility of exome copy number variations predictions

Background

Reproducibility is receiving increased attention across many domains of science and genomics is no exception. Efforts to identify copy number variations (CNVs) from exome sequence (ES) data have been increasing. Many algorithms have been published to discover CNVs from exomes and a major challenge is the reproducibility in other datasets. Here we test exome CNV calling reproducibility under three conditions: data generated by different sequencing centers; varying sample sizes; and varying capture methodology.

Methods

Four CNV tools were tested: eXome Hidden Markov Model (XHMM), Copy Number Inference From Exome Reads (CoNIFER), EXCAVATOR, and Copy Number Analysis for Targeted Resequencing (CONTRA). To examine the reproducibility, we ran the callers on four datasets, varying sample sizes of N = 10, 30, 75, 100, 300, and data with different capture methodology. We examined the false negative (FN) calls and false positive (FP) calls for potential limitations of the CNV callers. The positive predictive value (PPV) was measured by checking the CNV call concordance against single nucleotide polymorphism array.

Results

Using independently generated datasets, we examined the PPV for each dataset and observed wide range of PPVs. The PPV values were highly data dependent (p <0.001). For the sample sizes and capture method analyses, we tested the callers in triplicates. Both analyses resulted in wide ranges of PPVs, even for the same test. Interestingly, negative correlations between the PPV and the sample sizes were observed for CoNIFER (ρ = –0.80). Further examination of FN calls showed that 44 % of these were missed by all callers and were attributed to the CNV size (46 % spanned ≤3 exons). Overlap of the FP calls showed that FPs were unique to each caller, indicative of algorithm dependency.

Conclusions

Our results demonstrate that further improvements in CNV callers are necessary to improve reproducibility and to include wider spectrum of CNVs (including the small CNVs). These CNV callers should be evaluated on multiple independent, heterogeneously generated datasets of varying size to increase robustness and utility. These approaches to the evaluation of exome CNV are essential to support wide utility and applicability of CNV discovery in exome studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-016-0336-6) contains supplementary material, which is available to authorized users.

Vaccine-Induced Antibodies that Neutralize Group 1 and Group 2 Influenza A Viruses

Antibodies capable of neutralizing divergent influenza A viruses could form the basis of a universal vaccine. Here, from subjects enrolled in an H5N1 DNA/MIV-prime-boost influenza vaccine trial, we sorted hemagglutinin cross-reactive memory B cells and identified three antibody classes, each capable of neutralizing diverse subtypes of group 1 and group 2 influenza A viruses. Co-crystal structures with hemagglutinin revealed that each class utilized characteristic germline genes and convergent sequence motifs to recognize overlapping epitopes in the hemagglutinin stem. All six analyzed subjects had sequences from at least one multidonor class, and-in half the subjects-multidonor-class sequences were recovered from >40% of cross-reactive B cells. By contrast, these multidonor-class sequences were rare in published antibody datasets. Vaccination with a divergent hemagglutinin can thus increase the frequency of B cells encoding broad influenza A-neutralizing antibodies. We propose the sequence signature-quantified prevalence of these B cells as a metric to guide universal influenza A immunization strategies.

Developmental Pathway of the MPER-Directed HIV-1-Neutralizing Antibody 10E8

Antibody 10E8 targets the membrane-proximal external region (MPER) of HIV-1 gp41, neutralizes >97% of HIV-1 isolates, and lacks the auto-reactivity often associated with MPER-directed antibodies. The developmental pathway of 10E8 might therefore serve as a promising template for vaccine design, but samples from time-of-infection—often used to infer the B cell record—are unavailable. In this study, we used crystallography, next-generation sequencing (NGS), and functional assessments to infer the 10E8 developmental pathway from a single time point. Mutational analysis indicated somatic hypermutation of the 2^nd-heavy chain-complementarity determining region (CDR H2) to be critical for neutralization, and structures of 10E8 variants with V-gene regions reverted to genomic origin for heavy-and-light chains or heavy chain-only showed structural differences >2 Å relative to mature 10E8 in the CDR H2 and H3. To understand these developmental changes, we used bioinformatic sieving, maximum likelihood, and parsimony analyses of immunoglobulin transcripts to identify 10E8-lineage members, to infer the 10E8-unmutated common ancestor (UCA), and to calculate 10E8-developmental intermediates. We were assisted in this analysis by the preservation of a critical D-gene segment, which was unmutated in most 10E8-lineage sequences. UCA and early intermediates weakly bound a 26-residue-MPER peptide, whereas HIV-1 neutralization and epitope recognition in liposomes were only observed with late intermediates. Antibody 10E8 thus develops from a UCA with weak MPER affinity and substantial differences in CDR H2 and H3 from the mature 10E8; only after extensive somatic hypermutation do 10E8-lineage members gain recognition in the context of membrane and HIV-1 neutralization.

Detection and visualization of differential splicing in RNA-Seq data with JunctionSeq

Although RNA-Seq data provide unprecedented isoform-level expression information, detection of alternative isoform regulation (AIR) remains difficult, particularly when working with an incomplete transcript annotation. We introduce JunctionSeq, a new method that builds on the statistical techniques used by the well-established DEXSeq package to detect differential usage of both exonic regions and splice junctions. In particular, JunctionSeq is capable of detecting differential usage of novel splice junctions without the need for an additional isoform assembly step, greatly improving performance when the available transcript annotation is flawed or incomplete. JunctionSeq also provides a powerful and streamlined visualization toolset that allows bioinformaticians to quickly and intuitively interpret their results. We tested our method on publicly available data from several experiments performed on the rat pineal gland and Toxoplasma gondii, successfully detecting known and previously validated AIR genes in 19 out of 19 gene-level hypothesis tests. Due to its ability to query novel splice sites, JunctionSeq is still able to detect these differences even when all alternative isoforms for these genes were not included in the transcript annotation. JunctionSeq thus provides a powerful method for detecting alternative isoform regulation even with low-quality annotations. An implementation of JunctionSeq is available as an R/Bioconductor package.

Complex translocation disrupting TCF4 and altering TCF4 isoform expression segregates as mild autosomal dominant intellectual disability

BACKGROUND

Mutations of TCF4, which encodes a basic helix-loop-helix transcription factor, cause Pitt-Hopkins syndrome (PTHS) via multiple genetic mechanisms. TCF4 is a complex locus expressing multiple transcripts by alternative splicing and use of multiple promoters. To address the relationship between mutation of these transcripts and phenotype, we report a three-generation family segregating mild intellectual disability with a chromosomal translocation disrupting TCF4.

RESULTS

Using whole genome sequencing, we detected a complex unbalanced karyotype disrupting TCF4 (46,XY,del(14)(q23.3q23.3)del(18)(q21.2q21.2)del(18)(q21.2q21.2)inv(18)(q21.2q21.2)t(14;18)(q23.3;q21.2)(14pter®14q23.3::18q21.2®18q21.2::18q21.1®18qter;18pter®18q21.2::14q23.3®14qter). Subsequent transcriptome sequencing, qRT-PCR and nCounter analyses revealed that cultured skin fibroblasts and peripheral blood had normal expression of genes along chromosomes 14 or 18 and no marked changes in expression of genes other than TCF4. Affected individuals had 12-33 fold higher mRNA levels of TCF4 than did unaffected controls or individuals with PTHS. Although the derivative chromosome generated a PLEKHG3-TCF4 fusion transcript, the increased levels of TCF4 mRNA arose from transcript variants originating distal to the translocation breakpoint, not from the fusion transcript.

CONCLUSIONS

Although validation in additional patients is required, our findings suggest that the dysmorphic features and severe intellectual disability characteristic of PTHS are partially rescued by overexpression of those short TCF4 transcripts encoding a nuclear localization signal, a transcription activation domain, and the basic helix-loop-helix domain.

Maturation Pathway from Germline to Broad HIV-1 Neutralizer of a CD4-Mimic Antibody

Antibodies with ontogenies from VH1-2 or VH1-46-germline genes dominate the broadly neutralizing response against the CD4-binding site (CD4bs) on HIV-1. Here, we define with longitudinal sampling from time-of-infection the development of a VH1-46-derived antibody lineage that matured to neutralize 90% of HIV-1 isolates. Structures of lineage antibodies CH235 (week 41 from time-of-infection, 18% breadth), CH235.9 (week 152, 77%), and CH235.12 (week 323, 90%) demonstrated the maturing epitope to focus on the conformationally invariant portion of the CD4bs. Similarities between CH235 lineage and five unrelated CD4bs lineages in epitope focusing, length-of-time to develop breadth, and extraordinary level of somatic hypermutation suggested commonalities in maturation among all CD4bs antibodies. Fortunately, the required CH235-lineage hypermutation appeared substantially guided by the intrinsic mutability of the VH1-46 gene, which closely resembled VH1-2. We integrated our CH235-lineage findings with a second broadly neutralizing lineage and HIV-1 co-evolution to suggest a vaccination strategy for inducing both lineages.