The structure, function and evolution of a complete human chromosome 8.

The complete assembly of each human chromosome is essential for understanding human biology and evolution1,2. Here we use complementary long-read sequencing technologies to complete the linear assembly of human chromosome 8. Our assembly resolves the sequence of five previously long-standing gaps, including a 2.08-Mb centromeric α-satellite array, a 644-kb copy number polymorphism in the ß-defensin gene cluster that is important for disease risk, and an 863-kb variable number tandem repeat at chromosome 8q21.2 that can function as a neocentromere. We show that the centromeric α-satellite array is generally methylated except for a 73-kb hypomethylated region of diverse higher-order α-satellites enriched with CENP-A nucleosomes, consistent with the location of the kinetochore. In addition, we confirm the overall organization and methylation pattern of the centromere in a diploid human genome. Using a dual long-read sequencing approach, we complete high-quality draft assemblies of the orthologous centromere from chromosome 8 in chimpanzee, orangutan and macaque to reconstruct its evolutionary history. Comparative and phylogenetic analyses show that the higher-order α-satellite structure evolved in the great ape ancestor with a layered symmetry, in which more ancient higher-order repeats locate peripherally to monomeric α-satellites. We estimate that the mutation rate of centromeric satellite DNA is accelerated by more than 2.2-fold compared to the unique portions of the genome, and this acceleration extends into the flanking sequence.

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.

Reproductive outcomes in individuals with chromosomal reciprocal translocations.

PURPOSE: Patients with reciprocal balanced translocations (RBT) have a risk for recurrent pregnancy losses (RPL), affected child, and infertility. Currently, genetic counseling is based on karyotypes found among the products of conception (POC), although factors influencing the success of assisted reproductive technologies (ART) in RBT couples are not established. METHODS: Cytogenetic results from 261 POC and offspring of the parents (113 women and 90 men) with RBT were evaluated. Chromosome segregation modes and number of euploid embryos were assessed in couples undergoing in vitro fertilization. RESULTS: Patients with translocations involving an acrocentric chromosome have a higher risk of unbalanced gametes caused by a 3:1 segregation. Female RBT patients have a statistically higher risk of aneuploidy due to an interchromosomal effect. The rate of euploid embryos is low due to meiosis I malsegregation of RBT, meiosis II nondisjunction, additional whole chromosome or segmental aneusomies. RBT patients with RPL have a higher rate of miscarriage of euploid fetuses with RBT. CONCLUSION: Chromosome-specific factors, female gender, age, and history of RPL are the risk elements influencing pregnancy and in vitro fertilization success in RBT patients. Chromosomal microarray analysis of POC is necessary to provide an accurate and timely diagnosis for patients with adverse reproductive outcomes.

Resolving the complexity of the human genome using single-molecule sequencing.

The human genome is arguably the most complete mammalian reference assembly, yet more than 160 euchromatic gaps remain and aspects of its structural variation remain poorly understood ten years after its completion. To identify missing sequence and genetic variation, here we sequence and analyse a haploid human genome (CHM1) using single-molecule, real-time DNA sequencing. We close or extend 55% of the remaining interstitial gaps in the human GRCh37 reference genome--78% of which carried long runs of degenerate short tandem repeats, often several kilobases in length, embedded within (G+C)-rich genomic regions. We resolve the complete sequence of 26,079 euchromatic structural variants at the base-pair level, including inversions, complex insertions and long tracts of tandem repeats. Most have not been previously reported, with the greatest increases in sensitivity occurring for events less than 5 kilobases in size. Compared to the human reference, we find a significant insertional bias (3:1) in regions corresponding to complex insertions and long short tandem repeats. Our results suggest a greater complexity of the human genome in the form of variation of longer and more complex repetitive DNA that can now be largely resolved with the application of this longer-read sequencing technology.

Improved assembly and variant detection of a haploid human genome using single-molecule, high-fidelity long reads.

The sequence and assembly of human genomes using long-read sequencing technologies has revolutionized our understanding of structural variation and genome organization. We compared the accuracy, continuity, and gene annotation of genome assemblies generated from either high-fidelity (HiFi) or continuous long-read (CLR) datasets from the same complete hydatidiform mole human genome. We find that the HiFi sequence data assemble an additional 10% of duplicated regions and more accurately represent the structure of tandem repeats, as validated with orthogonal analyses. As a result, an additional 5 Mbp of pericentromeric sequences are recovered in the HiFi assembly, resulting in a 2.5-fold increase in the NG50 within 1 Mbp of the centromere (HiFi 480.6 kbp, CLR 191.5 kbp). Additionally, the HiFi genome assembly was generated in significantly less time with fewer computational resources than the CLR assembly. Although the HiFi assembly has significantly improved continuity and accuracy in many complex regions of the genome, it still falls short of the assembly of centromeric DNA and the largest regions of segmental duplication using existing assemblers. Despite these shortcomings, our results suggest that HiFi may be the most effective standalone technology for de novo assembly of human genomes.

