Integration of cytogenetic landmarks into the draft sequence of the human genome

Complete human recombination maps

Human recombination maps are a valuable resource for association and linkage studies and crucial for many inferences of population history and natural selection. Existing maps1-5 are based solely on cross-over (CO) recombination, omitting non-cross-overs (NCOs)-the more common form of recombination6-owing to the difficulty in detecting them. Using whole-genome sequence data in families, we estimate the number of NCOs transmitted from parent to offspring and derive complete, sex-specific recombination maps including both NCOs and COs. Mothers have fewer but longer NCOs than fathers, and oocytes accumulate NCOs in a non-regulated fashion with maternal age. Recombination, primarily NCO, is responsible for 1.8% (95% confidence interval: 1.3-2.3) and 11.3% (95% confidence interval: 9.0-13.6) of paternal and maternal de novo mutations, respectively, and may drive the increase in de novo mutations with maternal age. NCOs are substantially more prominent than COs in centromeres, possibly to avoid large-scale genomic changes that may cause aneuploidy. Our results demonstrate that NCOs highlight to a much greater extent than COs the differences in the meiotic process between the sexes, in which maternal NCOs may reflect the safeguarding of oocytes from infancy until ovulation.

SKI complex loss renders 9p21.3-deleted or MSI-H cancers dependent on PELO

Cancer genome alterations often lead to vulnerabilities that can be used to selectively target cancer cells. Various inhibitors of such synthetic lethal targets have been approved by the FDA or are in clinical trials, highlighting the potential of this approach1-3. Here we analysed large-scale CRISPR knockout screening data from the Cancer Dependency Map and identified a new synthetic lethal target, PELO, for two independent molecular subtypes of cancer: biallelic deletion of chromosomal region 9p21.3 or microsatellite instability-high (MSI-H). In 9p21.3-deleted cancers, PELO dependency emerges from biallelic deletion of the 9p21.3 gene FOCAD, a stabilizer of the superkiller complex (SKIc). In MSI-H cancers, PELO is required owing to MSI-H-associated mutations in TTC37 (also known as SKIC3), a critical component of the SKIc. We show that both cancer subtypes converge to destabilize the SKIc, which extracts mRNA from stalled ribosomes. In SKIc-deficient cells, PELO depletion induces the unfolded protein response, a stress response to accumulation of misfolded or unfolded nascent polypeptides. Together, our findings indicate PELO as a promising therapeutic target for a large patient population with cancers characterized as MSI-H with deleterious TTC37 mutations or with biallelic 9p21.3 deletions involving FOCAD.

The Theranostic Genome

Theranostic drugs represent an emerging path to deliver on the promise of precision medicine. However, bottlenecks remain in characterizing theranostic targets, identifying theranostic lead compounds, and tailoring theranostic drugs. To overcome these bottlenecks, we present the Theranostic Genome, the part of the human genome whose expression can be utilized to combine therapeutic and diagnostic applications. Using a deep learning-based hybrid human-AI pipeline that cross-references PubMed, the Gene Expression Omnibus, DisGeNET, The Cancer Genome Atlas and the NIH Molecular Imaging and Contrast Agent Database, we bridge individual genes in human cancers with respective theranostic compounds. Cross-referencing the Theranostic Genome with RNAseq data from over 17'000 human tissues identifies theranostic targets and lead compounds for various human cancers, and allows tailoring targeted theranostics to relevant cancer subpopulations. We expect the Theranostic Genome to facilitate the development of new targeted theranostics to better diagnose, understand, treat, and monitor a variety of human cancers.

A Novel Approach to Quantify Heterogeneity of Intrahepatic Cholangiocarcinoma: The Hidden-Genome Classifier

PURPOSE

Intrahepatic cholangiocarcinoma (IHC) is a heterogeneous tumor. The hidden-genome classifier, a supervised machine learning-based algorithm, was used to quantify tumor heterogeneity and improve classification.

EXPERIMENTAL DESIGN

A retrospective review of 1,370 patients with IHC, extrahepatic cholangiocarcinoma (EHC), gallbladder cancer (GBC), hepatocellular carcinoma (HCC), or biphenotypic tumors was conducted. A hidden-genome model classified 527 IHC based on genetic similarity to EHC/GBC or HCC. Genetic, histologic, and clinical data were correlated.

RESULTS

In this study, 410 IHC (78%) had >50% genetic homology with EHC/GBC; 122 (23%) had >90% homology ("biliary class"), characterized by alterations of KRAS, SMAD4, and CDKN2A loss; 117 IHC (22%) had >50% genetic homology with HCC; and 30 (5.7%) had >90% homology ("HCC class"), characterized by TERT alterations. Patients with biliary- versus non-biliary-class IHC had median overall survival (OS) of 1 year (95% CI, 0.77, 1.5) versus 1.8 years (95% CI, 1.6, 2.0) for unresectable disease and 2.4 years (95% CI, 2.1, NR) versus 5.1 years (95% CI, 4.8, 6.9) for resectable disease. Large-duct IHC (n = 28) was more common in the biliary class (n = 27); the HCC class was composed mostly of small-duct IHC (64%, P = 0.02). The hidden genomic classifier predicted OS independent of FGFR2 and IDH1 alterations. By contrast, the histology subtype did not predict OS.

CONCLUSIONS

IHC genetics form a spectrum with worse OS for tumors genetically aligned with EHC/GBC. The classifier proved superior to histologic subtypes for predicting OS independent of FGFR2 and IDH1 alterations. These results may explain the differential treatment responses seen in IHC and may direct therapy by helping stratify patients in future clinical trials.

Cell-specific regulation of the circadian clock by BMAL1 in the paraventricular nucleus: Implications for regulation of systemic biological rhythms

Circadian rhythms are internal biological rhythms driving temporal tissue-specific, metabolic programs. Loss of the circadian transcription factor BMAL1 in the paraventricular nucleus (PVN) of the hypothalamus reveals its importance in metabolic rhythms, but its functions in individual PVN cells are poorly understood. Here, loss of BMAL1 in the PVN results in arrhythmicity of processes controlling energy balance and alters peripheral diurnal gene expression. BMAL1 chromatin immunoprecipitation sequencing (ChIP-seq) and single-nucleus RNA sequencing (snRNA-seq) reveal its temporal regulation of target genes, including oxytocin (OXT), and restoring circulating OXT peaks in BMAL1-PVN knockout (KO) mice rescues absent activity rhythms. While glutamatergic neurons undergo day/night changes in expression of genes involved in cell morphogenesis, astrocytes and oligodendrocytes show gene expression changes in cytoskeletal organization and oxidative phosphorylation. Collectively, our findings show diurnal gene regulation in neuronal and non-neuronal PVN cells and that BMAL1 contributes to diurnal OXT secretion, which is important for systemic diurnal rhythms.

Genome-wide identification analysis in wild-type Solanum pinnatisectum reveals some genes defending against Phytophthora infestans

Solanum pinnatisectum exhibits strong resistance to late blight caused by Phytophthora infestans but only an incomplete genome assembly based on short Illumina reads has been published. In this study, we generated the first chromosome-level draft genome for the wild-type potato species S. pinnatisectum in China using Oxford Nanopore technology sequencing and Hi-C technology. The high-quality assembled genome size is 664 Mb with a scaffold N50 value of 49.17 Mb, of which 65.87% was occupied by repetitive sequences, and predominant long terminal repeats (42.51% of the entire genome). The genome of S. pinnatisectum was predicted to contain 34,245 genes, of which 99.34% were functionally annotated. Moreover, 303 NBS-coding disease resistance (R) genes were predicted in the S. pinnatisectum genome to investigate the potential mechanisms of resistance to late blight disease. The high-quality chromosome-level reference genome of S. pinnatisectum is expected to provide potential valuable resources for intensively and effectively investigating molecular breeding and genetic research in the future.

Genomic origin, fragmentomics, and transcriptional properties of long cell-free DNA molecules in human plasma

