Human Autosomes Research Articles

ASAR lncRNAs control DNA replication timing through interactions with multiple hnRNP/RNA binding proteins.

ASARs are a family of very-long noncoding RNAs that control replication timing on individual human autosomes, and are essential for chromosome stability. The eight known ASAR lncRNAs remain closely associated with their parent chromosomes. Analysis of RNA-protein interaction data (from ENCODE) revealed numerous RBPs with significant interactions with multiple ASAR lncRNAs, with several hnRNPs as abundant interactors. An ~7 kb domain within the ASAR6-141 lncRNA shows a striking density of RBP interaction sites. Genetic deletion and ectopic integration assays indicate that this ~7 kb RNA binding protein domain contains functional sequences for controlling replication timing of entire chromosomes in cis. shRNA-mediated depletion of 10 different RNA binding proteins, including HNRNPA1, HNRNPC, HNRNPL, HNRNPM, HNRNPU, or HNRNPUL1, results in dissociation of ASAR lncRNAs from their chromosome territories, and disrupts the synchronous replication that occurs on all autosome pairs, recapitulating the effect of individual ASAR knockouts on a genome-wide scale. Our results further demonstrate the role that ASARs play during the temporal order of genome-wide replication, and we propose that ASARs function as essential RNA scaffolds for the assembly of hnRNP complexes that help maintain the structural integrity of each mammalian chromosome.

ASAR lncRNAs control DNA replication timing through interactions with multiple hnRNP/RNA binding proteins

ASARs are a family of very-long noncoding RNAs that control replication timing on individual human autosomes, and are essential for chromosome stability. The eight known ASAR lncRNAs remain closely associated with their parent chromosomes. Analysis of RNA-protein interaction data (from ENCODE) revealed numerous RBPs with significant interactions with multiple ASAR lncRNAs, with several hnRNPs as abundant interactors. An ~7 kb domain within the ASAR6-141 lncRNA shows a striking density of RBP interaction sites. Genetic deletion and ectopic integration assays indicate that this ~7 kb RNA binding protein domain contains functional sequences for controlling replication timing of entire chromosomes in cis. shRNA-mediated depletion of 10 different RNA binding proteins, including HNRNPA1, HNRNPC, HNRNPL, HNRNPM, HNRNPU, or HNRNPUL1, results in dissociation of ASAR lncRNAs from their chromosome territories, and disrupts the synchronous replication that occurs on all autosome pairs, recapitulating the effect of individual ASAR knockouts on a genome-wide scale. Our results further demonstrate the role that ASARs play during the temporal order of genome-wide replication, and we propose that ASARs function as essential RNA scaffolds for the assembly of hnRNP complexes that help maintain the structural integrity of each mammalian chromosome.

On latent idealized models in symbolic datasets: unveiling signals in noisy sequencing data.

Data taking values on discrete sample spaces are the embodiment of modern biological research. "Omics" experiments based on high-throughput sequencing produce millions of symbolic outcomes in the form of reads (i.e., DNA sequences of a few dozens to a few hundred nucleotides). Unfortunately, these intrinsically non-numerical datasets often deviate dramatically from natural assumptions a practitioner might make, and the possible sources of this deviation are usually poorly characterized. This contrasts with numerical datasets where Gaussian-type errors are often well-justified. To overcome this hurdle, we introduce the notion of latent weight, which measures the largest expected fraction of samples from a probabilistic source that conform to a model in a class of idealized models. We examine various properties of latent weights, which we specialize to the class of exchangeable probability distributions. As proof of concept, we analyze DNA methylation data from the 22 human autosome pairs. Contrary to what is usually assumed in the literature, we provide strong evidence that highly specific methylation patterns are overrepresented at some genomic locations when latent weights are taken into account.

장애인 운동선수군에서 4종류의 후보유전자에 존재하는 유전자 다형성에 대한 분포 양상

