Accuracy in Copy Number Calling by qPCR and PRT : A Matter of DNA

7 p.

8p23 beta-defensin copy number determination by single-locus pseudogene-based paralog ratio tests risk bias due to low-frequency sequence variations.

BACKGROUND The copy number variation (CNV) in beta-defensin genes (DEFB) on human chromosome 8p23 has been proposed to contribute to the phenotypic differences in inflammatory diseases. However, determination of exact DEFB CN is a major challenge in association studies. Quantitative real-time PCR (qPCR), paralog ratio tests (PRT) and multiplex ligation-dependent probe amplification (MLPA) have been extensively used to determine DEFB CN in different laboratories, but inter-method inconsistencies were observed frequently. In this study we asked which one is superior among the three methods for DEFB CN determination. RESULTS We developed a clustering approach for MLPA and PRT to statistically correlate data from a single experiment. Then we compared qPCR, a newly designed PRT and MLPA for DEFB CN determination in 285 DNA samples. We found MLPA had the best convergence and clustering results of the raw data and the highest call rate. In addition, the concordance rates between MLPA or PRT and qPCR (32.12% and 37.99%, respectively) were unacceptably low with underestimated CN by qPCR. Concordance rate between MLPA and PRT (90.52%) was high but PRT systematically underestimated CN by one in a subset of samples. In these samples a sequence variant which caused complete PCR dropout of the respective DEFB cluster copies was found in one primer binding site of one of the targeted paralogous pseudogenes. CONCLUSION MLPA is superior to PRT and even more to qPCR for DEFB CN determination. Although the applied PRT provides in most cases reliable results, such a test is particularly sensitive to low-frequency sequence variations preferably accumulating in loci like pseudogenes which are most likely not under selective pressure. In the light of the superior performance of multiplex assays, the drawbacks of such single PRTs could be overcome by combining more test markers.

High-resolution SNP microarray investigation of copy number variations on chromosome 18 in a control cohort

Copy number variations (CNVs) as described in the healthy population are purported to contribute significantly to genetic heterogeneity. Recent studies have described CNVs using lymphoblastoid cell lines or by application of specifically developed algorithms to interrogate previously described data. However, the full extent of CNVs remains unclear. Using high-density SNP array, we have undertaken a comprehensive investigation of chromosome 18 for CNV discovery and characterisation of distribution and association with chromosome architecture. We identified 399 CNVs, of which loss represents 98%, 58% are less than 2.5 kb in size and 71% are intergenic. Intronic deletions account for the majority of copy number changes with gene involvement. Furthermore, one-third of CNVs do not have putative breakpoints within repetitive sequences. We conclude that replicative processes, mediated either by repetitive elements or microhomology, account for the majority of CNVs in the healthy population. Genomic instability involving the formation of a non-B structure is demonstrated in one region.

Absolute quantitation of Marek's disease virus and Herpesvirus of turkeys in chicken lymphocyte, feather tip and dust samples using real-time PCR

The further development of Taqman quantitative real-time PCR (qPCR) assays for the absolute quantitation of Marek's disease virus serotype 1 (MDV1) and Herpesvirus of turkeys (HVT) viruses is described and the sensitivity and reproducibility of each assay reported. Using plasmid DNA copies, the lower limit of detection was determined to be 5 copies for the MDV1 assay and 75 copies for the HVT assay. Both assays were found to be highly reproducible for Ct values and calculated copy numbers with mean intra- and inter-assay coefficients of variation being less than 5% for Ct and 20% for calculated copy number. The genome copy number of MDV1 and HVT viruses was quantified in PBL and feather tips from experimentally infected chickens, and field poultry dust samples. Parallelism was demonstrated between the plasmid-based standard curves, and standard curves derived from infected spleen material containing both viral and host DNA, allowing the latter to be used for absolute quantification. These methods should prove useful for the reliable differentiation and absolute quantitation of MDV1 and HVT viruses in a wide range of samples.

