p53-dependent change in replication timing of the human genome

The tumor suppressor gene p53 has roles in multiple cell-cycle checkpoints, including the G1/S transition, to prevent replication of cells with DNA damage. p53 is thought to be associated with regulation of replication timing during S-phase in the human genome. In the present study, we used p53-wild-type and p53-null HCT116 colon carcinoma cells to analyze p53-dependent changes in replication timing of the human genome. The percentage of HCT116 p53(−/−) cells in S-phase was higher than that of HCT116 p53(+/+) cells. We compared replication timing of human genes between the two cell lines using 25,000 human cDNA microarray. We identified genes that replicated earlier in HCT116 p53(−/−) cells than in HCT116 p53(+/+) cells. These genes included cell-cycle- and apoptosis-related genes. We propose that p53 plays a role in regulation of replication timing of the human genome through the control of cell-cycle checkpoints.     

On the origin of genomic adaptation at high temperature for prokaryotic organisms

For a long time, the central issue of evolutionary genomics was to find out the adaptive strategy of nucleic acid molecules of various microorganisms having different optimal growth temperatures (Topt). Long-standing controversies exist regarding the correlations between genomic G+C content and Topt, and this debate has not been yet settled. We address this problem by considering the fact that adaptation to growth at high temperature requires a coordinated set of evolutionary changes affecting: (i) nucleic acid thermostability and (ii) stability of codon-anticodon interactions. In the present study, we analyzed 16 prokaryotic genomes having intermediate G+C content and widely varying optimal growth temperatures. Results show that elevated growth temperature imposes selective constraints not only on nucleic acid level but also affects the stability of codon-anticodon interaction. We observed a decrease in the frequency of SSC and SSG codons with the increase in Topt to avoid the formation of side-by-side GC base pairs in the codon-anticodon interaction, thereby making it impossible for a genome to increase GC composition uniformly through the whole coding sequence. Thus, we suggest that any attempt to obtain a generalized relation between genomic GC composition and optimal growth temperature would hardly evolve any satisfactory result.     

Correlations between genomic GC levels and optimal growth temperatures: some comments

Regarding the existence of any specific correlation between optimal growth temperature and genomic GC levels, Musto et al. [FEBS Lett. 573 (2004) 73] have recently performed analysis on 20 prokaryotic families and showed that in most of the families there exists a positive correlation between these two parameters. On the basis of these results they claimed that optimal growth temperature is one of the factors that influence genomic GC composition in prokaryotes. In a subsequent article, Marashi and Ghalanbor [Biochem. Biophys. Res. Commun. 325 (2004) 381] have demonstrated that the correlation values change substantially when very few points in some of the families were excluded from the data set of Musto et al. [FEBS Lett. 573 (2004) 73]. But Marashi and Ghalanbor have not provided any reason behind this. The points excluded by Marashi and Ghalanbor are actually the outliers in the data set, which strongly affect the correlation coefficients. But the presence of outliers in large data set hardly had any effect on the correlation values. Marashi and Ghalanbor have excluded points from only those families that have small sample sizes and observed a substantial change in correlation coefficient values. Therefore, we argue that any conclusion drawn for a small sample size having outliers is always questionable. Although Musto's approach is a novel one, but to make any generalization one needs to be careful about the flawlessness in the data set.     

Biological Implications of Isochore Boundaries in the Human Genome

Abstract

The human genome is composed of large sequence segments with fairly homogeneous GC content, namely isochores, which have been linked to many important functions; biological implications of most isochore boundaries, however, remain elusive, partly due to the difficulty in determining these boundaries at high resolution. Using the segmentation algorithm based on the quadratic divergence, we re-determined all 79 boundaries of previously identified human isochores at single-nucleotide resolution, and then compared the boundary coordinates with other genome features. We found that 55.7% of isochore boundaries coincide with termini of repeat elements; 45.6% of isochore boundaries coincide with termini of highly conserved sequences based on alignment of 17 vertebrate genomes, i.e., the highly conserved genome sequence switches to a less or non-conserved one at the isochore boundary; some isochore boundaries coincide with abrupt change of CpG island distribution (note that one boundary can associate with more than one genome feature). In addition, sequences around isochore boundaries are highly conserved. It seems reasonable to deduce that the boundaries of all the isochores studied here would be replication timing sites in the human genome. These results suggest possible key roles of the isochore boundaries and may further our understanding of the human genome organization.     