Although multiple genetic polymorphisms are reported to affect the sports performance in athletes, but very little research has been performed on this field in the case of disability athletes. Thus, in this study, after selecting four kinds of candidate genes known to be associated with sports performances by previous studies, we analyzed whether the distributions of genetic polymorphisms within these candidate genes indicate the significant differences between controls and disability athletes, respectively. A total of 68 unrelated subjects were participated in this study, and they divided into controls and disability athletes showing the same number(34 people in both group) and gender distribution(each 17 male/17 female in both group), respectively. Also, the genetic polymorphisms in a kind of four candidate including AMPD1(C34T, rs17602729), MCT1(SLC16A1)(A1470T, rs1049434) and MMP3(G/A, rs679620) in human autosome, and AGTR2(C/A, rs11091046) in human X chromosome were analyzed by using SNaPshot method in both groups. Among genetic polymorphisms studied, in the case of rs679620 polymorphism within MMP3 gene associated with sports injury by previous studies, there was the marginally significant difference in allele frequences between both groups(x=df=1, p=.051), and especially, G allele in this genetic polymorphism was higher in disability athletes than in controls. Therefore, in the future. it seems necessary to investigate the distribution pattern of this genetic polymorphism in a disability group with sports injuries, or to investigate the distribution pattern of other genetic polymorphisms related to various sports injuries such as anterior cruciate ligamnet(ACL) rupture and Achilles tendon rupture.

Parent-of-origin detection and chromosome-scale haplotyping using long-read DNA methylation sequencing and Strand-seq

Hundreds of loci in human genomes have alleles that are methylated differentially according to their parent of origin. These imprinted loci generally show little variation across tissues, individuals, and populations. We show that such loci can be used to distinguish the maternal and paternal homologs for all human autosomes without the need for the parental DNA. We integrate methylation-detecting nanopore sequencing with the long-range phase information in Strand-seq data to determine the parent of origin of chromosome-length haplotypes for both DNA sequence and DNA methylation in five trios with diverse genetic backgrounds. The parent of origin was correctly inferred for all autosomes with an average mismatch error rate of 0.31% for SNVs and 1.89% for insertions or deletions (indels). Because our method can determine whether an inherited disease allele originated from the mother or the father, we predict that it will improve the diagnosis and management of many genetic diseases.

Knockdown of LINC01116 inhibits cell migration and invasion in head and neck squamous cell carcinoma through epithelial-mesenchymal transition pathway.

Long noncoding RNAs (lncRNAs) are linked to tumor development and progression. The aim of this study was to determine the prognostic significance and biological role of LINC01116 in head and neck squamous cell carcinoma (HNSCC). We identified 21 aberrantly expressed lncRNAs specific to HNSCC that were common in two microarray datasets. LINC01116 was highly overexpressed in HNSCC tissues and was correlated to shorter overall survival and relapse-free survival duration, as analyzed by the online Gene Expression Profiling Interactive Analysis platform. LINC01116 was also overexpressed in oral squamous cell carcinoma and nasopharyngeal carcinoma tissues, and LINC01116 silencing significantly inhibited the migration and invasion capacities of both cell lines by blocking the epithelial-mesenchymal transition process. In addition, 125 coexpressing genes were identified by circlncRNAnet, and were mainly located on human autosomes and enriched in transforming growth factor-β signaling pathway. These findings indicate that LINC01116 might be a potential therapeutic target for HNSCC.

Evaluation of a microhaplotypes panel for forensic genetics using massive parallel sequencing technology

Massive parallel DNA sequencing (MPS) makes it possible to explore a new type of genetic marker, known as microhaplotypes or microhaps. These loci were recently introduced in the landscape of forensic genetic and appear to be useful for identification purposes, reconstruction of family relationships, ancestry prediction and DNA mixtures deconvolution. Microhaplotypes loci, based on 89 loci in ALFRED, were selected and their genetic variations in 100 Italian individuals were evaluated by using MPS, in order to make inference about utility of a set of microhaps in forensic genetics. After MPS, the panel was reduced to 87 microhaps, comprised of 266 different SNPs and spread across 22 human autosomes. Genotype and haplotype frequencies were estimated, as well as the effective number of alleles at each locus (Ae), which relates to the usefulness of the locus in resolution of relationships and deconvolution of DNA mixtures. Overall, the Ae values for the 87 microhaps range from 1.010 to 8.344, with about 80% showing values greater than 2.0. Noteworthy, 32 microhaps display Ae values greater than 3.0 and 18 loci Ae above 4.0.To explore the suitability of microhaplotypes in mixture deconvolution, the probability of detecting a mixture, as a function of Ae, was inferred for different groups of loci. Considering the fourteen loci with Ae between 3.0 and 3.999 the probability of detecting a mixture was at least 0.99973, while considering the ten loci with Ae between 4.0 and 4.999 the probability was at least 0.99998. Moreover, when considering just the six loci with Ae between 5.0 and 5.999 the probability of detecting a mixture was at least 0.99984, while when considering just the two loci with Ae above 6 the probability was 0.97228. Combining these 32 MH loci, the theoretical probability of detecting a mixture was 0.999999999999973. These results make the subset of 32 loci with Ae above three informative for mixture resolution.The individual matching probabilities (PI) of the 87 microhaps ranged from 0.032 to 0.9802. Considering the 32 microhap loci with Ae greater than 3.0, the cumulative PI value was 1.6 × 10−33, while considering the 18 microhap loci with Ae above 4.0, the cumulative PI value was 2.34 × 10-21. Overall the results of this study confirmed the utility of microhaps in forensic genetics.