Correction: Comprehensive analysis of 204 sporadic hydatidiform moles: revisiting risk factors and their correlations with the molar genotypes.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Comprehensive analysis of 204 sporadic hydatidiform moles: revisiting risk factors and their correlations with the molar genotypes.

Hydatidiform mole (HM) is an aberrant human pregnancy characterized by excessive trophoblastic proliferation and abnormal embryonic development. HM has two morphological types, complete (CHM) and partial (PHM), and non-recurrent ones have three genotypic types, androgenetic monospermic, androgenetic dispermic, and triploid dispermic. Most available studies on risk factors predisposing to different types of HM and their malignant transformation mainly suffer from the lack of comprehensive genotypic analysis of large cohorts of molar tissues combined with accurate postmolar hCG follow-up. Moreover, 10-20% of patients with one HM have at least one non-molar miscarriage, which is higher than the frequency of two pregnancy losses in the general population (2-5%), suggesting a common genetic susceptibility to HM and miscarriages. However, the underlying causes of the miscarriages in these patients are unknown. Here, we comprehensively analyzed 204 HM, mostly from patients referred to the Quebec Registry of Trophoblastic Diseases and for which postmolar hCG monitoring is available, and 30 of their non-molar miscarriages. We revisited the risk of maternal age and neoplastic transformation across the different HM genotypic categories and investigated the presence of chromosomal abnormalities in their non-molar miscarriages. We confirm that androgenetic CHM is more prone to gestational trophoblastic neoplasia (GTN) than triploid dispermic PHM, and androgenetic dispermic CHM is more prone to high-risk GTN and choriocarcinoma (CC) than androgenetic monospermic CHM. We also confirm the association between increased maternal age and androgenetic CHM and their malignancies. Most importantly, we demonstrate for the first time that patients with an HM and miscarriages are at higher risk for aneuploid miscarriages [83.3%, 95% confidence interval (CI): 0.653-0.944] than women with sporadic (51.5%, 95% CI: 50.3-52.7%, p value = 0.0003828) or recurrent miscarriages (43.8%, 95% CI: 40.7-47.0%, p value = 0.00002). Our data suggest common genetic female germline defects predisposing to HM and aneuploid non-molar miscarriages in some patients.

Copy number alterations involving 59 ACMG-recommended secondary findings genes.

In clinical exome/genome sequencing, the American College of Medical Genetics and Genomics (ACMG) recommends reporting of secondary findings unrelated to a patient's phenotype when pathogenic single-nucleotide variants (SNVs) are observed in one of 59 genes associated with a life-threatening, medically actionable condition. Little is known about the incidence and sensitivity of chromosomal microarray analysis (CMA) for detection of pathogenic copy number variants (CNVs) comprising medically-actionable genes. Clinical CMA has been performed on 8865 individuals referred for molecular cytogenetic testing. We retrospectively reviewed the CMA results to identify patients with CNVs comprising genes included in the 59-ACMG list of secondary findings. We evaluated the clinical significance of these CNVs in respect to pathogenicity, phenotypic manifestations, and heritability. We identified 23 patients (0.26%) with relevant CNV either deletions comprising the entire gene or intragenic alterations involving one or more secondary findings genes. A number of patients and/or their family members with pathogenic CNVs manifest or expected to develop an anticipated clinical phenotype and would benefit from preventive management similar to the patients with pathogenic SNVs. To improve patients' care standardization should apply to reporting of both sequencing and CNVs obtained via clinical genome-wide analysis, including chromosomal microarray and exome/genome sequencing.

Four children with postnatally diagnosed mosaic trisomy 12: Clinical features, literature review, and current diagnostic capabilities of genetic testing.