Genes expressed in dental enamel development are associated with molar-incisor hypomineralization

Genetic disturbances during dental development influence variation of number and shape of the dentition. In this study, we tested if genetic variation in enamel formation genes is associated with molar-incisor hypomineralization (MIH), also taking into consideration caries experience. DNA samples from 163 cases with MIH and 82 unaffected controls from Turkey, and 71 cases with MIH and 89 unaffected controls from Brazil were studied. Eleven markers in five genes [ameloblastin (AMBN), amelogenin (AMELX), enamelin (ENAM), tuftelin (TUFT1), and tuftelin-interacting protein 11 (TFIP11)] were genotyped by the TaqMan method. Chi-square was used to compare allele and genotype frequencies between cases with MIH and controls. In the Brazilian data, distinct caries experience within the MIH group was also tested for association with genetic variation in enamel formation genes. The ENAM rs3796704 marker was associated with MIH in both populations (Brazil: p = 0.03; OR = 0.28; 95% C.I. = 0.06-1.0; Turkey: p = 1.22e-012; OR = 17.36; 95% C.I. = 5.98-56.78). Associations between TFIP11 (p = 0.02), ENAM (p = 0.00001), and AMELX (p = 0.01) could be seen with caries independent of having MIH or genomic DNA copies of Streptococcus mutans detected by real time PCR in the Brazilian sample. Several genes involved in enamel formation appear to contribute to MIH. © 2013.

Identificação de regiões com variação no número de cópias de segmentos de DNA em bovinos de raças autóctones espanholas

Pós-graduação em Genética e Melhoramento Animal - FCAV

A MULTILEVEL MODEL TO ADDRESS BATCH EFFECTS IN COPY NUMBER ESTIMATION USING SNP ARRAYS

Submicroscopic changes in chromosomal DNA copy number dosage are common and have been implicated in many heritable diseases and cancers. Recent high-throughput technologies have a resolution that permits the detection of segmental changes in DNA copy number that span thousands of basepairs across the genome. Genome-wide association studies (GWAS) may simultaneously screen for copy number-phenotype and SNP-phenotype associations as part of the analytic strategy. However, genome-wide array analyses are particularly susceptible to batch effects as the logistics of preparing DNA and processing thousands of arrays often involves multiple laboratories and technicians, or changes over calendar time to the reagents and laboratory equipment. Failure to adjust for batch effects can lead to incorrect inference and requires inefficient post-hoc quality control procedures that exclude regions that are associated with batch. Our work extends previous model-based approaches for copy number estimation by explicitly modeling batch effects and using shrinkage to improve locus-specific estimates of copy number uncertainty. Key features of this approach include the use of diallelic genotype calls from experimental data to estimate batch- and locus-specific parameters of background and signal without the requirement of training data. We illustrate these ideas using a study of bipolar disease and a study of chromosome 21 trisomy. The former has batch effects that dominate much of the observed variation in quantile-normalized intensities, while the latter illustrates the robustness of our approach to datasets where as many as 25% of the samples have altered copy number. Locus-specific estimates of copy number can be plotted on the copy-number scale to investigate mosaicism and guide the choice of appropriate downstream approaches for smoothing the copy number as a function of physical position. The software is open source and implemented in the R package CRLMM available at Bioconductor (http:www.bioconductor.org).

Transgenic DNA integrated into the oat genome is frequently interspersed by host DNA