IFISH is a publically available resource enabling versatile DNA FISH to study genome architecture

DNA fluorescence in situ hybridization (DNA FISH) is a powerful method to study chromosomal organization in single cells. At present, there is a lack of free resources of DNA FISH probes and probe design tools which can be readily applied. Here, we describe iFISH, an open-source repository currently comprising 380 DNA FISH probes targeting multiple loci on the human autosomes and chromosome X, as well as a genome-wide database of optimally designed oligonucleotides and a freely accessible web interface (http://ifish4u.org) that can be used to design DNA FISH probes. We individually validate 153 probes and take advantage of our probe repository to quantify the extent of intermingling between multiple heterologous chromosome pairs, showing a much higher extent of intermingling in human embryonic stem cells compared to fibroblasts. In conclusion, iFISH is a versatile and expandable resource, which can greatly facilitate the use of DNA FISH in research and diagnostics.

Defining the genetic control of human blood plasma N-glycome using genome-wide association study.

Glycosylation is a common post-translational modification of proteins. Glycosylation is associated with a number of human diseases. Defining genetic factors altering glycosylation may provide a basis for novel approaches to diagnostic and pharmaceutical applications. Here we report a genome-wide association study of the human blood plasma N-glycome composition in up to 3811 people measured by Ultra Performance Liquid Chromatography (UPLC) technology. Starting with the 36 original traits measured by UPLC, we computed an additional 77 derived traits leading to a total of 113 glycan traits. We studied associations between these traits and genetic polymorphisms located on human autosomes. We discovered and replicated 12 loci. This allowed us to demonstrate an overlap in genetic control between total plasma protein and IgG glycosylation. The majority of revealed loci contained genes that encode enzymes directly involved in glycosylation (FUT3/FUT6, FUT8, B3GAT1, ST6GAL1, B4GALT1, ST3GAL4, MGAT3 and MGAT5) and a known regulator of plasma protein fucosylation (HNF1A). However, we also found loci that could possibly reflect other more complex aspects of glycosylation process. Functional genomic annotation suggested the role of several genes including DERL3, CHCHD10, TMEM121, IGH and IKZF1. The hypotheses we generated may serve as a starting point for further functional studies in this research area.

Chromosomal localization of the adrenoleukodystrophy-related gene in man and mice.

We report here on the chromosomal mapping of the adrenoleukodystrophy-related (ALDR) gene on both the human and the mouse genomes. This gene encodes a peroxisomal ATP binding cassette transporter, closely related to the transporter identified as responsible for the adrenoleukodystrophy phenotype. ALDR maps on the syntenic region on murine and human autosomes. In addition, we could determine its position in relation to known microsatellite framework markers; this will allow to test its role in Zellweger syndrome and/ or related peroxisomal disorders.

Cost-effective high-throughput single-haplotype iterative mapping and sequencing for complex genomic structures.

The reference sequences of structurally complex regions can be obtained only through highly accurate clone-based approaches. We and others have successfully used single-haplotype iterative mapping and sequencing (SHIMS) 1.0 to assemble structurally complex regions across the sex chromosomes of several vertebrate species and to allow for targeted improvements to the reference sequences of human autosomes. However, SHIMS 1.0 is expensive and time consuming, requiring resources that only a genome center can provide. Here we introduce SHIMS 2.0, an improved SHIMS protocol that allows even a small laboratory to generate high-quality reference sequence from complex genomic regions. Using a streamlined and parallelized library-preparation protocol, and taking advantage of inexpensive high-throughput short-read-sequencing technologies, a small laboratory with both molecular biology and bioinformatics experience can sequence and assemble 192 large-insert bacterial artificial chromosome (BAC) or fosmid clones in 1 week. In SHIMS 2.0, in contrast to other pooling strategies, each clone is sequenced with a unique barcode, thus enabling clones containing nearly identical sequences to be multiplexed in a single sequencing run and assembled separately. Relative to SHIMS 1.0, SHIMS 2.0 decreases the required cost and time by two orders of magnitude while preserving high sequencing accuracy.

S-plot2: Rapid Visual and Statistical Analysis of Genomic Sequences.

With the daily release of data from whole genome sequencing projects, tools to facilitate comparative studies are hard-pressed to keep pace. Graphical software solutions can readily recognize synteny by measuring similarities between sequences. Nevertheless, regions of dissimilarity can prove to be equally informative; these regions may harbor genes acquired via lateral gene transfer (LGT), signify gene loss or gain, or include coding regions under strong selection. Previously, we developed the software S-plot. This tool employed an alignment-free approach for comparing bacterial genomes and generated a heatmap representing the genomes’ similarities and dissimilarities in nucleotide usage. In prior studies, this tool proved valuable in identifying genome rearrangements as well as exogenous sequences acquired via LGT in several bacterial species. Herein, we present the next generation of this tool, S-plot2. Similar to its predecessor, S-plot2 creates an interactive, 2-dimensional heatmap capturing the similarities and dissimilarities in nucleotide usage between genomic sequences (partial or complete). This new version, however, includes additional metrics for analysis, new reporting options, and integrated BLAST query functionality for the user to interrogate regions of interest. Furthermore, S-plot2 can evaluate larger sequences, including whole eukaryotic chromosomes. To illustrate some of the applications of the tool, 2 case studies are presented. The first examines strain-specific variation across the Pseudomonas aeruginosa genome and strain-specific LGT events. In the second case study, corresponding human, chimpanzee, and rhesus macaque autosomes were studied and lineage specific contributions to divergence were estimated. S-plot2 provides a means to both visually and quantitatively compare nucleotide sequences, from microbial genomes to eukaryotic chromosomes. The case studies presented illustrate just 2 potential applications of the tool, highlighting its capability to identify and investigate the variation in molecular divergence rates across sequences. S-plot2 is freely available through https://bitbucket.org/lkalesinskas/splot and is supported on the Linux and MS Windows operating systems.

Ribosomal DNA copy number amplification and loss in human cancers is linked to tumor genetic context, nucleolus activity, and proliferation.

Ribosomal RNAs (rRNAs) are transcribed from two multicopy DNA arrays: the 5S ribosomal DNA (rDNA) array residing in a single human autosome and the 45S rDNA array residing in five human autosomes. The arrays are among the most variable segments of the genome, exhibit concerted copy number variation (cCNV), encode essential components of the ribosome, and modulate global gene expression. Here we combined whole genome data from >700 tumors and paired normal tissues to provide a portrait of rDNA variation in human tissues and cancers of diverse mutational signatures, including stomach and lung adenocarcinomas, ovarian cancers, and others of the TCGA panel. We show that cancers undergo coupled 5S rDNA array expansion and 45S rDNA loss that is accompanied by increased estimates of proliferation rate and nucleolar activity. These somatic changes in rDNA CN occur in a background of over 10-fold naturally occurring rDNA CN variation across individuals and cCNV of 5S-45S arrays in some but not all tissues. Analysis of genetic context revealed associations between cancer rDNA CN amplification or loss and the presence of specific somatic alterations, including somatic SNPs and copy number gain/losses in protein coding genes across the cancer genome. For instance, somatic inactivation of the tumor suppressor gene TP53 emerged with a strong association with coupled 5S expansion / 45S loss in several cancers. Our results uncover frequent and contrasting changes in the 5S and 45S rDNA along rapidly proliferating cell lineages with high nucleolar activity. We suggest that 5S rDNA amplification facilitates increased proliferation, nucleolar activity, and ribosomal synthesis in cancer, whereas 45S rDNA loss emerges as a byproduct of transcription-replication conflict in rapidly replicating tumor cells. The observations raise the prospects of using the rDNA arrays as re-emerging targets for the design of novel strategies in cancer therapy.

Intricacies in arrangement of SNP haplotypes suggest \u201cGreat Admixture\u201d that created modern humans

BackgroundInferring history from genomic sequences is challenging and problematic because chromosomes are mosaics of thousands of small Identicalby-descent (IBD) fragments, each of them having their own unique story. However, the main events in recent evolution might be deciphered from comparative analysis of numerous loci. A paradox of why humans, whose effective population size is only 104, have nearly three million frequent SNPs is formulated and examined.ResultsWe studied 5398 loci evenly covering all human autosomes. Common haplotypes built from frequent SNPs that are present in people from various populations have been examined. We demonstrated highly non-random arrangement of alleles in common haplotypes. Abundance of mutually exclusive pairs of common haplotypes that have different alleles at every polymorphic position (so-called Yin/Yang haplotypes) was found in 56% of loci. A novel widely spread category of common haplotypes named Mosaic has been described. Mosaic consists of numerous pieces of Yin/Yang haplotypes and represents an ancestral stage of one of them. Scenarios of possible appearance of large number of frequent human SNPs and their habitual arrangement in Yin/Yang common haplotypes have been evaluated with an advanced genomic simulation algorithm.ConclusionsComputer modeling demonstrated that the observed arrangement of 2.9 million frequent SNPs could not originate from a sole stand-alone population. A “Great Admixture” event has been proposed that can explain peculiarities with frequent SNP distributions. This Great Admixture presumably occurred 100–300 thousand years ago between two ancestral populations that had been separated from each other about a million years ago. Our programs and algorithms can be applied to other species to perform evolutionary and comparative genomics.

Spatial Organization of Chromatin Domains and Compartments in Single Chromosomes

The spatial organization of chromatin critically affects genome function. Recent chromosome-conformation-capture studies have revealed topologically associating domains (TADs) as a conserved feature of chromatin organization, but how TADs are spatially organized in individual chromosomes remains unknown. Here, we developed an imaging method for mapping the spatial positions of numerous genomic regions along individual chromosomes and traced the positions of TADs in human interphase autosomes and X chromosomes. We observed that chromosome folding deviates from the ideal fractal-globule model at large length scales and that TADs are largely organized into two compartments spatially arranged in a polarized manner in individual chromosomes. Active and inactive X chromosomes adopt different folding and compartmentalization configurations. These results suggest that the spatial organization of chromatin domains can change in response to regulation.

Regional Differences in the Accumulation of SNPs on the Male-Specific Portion of the Human Y Chromosome Replicate Autosomal Patterns: Implications for Genetic Dating.

Factors affecting the rate and pattern of the mutational process are being identified for human autosomes, but the same relationships for the male specific portion of the Y chromosome (MSY) are not established. We considered 3,390 mutations occurring in 19 sequence bins identified by sequencing 1.5 Mb of the MSY from each of 104 present-day chromosomes. The occurrence of mutations was not proportional to the amount of sequenced bases in each bin, with a 2-fold variation. The regression of the number of mutations per unit sequence against a number of indicators of the genomic features of each bin, revealed the same fundamental patterns as in the autosomes. By considering the sequences of the same region from two precisely dated ancient specimens, we obtained a calibrated region-specific substitution rate of 0.716 × 10-9/site/year. Despite its lack of recombination and other peculiar features, the MSY then resembles the autosomes in displaying a marked regional heterogeneity of the mutation rate. An immediate implication is that a given figure for the substitution rate only makes sense if bound to a specific DNA region. By strictly applying this principle we obtained an unbiased estimate of the antiquity of lineages relevant to the genetic history of the human Y chromosome. In particular, the two deepest nodes of the tree highlight the survival, in Central-Western Africa, of lineages whose coalescence (291 ky, 95% C.I. 253–343) predates the emergence of anatomically modern features in the fossil record.

Characterization of Large Structural Genetic Mosaicism in Human Autosomes

Analyses of genome-wide association study (GWAS) data have revealed that detectable genetic mosaicism involving large (>2 Mb) structural autosomal alterations occurs in a fraction of individuals. We present results for a set of 24,849 genotyped individuals (total GWAS set II [TGSII]) in whom 341 large autosomal abnormalities were observed in 168 (0.68%) individuals. Merging data from the new TGSII set with data from two prior reports (the Gene-Environment Association Studies and the total GWAS set I) generated a large dataset of 127,179 individuals; we then conducted a meta-analysis to investigate the patterns of detectable autosomal mosaicism (n = 1,315 events in 925 [0.73%] individuals). Restricting to events >2 Mb in size, we observed an increase in event frequency as event size decreased. The combined results underscore that the rate of detectable mosaicism increases with age (p value = 5.5 × 10(-31)) and is higher in men (p value = 0.002) but lower in participants of African ancestry (p value = 0.003). In a subset of 47 individuals from whom serial samples were collected up to 6 years apart, complex changes were noted over time and showed an overall increase in the proportion of mosaic cells as age increased. Our large combined sample allowed for a unique ability to characterize detectable genetic mosaicism involving large structural events and strengthens the emerging evidence of non-random erosion of the genome in the aging population.

Multiple Break-Points Detection in Array CGH Data via the Cross-Entropy Method.

Array comparative genome hybridization (aCGH) is a widely used methodology to detect copy number variations of a genome in high resolution. Knowing the number of break-points and their corresponding locations in genomic sequences serves different biological needs. Primarily, it helps to identify disease-causing genes that have functional importance in characterizing genome wide diseases. For human autosomes the normal copy number is two, whereas at the sites of oncogenes it increases (gain of DNA) and at the tumour suppressor genes it decreases (loss of DNA). The majority of the current detection methods are deterministic in their set-up and use dynamic programming or different smoothing techniques to obtain the estimates of copy number variations. These approaches limit the search space of the problem due to different assumptions considered in the methods and do not represent the true nature of the uncertainty associated with the unknown break-points in genomic sequences. We propose the Cross-Entropy method, which is a model-based stochastic optimization technique as an exact search method, to estimate both the number and locations of the break-points in aCGH data. We model the continuous scale log-ratio data obtained by the aCGH technique as a multiple break-point problem. The proposed methodology is compared with well established publicly available methods using both artificially generated data and real data. Results show that the proposed procedure is an effective way of estimating number and especially the locations of break-points with high level of precision. Availability: The methods described in this article are implemented in the new R package breakpoint and it is available from the Comprehensive R Archive Network at http://CRAN.R-project.org/package=breakpoint.

An estimate of the average number of recessive lethal mutations carried by humans.

The effects of inbreeding on human health depend critically on the number and severity of recessive, deleterious mutations carried by individuals. In humans, existing estimates of these quantities are based on comparisons between consanguineous and nonconsanguineous couples, an approach that confounds socioeconomic and genetic effects of inbreeding. To overcome this limitation, we focused on a founder population that practices a communal lifestyle, for which there is almost complete Mendelian disease ascertainment and a known pedigree. Focusing on recessive lethal diseases and simulating allele transmissions, we estimated that each haploid set of human autosomes carries on average 0.29 (95% credible interval [0.10, 0.84]) recessive alleles that lead to complete sterility or death by reproductive age when homozygous. Comparison to existing estimates in humans suggests that a substantial fraction of the total burden imposed by recessive deleterious variants is due to single mutations that lead to sterility or death between birth and reproductive age. In turn, comparison to estimates from other eukaryotes points to a surprising constancy of the average number of recessive lethal mutations across organisms with markedly different genome sizes.

Mongolians in the Genetic Landscape of Central Asia: Exploring the Genetic Relations among Mongolians and Other World Populations.

Genetic data on North and Central Asian populations are underrepresented in the literature, especially for autosomal markers. In the present study we used 812 single nucleotide polymorphisms (SNPs) distributed across all the human autosomes and extensively studied at Yale to examine the affinities of two recently collected samples of populations: rural and cosmopolitan Mongolians from Ulaanbaatar and nomadic, Turkic-speaking Tsaatan from Mongolia near the Siberian border. We compare these two populations with each other and with a global set of populations and discuss their relationships to New World populations. Specifically, we analyze data on 521 autosomal loci (single SNPs and multi-SNP haplotypes) studied in 57 populations representing all the major geographical regions of the world. We conclude that these North and Central Asian populations are genetically distinct from all other populations in our study and may be close to the ancestral lineage leading to the New World populations.