人类基因组SNPs的研究现状及应用前景

基因组DNA是生物体各种生理、病理性状的物质基础,人类DNA序列变异约90%表现为单核苷酸多态性(singlenucleotidepolymorphisms,SNPs),这是一种常见的遗传变异类型,在人类基因组中广泛存在,被认为是人类疾病易感性和药物反应的决定性因素。本文主要介绍了SNPs的分类及特点、人类基因组SNPs的研究现状、SNPs在实践中的应用,以及SNPs在遗传作图、医药、遗传易感性、个体化医疗等方面的研究前景,并探讨了当前SNPs研究中存在的问题。     

Identification of isochore boundaries in the human genome using the technique of wavelet multiresolution analysis

Incorporated with the Z curve method, the technique of wavelet multiresolution (also known as multiscale) analysis has been proposed to identify the boundaries of isochores in the human genome. The human MHC sequence and the longest contigs of human chromosomes 21 and 22 are used as examples. The boundary between the isochores of Class III and Class II in the MHC sequence has been detected and found to be situated at the position 2,490,368bp. This result is in good agreement with the experimental evidence. An isochore with a length of about 7Mb in chromosome 21 has been identified and found to be gene- and Alu-poor. We have also found that the G+C content of chromosome 21 is more homogeneous than that of chromosome 22. Compared with the window-based methods, the present method has the highest resolution for identifying the boundaries of isochores, even at a scale of single base. Compared with the entropic segmentation method, the present method has the merits of more intuitiveness and less calculations. The important conclusion drawn in this study is that the segmentation points, at which the G+C content undergoes relatively dramatic changes, do exist in the human genome. These 'singularity' points may be considered to be candidates of isochore boundaries in the human genome. The method presented is a general one and can be used to analyze any other genomes.     

The effects of guanine and cytosine variation on dinucleotide frequency and amino acid composition in the human genome

Summary One hundred twelve human DNA sequences were analyzed with respect to dinucleotide frequency and amino acid composition. The variation in guanine and cytosine (G+C) content revealed: (1) at 2–3 and 3-1 doublet positions CG discrimination is attenuated at high G+C, but TA disfavor is enhanced, and (2) several amino acids are subject to G+C change. These findings have been reported in part for collections of sequences from various species. The present study confirms that in a single organism-the human-the G+C effects do exist. Aspects of the argument that connects G+C with protein thermal stability are also discussed.     

Prediction of whole-genome DNA G + C content within the genus Aeromonas based on housekeeping gene sequences

Different methods are available to determine the G + C content (e.g. thermal denaturation temperature or high performance liquid chromatography, HPLC), but obtained values may differ significantly between strains, as well as between laboratories. Recently, several authors have demonstrated that the genomic DNA G + C content of prokaryotes can be reliably estimated from one or several protein coding gene nucleotide sequences. Few G + C content values have been published for the Aeromonas species described and the data, when available, are often incomplete or provide only a range of values. Our aim in this current work was twofold. First, the genomic G + C content of the type or reference strains of all species and subspecies of the genus Aeromonas was determined with a traditional experimental method in the same laboratory. Second, we wanted to see if the sequence-based method to estimate the G + C content described by Fournier et al. [7] could be applied to determine the G + C content of the different species of Aeromonas from the sequences of the genes used in taxonomy or phylogeny for this genus.     

Validated hydrophilic interaction LC-MS/MS method for determination of N-methyl-2-pyrrolidinone residue: Applied to a depletion study of N-methyl-2-pyrrolidinone in swine liver following intramuscular administration of drug N-methyl-2-pyrrolidinone formulation