Recent studies have revealed an unexplored population of long cell-free DNA (cfDNA) molecules in human plasma using long-read sequencing technologies. However, the biological properties of long cfDNA molecules (>500 bp) remain largely unknown. To this end, we have investigated the origins of long cfDNA molecules from different genomic elements. Analysis of plasma cfDNA using long-read sequencing reveals an uneven distribution of long molecules from across the genome. Long cfDNA molecules show overrepresentation in euchromatic regions of the genome, in sharp contrast to short DNA molecules. We observe a stronger relationship between the abundance of long molecules and mRNA gene expression levels, compared with short molecules (Pearson's r = 0.71 vs. -0.14). Moreover, long and short molecules show distinct fragmentation patterns surrounding CpG sites. Leveraging the cleavage preferences surrounding CpG sites, the combined cleavage ratios of long and short molecules can differentiate patients with hepatocellular carcinoma (HCC) from non-HCC subjects (AUC = 0.87). We also investigated knockout mice in which selected nuclease genes had been inactivated in comparison with wild-type mice. The proportion of long molecules originating from transcription start sites are lower in Dffb-deficient mice but higher in Dnase1l3-deficient mice compared with that of wild-type mice. This work thus provides new insights into the biological properties and potential clinical applications of long cfDNA molecules.

ReConPlot: an R package for the visualization and interpretation of genomic rearrangements

MOTIVATION

Whole-genome sequencing studies of human tumours have revealed that complex forms of structural variation, collectively known as complex genome rearrangements (CGRs), are pervasive across diverse cancer types. Detection, classification, and mechanistic interpretation of CGRs requires the visualization of complex patterns of somatic copy number aberrations (SCNAs) and structural variants (SVs). However, there is a lack of tools specifically designed to facilitate the visualization and study of CGRs.

RESULTS

We present ReConPlot (REarrangement and COpy Number PLOT), an R package that provides functionalities for the joint visualization of SCNAs and SVs across one or multiple chromosomes. ReConPlot is based on the popular ggplot2 package, thus allowing customization of plots and the generation of publication-quality figures with minimal effort. Overall, ReConPlot facilitates the exploration, interpretation, and reporting of CGR patterns.

AVAILABILITY AND IMPLEMENTATION

The R package ReConPlot is available at https://github.com/cortes-ciriano-lab/ReConPlot. Detailed documentation and a tutorial with examples are provided with the package.

A cryptic microdeletion del(12)(p11.21p11.23) within an unbalanced translocation t(7;12)(q21.13;q23.1) implicates new candidate loci for intellectual disability and Kallmann syndrome

In a patient diagnosed with both Kallmann syndrome (KS) and intellectual disability (ID), who carried an apparently balanced translocation t(7;12)(q22;q24)dn, array comparative genomic hybridization (aCGH) disclosed a cryptic heterozygous 4.7 Mb deletion del(12)(p11.21p11.23), unrelated to the translocation breakpoint. This novel discovery prompted us to consider the possibility that the combination of KS and neurological disorder in this patient could be attributed to gene(s) within this specific deletion at 12p11.21-12p11.23, rather than disrupted or dysregulated genes at the translocation breakpoints. To further support this hypothesis, we expanded our study by screening five candidate genes at both breakpoints of the chromosomal translocation in a cohort of 48 KS patients. However, no mutations were found, thus reinforcing our supposition. In order to delve deeper into the characterization of the 12p11.21-12p11.23 region, we enlisted six additional patients with small copy number variations (CNVs) and analyzed eight individuals carrying small CNVs in this region from the DECIPHER database. Our investigation utilized a combination of complementary approaches. Firstly, we conducted a comprehensive phenotypic-genotypic comparison of reported CNV cases. Additionally, we reviewed knockout animal models that exhibit phenotypic similarities to human conditions. Moreover, we analyzed reported variants in candidate genes and explored their association with corresponding phenotypes. Lastly, we examined the interacting genes associated with these phenotypes to gain further insights. As a result, we identified a dozen candidate genes: TSPAN11 as a potential KS candidate gene, TM7SF3, STK38L, ARNTL2, ERGIC2, TMTC1, DENND5B, and ETFBKMT as candidate genes for the neurodevelopmental disorder, and INTS13, REP15, PPFIBP1, and FAR2 as candidate genes for KS with ID. Notably, the high-level expression pattern of these genes in relevant human tissues further supported their candidacy. Based on our findings, we propose that dosage alterations of these candidate genes may contribute to sexual and/or cognitive impairments observed in patients with KS and/or ID. However, the confirmation of their causal roles necessitates further identification of point mutations in these candidate genes through next-generation sequencing.

A microdeletion del(12)(p11.21p11.23) with a cryptic unbalanced translocation t(7;12)(q21.13;q23.1) implicates new candidate loci for intellectual disability and Kallmann syndrome

In an apparently balanced translocation t(7;12)(q22;q24)dn exhibiting both Kallmann syndrome (KS) and intellectual disability (ID), we detected a cryptic heterozygous 4.7 Mb del(12)(p11.21p11.23) unrelated to the translocation breakpoint. This new finding raised the possibility that KS combined with neurological disorder in this patient could be caused by gene(s) within this deletion at 12p11.21-12p11.23 instead of disrupted or dysregulated genes at the genomic breakpoints. Screening of five candidate genes at both breakpoints in 48 KS patients we recruited found no mutation, corroborating our supposition. To substantiate this hypothesis further, we recruited six additional subjects with small CNVs and analyzed eight individuals carrying small CNVs in this region from DECIPHER to dissect 12p11.21-12p11.23. We used multiple complementary approaches including a phenotypic-genotypic comparison of reported cases, a review of knockout animal models recapitulating the human phenotypes, and analyses of reported variants in the interacting genes with corresponding phenotypes. The results identified one potential KS candidate gene ( TSPAN11 ), seven candidate genes for the neurodevelopmental disorder ( TM7SF3 , STK38L , ARNTL2 , ERGIC2 , TMTC1 , DENND5B , and ETFBKMT ), and four candidate genes for KS with ID ( INTS13 , REP15 , PPFIBP1 , and FAR2 ). The high-level expression pattern in the relevant human tissues further suggested the candidacy of these genes. We propose that the dosage alterations of the candidate genes may contribute to sexual and/or cognitive impairment in patients with KS and/or ID. Further identification of point mutations through next generation sequencing will be necessary to confirm their causal roles.

Combined Focused Next-Generation Sequencing Assays to Guide Precision Oncology in Solid Tumors: A Retrospective Analysis from an Institutional Molecular Tumor Board

BACKGROUND

Increasing knowledge of cancer biology and an expanding spectrum of molecularly targeted therapies provide the basis for precision oncology. Despite extensive gene diagnostics, previous reports indicate that less than 10% of patients benefit from this concept.

METHODS

We retrospectively analyzed all patients referred to our center's Molecular Tumor Board (MTB) from 2018 to 2021. Molecular testing by next-generation sequencing (NGS) included a 67-gene panel for the detection of short-sequence variants and copy-number alterations, a 53- or 137-gene fusion panel and an ultra-low-coverage whole-genome sequencing for the detection of additional copy-number alterations outside the panel's target regions. Immunohistochemistry for microsatellite instability and PD-L1 expression complemented NGS.

RESULTS

A total of 109 patients were referred to the MTB. In all, 78 patients received therapeutic proposals (70 based on NGS) and 33 were treated accordingly. Evaluable patients treated with MTB-recommended therapy (n = 30) had significantly longer progression-free survival than patients treated with other therapies (n = 17) (4.3 vs. 1.9 months, p = 0.0094). Seven patients treated with off-label regimens experienced major clinical benefits.

CONCLUSION

The combined focused sequencing assays detected targetable alterations in the majority of patients. Patient benefits appeared to lie in the same range as with large-scale sequencing approaches.

The 'Alu-ome' shapes the epigenetic environment of regulatory elements controlling cellular defense

Promoters and enhancers are sites of transcription initiation (TSSs) and carry specific histone modifications, including H3K4me1, H3K4me3, and H3K27ac. Yet, the principles governing the boundaries of such regulatory elements are still poorly characterized. Alu elements are good candidates for a boundary function, being highly abundant in gene-rich regions, while essentially excluded from regulatory elements. Here, we show that the interval ranging from TSS to first upstream Alu, accommodates all H3K4me3 and most H3K27ac marks, while excluding DNA methylation. Remarkably, the average length of these intervals greatly varies in-between tissues, being longer in stem- and shorter in immune-cells. The very shortest TSS-to-first-Alu intervals were observed at promoters active in T-cells, particularly at immune genes, where first-Alus were traversed by RNA polymerase II transcription, while accumulating H3K4me1 signal. Finally, DNA methylation at first-Alus was found to evolve with age, regressing from young to middle-aged, then recovering later in life. Thus, the first-Alus upstream of TSSs appear as dynamic boundaries marking the transition from DNA methylation to active histone modifications at regulatory elements, while also participating in the recording of immune gene transcriptional events by positioning H3K4me1-modified nucleosomes.

Genome-Derived Classification Signature for Ampullary Adenocarcinoma to Improve Clinical Cancer Care

PURPOSE

The clinical behavior of ampullary adenocarcinoma varies widely. Targeted tumor sequencing may better define biologically distinct subtypes to improve diagnosis and management.

EXPERIMENTAL DESIGN

The hidden-genome algorithm, a multilevel meta-feature regression model, was trained on a prospectively sequenced cohort of 3,411 patients (1,001 pancreatic adenocarcinoma, 165 distal bile-duct adenocarcinoma, 2,245 colorectal adenocarcinoma) and subsequently applied to targeted panel DNA-sequencing data from ampullary adenocarcinomas. Genomic classification (i.e., colorectal vs. pancreatic) was correlated with standard histologic classification [i.e., intestinal (INT) vs. pancreatobiliary (PB)] and clinical outcome.

RESULTS

Colorectal genomic subtype prediction was primarily influenced by mutations in APC and PIK3CA, tumor mutational burden, and DNA mismatch repair (MMR)-deficiency signature. Pancreatic genomic-subtype prediction was dictated by KRAS gene alterations, particularly KRAS G12D, KRAS G12R, and KRAS G12V. Distal bile-duct adenocarcinoma genomic subtype was most influenced by copy-number gains in the MDM2 gene. Despite high (73%) concordance between immunomorphologic subtype and genomic category, there was significant genomic heterogeneity within both histologic subtypes. Genomic scores with higher colorectal probability were associated with greater survival compared with those with a higher pancreatic probability.

CONCLUSIONS

The genomic classifier provides insight into the heterogeneity of ampullary adenocarcinoma and improves stratification, which is dictated by the proportion of colorectal and pancreatic genomic alterations. This approach is reproducible with available molecular testing and obviates subjective histologic interpretation.

Anchoring the CerEla1.0 Genome Assembly to Red Deer (Cervus elaphus) and Cattle (Bos taurus) Chromosomes and Specification of Evolutionary Chromosome Rearrangements in Cervidae

Simple Summary

The red deer (Cervus elaphus) de novo genome assembly (CerEla1.0) has provided a great resource for genetic studies in various deer species. In this study, we used gene order comparisons between C. elaphus CerEla1.0 and B. taurus ARS-UCD1.2 genome assemblies and fluorescence in situ hybridization (FISH) with bovine BAC probes to verify the red deer-bovine chromosome relationships and anchor the CerEla1.0 C-scaffolds to karyotypes of both species. We showed the homology between bovine and deer chromosomes and determined the centromere-telomere orientation of the CerEla1.0 C-scaffolds. Using a set of BAC probes, we were able to narrow the positions of evolutionary chromosome breakpoints defining the family Cervidae. In addition, we revealed several errors in the current CerEla1.0 genome assembly. Finally, we expanded our analysis to other Cervidae and confirmed the locations of the cervid evolutionary fissions and orientation of the fused chromosomes in eight cervid species. Our results can serve as a basis for necessary improvements of the red deer genome assembly and provide support to other genetic studies in Cervidae.

Abstract

The family Cervidae groups a range of species with an increasing economic significance. Their karyotypes share 35 evolutionary conserved chromosomal segments with cattle (Bos taurus). Recent publication of the annotated red deer (Cervus elaphus) whole genome assembly (CerEla1.0) has provided a basis for advanced genetic studies. In this study, we compared the red deer CerEla1.0 and bovine ARS-UCD1.2 genome assembly and used fluorescence in situ hybridization with bovine BAC probes to verify the homology between bovine and deer chromosomes, determined the centromere-telomere orientation of the CerEla1.0 C-scaffolds and specified positions of the cervid evolutionary chromosome breakpoints. In addition, we revealed several incongruences between the current deer and bovine genome assemblies that were shown to be caused by errors in the CerEla1.0 assembly. Finally, we verified the centromere-to-centromere orientation of evolutionarily fused chromosomes in seven additional deer species, giving a support to previous studies on their chromosome evolution.

Copy-number variation of the NPHP1 gene in patients with juvenile Nephronophthisis

Objective: Juvenile nephronophthisis (NPHP) is an autosomal recessive cystic disease of the kidney. It represents the most frequent genetic cause of chronic renal failure in children. Methods: we investigated clinical and molecular features in two children with Juvenile nephronophthisis using firstly Multiplex ligation-dependent probe amplification (MLPA) and secondly multiplex PCR. Results: we report a homozygous NPHP1 deletion in two children. Conclusion: NPHP1 deletion analysis using diagnostic methods (e.g. MLPA, Multiplex PCR) should always be considered in patients with nephronophthisis, especially from consanguineous families. Our results provide insights into genotype-phenotype correlations in juvenile nephronophthisis that can be utilized in genetic counseling.

Characterization of long G4-rich enhancer-associated genomic regions engaging in a novel loop:loop 'G4 Kissing' interaction

Mammalian antibody switch regions (∼1500 bp) are composed of a series of closely neighboring G4-capable sequences. Whereas numerous structural and genome-wide analyses of roles for minimal G4s in transcriptional regulation have been reported, Long G4-capable regions (LG4s)-like those at antibody switch regions-remain virtually unexplored. Using a novel computational approach we have identified 301 LG4s in the human genome and find LG4s prone to mutation and significantly associated with chromosomal rearrangements in malignancy. Strikingly, 217 LG4s overlap annotated enhancers, and we find the promoters regulated by these enhancers markedly enriched in G4-capable sequences suggesting G4s facilitate promoter-enhancer interactions. Finally, and much to our surprise, we also find single-stranded loops of minimal G4s within individual LG4 loci are frequently highly complementary to one another with 178 LG4 loci averaging >35 internal loop:loop complements of >8 bp. As such, we hypothesized (then experimentally confirmed) that G4 loops within individual LG4 loci directly basepair with one another (similar to characterized stem-loop kissing interactions) forming a hitherto undescribed, higher-order, G4-based secondary structure we term a 'G4 Kiss or G4K'. In conclusion, LG4s adopt novel, higher-order, composite G4 structures directly contributing to the inherent instability, regulatory capacity, and maintenance of these conspicuous genomic regions.

Genome-wide unique insertion sequences among five Brucella species and demonstration of differential identification of Brucella by multiplex PCR assay

Brucellosis is a neglected zoonotic disease caused by alpha proteobacterial genus Brucella comprising of facultative intracellular pathogenic species that can infect both animals and humans. In this study, we aimed to identify genome-wide unique insertion sequence (IS) elements among Brucella abortus, B. melitensis, B. ovis, B. suis and B. canis for use in species differentiation by conducting an intensive in silico-based comparative genomic analysis. As a result, 25, 27, 37, 86 and 3 unique ISs were identified respectively and they had a striking pattern of distribution among them. To explain, a particular IS would be present in four species with 100% identity whereas completely absent in the fifth species. However, flanking regions of that IS element would be highly identical and conserved in all five species. Species-specific primers designed on these flanking conserved regions resulted in two different amplicons grouping the species into two: one that possesses IS and the other that lacks it. Seeking for species-specific amplicon size for particular species was sufficient to identify it irrespective of biovar. A multiplex PCR developed using these primers resulted in successful differentiation of the five species irrespective of biovars with significant specificity and sensitivity when examined on clinical samples.

Cellular and genomic approaches for exploring structural chromosomal rearrangements

Human chromosomes are arranged in a linear and conserved sequence order that undergoes further spatial folding within the three-dimensional space of the nucleus. Although structural variations in this organization are an important source of natural genetic diversity, cytogenetic aberrations can also underlie a number of human diseases and disorders. Approaches for studying chromosome structure began half a century ago with karyotyping of Giemsa-banded chromosomes and has now evolved to encompass high-resolution fluorescence microscopy, reporter-based assays, and next-generation DNA sequencing technologies. Here, we provide a general overview of experimental methods at different resolution and sensitivity scales and discuss how they can be complemented to provide synergistic insight into the study of human chromosome structural rearrangements. These approaches range from kilobase-level resolution DNA fluorescence in situ hybridization (FISH)-based imaging approaches of individual cells to genome-wide sequencing strategies that can capture nucleotide-level information from diverse sample types. Technological advances coupled to the combinatorial use of multiple methods have resulted in the discovery of new rearrangement classes along with mechanistic insights into the processes that drive structural alterations in the human genome.

A Multireference-Based Whole Genome Assembly for the Obligate Ant-Following Antbird, Rhegmatorhina melanosticta (Thamnophilidae)

Current generation high-throughput sequencing technology has facilitated the generation of more genomic-scale data than ever before, thus greatly improving our understanding of avian biology across a range of disciplines. Recent developments in linked-read sequencing (Chromium 10×) and reference-based whole-genome assembly offer an exciting prospect of more accessible chromosome-level genome sequencing in the near future. We sequenced and assembled a genome of the Hairy-crested Antbird (Rhegmatorhina melanosticta), which represents the first publicly available genome for any antbird (Thamnophilidae). Our objectives were to (1) assemble scaffolds to chromosome level based on multiple reference genomes, and report on differences relative to other genomes, (2) assess genome completeness and compare content to other related genomes, and (3) assess the suitability of linked-read sequencing technology for future studies in comparative phylogenomics and population genomics studies. Our R. melanosticta assembly was both highly contiguous (de novo scaffold N50 = 3.3 Mb, reference based N50 = 53.3 Mb) and relatively complete (contained close to 90% of evolutionarily conserved single-copy avian genes and known tetrapod ultraconserved elements). The high contiguity and completeness of this assembly enabled the genome to be successfully mapped to the chromosome level, which uncovered a consistent structural difference between R. melanosticta and other avian genomes. Our results are consistent with the observation that avian genomes are structurally conserved. Additionally, our results demonstrate the utility of linked-read sequencing for non-model genomics. Finally, we demonstrate the value of our R. melanosticta genome for future researchers by mapping reduced representation sequencing data, and by accurately reconstructing the phylogenetic relationships among a sample of thamnophilid species.

Genome-wide mapping and profiling of γH2AX binding hotspots in response to different replication stress inducers

BACKGROUND

Replication stress (RS) gives rise to DNA damage that threatens genome stability. RS can originate from different sources that stall replication by diverse mechanisms. However, the mechanism underlying how different types of RS contribute to genome instability is unclear, in part due to the poor understanding of the distribution and characteristics of damage sites induced by different RS mechanisms.

RESULTS

We use ChIP-seq to map γH2AX binding sites genome-wide caused by aphidicolin (APH), hydroxyurea (HU), and methyl methanesulfonate (MMS) treatments in human lymphocyte cells. Mapping of γH2AX ChIP-seq reveals that APH, HU, and MMS treatments induce non-random γH2AX chromatin binding at discrete regions, suggesting that there are γH2AX binding hotspots in the genome. Characterization of the distribution and sequence/epigenetic features of γH2AX binding sites reveals that the three treatments induce γH2AX binding at largely non-overlapping regions, suggesting that RS may cause damage at specific genomic loci in a manner dependent on the fork stalling mechanism. Nonetheless, γH2AX binding sites induced by the three treatments share common features including compact chromatin, coinciding with larger-than-average genes, and depletion of CpG islands and transcription start sites. Moreover, we observe significant enrichment of SINEs in γH2AX sites in all treatments, indicating that SINEs may be a common barrier for replication polymerases.

CONCLUSIONS

Our results identify the location and common features of genome instability hotspots induced by different types of RS, and help in deciphering the mechanisms underlying RS-induced genetic diseases and carcinogenesis.

Next Generation Sequencing of Prenatal Structural Chromosomal Rearrangements Using Large-Insert Libraries

Precise tests for genomic structural variation (SV) are essential for accurate diagnosis of prenatal genome abnormalities. The two most ubiquitous traditional methods for prenatal SV assessment, karyotyping and chromosomal microarrays, do not provide sufficient resolution for some clinically actionable SVs. Standard whole-genome sequencing (WGS) overcomes shortcomings of traditional techniques by providing base-pair resolution of the entire accessible genome. However, while sequencing costs have continued to decline in recent years, conventional WGS costs remain high for most routine clinical applications. Here, we describe a specialized WGS technique using large inserts (liWGS; also known as "jumping libraries") to resolve large (>5000-10,000 nucleotides) SVs at kilobase-resolution in prenatal samples, and at a fraction of the cost of standard WGS. We explicate the protocols for generating liWGS libraries and supplement with an overview for processing and analyzing liWGS data.

Exploring chromatin hierarchical organization via Markov State Modelling

We propose a new computational method for exploring chromatin structural organization based on Markov State Modelling of Hi-C data represented as an interaction network between genomic loci. A Markov process describes the random walk of a traveling probe in the corresponding energy landscape, mimicking the motion of a biomolecule involved in chromatin function. By studying the metastability of the associated Markov State Model upon annealing, the hierarchical structure of individual chromosomes is observed, and corresponding set of structural partitions is identified at each level of hierarchy. Then, the notion of effective interaction between partitions is derived, delineating the overall topology and architecture of chromosomes. Mapping epigenetic data on the graphs of intra-chromosomal effective interactions helps in understanding how chromosome organization facilitates its function. A sketch of whole-genome interactions obtained from the analysis of 539 partitions from all 23 chromosomes, complemented by distributions of gene expression regulators and epigenetic factors, sheds light on the structure-function relationships in chromatin, delineating chromosomal territories, as well as structural partitions analogous to topologically associating domains and active / passive epigenomic compartments. In addition to the overall genome architecture shown by effective interactions, the affinity between partitions of different chromosomes was analyzed as an indicator of the degree of association between partitions in functionally relevant genomic interactions. The overall static picture of whole-genome interactions obtained with the method presented in this work provides a foundation for chromatin structural reconstruction, for the modelling of chromatin dynamics, and for exploring the regulation of genome function. The algorithms used in this study are implemented in a freely available Python package ChromaWalker (https://bitbucket.org/ZhenWahTan/chromawalker).

Patterns and mechanisms of structural variations in human cancer

Next-generation sequencing technology has enabled the comprehensive detection of genomic alterations in human somatic cells, including point mutations, chromosomal rearrangements, and structural variations (SVs). Using sophisticated bioinformatics algorithms, unbiased catalogs of SVs are emerging from thousands of human cancer genomes for the first time. Via careful examination of SV breakpoints at single-nucleotide resolution as well as local DNA copy number changes, diverse patterns of genomic rearrangements are being revealed. These "SV signatures" provide deep insight into the mutational processes that have shaped genome changes in human somatic cells. This review summarizes the characteristics of recently identified complex SVs, including chromothripsis, chromoplexy, microhomology-mediated breakage-induced replication (MMBIR), and others, to provide a holistic snapshot of the current knowledge on genomic rearrangements in somatic cells.

A new targeted capture method using bacterial artificial chromosome (BAC) libraries as baits for sequencing relatively large genes

To analyze a specific genome region using next-generation sequencing technologies, the enrichment of DNA libraries with targeted capture methods has been standardized. For enrichment of mitochondrial genome, a previous study developed an original targeted capture method that use baits constructed from long-range polymerase chain reaction (PCR) amplicons, common laboratory reagents, and equipment. In this study, a new targeted capture method is presented, that of bacterial artificial chromosome (BAC) double capture (BDC), modifying the previous method, but using BAC libraries as baits for sequencing a relatively large gene. We applied the BDC approach for the 214 kb autosomal region, ring finger protein 213, which is the susceptibility gene of moyamoya disease (MMD). To evaluate the reliability of BDC, cost and data quality were compared with those of a commercial kit. While the ratio of duplicate reads was higher, the cost was less than that of the commercial kit. The data quality was sufficiently the same as that of the kit. Thus, BDC can be an easy, low-cost, and useful method for analyzing individual genome regions with substantial length.

Evaluation of preimplantation genetic testing for chromosomal structural rearrangement by a commonly used next generation sequencing workflow

OBJECTIVES

To evaluate the applicability of a commonly used next generation sequencing workflow in detecting unbalanced meiotic segregation products for reciprocal translocation and inversion carriers.

STUDY DESIGN

All preimplantation genetic testing treatment cycles performed for reciprocal translocation or inversion carriers from 2012 to April 2017 were included. Three hundreds and forty-two archived whole genome amplified DNA, which had previously analyzed by array comparative genomic hybridization (aCGH), were retrospectively analyzed by next generation sequencing (NGS). Concordance on overall diagnosis and segmental aneuploidies related to the translocation/inversion breakpoints between aCGH and NGS were determined.

RESULTS

Retrospective analysis of 287 blastomere biopsies and 55 trophectoderm (TE) biopsies showed that the concordance rate on the overall diagnosis between aCGH and NGS on abnormal samples was 100% (266/266), irrespective to the type of biopsy. The concordance rates of normal biopsies were 98.4% (61/62) on blastomere and 78.6% (11/14) on TE biopsies. NGS detected a de novo segmental aneuploidy on one blastomere biopsy and three possible low level mosaic aneuploidies on 3 TE biopsies, which were previously concluded as euploid by aCGH. Using the karyotype of reciprocal translocation/inversion carriers, size of anticipated segmental aneuploidies could be calculated and be used to predict the applicability of NGS before proceeding to treatment.

CONCLUSION

This is the first report to evaluate the applicability of a commercial NGS-based workflow for preimplantation testing for reciprocal translocations/inversions. Our study demonstrated that NGS can diagnose unbalanced translocation/inversion products with the same efficiency as aCGH. The applicability of NGS, with respect to individual karyotype, can be predicted before proceeding to treatment.

Next generation sequencing: implications in personalized medicine and pharmacogenomics

A breakthrough in next generation sequencing (NGS) in the last decade provided an unprecedented opportunity to investigate genetic variations in humans and their roles in health and disease. NGS offers regional genomic sequencing such as whole exome sequencing of coding regions of all genes, as well as whole genome sequencing. RNA-seq offers sequencing of the entire transcriptome and ChIP-seq allows for sequencing the epigenetic architecture of the genome. Identifying genetic variations in individuals can be used to predict disease risk, with the potential to halt or retard disease progression. NGS can also be used to predict the response to or adverse effects of drugs or to calculate appropriate drug dosage. Such a personalized medicine also provides the possibility to treat diseases based on the genetic makeup of the patient. Here, we review the basics of NGS technologies and their application in human diseases to foster human healthcare and personalized medicine.

In vivo evasion of MxA by avian influenza viruses requires human signature in the viral nucleoprotein

Deeg et al. show a novel line of transgenic mice expressing restriction factor MxA exhibits robust resistance to influenza viruses of avian but not human origin. In vivo evasion of MxA is mediated by distinct amino acids in the nucleoprotein of human influenza viruses.

Zoonotic transmission of influenza A viruses can give rise to devastating pandemics, but currently it is impossible to predict the pandemic potential of circulating avian influenza viruses. Here, we describe a new mouse model suitable for such risk assessment, based on the observation that the innate restriction factor MxA represents an effective species barrier that must be overcome by zoonotic viruses. Our mouse lacks functional endogenous Mx genes but instead carries the human MX1 locus as a transgene. Such transgenic mice were largely resistant to highly pathogenic avian H5 and H7 influenza A viruses, but were almost as susceptible to infection with influenza viruses of human origin as nontransgenic littermates. Influenza A viruses that successfully established stable lineages in humans have acquired adaptive mutations which allow partial MxA escape. Accordingly, an engineered avian H7N7 influenza virus carrying a nucleoprotein with signature mutations typically found in human virus isolates was more virulent in transgenic mice than parental virus, demonstrating that a few amino acid changes in the viral target protein can mediate escape from MxA restriction in vivo. Similar mutations probably need to be acquired by emerging influenza A viruses before they can spread in the human population.

The UCSC Genome Browser database: 2017 update

Since its 2001 debut, the University of California, Santa Cruz (UCSC) Genome Browser (http://genome.ucsc.edu/) team has provided continuous support to the international genomics and biomedical communities through a web-based, open source platform designed for the fast, scalable display of sequence alignments and annotations landscaped against a vast collection of quality reference genome assemblies. The browser's publicly accessible databases are the backbone of a rich, integrated bioinformatics tool suite that includes a graphical interface for data queries and downloads, alignment programs, command-line utilities and more. This year's highlights include newly designed home and gateway pages; a new 'multi-region' track display configuration for exon-only, gene-only and custom regions visualization; new genome browsers for three species (brown kiwi, crab-eating macaque and Malayan flying lemur); eight updated genome assemblies; extended support for new data types such as CRAM, RNA-seq expression data and long-range chromatin interaction pairs; and the unveiling of a new supported mirror site in Japan.

Stacked thin layers of metaphase chromatin explain the geometry of chromosome rearrangements and banding

The three-dimensional organization of tightly condensed chromatin within metaphase chromosomes has been one of the most challenging problems in structural biology since the discovery of the nucleosome. This study shows that chromosome images obtained from typical banded karyotypes and from different multicolour cytogenetic analyses can be used to gain information about the internal structure of chromosomes. Chromatin bands and the connection surfaces in sister chromatid exchanges and in cancer translocations are planar and orthogonal to the chromosome axis. Chromosome stretching produces band splitting and even the thinnest bands are orthogonal and well defined, indicating that short stretches of DNA can occupy completely the chromosome cross-section. These observations impose strong physical constraints on models that attempt to explain chromatin folding in chromosomes. The thin-plate model, which consists of many stacked layers of planar chromatin perpendicular to the chromosome axis, is compatible with the observed orientation of bands, with the existence of thin bands, and with band splitting; it is also compatible with the orthogonal orientation and planar geometry of the connection surfaces in chromosome rearrangements. The results obtained provide a consistent interpretation of the chromosome structural properties that are used in clinical cytogenetics for the diagnosis of hereditary diseases and cancers.

Evolution of genetic techniques: past, present, and beyond

Genetics is the study of heredity, which means the study of genes and factors related to all aspects of genes. The scientific history of genetics began with the works of Gregor Mendel in the mid-19th century. Prior to Mendel, genetics was primarily theoretical whilst, after Mendel, the science of genetics was broadened to include experimental genetics. Developments in all fields of genetics and genetic technology in the first half of the 20th century provided a basis for the later developments. In the second half of the 20th century, the molecular background of genetics has become more understandable. Rapid technological advancements, followed by the completion of Human Genome Project, have contributed a great deal to the knowledge of genetic factors and their impact on human life and diseases. Currently, more than 1800 disease genes have been identified, more than 2000 genetic tests have become available, and in conjunction with this at least 350 biotechnology-based products have been released onto the market. Novel technologies, particularly next generation sequencing, have dramatically accelerated the pace of biological research, while at the same time increasing expectations. In this paper, a brief summary of genetic history with short explanations of most popular genetic techniques is given.

DNA methylation at hepatitis B viral integrants is associated with methylation at flanking human genomic sequences

Integration of DNA viruses into the human genome plays an important role in various types of tumors, including hepatitis B virus (HBV)–related hepatocellular carcinoma. However, the molecular details and clinical impact of HBV integration on either human or HBV epigenomes are unknown. Here, we show that methylation of the integrated HBV DNA is related to the methylation status of the flanking human genome. We developed a next-generation sequencing-based method for structural methylation analysis of integrated viral genomes (denoted G-NaVI). This method is a novel approach that enables enrichment of viral fragments for sequencing using unique baits based on the sequence of the HBV genome. We detected integrated HBV sequences in the genome of the PLC/PRF/5 cell line and found variable levels of methylation within the integrated HBV genomes. Allele-specific methylation analysis revealed that the HBV genome often became significantly methylated when integrated into highly methylated host sites. After integration into unmethylated human genome regions such as promoters, however, the HBV DNA remains unmethylated and may eventually play an important role in tumorigenesis. The observed dynamic changes in DNA methylation of the host and viral genomes may functionally affect the biological behavior of HBV. These findings may impact public health given that millions of people worldwide are carriers of HBV. We also believe our assay will be a powerful tool to increase our understanding of the various types of DNA virus-associated tumorigenesis.

A recurrent deletion syndrome at chromosome bands 2p11.2-2p12 flanked by segmental duplications at the breakpoints and including REEP1

We identified an identical and recurrent 9.4-Mbp deletion at chromosome bands 2p11.2-2p12, which occurred de novo in two unrelated patients. It is flanked at the distal and proximal breakpoints by two homologous segmental duplications consisting of low copy repeat (LCR) blocks in direct orientation, which have >99% sequence identity. Despite the fact that the deletion was almost 10 Mbp in size, the patients showed a relatively mild clinical phenotype, that is, mild-to-moderate intellectual disability, a happy disposition, speech delay and delayed motor development. Their phenotype matches with that of previously described patients. The 2p11.2-2p12 deletion includes the REEP1 gene that is associated with spastic paraplegia and phenotypic features related to this are apparent in most 2p11.2-2p12 deletion patients, but not in all. Other hemizygous genes that may contribute to the clinical phenotype include LRRTM1 and CTNNA2. We propose a recurrent but rare 2p11.2-2p12 deletion syndrome based on (1) the identical, non-random localisation of the de novo deletion breakpoints in two unrelated patients and a patient from literature, (2) the patients' phenotypic similarity and their phenotypic overlap with other 2p deletions and (3) the presence of highly identical LCR blocks flanking both breakpoints, consistent with a non-allelic homologous recombination (NAHR)-mediated rearrangement.

Genomic characterization of large heterochromatic gaps in the human genome assembly

The largest gaps in the human genome assembly correspond to multi-megabase heterochromatic regions composed primarily of two related families of tandem repeats, Human Satellites 2 and 3 (HSat2,3). The abundance of repetitive DNA in these regions challenges standard mapping and assembly algorithms, and as a result, the sequence composition and potential biological functions of these regions remain largely unexplored. Furthermore, existing genomic tools designed to predict consensus-based descriptions of repeat families cannot be readily applied to complex satellite repeats such as HSat2,3, which lack a consistent repeat unit reference sequence. Here we present an alignment-free method to characterize complex satellites using whole-genome shotgun read datasets. Utilizing this approach, we classify HSat2,3 sequences into fourteen subfamilies and predict their chromosomal distributions, resulting in a comprehensive satellite reference database to further enable genomic studies of heterochromatic regions. We also identify 1.3 Mb of non-repetitive sequence interspersed with HSat2,3 across 17 unmapped assembly scaffolds, including eight annotated gene predictions. Finally, we apply our satellite reference database to high-throughput sequence data from 396 males to estimate array size variation of the predominant HSat3 array on the Y chromosome, confirming that satellite array sizes can vary between individuals over an order of magnitude (7 to 98 Mb) and further demonstrating that array sizes are distributed differently within distinct Y haplogroups. In summary, we present a novel framework for generating initial reference databases for unassembled genomic regions enriched with complex satellite DNA, and we further demonstrate the utility of these reference databases for studying patterns of sequence variation within human populations.

At least 5–10% of the human genome remains unassembled, unmapped, and poorly characterized. The reference assembly annotates these missing regions as multi-megabase heterochromatic gaps, found primarily near centromeres and on the short arms of the acrocentric chromosomes. This missing fraction of the genome consists predominantly of long arrays of near-identical tandem repeats called satellite DNA. Due to the repetitive nature of satellite DNA, sequence assembly algorithms cannot uniquely align overlapping sequence reads, and thus satellite-rich domains have been omitted from the reference assembly and from most genome-wide studies of variation and function. Existing methods for analyzing some satellite DNAs cannot be easily extended to a large portion of satellites whose repeat structures are complex and largely uncharacterized, such as Human Satellites 2 and 3 (HSat2,3). Here we characterize HSat2,3 using a novel approach that does not depend on having a well-defined repeat structure. By classifying genome-wide HSat2,3 sequences into subfamilies and localizing them to chromosomes, we have generated an initial HSat2,3 genomic reference, which serves as a critical foundation for future studies of variation and function in these regions. This approach should be generally applicable to other classes of satellite DNA, in both the human genome and other complex genomes.

The past, present, and future of human centromere genomics

The centromere is the chromosomal locus essential for chromosome inheritance and genome stability. Human centromeres are located at repetitive alpha satellite DNA arrays that compose approximately 5% of the genome. Contiguous alpha satellite DNA sequence is absent from the assembled reference genome, limiting current understanding of centromere organization and function. Here, we review the progress in centromere genomics spanning the discovery of the sequence to its molecular characterization and the work done during the Human Genome Project era to elucidate alpha satellite structure and sequence variation. We discuss exciting recent advances in alpha satellite sequence assembly that have provided important insight into the abundance and complex organization of this sequence on human chromosomes. In light of these new findings, we offer perspectives for future studies of human centromere assembly and function.

High-throughput genome scaffolding from in vivo DNA interaction frequency

Despite advances in DNA-sequencing technology, assembly of complex genomes remains a major challenge, particularly for genomes sequenced using short reads, which yield highly fragmented assemblies. Here we show that genome-wide in vivo chromatin interaction frequency data, which are measurable with chromosome conformation capture–based experiments, can be used as genomic distance proxies to accurately position individual contigs without requiring any sequence overlap. We also use these data to construct approximate genome scaffolds de novo. Applying our approach to incomplete regions of the human genome, we predict the positions of 65 previously unplaced contigs, in agreement with alternative methods in 26/31 cases attempted in common. Our approach can theoretically bridge any gap size and should be applicable to any species for which global chromatin interaction data can be generated.

Hominoid fission of chromosome 14/15 and the role of segmental duplications

Ape chromosomes homologous to human chromosomes 14 and 15 were generated by a fission event of an ancestral submetacentric chromosome, where the two chromosomes were joined head-to-tail. The hominoid ancestral chromosome most closely resembles the macaque chromosome 7. In this work, we provide insights into the evolution of human chromosomes 14 and 15, performing a comparative study between macaque boundary region 14/15 and the orthologous human regions. We construct a 1.6-Mb contig of macaque BAC clones in the region orthologous to the ancestral hominoid fission site and use it to define the structural changes that occurred on human 14q pericentromeric and 15q subtelomeric regions. We characterize the novel euchromatin–heterochromatin transition region (∼20 Mb) acquired during the neocentromere establishment on chromosome 14, and find it was mainly derived through pericentromeric duplications from ancestral hominoid chromosomes homologous to human 2q14–qter and 10. Further, we show a relationship between evolutionary hotspots and low-copy repeat loci for chromosome 15, revealing a possible role of segmental duplications not only in mediating but also in “stitching” together rearrangement breakpoints.

Segmental duplications in the silkworm genome

Background

Segmental duplications (SDs) or low-copy repeats play important roles in both gene and genome evolution. SDs have been extensively investigated in many organisms, however, there is no information about SDs in the silkworm, Bombyx mori.

Result

In this study, we identified and annotated the SDs in the silkworm genome. Our results suggested that SDs constitute ~1.4% of the silkworm genome sequence (≥1 kb in length and ≥90% in the identity of sequence); the number is similar to that in Drosophila melanogaster but smaller than mammalian organisms. Almost half (42%) of the SD sequences are not assigned to chromosomes, indicating that the SDs are challenges for the assembling of genome sequences. We also provided experimental validation of large duplications using qPCR. The analysis of SD content indicated that the genes related to immunity, detoxification, reproduction, and environmental signal recognition are significantly enriched in the silkworm SDs.

Conclusion

Our results suggested that segmental duplications have been problematic for sequencing and assembling of the silkworm genome. SDs may have important biological significances in immunity, detoxification, reproduction, and environmental signal recognition in the silkworm. This study provides insight into the evolution of the silkworm genome and an invaluable resource for insect genomics research.

Identifying gene locus associations with promyelocytic leukemia nuclear bodies using immuno-TRAP

Important insights into nuclear function would arise if gene loci physically interacting with particular subnuclear domains could be readily identified. Immunofluorescence microscopy combined with fluorescence in situ hybridization (immuno-FISH), the method that would typically be used in such a study, is limited by spatial resolution and requires prior assumptions for selecting genes to probe. Our new technique, immuno-TRAP, overcomes these limitations. Using promyelocytic leukemia nuclear bodies (PML NBs) as a model, we used immuno-TRAP to determine if specific genes localize within molecular dimensions with these bodies. Although we confirmed a TP53 gene-PML NB association, immuno-TRAP allowed us to uncover novel locus-PML NB associations, including the ABCA7 and TFF1 loci and, most surprisingly, the PML locus itself. These associations were cell type specific and reflected the cell's physiological state. Combined with microarrays or deep sequencing, immuno-TRAP provides powerful opportunities for identifying gene locus associations with potentially any nuclear subcompartment.

Megalencephaly syndromes and activating mutations in the PI3K-AKT pathway: MPPH and MCAP

The megalencephaly-polymicrogyria-polydactyly-hydrocephalus (MPPH) and megalencephaly-capillary malformation (MCAP) syndromes are highly recognizable and partly overlapping disorders of brain overgrowth (megalencephaly). Both syndromes are characterized by congenital or early postnatal megalencephaly, with a high risk for progressive ventriculomegaly leading to hydrocephalus and cerebellar tonsillar ectopia leading to Chiari malformation, and cortical brain abnormalities, specifically polymicrogyria. MCAP is further characterized by distinct cutaneous capillary malformations, finger or toe syndactyly, postaxial polydactyly, variable connective tissue dysplasia and mild focal or segmental body overgrowth, among other features. MPPH, on the other hand, lacks consistent vascular or somatic manifestations besides postaxial polydactyly in almost half of reported individuals. We identified de novo germline mutations in PIK3R2 and AKT3 in individuals with MPPH, and both postzygotic, mosaic and rare germline mutations in PIK3CA in individuals with MCAP. PIK3R2, AKT3, and PIK3CA are members of the critical phosphatidylinositol-3-kinase (PI3K)-vakt murine thymoma viral oncogene homolog (AKT) pathway that is well implicated in cell growth, proliferation, survival, apoptosis, among other diverse cellular functions. The identified mutations in these three genes have been shown to lead to gain of function and activation of the PI3K-AKT pathway. Germline and postzygotic mutations of PIK3CA and other PI3K-AKT-mTOR pathway genes have also been identified in several other overgrowth syndromes, highlighting the key role of this signaling pathway in normal development and pathophysiology of a large group of congenital anomalies.

Using population admixture to help complete maps of the human genome

Tens of millions of base pairs of euchromatic human genome sequence, including many protein-coding genes, have no known location in the human genome. We describe an approach for localizing the human genome's missing pieces using the patterns of genome sequence variation created by population admixture. We mapped the locations of 70 scaffolds spanning 4 million base pairs of the human genome's unplaced euchromatic sequence, including more than a dozen protein-coding genes, and identified 8 new large interchromosomal segmental duplications. We find that most of these sequences are hidden in the genome's heterochromatin, particularly its pericentromeric regions. Many cryptic, pericentromeric genes are expressed at the RNA level and have been maintained intact for millions of years while their expression patterns diverged from those of paralogous genes elsewhere in the genome. We describe how knowledge of the locations of these sequences can inform disease association and genome biology studies.

The UCSC Genome Browser

The University of California Santa Cruz (UCSC) Genome Browser is a popular Web-based tool for quickly displaying a requested portion of a genome at any scale, accompanied by a series of aligned annotation "tracks." The annotations generated by the UCSC Genome Bioinformatics Group and external collaborators include gene predictions, mRNA and expressed sequence tag alignments, simple nucleotide polymorphisms, expression and regulatory data, phenotype and variation data, and pairwise and multiple-species comparative genomics data. All information relevant to a region is presented in one window, facilitating biological analysis and interpretation. The database tables underlying the Genome Browser tracks can be viewed, downloaded, and manipulated using another Web-based application, the UCSC Table Browser. Users can upload personal datasets in a wide variety of formats as custom annotation tracks in both browsers for research or educational purposes. This unit describes how to use the Genome Browser and Table Browser for genome analysis, download the underlying database tables, and create and display custom annotation tracks.

Complexity in cancer biology: is systems biology the answer?

Complex phenotypes emerge from the interactions of thousands of macromolecules that are organized in multimolecular complexes and interacting functional modules. In turn, modules form functional networks in health and disease. Omics approaches collect data on changes for all genes and proteins and statistical analysis attempts to uncover the functional modules that perform the functions that characterize higher levels of biological organization. Systems biology attempts to transcend the study of individual genes/proteins and to integrate them into higher order information. Cancer cells exhibit defective genetic and epigenetic networks formed by altered complexes and network modules arising in different parts of tumor tissues that sustain autonomous cell behavior which ultimately lead tumor growth. We suggest that an understanding of tumor behavior must address not only molecular but also, and more importantly, tumor cell heterogeneity, by considering cancer tissue genetic and epigenetic networks, by characterizing changes in the types, composition, and interactions of complexes and networks in the different parts of tumor tissues, and by identifying critical hubs that connect them in time and space.

Self-assembly of thin plates from micrococcal nuclease-digested chromatin of metaphase chromosomes

The three-dimensional organization of the enormously long DNA molecules packaged within metaphase chromosomes has been one of the most elusive problems in structural biology. Chromosomal DNA is associated with histones and different structural models consider that the resulting long chromatin fibers are folded forming loops or more irregular three-dimensional networks. Here, we report that fragments of chromatin fibers obtained from human metaphase chromosomes digested with micrococcal nuclease associate spontaneously forming multilaminar platelike structures. These self-assembled structures are identical to the thin plates found previously in partially denatured chromosomes. Under metaphase ionic conditions, the fragments that are initially folded forming the typical 30-nm chromatin fibers are untwisted and incorporated into growing plates. Large plates can be self-assembled from very short chromatin fragments, indicating that metaphase chromatin has a high tendency to generate plates even when there are many discontinuities in the DNA chain. Self-assembly at 37°C favors the formation of thick plates having many layers. All these results demonstrate conclusively that metaphase chromatin has the intrinsic capacity to self-organize as a multilayered planar structure. A chromosome structure consistent of many stacked layers of planar chromatin avoids random entanglement of DNA, and gives compactness and a high physical consistency to chromatids.

Hyperspectral fluorescence microfluidic (HFM) microscopy

We present an imaging system that collects hyperspectral images of cells travelling through a microfluidic channel. Using a single monochrome camera and a linear variable bandpass filter (LVF), the system captures a bright field image and a set of hyperspectral fluorescence images for each cell. While the bandwidth of the LVF is 20 nm, we have demonstrated that we can determine the peak wavelength of a fluorescent object's emission spectrum with an accuracy of below 3 nm. In addition, we have used this system to capture fluorescence spectra of individual spatially resolved cellular organelles and to spectrally resolve multiple fluorophores in individual cells.

Clone DB: an integrated NCBI resource for clone-associated data

The National Center for Biotechnology Information (NCBI) Clone DB (http://www.ncbi.nlm.nih.gov/clone/) is an integrated resource providing information about and facilitating access to clones, which serve as valuable research reagents in many fields, including genome sequencing and variation analysis. Clone DB represents an expansion and replacement of the former NCBI Clone Registry and has records for genomic and cell-based libraries and clones representing more than 100 different eukaryotic taxa. Records provide details of library construction, associated sequences, map positions and information about resource distribution. Clone DB is indexed in the NCBI Entrez system and can be queried by fields that include organism, clone name, gene name and sequence identifier. Whenever possible, genomic clones are mapped to reference assemblies and their map positions provided in clone records. Clones mapping to specific genomic regions can also be searched for using the NCBI Clone Finder tool, which accepts queries based on sequence coordinates or features such as gene or transcript names. Clone DB makes reports of library, clone and placement data on its FTP site available for download. With Clone DB, users now have available to them a centralized resource that provides them with the tools they will need to make use of these important research reagents.

Identification of nickel response genes in abnormal early developments of sea urchin by differential display polymerase chain reaction

Bioassays and biomarkers have been previously developed to assess the effects of heavy metal contaminants on the early life stages of the sea urchin. In this study, malformation in the early developmental processes was observed in sea urchin (Strongylocentrotus intermedius) larvae exposed to 10 ppm Ni for over 30 h. The most critical stage at which the triggering of nickel effects takes place is thought to be the blastula stage, which occurs after fertilization in larval development. To investigate the molecular-level responses of sea urchin exposed to heavy metal stress and to explore the differentially expressed genes that are induced or repressed by nickel, differential display polymerase chain reaction (DD-PCR) was used with sea urchin mRNAs. The malformation-related genes expressed in the early life stages of the sea urchin were cloned from larvae exposed to 10 ppm of nickel for 15 h, and accessed via DD-PCR. Sequence analysis results revealed that each of the genes evidenced high homology with EGF2, PCSK9, serine/threonine protein kinase, apolipophorin precursor protein, and MGC80921 protein/transcript variant 2. This result may prove useful in the development of novel biomarkers for the assessment of heavy metal stresses on sea urchin embryos.

Translocations disrupting PHF21A in the Potocki-Shaffer-syndrome region are associated with intellectual disability and craniofacial anomalies

Potocki-Shaffer syndrome (PSS) is a contiguous gene disorder due to the interstitial deletion of band p11.2 of chromosome 11 and is characterized by multiple exostoses, parietal foramina, intellectual disability (ID), and craniofacial anomalies (CFAs). Despite the identification of individual genes responsible for multiple exostoses and parietal foramina in PSS, the identity of the gene(s) associated with the ID and CFA phenotypes has remained elusive. Through characterization of independent subjects with balanced translocations and supportive comparative deletion mapping of PSS subjects, we have uncovered evidence that the ID and CFA phenotypes are both caused by haploinsufficiency of a single gene, PHF21A, at 11p11.2. PHF21A encodes a plant homeodomain finger protein whose murine and zebrafish orthologs are both expressed in a manner consistent with a function in neurofacial and craniofacial development, and suppression of the latter led to both craniofacial abnormalities and neuronal apoptosis. Along with lysine-specific demethylase 1 (LSD1), PHF21A, also known as BHC80, is a component of the BRAF-histone deacetylase complex that represses target-gene transcription. In lymphoblastoid cell lines from two translocation subjects in whom PHF21A was directly disrupted by the respective breakpoints, we observed derepression of the neuronal gene SCN3A and reduced LSD1 occupancy at the SCN3A promoter, supporting a direct functional consequence of PHF21A haploinsufficiency on transcriptional regulation. Our finding that disruption of PHF21A by translocations in the PSS region is associated with ID adds to the growing list of ID-associated genes that emphasize the critical role of transcriptional regulation and chromatin remodeling in normal brain development and cognitive function.

Evaluation of 22q11.2 deletion in Cleft Palate patients

BACKGROUND

Cleft palate is the commonest multifactorial epigenetic disorder with a prevalence of 0.43-2.45 per 1000. The objectives of this study were to evaluate the clinical features and identify the 22q11.2 deletion in patients with cleft palate in Sri Lanka.

MATERIALS AND METHODS

Cleft patients attending a Teaching Hospital in Sri Lanka were recruited for this study. The relevant data were obtained from review of case notes, interviews, and examination of patients according to a standard evaluation sheet. Quantitative multiplex polymerase chain reaction (PCR) was performed to identify the 22q11.2 deletion. A gel documentation system (Bio-Doc) was used to quantify the PCR product following electrophoresis on 0.8% agarose gel.

RESULTS AND CONCLUSION

There were 162 cleft palate patients of whom 59% were females. A total of 92 cleft palate subjects (56.2%) had other associated clinical features. Dysmorphic features (25.27%) and developmental delays (25.27%) were the commonest medical problems encountered. The cleft was limited to the soft palate in 125 patients, while in 25 patients it involved both the hard and the soft palate. There were seven subjects with bifid uvula and five subjects with submucous cleft palate. None of the patients had 22q11.2 deletion in this study population. A multicentered large population-based study is needed to confirm the results of this study and to develop guidelines on the appropriate use of 22q11.2 deletion testing, which are valid for cleft palate patients in Sri Lanka.

Expression profiling of uterine leiomyomata cytogenetic subgroups reveals distinct signatures in matched myometrium: transcriptional profilingof the t(12;14) and evidence in support of predisposing genetic heterogeneity

Uterine leiomyomata (UL), the most common neoplasm in reproductive-age women, are classified into distinct genetic subgroups based on recurrent chromosome abnormalities. To develop a molecular signature of UL with t(12;14)(q14-q15;q23-q24), we took advantage of the multiple UL arising as independent clonal lesions within a single uterus. We compared genome-wide expression levels of t(12;14) UL to non-t(12;14) UL from each of nine women in a paired analysis, with each sample weighted for the percentage of t(12;14) cells to adjust for mosaicism with normal cells. This resulted in a transcriptional profile that confirmed HMGA2, known to be overexpressed in t(12;14) UL, as the most significantly altered gene. Pathway analysis of the differentially expressed genes showed significant association with cell proliferation, particularly G1/S checkpoint regulation. This is consistent with the known larger size of t(12;14) UL relative to karyotypically normal UL or to UL in the deletion 7q22 subgroup. Unsupervised hierarchical clustering demonstrated that patient variability is relatively dominant to the distinction of t(12;14) UL compared with non-t(12;14) UL or of t(12;14) UL compared with del(7q) UL. The paired design we employed is therefore important to produce an accurate t(12;14) UL-specific gene list by removing the confounding effects of genotype and environment. Interestingly, myometrium not only clustered away from the tumors, but generally separated based on associated t(12;14) versus del(7q) status. Nine genes were identified whose expression can distinguish the myometrium origin. This suggests an underlying constitutional genetic predisposition to these somatic changes which could potentially lead to improved personalized management and treatment.

INCONSISTENCIES BETWEEN MAPS OF HUMAN CHROMOSOME 22 CORRELATE WITH INCREASED FREQUENCY OF DISEASE-RELATED LOCI

The relationships between genetic or radiation hybrid (RH) and sequence maps of chromosome 22 have been reconsidered based on the sequence map. Integrated maps have been constructed by retaining only common markers between genetic or RH maps and the sequence map. Local inversions of markers have been detected. Ratios between either genetic or RH distances and sequence-based distances have been calculated for each map interval. Hot zones for recombination or radiation breakage have been delineated by merging together intervals displaying high distance ratios and located close to each other for sequence-constrained RH maps, and for female and male genetic maps. A statistically significant positive correlation was found between the distribution of disease-related genes and the hot zones for recombination or radiation breakage on both female genetic and Stanford-G3 RH maps. This observation indicates that investigation of chromosomal regions displaying inconsistencies between RH or genetic linkage and sequence-based maps can accelerate the initial phase of identification of disease-associated genes.