Children or adults with mosaic trisomy 12 diagnosed postnatally are extremely rare. Only a small number of patients with this mosaicism have been reported in the literature. The clinical manifestation of mosaic trisomy 12 is variable, ranging from mild developmental delay to severe congenital anomaly and neonatal death. The trisomy 12 cells are not usually able to be detected by phytohemagglutinin stimulated peripheral blood chromosome analysis. The variability of phenotypes and the limited number of patients with this anomaly pose a challenge to predict the clinical outcomes. In this study, we present the phenotypes and laboratory findings in four patients and review the 11 previously reported patients with mosaic trisomy 12 diagnosed postnatally, as well as 11 patients with mosaic trisomy 12 diagnosed prenatally. The findings of this study provide useful information for laboratory diagnosis and clinical management of these patients.

DNA Methylation Profiling of Uniparental Disomy Subjects Provides a Map of Parental Epigenetic Bias in the Human Genome.

Genomic imprinting is a mechanism in which gene expression varies depending on parental origin. Imprinting occurs through differential epigenetic marks on the two parental alleles, with most imprinted loci marked by the presence of differentially methylated regions (DMRs). To identify sites of parental epigenetic bias, here we have profiled DNA methylation patterns in a cohort of 57 individuals with uniparental disomy (UPD) for 19 different chromosomes, defining imprinted DMRs as sites where the maternal and paternal methylation levels diverge significantly from the biparental mean. Using this approach we identified 77 DMRs, including nearly all those described in previous studies, in addition to 34 DMRs not previously reported. These include a DMR at TUBGCP5 within the recurrent 15q11.2 microdeletion region, suggesting potential parent-of-origin effects associated with this genomic disorder. We also observed a modest parental bias in DNA methylation levels at every CpG analyzed across ∼1.9 Mb of the 15q11-q13 Prader-Willi/Angelman syndrome region, demonstrating that the influence of imprinting is not limited to individual regulatory elements such as CpG islands, but can extend across entire chromosomal domains. Using RNA-seq data, we detected signatures consistent with imprinted expression associated with nine novel DMRs. Finally, using a population sample of 4,004 blood methylomes, we define patterns of epigenetic variation at DMRs, identifying rare individuals with global gain or loss of methylation across multiple imprinted loci. Our data provide a detailed map of parental epigenetic bias in the human genome, providing insights into potential parent-of-origin effects.

Phenotypic association of 15q11.2 CNVs of the region of breakpoints 1-2 (BP1-BP2) in a large cohort of samples referred for genetic diagnosis.

In view of conflicting reports on the pathogenicity of 15q11.2 CNVs of the breakpoints 1-2 (BP1-BP2) region and lack of association with a specific phenotype, we collected phenotypic data on 51,462 patients referred for genetic testing at two centers (Magee-Womens Hospital of UPMC and Baylor Genetics Laboratories, Baylor College of Medicine). Using array CGH, 262 patients with deletions and 215 with duplications were identified and tested for their association with four phenotypes (developmental delay, dysmorphic features, autism group of disorders, and epilepsy/seizures). Only association of deletions with dysmorphic features was observed (P = 0.013) with low penetrance (3.8%). Our results, viewed in the context of other reports suggesting the lack of a clear phenotypic outcome, underscore the need for detailed phenotypic studies to better understand the pathogenicity of 15q11.2 (BP1-BP2) CNVs.

Partial 5p Deletion and Partial 5q Duplication in a Patient with Multiple Congenital Anomalies: A Two-Step Mechanism in Chromosomal Rearrangement Mediated by Non-Allelic Homologous Recombination.

We describe a 5-month-old female who presented with clinical features of 5p deletion syndrome, including high-pitched cry, microcephaly, micrognathia, bilateral preauricular tags, bifid uvula, abnormal palmar creases, bilateral hypoplastic nipples, feeding difficulties, and developmental delay. In addition, the patient also had a cardiac defect, proximal esophageal atresia, and distal tracheoesophageal fistula. aCGH of the patient revealed a 22.9-Mb deletion of chromosome 5p15.33p14.3 and an 8.28-Mb duplication of chromosome 5q12.1q13.2. Parental chromosome analysis indicated that these alterations are de novo. Chromosome and FISH analysis demonstrated that the 5q12.1q13.2 duplicated segment was attached to the 5p14.3 region with the band 5q12.1 more distal to the centromere than the band 5q13.2. Based on the bioinformatic analysis, we postulate a mechanism for the formation of this complex rearrangement of chromosome 5 by 2-step-wise events mediate by nonallelic homologous recombination between low copy repeats. To the best of our knowledge this rearrangement found in our patient has not been reported in the literature. This report demonstrates the value of chromosome analysis in conjunction with FISH and aCGH for identification of complex rearrangements which cannot be revealed by array analysis alone.

Expression and trafficking of placental microRNAs at the feto-maternal interface.

