Search

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.

Most colonial marine invertebrates are capable of allorecognition, the ability to distinguish between themselves and conspecifics. One long-standing question is whether invertebrate allorecognition genes are homologous to vertebrate histocompatibility genes. In the cnidarian Hydractinia symbiolongicarpus, allorecognition is controlled by at least two genes, Allorecognition 1 (Alr1) and Allorecognition 2 (Alr2), which encode highly polymorphic cell surface proteins that serve as markers of self. Here, we show that Alr1 and Alr2 are part of a family of 41 Alr genes, all of which reside a single genomic interval called the Allorecognition Complex (ARC). Using sensitive homology searches and highly accurate structural predictions, we demonstrate that the Alr proteins are members of the immunoglobulin superfamily (IgSF) with V-set and I-set Ig domains unlike any previously identified in animals. Specifically, their primary amino acid sequences lack many of the motifs considered diagnostic for V-set and I-set domains, yet they adopt secondary and tertiary structures nearly identical to canonical Ig domains. Thus, the V-set domain, which played a central role in the evolution of vertebrate adaptive immunity, was present in the last common ancestor of cnidarians and bilaterians. Unexpectedly, several Alr proteins also have immunoreceptor tyrosine-based activation motifs (ITAMs) and immunoreceptor tyrosine-based inhibitory motifs (ITIMs) in their cytoplasmic tails, suggesting they could participate in pathways homologous to those that regulate immunity in humans and flies. This work expands our definition of the IgSF with the addition of a family of unusual members, several of which play a role in invertebrate histocompatibility.Significance StatementThe immunoglobulin superfamily (IgSF) is one of the largest and most functionally versatile domain families in animal genomes. Although their amino acid sequences can vary considerably, IgSF domains have been traditionally defined by conserved residues at several key positions in their fold. Here, we sequenced an invertebrate histocompatibility complex and discovered a family of IgSF genes with amino acid sequences that lack most of these residues yet are predicted to adopt folds virtually identical to canonical V-set and I-set IgSF domains. This work broadens the definition of the IgSF and shows that the V-set domain was present earlier in animal evolution than previously appreciated.

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.

Bladder exstrophy (BE) is a rare, lower ventral midline defect with the bladder and part of the urethra exposed. The etiology of BE is unknown but thought to be influenced by genetic variation with more recent studies suggesting a role for rare variants. As such, we conducted paired-end exome sequencing in 26 child/mother/father trios. Three children had rare (allele frequency ≤ 0.0001 in several public databases) inherited variants in TSPAN4, one with a loss-of-function variant and two with missense variants. Two children had loss-of-function variants in TUBE1. Four children had rare missense or nonsense variants (one per child) in WNT3, CRKL, MYH9, or LZTR1, genes previously associated with BE. We detected 17 de novo missense variants in 13 children and three de novo loss-of-function variants (AKR1C2, PRRX1, PPM1D) in three children (one per child). We also detected rare compound heterozygous loss-of-function variants in PLCH2 and CLEC4M and rare inherited missense or loss-of-function variants in additional genes applying autosomal recessive (three genes) and X-linked recessive inheritance models (13 genes). Variants in two genes identified may implicate disruption in cell migration (TUBE1) and adhesion (TSPAN4) processes, mechanisms proposed for BE, and provide additional evidence for rare variants in the development of this defect.

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.

Background Adult-onset inflammatory syndromes often manifest with overlapping clinical features. Variants in ubiquitin-related genes, previously implicated in autoinflammatory disease, may...