Integration of transgenic DNA into the plant genome was investigated in 13 transgenic oat (Avena sativa L.) lines produced using microprojectile bombardment with one or two cotransformed plasmids. In all transformation events, the transgenic DNA integrated into the plant genome consisted of intact transgene copies that were accompanied by multiple, rearranged, and/or truncated transgene fragments. All fragments of transgenic DNA cosegregated, indicating that they were integrated at single gene loci. Analysis of the structure of the transgenic loci indicated that the transgenic DNA was interspersed by the host genomic DNA. The number of insertions of transgenic DNA within the transgene loci varied from 2 to 12 among the 13 lines. Restriction endonucleases that do not cleave the introduced plasmids produced restriction fragments ranging from 3.6 to about 60 kb in length hybridizing to a probe comprising the introduced plasmids. Although the size of the interspersing host DNA within the transgene locus is unknown, the sizes of the transgene-hybridizing restriction fragments indicated that the entire transgene locus must be at least from 35–280 kb. The observation that all transgenic lines analyzed exhibited genomic interspersion of multiple clustered transgenes suggests a predominating integration mechanism. We propose that transgene integration at multiple clustered DNA replication forks could account for the observed interspersion of transgenic DNA with host genomic DNA within transgenic loci.

Chloroplast DNA variation and the recent radiation of the giant senecios (Asteraceae) on the tall mountains of eastern Africa.

Chloroplast DNA restriction-site variation was surveyed among 40 accessions representing all 11 species of giant senecios (Dendrosenecio, Asteraceae) at all but one known location, plus three outgroup species. Remarkably little variation (only 9 variable sites out of roughly 1000 sites examined) was found among the 40 giant senecio accessions, yet as a group they differ significantly (at 18 sites) from Cineraria deltoidea, the closest known relative. This pattern indicates that the giant senecios underwent a recent dramatic radiation in eastern Africa and evolved from a relatively isolated lineage within the Senecioneae. Biogeographic interpretation of the molecular phylogeny suggests that the giant senecios originated high on Mt. Kilimanjaro, with subsequent dispersion to the Aberdares, Mt. Kenya, and the Cherangani Hills, followed by dispersion westward to the Ruwenzori Mountains, and then south to the Virunga Mountains, Mt. Kahuzi, and Mt. Muhi, but with dispersion back to Mt. Elgon. Geographic radiation was an important antecedent to the diversification in eastern Africa, which primarily involved repeated altitudinal radiation, both up and down the mountains, leading to morphological parallelism in both directions. In general, the plants on a given mountain are more closely related to each other than they are to plants on other mountains, and plants on nearby mountains are more closely related to each other than they are to plants on more distant mountains. The individual steps of the geographic radiation have occurred at various altitudes, some clearly the result of intermountain dispersal. The molecular evidence suggests that two species are extant ancestors to other species on the same or nearby mountains.

No mutations detected in the INSR gene in a chromosome 19p13 linked migraine pedigree

The aim of this study was to investigate through direct sequencing the insulin receptor (INSR) gene in DNA samples from a migraine affected family previously showing linkage to chromosome 19p13 in an attempt to detect disease associated mutations. Migraine is a common debilitating disorder with a significant genetic component. At present, the number and type of genes involved in the common forms of migraine are not clear. The INSR gene on chromosome 19p13.3-13.2 is a gene of interest since a number of single nucleotide polymorphisms (SNPs) located within the gene have been implicated in migraine with (MA) and without aura (MO). Six DNA samples obtained from non-founding migraine affected members of migraine family 1 (MF1) were used in this study. Genomic DNA was sequenced for the INSR gene in exons 1-22 and the promoter region. In the six migraine family member samples, previously reported SNPs were detected within two exonic DNA coding regions of the INSR gene. These SNPs, in exons 13 and 17, do not alter the normal INSR polypeptide sequence. In addition, intron 7 also revealed a DNA base sequence variation. For the 5' untranslated promoter region of the gene, no mutations or polymorphisms were detected. In conclusion, this study detected no INSR mutations in affected members of a chromosome 19 linked migraine pedigree. Hence, migraine linkage to this chromosomal region may involve other candidate genes.

Molecular biomarkers in non-small-cell lung cancer : a retrospective analysis of data from the phase 3 FLEX study

