A Scalable Computational Approach for Simulating Complexes of Multiple Chromosomes.

Significant efforts have been recently made to obtain the three-dimensional (3D) structure of the genome with the goal of understanding how structures may affect gene regulation and expression. Chromosome conformational capture techniques such as Hi-C, have been key in uncovering the quantitative information needed to determine chromatin organization. Complementing these experimental tools, co-polymers theoretical methods are necessary to determine the ensemble of three-dimensional structures associated to the experimental data provided by Hi-C maps. Going beyond just structural information, these theoretical advances also start to provide an understanding of the underlying mechanisms governing genome assembly and function. Recent theoretical work, however, has been focused on single chromosome structures, missing the fact that, in the full nucleus, interactions between chromosomes play a central role in their organization. To overcome this limitation, MiChroM (Minimal Chromatin Model) has been modified to become capable of performing these multi-chromosome simulations. It has been upgraded into a fast and scalable software version, which is able to perform chromosome simulations using GPUs via OpenMM Python API, called Open-MiChroM. To validate the efficiency of this new version, analyses for GM12878 individual autosomes were performed and compared to earlier studies. This validation was followed by multi-chain simulations including the four largest human chromosomes (C1-C4). These simulations demonstrated the full power of this new approach. Comparison to Hi-C data shows that these multiple chromosome interactions are essential for a more accurate agreement with experimental results. Without any changes to the original MiChroM potential, it is now possible to predict experimentally observed inter-chromosome contacts. This scalability of Open-MiChroM allow for more audacious investigations, looking at interactions of multiple chains as well as moving towards higher resolution chromosomes models.

Mutations in Smad-interacting protein 1 gene are responsible for absence of its expression in Hirschsprung's disease.

Hirschsprung's disease (HSCR) is a common congenital malformation of the enteric nervous system. The pathophysiological basis remains unclear. Recently, the SIP1 gene has been recognized as being involved in the pathogenesis of symptomatic HSCR patients with 2q22 chromosomal rearrangement. In this study, mutations in SIP1 were analyzed to explore the relationship between SIP1 and HSCR. All exons of SIP1 were amplified and then analyzed by PCR-restriction fragment length polymorphism (PCR-RFLP) and DNA sequencing. SIP1 expression was determined by immunohistochemistry, Western blot and real-time quantitative PCR. By PCR-RFLP, three different electrophoretic bands of 536, 428 and 257 bp representing different genotypes were demonstrated accordingly. DNA sequencing revealed a heterozygous absence of codon 157 GTG → GTA exchange at exon 7. Simultaneously, exchanges of GCC → ACC at codon 351 and ACC → GCC at codon 395 were also observed in exon 8. All the exchanges caused a missense mutation. By immunohistochemistry, SIP1 was ectopically expressed in the aganglionic segment of HSCR without mutation. For comparison, in HSCR with mutation either in exon 7 or exon 8, SIP1 immunoreactivity disappeared in all structures. The protein and mRNA levels determined by Western blot and real-time quantitative PCR were consistent with that of immunohistochemistry. In summary, mutations of the SIP1 gene were detected in HSCR. These mutations in SIP1 were responsible for the absence of its expression in HSCR and contributed to the pathogenesis of this disease.

Neurological manifestations of 2q31 microdeletion syndrome.

Microdeletion of 2q31 involving the HOXD gene cluster is a rare syndrome. The deletion of the HOXD gene cluster is thought to result in skeletal anomalies in these patients. HOX genes encode highly conserved transcription factors that control cell fate and the regional identities along the primary body and limb axes. We experienced a new patient with 2q31 microdeletion encompassing the HOXD gene cluster and some neighboring genes including the ZNF385B. The patient showed digital anomalies, growth failure, epileptic seizures, and intellectual disability. Magnetic resonance imaging showed delayed myelination and low signal intensity in the basal ganglia. The ZNF385B is a zinc finger protein expressed in brain. Disruption of ZNF385B was suspected to be responsible for the neurological features of this syndrome.

