Systematic biases in DNA copy number originate from isolation procedures

Background

The ability to accurately detect DNA copy number variation in both a sensitive and quantitative manner is important in many research areas. However, genome-wide DNA copy number analyses are complicated by variations in detection signal.

Results

While GC content has been used to correct for this, here we show that coverage biases are tissue-specific and independent of the detection method as demonstrated by next-generation sequencing and array CGH. Moreover, we show that DNA isolation stringency affects the degree of equimolar coverage and that the observed biases coincide with chromatin characteristics like gene expression, genomic isochores, and replication timing.

Conclusion

These results indicate that chromatin organization is a main determinant for differential DNA retrieval. These findings are highly relevant for germline and somatic DNA copy number variation analyses.     

Codon usage of human hepatitis C virus clearance genes in relation to its expression

Hepatitis C virus (HCV) infection is among the leading causes of hepatocellular carcinoma and liver cirrhosis globally, with a high economic burden. The disease progression is well established, but less is known about the spontaneous HCV infection clearance. This study tries to establish the relationship between codon biasness and expression of HCV clearance candidate genes in normal and HCV infected liver tissues. A total of 112 coding sequences comprising 151 679 codons were subjected to the computation of codon indices, namely relative synonymous codon usage, an effective number of codon (Nc), frequency of optimal codon, codon adaptation index, codon bias index, and base compositions. Codon indices report of GC3s, GC12, hydropathicity, and aromaticity implicates both mutational and translational selection in the candidate gene set. This was further correlated with the differentially expressed genes among the selected genes using BioGPS. A significant correlation is observed between the gene expression of normal liver and cancerous liver tissues with codon bias (Nc). Gene expression is also correlated with relative codon bias values, indicating that CCL5, APOA2, CD28, IFITM1, and TNFSF4 genes have higher expression. These results are quite encouraging in selecting the high responsive genes in HCV clearance. However, there could be additional genes which could also orchestrate the clearance role with the above mentioned first line of defensive genes.     

Approximation of bias and mean-squared error in two-sample Mendelian randomization analyses

Mendelian randomization (MR) is a type of instrumental variable (IV) analysis that uses genetic variants as IVs for a risk factor to study its causal effect on an outcome. Extensive investigations on the performance of IV analysis procedures, such as the one based on the two-stage least squares (2SLS) procedure, have been conducted under the one-sample scenario, where measures on IVs, the risk factor, and the outcome are assumed to be available for each study participant. Recent MR analysis usually is performed with data from two independent or partially overlapping genetic association studies (two-sample setting), with one providing information on the association between the IVs and the outcome, and the other on the association between the IVs and the risk factor. We investigate the performance of 2SLS in the two-sample–based MR when the IVs are weakly associated with the risk factor. We derive closed form formulas for the bias and mean squared error of the 2SLS estimate and verify them with numeric simulations under realistic circumstances. Using these analytic formulas, we can study the pros and cons of conducting MR analysis under one-sample and two-sample settings and assess the impact of having overlapping samples. We also propose and validate a bias-corrected estimator for the causal effect.     

Testing small study effects in multivariate meta-analysis

Small study effects occur when smaller studies show different, often larger, treatment effects than large ones, which may threaten the validity of systematic reviews and meta-analyses. The most well-known reasons for small study effects include publication bias, outcome reporting bias, and clinical heterogeneity. Methods to account for small study effects in univariate meta-analysis have been extensively studied. However, detecting small study effects in a multivariate meta-analysis setting remains an untouched research area. One of the complications is that different types of selection processes can be involved in the reporting of multivariate outcomes. For example, some studies may be completely unpublished while others may selectively report multiple outcomes. In this paper, we propose a score test as an overall test of small study effects in multivariate meta-analysis. Two detailed case studies are given to demonstrate the advantage of the proposed test over various naive applications of univariate tests in practice. Through simulation studies, the proposed test is found to retain nominal Type I error rates with considerable power in moderate sample size settings. Finally, we also evaluate the concordance between the proposed tests with the naive application of univariate tests by evaluating 44 systematic reviews with multiple outcomes from the Cochrane Database.     

Female-specific resource limitation does not make the opportunity for selection more female biased