Background: Findings from the phase 3 FLEX study showed that the addition of cetuximab to cisplatin and vinorelbine significantly improved overall survival, compared with cisplatin and vinorelbine alone, in the first-line treatment of EGFR-expressing, advanced non-small-cell lung cancer (NSCLC). We investigated whether candidate biomarkers were predictive for the efficacy of chemotherapy plus cetuximab in this setting. Methods: Genomic DNA extracted from formalin-fixed paraffin-embedded (FFPE) tumour tissue of patients enrolled in the FLEX study was screened for KRAS codon 12 and 13 and EGFR kinase domain mutations with PCR-based assays. In FFPE tissue sections, EGFR copy number was assessed by dual-colour fluorescence in-situ hybridisation and PTEN expression by immunohistochemistry. Treatment outcome was investigated according to biomarker status in all available samples from patients in the intention-to-treat population. The primary endpoint in the FLEX study was overall survival. The FLEX study, which is ongoing but not recruiting participants, is registered with ClinicalTrials.gov, number NCT00148798. Findings: KRAS mutations were detected in 75 of 395 (19%) tumours and activating EGFR mutations in 64 of 436 (15%). EGFR copy number was scored as increased in 102 of 279 (37%) tumours and PTEN expression as negative in 107 of 303 (35%). Comparisons of treatment outcome between the two groups (chemotherapy plus cetuximab vs chemotherapy alone) according to biomarker status provided no indication that these biomarkers were of predictive value. Activating EGFR mutations were identified as indicators of good prognosis, with patients in both treatment groups whose tumours carried such mutations having improved survival compared with those whose tumours did not (chemotherapy plus cetuximab: median 17·5 months [95% CI 11·7-23·4] vs 8·5 months [7·1-10·8], hazard ratio [HR] 0·52 [0·32-0·84], p=0·0063; chemotherapy alone: 23·8 months [15·2-not reached] vs 10·0 months [8·7-11·0], HR 0·35 [0·21-0·59], p<0·0001). Expression of PTEN seemed to be a potential indicator of good prognosis, with patients whose tumours expressed PTEN having improved survival compared with those whose tumours did not, although this finding was not significant (chemotherapy plus cetuximab: median 11·4 months [8·6-13·6] vs 6·8 months [5·9-12·7], HR 0·80 [0·55-1·16], p=0·24; chemotherapy alone: 11·0 months [9·2-12·6] vs 9·3 months [7·6-11·9], HR 0·77 [0·54-1·10], p=0·16). Interpretation: The efficacy of chemotherapy plus cetuximab in the first-line treatment of advanced NSCLC seems to be independent of each of the biomarkers assessed. Funding: Merck KGaA. © 2011 Elsevier Ltd.

JK1 (FAM134B) gene and colorectal cancer : a pilot study on the gene copy number alterations and correlations with clinicopathological parameters

AIMS The aims of the study are to characterize changes in JK-1 (FAM134B) at the DNA level in colorectal adenocarcinoma and adenoma and exploring the possible correlations with clinical and pathological features. METHOD JK-1 gene DNA copy number changes were studied in 211 colorectal carcinomas, 32 colorectal adenoma and 20 colorectal non-cancer colorectal tissue samples by real-time quantitative polymerase chain reaction. The results were correlated with clinical and pathological parameters. RESULTS Colorectal adenomas were more likely to be amplified than deleted with regard to JK-1 (FAM134B) DNA copy number change. The copy number level of JK-1 (FAM134B) DNA in colorectal adenocarcinomas was significantly lower in comparison to colorectal adenomas. Changes in JK-1 (FAM134B) DNA copy number were associated with histological subtypes, and cancer stage. Lower copy numbers were associated with higher tumor stage, lymph node stage and overall pathological stage of cancer. Conversely, higher DNA copy numbers were detected more often in the mucinous adenocarcinoma. CONCLUSIONS This is the first study showing significant correlations of the JK-1 (FAM134B) gene copy number alterations with clinical and pathological features in a large cohort of pre-invasive and invasive colorectal malignancies. The changes in DNA copy number associated with progression of colorectal malignancies reflect that JK-1 (FAM134B) gene could play a role in controlling some steps in development of the invasive phenotypes.