Possible genes responsible for developmental delay observed in patients with rare 2q23q24 microdeletion syndrome: Literature review and description of an additional patient.

Cases of 2q23q24 microdeletion syndrome are rare. Patients with chromosomal deletions in this region often show language impairment and/or developmental delay of variable severity. Previous genotype-phenotype correlation study suggested GALNT13 and KCNJ3 as possible candidate genes for such phenotypes. We identified a new overlapping deletion in a patient with severe developmental delay. The identified deletion extended toward the distal 2q24.1 region, and more severe phenotypes in the present patient were considered to be related to the additionally deleted genes including NR4A2 and GPD2. Previously reported chromosomal translocation and the mutation identified in GPD2 suggested that this gene would be responsible for the developmental delay. Re-evaluation for the critical region for behavior abnormalities commonly observed in the patients with overlapping deletions of this region suggested that KCNJ3 rather than GALNT13 may be responsible for abnormal behaviors, although there was phenotypic variability. Combinatory deletions involving KCNJ3 and GPD2 may lead to more severe developmental delay. Further studies would be necessary to establish clearer genotype-phenotype correlation in patients with 2q23q24 microdeletion syndrome.

Deletion of protein tyrosine phosphatase, non-receptor type 4 (PTPN4) in twins with a Rett syndrome-like phenotype.

Rett syndrome (RTT), a neurodevelopmental disorder that predominantly affects females, is primarily caused by variants in MECP2. Variants in other genes such as CDKL5 and FOXG1 are usually associated with individuals who manifest distinct phenotypes that may overlap with RTT. Individuals with phenotypes suggestive of RTT are typically screened for variants in MECP2 and then subsequently the other genes dependent on the specific phenotype. Even with this screening strategy, there are individuals in whom no causative variant can be identified, suggesting that there are other novel genes that contribute to the RTT phenotype. Here we report a de novo deletion of protein tyrosine phosphatase, non-receptor type 4 (PTPN4) in identical twins with a RTT-like phenotype. We also demonstrate the reduced expression of Ptpn4 in a Mecp2 null mouse model of RTT, as well as the activation of the PTPN4 promoter by MeCP2. Our findings suggest that PTPN4 should be considered for addition to the growing list of genes that warrant screening in individuals with a RTT-like phenotype.

Genome-wide association study identifies variants in PMS1 associated with serum ferritin in a Chinese population.

Only a small proportion of genetic variation in serum ferritin has been explained by variant genetic studies, and genome-wide association study (GWAS) for serum ferritin has not been investigated widely in Chinese population. We aimed at exploring the novel genetic susceptibility to serum ferritin, and performed this two stage GWAS in a healthy Chinese population of 3,495 men aged 20-69 y, including 1,999 unrelated subjects in the first stage and 1,496 independent individuals in the second stage. Serum ferritin was measured with electrochemiluminescence immunoassay, and DNA samples were collected for genotyping. A total of 1,940,243 SNPs were tested by using multivariate linear regression analysis. After adjusting for population stratification, age and BMI, the rs5742933 located in the 5'UTR region of PMS1 gene on chromosome 2 was the most significantly associated with ferritin concentrations (P-combined  = 2.329×10(-10)) (ß â€Š=  -0.11, 95% CI: -0.14, -0.07). Moreover, this marker was about 200 kb away from the candidate gene SLC40A1 which is responsible for iron export. PMS1 gene was the novel genetic susceptibility to serum ferritin in Chinese males and its relation to SLC40A1 needs further study.

Identification of the genomic insertion site of Pmel-1 TCR α and ß transgenes by next-generation sequencing.