Competition for limiting resources and stress can magnify variance in fitness and therefore selection. But even in a common environment, the strength of selection can differ across the sexes, as their fitness is often limited by different factors. Indeed, most taxa show stronger selection in males, a bias often ascribed to intense competition for access to mating partners. This sex bias could reverberate on many aspects of evolution, from speed of adaptation to genome evolution. It is unclear, however, whether stronger opportunity for selection in males is a pattern robust to sex-specific stress or resource limitation. We test this in the model species Callosobruchus maculatus by comparing female and male opportunity for selection (i) with and without limitation of quality oviposition sites, and (ii) under delayed age at oviposition. Decreasing the abundance of the resource key to females or increasing their reproductive age was challenging, as shown by a reduction in mean fitness, but opportunity for selection remained stronger in males across all treatments, and even more so when oviposition sites were limiting. This suggests that males remain the more variable sex independent of context, and that the opportunity for selection through males is indirectly affected by female-specific resource limitation.     

Compensatory nature of Chargaff’s second parity rule

The second parity rule of Chargaff (A≈T and G≈C within one strand) holds all over the living world with minor exceptions. It is maintained with higher accuracy for long sequences. The question addressed in the article is how different sequence types, with different biases from the parity, contribute to the general effect. It appears that the sequence segments with biases of opposite sign are intermingled, so that with sufficient sequence lengths the parity is established. The parity rule seems to be a cumulative result of a number of independent processes in the genome evolution, with the parity as their intrinsic property. Symmetrical appearance of simple repeats and of Alu sequences in the human DNA strands, and other contributions to the Chargaff parity II rule are discussed.     

Molecular characterization influencing metal resistance in the Cupriavidus/Ralstonia genomes

Our environment is stressed with a load of heavy and toxic metals. Microbes, abundant in our environment, are found to adapt well to this metal-stressed condition. A comparative study among five Cupriavidus/Ralstonia genomes can offer a better perception of their evolutionary mechanisms to adapt to these conditions. We have studied codon usage among 1051 genes common to all these organisms and identified 15 optimal codons frequently used in highly expressed genes present within 1051 genes. We found the core genes of Cupriavidus metallidurans CH34 have a different optimal codon choice for arginine, glycine and alanine in comparison with the other four bacteria. We also found that the synonymous codon usage bias within these 1051 core genes is highly correlated with their gene expression. This supports that translational selection drives synonymous codon usage in the core genes of these genomes. Synonymous codon usage is highly conserved in the core genes of these five genomes. The only exception among them is C. metallidurans CH34. This genomewide shift in synonymous codon choice in C. metallidurans CH34 may have taken place due to the insertion of new genes in its genomes facilitating them to survive in heavy metal containing environment and the co-evolution of the other genes in its genome to achieve a balance in gene expression. Structural studies indicated the presence of a longer N-terminal region containing a copper-binding domain in the cupC proteins of C. metallidurans CH3 that helps it to attain higher binding efficacy with copper in comparison with its orthologs.     

Yeast Pif1 Helicase Exhibits a One-base-pair Stepping Mechanism for Unwinding Duplex DNA

Kinetic analysis of the DNA unwinding and translocation activities of helicases is necessary for characterization of the biochemical mechanism(s) for this class of enzymes. Saccharomyces cerevisiae Pif1 helicase was characterized using presteady state kinetics to determine rates of DNA unwinding, displacement of streptavidin from biotinylated DNA, translocation on single-stranded DNA (ssDNA), and ATP hydrolysis activities. Unwinding of substrates containing varying duplex lengths was fit globally to a model for stepwise unwinding and resulted in an unwinding rate of ∼75 bp/s and a kinetic step size of 1 base pair. Pif1 is capable of displacing streptavidin from biotinylated oligonucleotides with a linear increase in the rates as the length of the oligonucleotides increased. The rate of translocation on ssDNA was determined by measuring dissociation from varying lengths of ssDNA and is essentially the same as the rate of unwinding of dsDNA, making Pif1 an active helicase. The ATPase activity of Pif1 on ssDNA was determined using fluorescently labeled phosphate-binding protein to measure the rate of phosphate release. The quantity of phosphate released corresponds to a chemical efficiency of 0.84 ATP/nucleotides translocated. Hence, when all of the kinetic data are considered, Pif1 appears to move along DNA in single nucleotide or base pair steps, powered by hydrolysis of 1 molecule of ATP.