Comparative analysis of genome maintenance genes in naked mole rat, mouse, and human

Genome maintenance (GM) is an essential defense system against aging and cancer, as both are characterized by increased genome instability. Here, we compared the copy number variation and mutation rate of 518 GM-associated genes in the naked mole rat (NMR), mouse, and human genomes. GM genes appeared to be strongly conserved, with copy number variation in only four genes. Interestingly, we found NMR to have a higher copy number of CEBPG, a regulator of DNA repair, and TINF2, a protector of telomere integrity. NMR, as well as human, was also found to have a lower rate of germline nucleotide substitution than the mouse. Together, the data suggest that the long-lived NMR, as well as human, has more robust GM than mouse and identifies new targets for the analysis of the exceptional longevity of the NMR.

Schizophrenia Genetic Variants are not Associated with Intelligence

Background Schizophrenia is associated with lower pre-morbid intelligence (IQ) in addition to (pre-morbid) cognitive decline. Both schizophrenia and IQ are highly heritable traits. Therefore, we hypothesized that genetic variants associated with schizophrenia, including copy number variants (CNVs) and a polygenic schizophrenia (risk) score (PSS), may influence intelligence. Method IQ was estimated with the Wechsler Adult Intelligence Scale (WAIS). CNVs were determined from single nucleotide polymorphism (SNP) data using the QuantiSNP and PennCNV algorithms. For the PSS, odds ratios for genome-wide SNP data were calculated in a sample collected by the Psychiatric Genome-Wide Association Study (GWAS) Consortium (8690 schizophrenia patients and 11 831 controls). These were used to calculate individual PSSs in our independent sample of 350 schizophrenia patients and 322 healthy controls. Results Although significantly more genes were disrupted by deletions in schizophrenia patients compared to controls (p = 0.009), there was no effect of CNV measures on IQ. The PSS was associated with disease status (R 2 = 0.055, p = 2.1 × 10 -7) and with IQ in the entire sample (R 2 = 0.018, p = 0.0008) but the effect on IQ disappeared after correction for disease status. Conclusions Our data suggest that rare and common schizophrenia-associated variants do not explain the variation in IQ in healthy subjects or in schizophrenia patients. Thus, reductions in IQ in schizophrenia patients may be secondary to other processes related to schizophrenia risk. © Cambridge University Press 2013.

I. Studies on bacteriophage DNA's and minicircular DNA's in M. lysodeikticus and E. coli 15. II. Flow dichroism of DNA solutions

Part I

Chapter 1.....A physicochemical study of the DNA molecules from the three bacteriophages, N1, N5, and N6, which infect the bacterium, M. lysodeikticus, has been made. The molecular weights, as measured by both electron microscopy and sedimentation velocity, are 23 x 10⁶ for N5 DNA and 31 x 10⁶ for N1 and N6 DNA's. All three DNA's are capable of thermally reversible cyclization. N1 and N6 DNA's have identical or very similar base sequences as judged by membrane filter hybridization and by electron microscope heteroduplex studies. They have identical or similar cohesive ends. These results are in accord with the close biological relation between N1 and N6 phages. N5 DNA is not closely related to N1 or N6 DNA. The denaturation T_m of all three DNA's is the same and corresponds to a (GC) content of 70%. However, the buoyant densities in CsCl of Nl and N6 DNA's are lower than expected, corresponding to predicted GC contents of 64 and 67%. The buoyant densities in Cs₂SO₄ are also somewhat anomalous. The buoyant density anomalies are probably due to the presence of odd bases. However, direct base composition analysis of N1 DNA by anion exchange chromatography confirms a GC content of 70%, and, in the elution system used, no peaks due to odd bases are present.