Significance HLA alleles have previously been implicated with scleroderma risk, but, in this study, using a European American ancestral cohort and a newly recruited large cohort of African Americans, we comprehensively define the HLA alleles and amino acid residues associated with scleroderma. Scleroderma is characterized by mutually exclusive and specific autoantibodies. We demonstrated ancestry-predominant HLA alleles that were much more strongly associated with autoantibody subsets of scleroderma than with the overall risk of SSc. We bioinformatically predicted immunodominant peptides of self-antigens and demonstrated homology of these peptides with viral protein sequences from Mimiviridae and Phycodnaviridae families. Our findings suggest the hypothesis that scleroderma-specific autoantibodies may arise through molecular mimicry, driven by the interaction of specific viral antigens with corresponding HLA α/β heterodimers., Systemic sclerosis (SSc) is a clinically heterogeneous autoimmune disease characterized by mutually exclusive autoantibodies directed against distinct nuclear antigens. We examined HLA associations in SSc and its autoantibody subsets in a large, newly recruited African American (AA) cohort and among European Americans (EA). In the AA population, the African ancestry-predominant HLA-DRB1*08:04 and HLA-DRB1*11:02 alleles were associated with overall SSc risk, and the HLA-DRB1*08:04 allele was strongly associated with the severe antifibrillarin (AFA) antibody subset of SSc (odds ratio = 7.4). These African ancestry-predominant alleles may help explain the increased frequency and severity of SSc among the AA population. In the EA population, the HLA-DPB1*13:01 and HLA-DRB1*07:01 alleles were more strongly associated with antitopoisomerase (ATA) and anticentromere antibody-positive subsets of SSc, respectively, than with overall SSc risk, emphasizing the importance of HLA in defining autoantibody subtypes. The association of the HLA-DPB1*13:01 allele with the ATA+ subset of SSc in both AA and EA patients demonstrated a transancestry effect. A direct correlation between SSc prevalence and HLA-DPB1*13:01 allele frequency in multiple populations was observed (r = 0.98, P = 3 × 10−6). Conditional analysis in the autoantibody subsets of SSc revealed several associated amino acid residues, mostly in the peptide-binding groove of the class II HLA molecules. Using HLA α/β allelic heterodimers, we bioinformatically predicted immunodominant peptides of topoisomerase 1, fibrillarin, and centromere protein A and discovered that they are homologous to viral protein sequences from the Mimiviridae and Phycodnaviridae families. Taken together, these data suggest a possible link between HLA alleles, autoantibodies, and environmental triggers in the pathogenesis of SSc.

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 cleavage1–7. The physiological relevance of this cleavage remains unclear, though it is believed to inhibit activation of RIPK3 and necroptosis8. 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 Ripk1D325A mutant mouse strain. While Ripk1-/- mice die postnatally from systemic inflammation, Ripk1D325A/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 Ripk1D325A/D325A embryonic lethality, however the mice died before weaning from multi organ inflammation in a RIPK3 dependent manner. Consistently, Ripk1D325A/D325A and Ripk1D325A/+ cells were hypersensitive to RIPK3 dependent TNF-induced apoptosis and necroptosis. Heterozygous Ripk1D325A/+ 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.

In 2001, Celera Genomics and the International Human Genome Sequencing Consortium published their initial drafts of the human genome, which revolutionized the field of genomics. While these drafts and the updates that followed effectively covered the euchromatic fraction of the genome, the heterochromatin and many other complex regions were left unfinished or erroneous. Addressing this remaining 8% of the genome, the Telomere-to-Telomere (T2T) Consortium has finished the first truly complete 3.055 billion base pair (bp) sequence of a human genome, representing the largest improvement to the human reference genome since its initial release. The new T2T-CHM13 reference includes gapless assemblies for all 22 autosomes plus Chromosome X, corrects numerous errors, and introduces nearly 200 million bp of novel sequence containing 2,226 paralogous gene copies, 115 of which are predicted to be protein coding. The newly completed regions include all centromeric satellite arrays and the short arms of all five acrocentric chromosomes, unlocking these complex regions of the genome to variational and functional studies for the first time.

Identifying the molecular mechanisms by which genome-wide association study (GWAS) loci influence traits remains challenging. Chromatin accessibility quantitative trait loci (caQTLs) help identify GWAS loci that may alter GWAS traits by modulating chromatin structure, but caQTLs have been identified in a limited set of human tissues. Here we mapped caQTLs in human liver tissue in 20 liver samples and identified 3,123 caQTLs. The caQTL variants are enriched in liver tissue promoter and enhancer states and frequently disrupt binding motifs of transcription factors expressed in liver. We predicted target genes for 861 caQTL peaks using proximity, chromatin interactions, correlation with promoter accessibility or gene expression, and colocalization with expression QTLs. Using GWAS signals for 19 liver function and/or cardiometabolic traits, we identified 110 colocalized caQTLs and GWAS signals, 56 of which contained a predicted caPeak target gene. At the LITAF LDL-cholesterol GWAS locus, we validated that a caQTL variant showed allelic differences in protein binding and transcriptional activity. These caQTLs contribute to the epigenomic characterization of human liver and help identify molecular mechanisms and genes at GWAS loci.

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., High-coverage, ultra-long-read nanopore sequencing is used to create a new human genome assembly that improves on the coverage and accuracy of the current reference (GRCh38) and includes the gap-free, telomere-to-telomere sequence of the X chromosome.

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.

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.