A hydrophilic interaction high-performance liquid chromatography coupled with tandem mass spectrometry method for determination of N-methyl-2-pyrrolidinone in swine liver was developed and validated. After the fortification of N-methyl-d(3)-2-pyrrolidinone-d(6) as the deuterium-labeled internal standard, N-methyl-2-pyrrolidinone in swine liver was extracted by acetonitrile and the supernatant was led through a C18+WAX mixed-mode SPE cartridge for removal of the matrix interferences. The final eluate was acidified by formic acid and then injected onto a 3μm 15cm×2.1mm TX column for hydrophilic interaction chromatographic analysis. Mass spectrometry detection was carried on a PE Sciex API 4000 triple quadrupole mass spectrometer using positive turbo-ion spray ionization mode. The MRM transitions were 100→58 for N-methyl-2-pyrrolidinone and 109→62 for N-methyl-d(3)-2-pyrrolidinone-d(6). Solvent calibration standards could be readily used for quantitative analysis of N-methyl-2-pyrrolidinone with excellent precision and accuracy, although there are endogenous levels of N-methyl-2-pyrrolidinone in many blank matrices. The true recovery was nearly 100% and the MRM signal of N-methyl-2-pyrrolidinone was suppressed about 30% because of the matrix effect. Nevertheless, N-methyl-d(3)-2-pyrrolidinone-d(6) completely compensated the ion-suppression effect and the injection-to-injection variation. The detection limit was 5ngg(-1) swine liver. The validated method was applied to a depletion study of N-methyl-2-pyrrolidinone in swine liver following intramuscular administration of a drug N-methyl-2-pyrrolidinone formulation.     

Ribosomal DNA replication time coordinates completion of genome replication and anaphase in yeast

1. Download : Download high-res image (134KB)
2. Download : Download full-size image

     

A survey of single nucleotide polymorphisms identified from whole‐genome sequencing and their functional effect in the porcine genome,

Genetic variants detected from sequence have been used to successfully identify causal variants and map complex traits in several organisms. High and moderate impact variants, those expected to alter or disrupt the protein coded by a gene and those that regulate protein production, likely have a more significant effect on phenotypic variation than do other types of genetic variants. Hence, a comprehensive list of these functional variants would be of considerable interest in swine genomic studies, particularly those targeting fertility and production traits. Whole‐genome sequence was obtained from 72 of the founders of an intensely phenotyped experimental swine herd at the U.S. Meat Animal Research Center (USMARC). These animals included all 24 of the founding boars (12 Duroc and 12 Landrace) and 48 Yorkshire–Landrace composite sows. Sequence reads were mapped to the Sscrofa10.2 genome build, resulting in a mean of 6.1 fold (×) coverage per genome. A total of 22 342 915 high confidence SNPs were identified from the sequenced genomes. These included 21 million previously reported SNPs and 79% of the 62 163 SNPs on the PorcineSNP60 BeadChip assay. Variation was detected in the coding sequence or untranslated regions (UTRs) of 87.8% of the genes in the porcine genome: loss‐of‐function variants were predicted in 504 genes, 10 202 genes contained nonsynonymous variants, 10 773 had variation in UTRs and 13 010 genes contained synonymous variants. Approximately 139 000 SNPs were classified as loss‐of‐function, nonsynonymous or regulatory, which suggests that over 99% of the variation detected in our pigs could potentially be ignored, allowing us to focus on a much smaller number of functional SNPs during future analyses.     

Scanning the genome for gene single nucleotide polymorphisms involved in adaptive population differentiation in white spruce

Conifers are characterized by a large genome size and a rapid decay of linkage disequilibrium, most often within gene limits. Genome scans based on noncoding markers are less likely to detect molecular adaptation linked to genes in these species. In this study, we assessed the effectiveness of a genome-wide single nucleotide polymorphism (SNP) scan focused on expressed genes in detecting local adaptation in a conifer species. Samples were collected from six natural populations of white spruce ( Picea glauca ) moderately differentiated for several quantitative characters. A total of 534 SNPs representing 345 expressed genes were analysed. Genes potentially under natural selection were identified by estimating the differentiation in SNP frequencies among populations ( F _ST) and identifying outliers, and by estimating local differentiation using a Bayesian approach. Both average expected heterozygosity and population differentiation estimates ( H _E = 0.270 and F _ST = 0.006) were comparable to those obtained with other genetic markers. Of all genes, 5.5% were identified as outliers with F _ST at the 95% confidence level, while 14% were identified as candidates for local adaptation with the Bayesian method. There was some overlap between the two gene sets. More than half of the candidate genes for local adaptation were specific to the warmest population, about 20% to the most arid population, and 15% to the coldest and most humid higher altitude population. These adaptive trends were consistent with the genes' putative functions and the divergence in quantitative traits noted among the populations. The results suggest that an approach separating the locus and population effects is useful to identify genes potentially under selection. These candidates are worth exploring in more details at the physiological and ecological levels.     