The pmel-1 T cell receptor transgenic mouse has been extensively employed as an ideal model system to study the mechanisms of tumor immunology, CD8+ T cell differentiation, autoimmunity and adoptive immunotherapy. The 'zygosity' of the transgene affects the transgene expression levels and may compromise optimal breeding scheme design. However, the integration sites for the pmel-1 mouse have remained uncharacterized. This is also true for many other commonly used transgenic mice created before the modern era of rapid and inexpensive next-generation sequencing. Here, we show that whole genome sequencing can be used to determine the exact pmel-1 genomic integration site, even with relatively 'shallow' (8X) coverage. The results were used to develop a validated polymerase chain reaction-based genotyping assay. For the first time, we provide a quick and convenient polymerase chain reaction method to determine the dosage of pmel-1 transgene for this freely and publically available mouse resource. We also demonstrate that next-generation sequencing provides a feasible approach for mapping foreign DNA integration sites, even when information of the original vector sequences is only partially known.

Heterogeneity of primary glioblastoma cells in the expression of caspase-8 and the response to TRAIL-induced apoptosis.

Recent studies suggest that cancer stem cells (CSCs) are responsible for cancer resistance to therapies. We therefore investigated how glioblastoma-derived CSCs respond to the treatment of tumor necrosis factor-related apoptosis-inducing ligand (TRAIL). Neurospheres were generated from glioblastomas, characterized for CSC properties including self-renewal, cell differentiation and xenograft formation capacity, and analyzed for TRAIL-induced apoptosis, CASP8 genomic status, and caspase-8 protein expression. The neurosphere NSC326 was sensitive to TRAIL-induced apoptosis as evidenced by cell death and caspase-8, -3, and -7 enzymatic activities. In contrast, however, the neurosphere NSC189 was TRAIL-resistant. G-banding analysis identified five chromosomally distinguishable cell populations in the neurospheres. Fluorescence in situ hybridization revealed the variation of chromosome 2 copy number in these populations and the loss of CASP8 locus in 2q33-34 region in a small set of cell populations in the neurosphere. Immunohistochemistry of NSC189 cell blocks revealed the lack of caspase-8 protein in a subset of neurosphere cells. Western blotting and immunohistochemistry of human glioblastoma tumors demonstrated the expression of caspase-8 protein in the vast majority of the tumors as compared to normal human brain tissues that lack the caspase-8 expression. This study shows heterogeneity of glioblastomas and derived CSCs in the genomic status of CASP8, expression of caspase-8, and thus responsiveness to TRAIL-induced apoptosis. Clinic trials may consider genomic analysis of the cancer tissue to identify the genomic loss of CASP8 and use it as a genomic marker to predict the resistance of glioblastomas to TRAIL apoptosis pathway-targeted therapies.

Dynamics of bivalent chromatin domains upon drug induced reactivation and resilencing in cancer cells.

Epigenetic deregulation revealed by altered profiles of DNA methylation and histone modifications is a frequent event in cancer cells and results in abnormal patterns of gene expression. Cancer silenced genes constitute prime therapeutic targets and considerable progress has been made in the epigenetic characterization of the chromatin scenarios associated with their inactivation and drug induced reactivation. Despite these advances, the mechanisms involved in the maintenance or resetting of epigenetic states in both physiological and pharmacological situations are poorly known. To get insights into the dynamics of chromatin regulation upon drug-induced reactivation, we have investigated the epigenetic profiles of two chromosomal regions undergoing long range epigenetic silencing in colon cancer cells in time-course settings after exposure of cells to chromatin reactivating agents. The DNA methylation states and the balance between histone H3K4 methylation and H3K27 methylation marks clearly define groups of genes with alternative responses to therapy. We show that the expected epigenetic remodeling induced by the reactivating drugs, just achieves a transient disruption of the bivalent states, which overcome the treatment and restore the transcriptional silencing approximately four weeks after drug exposure. The interplay between DNA methylation and bivalent histone marks appears to configure a plastic but stable chromatin scenario that is fully restored in silenced genes after drug withdrawal. These data suggest that improvement of epigenetic therapies may be achieved by designing strategies with long lasting effects.

A novel MERTK deletion is a common founder mutation in the Faroe Islands and is responsible for a high proportion of retinitis pigmentosa cases.