During pregnancy, placental trophoblasts at the feto-maternal interface produce a broad repertoire of microRNA (miRNA) species. These species include miRNA from the primate-specific chromosome 19 miRNA cluster (C19MC), which is expressed nearly exclusively in the placenta. Trafficking of these miRNAs among the maternal, placental, and fetal compartments is unknown. To determine miRNA expression and trafficking patterns during pregnancy, we sequenced miRNAs in triads of human placenta and of maternal and fetal blood and found large subject-to-subject variability, with C19MC exhibiting compartment-specific expression. We therefore created humanized mice that transgenically express the entire 160-kb human C19MC locus or lentivirally express C19MC miRNA members selectively in the placenta. C19MC transgenic mice expressed a low level of C19MC miRNAs in diverse organs. When pregnant, female C19MC mice exhibited a strikingly elevated (>40-fold) expression of C19MC miRNA in the placenta, compared with other organs, that resembled C19MC miRNAs patterns in humans. Our mouse models showed that placental miRNA traffic primarily to the maternal circulation and that maternal miRNA can traffic to the placenta and even into the fetal compartment. These findings define an extraordinary means of nonhormonal, miRNA-based communication between the placenta and feto-maternal compartments.-Chang, G., Mouillet, J.-F., Mishima, T., Chu, T., Sadovsky, E., Coyne, C. B., Parks, W. T., Surti, U., Sadovsky, Y. Expression and trafficking of placental microRNAs at the feto-maternal interface.

Diploid/triploid mixoploidy: A consequence of asymmetric zygotic segregation of parental genomes.

Triploidy is the presence of an extra haploid set of chromosomes and can exist in complete or mosaic form. The extra haploid set of chromosomes in triploid cells can be of maternal or paternal origin. Diploid/triploid mixoploidy is a unique form of triploid mosaicism that requires the aberrant segregation of entire parental genomes into distinct blastomere lineages (heterogoneic cell division) at the earliest zygotic divisions. Here we report on eight cases of diploid/triploid mixoploidy from our institution and conduct a comprehensive review of the literature. The parental origin of the extra set of chromosomes was determined in two cases; and, based on phenotypic evidence we propose the parental origin in the other cases. One case with complex mixoploidy appears to have a digynic origin in addition to the involvement of two different sperm. Of our eight cases, only one resulted in the birth of a live healthy child. The other pregnancies ended in miscarriage, elective termination of pregnancy, intrauterine fetal demise or neonatal death. A review of the literature and the results of our cases show that a preponderance of recognized cases of diploid/triploid mixoploidy has a digynic origin.

MCM9 mutations are associated with ovarian failure, short stature, and chromosomal instability.

Premature ovarian failure (POF) is genetically heterogeneous and manifests as hypergonadotropic hypogonadism either as part of a syndrome or in isolation. We studied two unrelated consanguineous families with daughters exhibiting primary amenorrhea, short stature, and a 46,XX karyotype. A combination of SNP arrays, comparative genomic hybridization arrays, and whole-exome sequencing analyses identified homozygous pathogenic variants in MCM9, a gene implicated in homologous recombination and repair of double-stranded DNA breaks. In one family, the MCM9 c.1732+2T>C variant alters a splice donor site, resulting in abnormal alternative splicing and truncated forms of MCM9 that are unable to be recruited to sites of DNA damage. In the second family, MCM9 c.394C>T (p.Arg132(∗)) results in a predicted loss of functional MCM9. Repair of chromosome breaks was impaired in lymphocytes from affected, but not unaffected, females in both families, consistent with MCM9 function in homologous recombination. Autosomal-recessive variants in MCM9 cause a genomic-instability syndrome associated with hypergonadotropic hypogonadism and short stature. Preferential sensitivity of germline meiosis to MCM9 functional deficiency and compromised DNA repair in the somatic component most likely account for the ovarian failure and short stature.

Single haplotype assembly of the human genome from a hydatidiform mole.