New technologies and analysis methods are enabling genomic structural variants (SVs) to be detected with ever-increasing accuracy, resolution and comprehensiveness. To help translate these methods to routine research and clinical practice, we developed a sequence-resolved benchmark set for identification of both false-negative and false-positive germline large insertions and deletions. To create this benchmark for a broadly consented son in a Personal Genome Project trio with broadly available cells and DNA, the Genome in a Bottle Consortium integrated 19 sequence-resolved variant calling methods from diverse technologies. The final benchmark set contains 12,745 isolated, sequence-resolved insertion (7,281) and deletion (5,464) calls ≥50 base pairs (bp). The Tier 1 benchmark regions, for which any extra calls are putative false positives, cover 2.51 Gbp and 5,262 insertions and 4,095 deletions supported by ≥1 diploid assembly. We demonstrate that the benchmark set reliably identifies false negatives and false positives in high-quality SV callsets from short-, linked- and long-read sequencing and optical mapping.

Background Dermatophagoides pteronyssinus (DP) and Dermatophagoides farinae (DF) are highly similar disease-associated mites with frequently overlapping geographic distributions. A draft genome of DP was assembled to identify candidate allergens in DP homologous to those in DF, investigate allergen isoforms, and facilitate comparisons with related Acari. Methods PacBio and Illumina whole genome sequencing was performed on DP. Assembly and reconstruction of the genomes were optimized for isoform identification in a heterogeneous population. Bioinformatic analyses of Acari genomes were peformed. Results The predicted size of the DP nuclear genome is 52.5 Mb. A predicted protein set of 19,368 proteins was identified, including all 19 currently recognized allergens from this species. Orthologs for 12 allergens established for DF were found. The population of DP mites showed a high level of heterozygosity that allowed the identification of 43 new isoforms for both established and candidate allergens in DP, including a new isoform for the major allergen Der p 23. Reanalyzing the previous DF data assuming heterozygosity, 14 new allergen isoforms could be indentified. Some new isoforms were observed in both species suggesting that these isoforms pre-dated speciation. The high quality of both genomes allowed an examination of synteny which showed many allergen orthologs are physically clustered but with species specific exon/intron structures. Comparative genomic analyses with other Acariformes mites showed that most of the allergen homologs are widely conserved within this Superorder. Conclusions Candidate allergens in DP were identified to facilitate future serological studies. While DP and DF are highly similar genetically, species-specific allergen isoforms exist to facilitate molecular differentiation.

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.

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 2017. © 2017 American Cancer Society.

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

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. © 2017 American Cancer Society.

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.

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.

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.

New technologies and analysis methods are enabling genomic structural variants (SVs) to be detected with ever-increasing accuracy, resolution, and comprehensiveness. To help translate these methods to routine research and clinical practice, we developed the first sequence-resolved benchmark set for identification of both false negative and false positive germline large insertions and deletions. To create this benchmark for a broadly consented son in a Personal Genome Project trio with broadly available cells and DNA, the Genome in a Bottle (GIAB) Consortium integrated 19 sequence-resolved variant calling methods from diverse technologies. The final benchmark set contains 12745 isolated, sequence-resolved insertion (7281) and deletion (5464) calls ≥50 base pairs (bp). The Tier 1 benchmark regions, for which any extra calls are putative false positives, cover 2.51 Gbp and 5262 insertions and 4095 deletions supported by ≥1 diploid assembly. We demonstrate the benchmark set reliably identifies false negatives and false positives in high-quality SV callsets from short-, linked-, and long-read sequencing and optical mapping.

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.

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. This article is protected by copyright. All rights reserved

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.