PURPOSE: The aim of the study was to elucidate the genetic background of retinitis pigmentosa (RP) in a Faroe Islands population, a genetic isolate in the North Atlantic Ocean. METHODS: Blood samples were collected from subjects diagnosed with RP and their families. DNA from affected individuals underwent single nucleotide polymorphism microarray analysis and homozygosity mapping followed by sequence analysis of candidate genes. RESULTS: We identified 25 cases of nonsyndromic RP corresponding to a prevalence of 1 in 1,900. Single nucleotide polymorphism analysis revealed a homozygous region on chromosome 2q, common to patients in four families, which harbored the RP gene MER tyrosine kinase protooncogene (MERTK). A deletion of 91 kb was identified in seven patients, representing 30% of the analyzed Faroese cases of nonsyndromic RP. The clinical course of six patients who were homozygous for the deletion showed onset in the first decade followed by a rapid deterioration of both rod and cone photoreceptor function. Early macular involvement was present, in accordance with that of other reported patients with MERTK mutations. CONCLUSIONS: Previous studies have shown a frequency of less than 1% of MERTK mutations in RP patients. The 91-kb deletion encompassing exons 1-7 of MERTK is a common founder mutation in the Faroe Islands, responsible for around 30% of RP, and together with mutations in protocadherin 21 (PCDH21) accounts for more than half of the retinal dystrophy cases.

A recurrent mutation in CRYGD is associated with autosomal dominant congenital coralliform cataract in two unrelated Chinese families.

PURPOSE: Congenital cataract is a clinically and genetically heterogeneous lens disorder. The purpose of this study was to identify the mutation responsible for autosomal dominant congenital coralliform cataracts in two Chinese families and to investigate the relationship between virulence genes and lens morphology. METHODS: Patients received a physical examination, and blood samples were collected for DNA extraction. Mutation analysis was performed by direct sequencing of the candidate genes: gammaC-crystallin (CRYGC), gammaD-crystallin (CRYGD), gammaS-crystallin (CRYGS), gap-junction protein, alpha 8 (GJA8), gap-junction protein, alpha 3 (GJA3), and alphaA-crystallin (CRYAA). RESULTS: The affected individuals in two families had congenital coralliform cataracts. Mutational analysis of the CRYGD identified a C→A transversion at nucleotide position c.70 in exon 2, which resulted in a threonine substitution for proline at amino-acid residue 24 (P24T). This mutation was identified in all affected individuals but was not found in healthy relatives or 100 control chromosomes from the same ethnic background. CONCLUSIONS: Our results indicated that the P24T mutation of CRYGD was responsible for two Chinese pedigrees with congenital coralliform cataracts. CRYGD and coralliform cataracts are highly related, and P24T may be a hot-point mutation for this disorder.

Gene content and function of the ancestral chromosome fusion site in human chromosome 2q13-2q14.1 and paralogous regions.

Various portions of the region surrounding the site where two ancestral chromosomes fused to form human chromosome 2 are duplicated elsewhere in the human genome, primarily in subtelomeric and pericentromeric locations. At least 24 potentially functional genes and 16 pseudogenes reside in the 614-kb of sequence surrounding the fusion site and paralogous segments on other chromosomes. By comparing the sequences of genomic copies and transcripts, we show that at least 18 of the genes in these paralogous regions are transcriptionally active. Among these genes are new members of the cobalamin synthetase W domain (CBWD) and forkhead domain FOXD4 gene families. Copies of RPL23A and SNRPA1 on chromosome 2 are retrotransposed-processed pseudogenes that were included in segmental duplications; we find 53 RPL23A pseudogenes in the human genome and map the functional copy of SNRPA1 to 15qter. The draft sequence of the human genome also provides new information on the location and intron-exon structure of functional copies of other 2q-fusion genes (PGM5, retina-specific F379, helicase CHLR1, and acrosin). This study illustrates that the duplication and rearrangement of subtelomeric and pericentromeric regions have functional relevance to human biology; these processes can change gene dosage and/or generate genes with new functions.