Sequence context analysis of 8.2 million single nucleotide polymorphisms in the human genome

We analyzed n-mers (n=3-8) in the local environment of 8,249,446 human SNPs and compared their distribution with that in the genome reference sequences. The results revealed that the short sequences, which contained at least one CpG dinucleotide, occurred more frequently in the local SNP sequences than in the genome sequences. To exclude the hypermutability effect of the methylated CpG dinucleotides on the sequence context of SNPs, we examined the distribution patterns for each of the six categories of substitution. We observed the similar pattern (i.e., CpG-containing n-mers vs. non-CpG-containing n-mers) in SNP categories A/G, C/T and C/G but the opposite pattern in category A/T. We next identified 34,928 putative CpG islands in the human genome and located 133,591 SNPs within these islands. In the CpG islands, CpG SNPs were 3.92-fold less prevalent relative to the presence of CpG dinucleotides. Conversely, in the human genome, the frequency of CpG dinucleotides at the polymorphic sites was 6.09 times that in the genome reference sequences. These results support the previous views of mutational suppression at the CpG sites in the CpG islands and hypermutability of the methylated CpG dinucleotides that are prevalent in the non-CpG island sequences in the human genome. Our study represents a comprehensive investigation of the sequence context of SNPs in the human genome and in human CpG islands.     

Gap-filling and bypass at the replication fork are both active mechanisms for tolerance of low-dose ultraviolet-induced DNA damage in the human genome

Ultraviolet (UV)-induced DNA damage are removed by nucleotide excision repair (NER) or can be tolerated by specialized translesion synthesis (TLS) polymerases, such as Polη. TLS may act at stalled replication forks or through an S-phase independent gap-filling mechanism. After UVC irradiation, Polη-deficient (XP-V) human cells were arrested in early S-phase and exhibited both single-strand DNA (ssDNA) and prolonged replication fork stalling, as detected by DNA fiber assay. In contrast, NER deficiency in XP-C cells caused no apparent defect in S-phase progression despite the accumulation of ssDNA and a G2-phase arrest. These data indicate that while Polη is essential for DNA synthesis at ongoing damaged replication forks, NER deficiency might unmask the involvement of tolerance pathway through a gap-filling mechanism. ATR knock down by siRNA or caffeine addition provoked increased cell death in both XP-V and XP-C cells exposed to low-dose of UVC, underscoring the involvement of ATR/Chk1 pathway in both DNA damage tolerance mechanisms. We generated a unique human cell line deficient in XPC and Polη proteins, which exhibited both S- and G2-phase arrest after UVC irradiation, consistent with both single deficiencies. In these XP-C/Polη^KD cells, UVC-induced replicative intermediates may collapse into double-strand breaks, leading to cell death. In conclusion, both TLS at stalled replication forks and gap-filling are active mechanisms for the tolerance of UVC-induced DNA damage in human cells and the preference for one or another pathway depends on the cellular genotype.     

Detection of single nucleotide variations in expressed exons of the human genome using RNA-Seq

Whole-genome resequencing is still a costly method to detect genetic mutations that lead to altered forms of proteins and may be associated with disease development. Since the majority of disease-related single nucleotide variations (SNVs) are found in protein-coding regions, we propose to identify SNVs in expressed exons of the human genome using the recently developed RNA-Seq technique. We identify 12 176 and 10 621 SNVs, respectively, in Jurkat T cells and CD4⁺ T cells from a healthy donor. Interestingly, our data show that one copy of the TAL-1 proto-oncogene has a point mutation in 3′ UTR and only the mutant allele is expressed in Jurkat cells. We provide a comprehensive dataset for further understanding the cancer biology of Jurkat cells. Our results indicate that this is a cost-effective and efficient strategy to systematically identify SNVs in the expressed regions of the human genome.