BACKGROUND Next-generation sequencing (NGS) data are used for both clinical care and clinical research. DNA sequence variants identified using NGS are often returned to patients/participants as part of clinical or research protocols. The current standard of care is to validate NGS variants using Sanger sequencing, which is costly and time-consuming. METHODS We performed a large-scale, systematic evaluation of Sanger-based validation of NGS variants using data from the ClinSeq® project. We first used NGS data from 19 genes in 5 participants, comparing them to high-throughput Sanger sequencing results on the same samples, and found no discrepancies among 234 NGS variants. We then compared NGS variants in 5 genes from 684 participants against data from Sanger sequencing. RESULTS Of over 5800 NGS-derived variants, 19 were not validated by Sanger data. Using newly designed sequencing primers, Sanger sequencing confirmed 17 of the NGS variants, and the remaining 2 variants had low quality scores from exome sequencing. Overall, we measured a validation rate of 99.965% for NGS variants using Sanger sequencing, which was higher than many existing medical tests that do not necessitate orthogonal validation. CONCLUSIONS A single round of Sanger sequencing is more likely to incorrectly refute a true-positive variant from NGS than to correctly identify a false-positive variant from NGS. Validation of NGS-derived variants using Sanger sequencing has limited utility, and best practice standards should not include routine orthogonal Sanger validation of NGS variants.

Patients with autosomal dominant vibratory urticaria have localized hives and systemic manifestations in response to dermal vibration, with coincident degranulation of mast cells and increased histamine levels in serum. We identified a previously unknown missense substitution in ADGRE2 (also known as EMR2), which was predicted to result in the replacement of cysteine with tyrosine at amino acid position 492 (p.C492Y), as the only nonsynonymous variant cosegregating with vibratory urticaria in two large kindreds. The ADGRE2 receptor undergoes autocatalytic cleavage, producing an extracellular subunit that noncovalently binds a transmembrane subunit. We showed that the variant probably destabilizes an autoinhibitory subunit interaction, sensitizing mast cells to IgE-independent vibration-induced degranulation. (Funded by the National Institutes of Health.)

Broadly neutralizing antibodies (bNAbs) against HIV-1 Env V1V2 arise in multiple donors. However, atomic-level interactions had previously been determined only with antibodies from a single donor, thus making commonalities in recognition uncertain. Here we report the cocrystal structure of V1V2 with antibody CH03 from a second donor and model Env interactions of antibody CAP256-VRC26 from a third donor. These V1V2-directed bNAbs used strand-strand interactions between a protruding antibody loop and a V1V2 strand but differed in their N-glycan recognition. Ontogeny analysis indicated that protruding loops develop early, and glycan interactions mature over time. Altogether, the multidonor information suggested that V1V2-directed bNAbs form an 'extended class', for which we engineered ontogeny-specific antigens: Env trimers with chimeric V1V2s that interacted with inferred ancestor and intermediate antibodies. The ontogeny-based design of vaccine antigens described here may provide a general means for eliciting antibodies of a desired class.

© 2018 While consensus regarding the return of secondary genomic findings in the clinical setting has been reached, debate about such findings in the research setting remains. We developed a hybrid, research-clinical translational genomics process for research exome data coupled with a CLIA-validated secondary findings analysis. Eleven intramural investigators from ten institutes at the National Institutes of Health piloted this process. Nearly 1,200 individuals were sequenced and 14 secondary findings were identified in 18 participants. Positive secondary findings were returned by a genetic counselor following a standardized protocol, including referrals for specialty follow-up care for the secondary finding local to the participants. Interviews were undertaken with 13 participants 4 months after receipt of a positive report. These participants reported minimal psychologic distress within a process to assimilate their results. Of the 13, 9 reported accessing the recommended health care services. A sample of 107 participants who received a negative findings report were surveyed 4 months after receiving it. They demonstrated good understanding of the negative secondary findings result and most expressed reassurance (64%) from that report. However, a notable minority (up to 17%) expressed confusion regarding the distinction of primary from secondary findings. This pilot shows it is feasible to couple CLIA-compliant secondary findings to research sequencing with minimal harms. Participants managed the surprise of a secondary finding with most following recommended follow up, yet some with negative findings conflated secondary and primary findings. Additional work is needed to understand barriers to follow-up care and help participants distinguish secondary from primary findings.

A genomic assembly of the goldfish (Carrassius auratus) increases the utility of goldfish in research and agriculture., 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.