Identification and characterization of a novel ABCA subfamily member, ABCA12, located in the lamellar ichthyosis region on 2q34.

The ABCA subfamily of ABC transporters includes ten members to date. In this study, we describe an additional gene, ABCA12. Four full-length cDNA sequences have been obtained from human placenta that contain two different polyadenylation sites and two splicing forms, coding for ABCA12 isoforms of 2,595 and 2,516 amino acid residues. Both isoforms are predicted to have two ATP-binding domains (nucleotide binding domain, NBD) and two transmembrane (TM) domains, features shared by all other ABCA subfamily proteins. ABCA12 is most closely related to ABCA1, with an amino acid similarity of 47%. Northern blot analysis demonstrates that a 9.5-kb transcript is mainly expressed in the stom- ach. ABCA12 was mapped to human chromosome 2q34. Two other genes from ABCA subfamily are associated with human inherited diseases, ABCA1 with the cholesterol transport disorders Tangier disease and familial hypoalphalipoproteinemia, and ABCA4 with several retinal degeneration disorders. The ABCA12 gene is located in a region of chromosome 2q34 that harbors the genes for lamellar ichthyosis, polymorphic congenital cataract, and insulin-dependent diabetes mellitus (IDDM13), and therefore is a positional candidate for these pathologies.

Mapping of human DNA-binding nuclear protein (NP220) to chromosome band 2p13.1-p13.2 and its relation to matrin 3.

Human NP220 (hNP220) is a novel DNA-binding nuclear protein, which has an arginine/serine-rich motif and polypyrimidine tract-binding motif, and NP220s and matrin 3 are thought to form a novel family of nuclear proteins. We have determined a chromosomal localization of the cDNA encoding human NP220 to 2p13.1-p13.2 by using fluorescence in situ hybridization. Human matrin 3 cDNA was mapped to chromosomes 1p13.1-p21.1 and 5q31.3, demonstrating that these novel nuclear proteins with similar functions are on different chromosomes.

The atomic force microscope as a new microdissecting tool for the generation of genetic probes.

The atomic force microscope (AFM) can be used to visualize and to manipulate biological material with relative case and high resolution. This study was carried out to investigate whether probe sets, specific for subregions of the human genome and useful for the painting of chromosome bands, can be established by PCR amplification of AFM-dissected chromosome regions. Compared to standard microdissection techniques, the AFM can be used with much higher precision for the dissection of the region of interest and subsequent nanoextraction of DNA material. After scanning the area of interest in noncontact mode AFM, chromosome bands were cut by the AFM tip at high force. The genetic material of a single cut attached itself to the tip and was extracted and amplified using degenerate oligonucleotide-primed-PCR. Subsequent to hapten labeling, fluorescence in situ hybridization was performed and chromosome band-specific probes were visualized by standard fluorescence microscopy.

Comparative fluorescence in situ hybridization mapping of primate chromosomes with Alu polymerase chain reaction generated probes from human/rodent somatic cell hybrids.

We have used Alu polymerase chain reaction generated probes from rearranged human/rodent somatic cell hybrids for fluorescence in situ hybridization and comparative mapping of some intrachromosomal changes in the karyotypes of great apes (Pan troglodytes, P. paniscus, Gorilla gorilla, Pongo pygmaeus), a gibbon (Hylobates lar), and an Old World monkey (Macaca fuscata). Probes containing chromosomes 2 and 18 fragments confirmed inversions already suggested by the banding pattern of great ape homologues. However, a chromosome 3 fragment showed complex rearrangements in the gibbon and macaque karyotype which were previously not well defined from banding. 'Subchromosomal painting' will allow the identification of intrachromosomal changes on the basis of DNA homology and provides a powerful method to study karyological and genomic evolution.

A region-specific microdissection library for human chromosome 2p23-p25 and the analysis of an interstitial deletion of 2p23.3-p25.1.