A complete reference assembly is essential for accurately interpreting individual genomes and associating variation with phenotypes. While the current human reference genome sequence is of very high quality, gaps and misassemblies remain due to biological and technical complexities. Large repetitive sequences and complex allelic diversity are the two main drivers of assembly error. Although increasing the length of sequence reads and library fragments can improve assembly, even the longest available reads do not resolve all regions. In order to overcome the issue of allelic diversity, we used genomic DNA from an essentially haploid hydatidiform mole, CHM1. We utilized several resources from this DNA including a set of end-sequenced and indexed BAC clones and 100× Illumina whole-genome shotgun (WGS) sequence coverage. We used the WGS sequence and the GRCh37 reference assembly to create an assembly of the CHM1 genome. We subsequently incorporated 382 finished BAC clone sequences to generate a draft assembly, CHM1_1.1 (NCBI AssemblyDB GCA_000306695.2). Analysis of gene, repetitive element, and segmental duplication content show this assembly to be of excellent quality and contiguity. However, comparison to assembly-independent resources, such as BAC clone end sequences and PacBio long reads, indicate misassembled regions. Most of these regions are enriched for structural variation and segmental duplication, and can be resolved in the future. This publicly available assembly will be integrated into the Genome Reference Consortium curation framework for further improvement, with the ultimate goal being a completely finished gap-free assembly.

Beyond Down syndrome phenotype: Paternally derived isodicentric chromosome 21 with partial monosomy 21q22.3.

Inverted isodicentric chromosome 21 is a rare form of chromosomal rearrangement that may result in trisomy 21; sometimes this rearrangement may also lead to segmental monosomy of the terminal long arm of chromosome 21. In this report, we describe the prenatal diagnosis and neonatal follow-up of a child with a paternally derived, de novo isodicentric chromosome 21 and a concurrent ∼1.2 Mb deletion of the 21q22.3 region [46,XX,idic(21)(q22.3)]. This child presented with unusual phenotype of Down syndrome and additional defects including esophageal atresia and tethered cord syndrome. The resulting phenotype in this infant might be a coalescence of the partial trisomy and monosomy 21, as well as homozygosity for idic (21). The utilization of chromosomal microarray in this case enabled accurate characterization of a rare chromosome abnormality, potentially contributes to future phenotype-genotype correlation and produced evidence for a molecular mechanism underlying this rearrangement.

Pathogenetics of alveolar capillary dysplasia with misalignment of pulmonary veins.

Alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV) is a lethal lung developmental disorder caused by heterozygous point mutations or genomic deletion copy-number variants (CNVs) of FOXF1 or its upstream enhancer involving fetal lung-expressed long noncoding RNA genes LINC01081 and LINC01082. Using custom-designed array comparative genomic hybridization, Sanger sequencing, whole exome sequencing (WES), and bioinformatic analyses, we studied 22 new unrelated families (20 postnatal and two prenatal) with clinically diagnosed ACDMPV. We describe novel deletion CNVs at the FOXF1 locus in 13 unrelated ACDMPV patients. Together with the previously reported cases, all 31 genomic deletions in 16q24.1, pathogenic for ACDMPV, for which parental origin was determined, arose de novo with 30 of them occurring on the maternally inherited chromosome 16, strongly implicating genomic imprinting of the FOXF1 locus in human lungs. Surprisingly, we have also identified four ACDMPV families with the pathogenic variants in the FOXF1 locus that arose on paternal chromosome 16. Interestingly, a combination of the severe cardiac defects, including hypoplastic left heart, and single umbilical artery were observed only in children with deletion CNVs involving FOXF1 and its upstream enhancer. Our data demonstrate that genomic imprinting at 16q24.1 plays an important role in variable ACDMPV manifestation likely through long-range regulation of FOXF1 expression, and may be also responsible for key phenotypic features of maternal uniparental disomy 16. Moreover, in one family, WES revealed a de novo missense variant in ESRP1, potentially implicating FGF signaling in the etiology of ACDMPV.

Genomic Characterization of a Metastatic Alveolar Rhabdomyosarcoma Case Using FISH Studies and CGH+SNP Microarray Revealing FOXO1-PAX7 Rearrangement with MYCN and MDM2 Amplification and RB1 Region Loss.

Rhabdomyosarcomas (RMS) are rare, heterogeneous, soft tissue sarcomas and a common type of childhood malignancy with a distinct histomorphology. At the molecular level, alveolar rhabdomyosarcoma (ARMS), a subtype of RMS, harbors a signature genetic makeup characterized by specific translocations. The type of translocation and associated genetic aberrations correlate with disease progression, hence we used multiple molecular modalities including high-resolution array comparative genomic hybridization to explore the oncogenic gene fusion and associated copy number variations in a case of metastatic ARMS. We describe a case where traditional cytogenetic and molecular methods yielded inconclusive results in detecting the FOXO1 gene rearrangement. However, microarray analysis identified the essential FOXO1-PAX7 aberration and additional submicroscopic genomic alterations, including amplification of MYCN and MDM2 and deletion of RB1.