SummaryFor over a thousand years throughout Asia, the common goldfish (Carassius auratus) was raised for both food and as an ornamental pet. Selective breeding over more than 500 years has created a wide array of body and pigmentation variation particularly valued by ornamental fish enthusiasts. As a very close relative of the common carp (Cyprinus carpio), goldfish shares the recent genome duplication that occurred approximately 14-16 million years ago (mya) in their common ancestor. The combination of centuries of breeding and a wide array of interesting body morphologies is an exciting opportunity to link genotype to phenotype as well as understanding the dynamics of genome evolution and speciation. Here we generated a high-quality draft sequence of a “Wakin” goldfish using 71X PacBio long-reads. We identified 70,324 coding genes and more than 11,000 non-coding transcripts. We found that the two sub-genomes in goldfish retained extensive synteny and collinearity between goldfish and zebrafish. However, “ohnologous” genes were lost quickly after the carp whole-genome duplication, and the expression of 30% of the retained duplicated gene diverged significantly across seven tissues sampled. Loss of sequence identity and/or exons determined the divergence of the expression across all tissues, while loss of conserved, non-coding elements determined expression variance between different tissues. This draft assembly also provides an important resource for comparative genomics with the very commonly used zebrafish model (Danio rerio), and for understanding the underlying genetic causes of goldfish variants.

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.

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.

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. Electronic supplementary material The online version of this article (doi:10.1186/s13023-016-0439-6) contains supplementary material, which is available to authorized users.

Somatic mutational landscape of AML with inv(16) or t(8;21) identifies patterns of clonal evolution in relapse leukemia

// Bo Yu 1 , Valya R. Russanova 2 , Silvia Gravina 3 , Stephen Hartley 4 , James C. Mullikin 4, 5 , Alice Ignezweski 6 , James Graham 6 , James H. Segars 7 , Alan H. DeCherney 7 , Bruce H. Howard 2 1 Department of Obstetrics and Gynecology & Women's Health, Albert Einstein College of Medicine, Bronx, New York 10461, USA 2 Program in Genomics of Differentiation, Eunice Kennedy Shriver National Institute for Child Health and Human Development, National Institutes of Health, Bethesda, Maryland 20892, USA 3 Department of Genetics, Albert Einstein College of Medicine, Bronx, New York 10461, USA 4 Comparative Genomics Unit, Genome Technology Branch, National Human Genome Research Institute, National Institutes of Health, Bethesda, Maryland 20892, USA 5 NIH Intramural Sequencing Center, National Human Genome Research Institute, National Institutes of Health, Bethesda, Maryland 20892, USA 6 Shady Grove Fertility Reproductive Science Center, Rockville, Maryland 20850, USA 7 Program in Reproductive and Adult Endocrinology, Eunice Kennedy Shriver National Institute for Child Health and Human Development, National Institutes of Health, Bethesda, Maryland 20892, USA Correspondence to: Bo Yu, e-mail: boyumich2@gmail.com Keywords: DNA methylation, transcription end site, fertility, ovarian granulosa cell, transcriptome Received: November 14, 2014 Accepted: December 08, 2014 Published: February 17, 2015 ABSTRACT Diminished ovarian function occurs early and is a primary cause for age-related decline in female fertility; however, its underlying mechanism remains unclear. This study investigated the roles that genome and epigenome structure play in age-related changes in gene expression and ovarian function, using human ovarian granulosa cells as an experimental system. DNA methylomes were compared between two groups of women with distinct age-related differences in ovarian functions, using both Methylated DNA Capture followed by Next Generation Sequencing (MethylCap-seq) and Reduced Representation Bisulfite Sequencing (RRBS); their transcriptomes were investigated using mRNA-seq. Significant, non-random changes in transcriptome and DNA methylome features are observed in human ovarian granulosa cells as women age and their ovarian functions deteriorate. The strongest correlations between methylation and the age-related changes in gene expression are not confined to the promoter region; rather, high densities of hypomethylated CpG-rich regions spanning the gene body are preferentially associated with gene down-regulation. This association is further enhanced where CpG regions are localized near the 3'-end of the gene. Such features characterize several genes crucial in age-related decline in ovarian function, most notably the AMH (Anti-Mullerian Hormone) gene. The genome-wide correlation between the density of hypomethylated intragenic and 3'-end regions and gene expression suggests previously unexplored mechanisms linking epigenome structure to age-related physiology and pathology.