A region-specific library for human chromosome 2p23-p25 was constructed using microdissection and polymerase chain reaction (PCR)-mediated microcloning techniques. This library is large, comprising 300,000 recombinant microclones. The insert sizes range between 50-600 base pairs (bp) with a mean of 200 bp. About 50%-60% of the clones contain unique or very low copy number sequence inserts as determined by their weak or no hybridization to total human DNA. A subset of 48 microclones that did not hybridize to total human DNA after colony hybridization was analyzed, and 26 (54%) clones were shown to contain single-copy inserts and hybridize to human chromosome 2 DNAs, indicating that they are human chromosome 2 specific. The human genomic fragments identified by these clones after cleavage with HindIII have also been characterized. The single-copy microclones were used to analyze an interstitial deletion in the 2p23.3-p25.1 region--46,XY, del(2) (pter-->p25.1::p23.3-->qter)--previously reported in a patient with severe growth and mental retardation and multiple anomalies. Of the 26 microclones analyzed, 14 clones were mapped to the deletion region. The availability of the 2p23-p25 region-specific library and the probes derived from the library should be valuable for fine structure physical mapping analysis and the cloning of disease-related genes localized to the region. These studies also demonstrate the efficiency with which useful probes can be quickly generated for genome studies and for positional cloning.

Chromosomal mapping of the human interleukin-1 receptor antagonist gene (IL-1RN) and isolation of specific YAC clones.

Using a panel of somatic rodent-human cell hybrids, we show that the interleukin-1 receptor antagonist gene (IL-1RN) maps to the long arm of human chromosome 2. Linkage studies permitted the regional localization of this gene to band q14-21. This is the same region in which the IL-1 alpha and IL-1 beta genes are localized. Three yeast artificial chromosome (YAC) clones containing the IL-1RN gene were isolated, and these will be used for further characterization of this chromosome 2 region.

Human immunoglobulin genes of the kappa type. The long-range map of an orphon V kappa gene region.

As was previously shown by Zimmer et al. (EMBO J. 9, 1535-1542, 1990 and Biol. Chem. Hoppe-Seyler 371, 939-951, 1990), the so-called W regions comprising 11 V kappa pseudogenes are located on the long arm of chromosome 2, very closely to the centromere. They are probably derived by a pericentric inversion and amplification events from gene regions of the kappa locus, which is located on the short arm of chromosome 2 also very closely to the centromere. The restriction map of the W regions was now extended from the previous 1.2 Mb to 4.3 Mb and, at the same time, revised with respect to certain features. This was made possible by a new hybridization probe specific for the Wc region and by the improved resolution and extended range of pulsed field gel electrophoresis. On the basis of the long-range maps of the W regions and the kappa locus the V kappa genes of the kappa locus have to be at least 2.5 Mb apart. This distance can be taken also as a minimal estimate for the size of the centromere DNA of chromosome 2; it is quite possible that the size is much larger.

A simple and efficient amplification method of DNA with unknown sequences and its application to microdissection/microcloning.

An alternative method for amplification of DNA with unknown sequences was developed. This involves the direct ligation of a primer oligodeoxyribonucleotide itself to restricted DNA fragments with unknown sequences to be amplified by PCR. The oligonucleotide need not be phosphorylated and need not be annealed with its complementary oligonucleotide in advance for ligation. The ligation reaction seems to be independent of the concentration of unknown DNA, proceeds in short time, and is efficient. The ligation efficiency was more than 30% at a low concentration, 10 fg/microliters, of DNA. This method was applied to a microdissection/microcloning of the short arm of human chromosome 2. Of 65 clones screened for the highly repetitive sequences with total human genomic DNA, eleven (17%) were positive. Their inserts ranged in size from 150 to 1,200 bp (average, 460 bp). In Southern blot analysis, thirty consecutive clones all detected signals common to both total human genomic DNA and mouse-human hybrid cell DNA containing only chromosome 2 of human origin. Among them, 24 (80%) were unique sequences, and 6 (20%) were multi-copy (or intermediate-repeat) sequences. Thus, this method is simple and efficient, and provides an alternative way to amplify unknown DNA.