Chapter 2.....A covalently closed circular DNA form has been observed as an intracellular form during both productive and abortive infection processes in M. lysodeikticus. This species has been isolated by the method of CsC1-ethidium bromide centrifugation and examined with an electron microscope.

Chapter 3.....A minute circular DNA has been discovered as a homogeneous population in M. lysodeikticus. Its length and molecular weight as determined by electron microscopy are 0.445 μ and 0.88 x 10⁶ daltons respectively. There is about one minicircle per bacterium.

Chapter 4.....Several strains of E. coli 15 harbor a prophage. Viral growth can be induced by exposing the host to mitomycin C or to uv irradiation. The coliphage 15 particles from E. coli 15 and E, coli 15 T^- appear as normal phage with head and tail structure; the particles from E. coli 15 TAU are tailless. The complete particles exert a colicinogenic activity on E.coli 15 and 15 T^-, the tailless particles do not. No host for a productive viral infection has been found and the phage may be defective. The properties of the DNA of the virus have been studied, mainly by electron microscopy. After induction but before lysis, a closed circular DNA with a contour length of about 11.9 μ is found in the bacterium; the mature phage DNA is a linear duplex and 7.5% longer than the intracellular circular form. This suggests the hypothesis that the mature phage DNA is terminally repetitious and circularly permuted. The hypothesis was confirmed by observing that denaturation and renaturation of the mature phage DNA produce circular duplexes with two single-stranded branches corresponding to the terminal repetition. The contour length of the mature phage DNA was measured relative to φX RFII DNA and λ DNA; the calculated molecular weight is 27 x 10⁶. The length of the single-stranded terminal repetition was compared to the length of φX 174 DNA under conditions where single-stranded DNA is seen in an extended form in electron micrographs. The length of the terminal repetition is found to be 7.4% of the length of the nonrepetitious part of the coliphage 15 DNA. The number of base pairs in the terminal repetition is variable in different molecules, with a fractional standard deviation of 0.18 of the average number in the terminal repetition. A new phenomenon termed "branch migration" has been discovered in renatured circular molecules; it results in forked branches, with two emerging single strands, at the position of the terminal repetition. The distribution of branch separations between the two terminal repetitions in the population of renatured circular molecules was studied. The observed distribution suggests that there is an excluded volume effect in the renaturation of a population of circularly permuted molecules such that strands with close beginning points preferentially renature with each other. This selective renaturation and the phenomenon of branch migration both affect the distribution of branch separations; the observed distribution does not contradict the hypothesis of a random distribution of beginning points around the chromosome.

Chapter 5....Some physicochemical studies on the minicircular DNA species in E. coli 15 (0.670 μ, 1.47 x 10⁶ daltons) have been made. Electron microscopic observations showed multimeric forms of the minicircle which amount to 5% of total DNA species and also showed presumably replicating forms of the minicircle. A renaturation kinetic study showed that the minicircle is a unique DNA species in its size and base sequence. A study on the minicircle replication has been made under condition in which host DNA synthesis is synchronized. Despite experimental uncertainties involved, it seems that the minicircle replication is random and the number of the minicircles increases continuously throughout a generation of the host, regardless of host DNA synchronization.

Part II

The flow dichroism of dilute DNA solutions (A₂₆₀≈0.1) has been studied in a Couette-type apparatus with the outer cylinder rotating and with the light path parallel to the cylinder axis. Shear gradients in the range of 5-160 sec.^-1 were studied. The DNA samples were whole, "half," and "quarter" molecules of T4 bacteriophage DNA, and linear and circular λb₂b₅c DNA. For the linear molecules, the fractional flow dichroism is a linear function of molecular weight. The dichroism for linear A DNA is about 1.8 that of the circular molecule. For a given DNA, the dichroism is an approximately linear function of shear gradient, but with a slight upward curvature at low values of G, and some trend toward saturation at larger values of G. The fractional dichroism increases as the supporting electrolyte concentration decreases.