Major influence of repetitive elements on disease-associated copy number variants (CNVs)

Copy number variants (CNVs) are important contributors to the human pathogenic genetic diversity as demonstrated by a number of cases reported in the literature. The high homology between repetitive elements may guide genomic stability which will give rise to CNVs either by non-allelic homologous recombination (NAHR) or non-homologous end joining (NHEJ). Here, we present a short guide based on previously documented cases of disease-associated CNVs in order to provide a general view on the impact of repeated elements on the stability of the genomic sequence and consequently in the origin of the human pathogenic variome.     

Gene copy number variations in adaptive evolution: The genomic distribution of gene copy number variations revealed by genetic mapping and their adaptive role in an undomesticated species,white spruce (Picea glauca)

Gene copy number variation (CNV) has been associated with phenotypic variability in animals and plants, but a genomewide understanding of their impacts on phenotypes is largely restricted to human and agricultural systems. As such, CNVs have rarely been considered in investigations of the genomic architecture of adaptation in wild species. Here, we report on the genetic mapping of gene CNVs in white spruce, which lacks a contiguous assembly of its large genome (~20 Gb), and their relationships with adaptive phenotypic variation. We detected 3,911 gene CNVs including de novo structural variations using comparative genome hybridization on arrays (aCGH) in a large progeny set. We inferred the heterozygosity at CNV loci within parents by comparing haploid and diploid tissues and genetically mapped 82 gene CNVs. Our analysis showed that CNVs were distributed over 10 linkage groups and identified four CNV hotspots that we predict to occur in other species of the Pinaceae. Significant relationships were found between 29 of the gene CNVs and adaptive traits based on regression analyses with timings of bud set and bud flush, and height growth, suggesting a role for CNVs in climate adaptation. The importance of CNVs in adaptive evolution of white spruce was also indicated by functional gene annotations and the clustering of 31% of the mapped adaptive gene CNVs in CNV hotspots. Taken together, these results illustrate the feasibility of studying CNVs in undomesticated species and represent a major step towards a better understanding of the roles of CNVs in adaptive evolution.     

The genetic effect of copy number variations on the risk of type 2 diabetes in a Korean population

Background

Unlike Caucasian populations, genetic factors contributing to the risk of type 2 diabetes mellitus (T2DM) are not well studied in Asian populations. In light of this, and the fact that copy number variation (CNV) is emerging as a new way to understand human genomic variation, the objective of this study was to identify type 2 diabetes–associated CNV in a Korean cohort.

Methodology/Principal Findings

Using the Illumina HumanHap300 BeadChip (317,503 markers), genome-wide genotyping was performed to obtain signal and allelic intensities from 275 patients with type 2 diabetes mellitus (T2DM) and 496 nondiabetic subjects (Total n = 771). To increase the sensitivity of CNV identification, we incorporated multiple factors using PennCNV, a program that is based on the hidden Markov model (HMM). To assess the genetic effect of CNV on T2DM, a multivariate logistic regression model controlling for age and gender was used. We identified a total of 7,478 CNVs (average of 9.7 CNVs per individual) and 2,554 CNV regions (CNVRs; 164 common CNVRs for frequency>1%) in this study. Although we failed to demonstrate robust associations between CNVs and the risk of T2DM, our results revealed a putative association between several CNVRs including chr15:45994758–45999227 (P = 8.6E-04, P^corr = 0.01) and the risk of T2DM. The identified CNVs in this study were validated using overlapping analysis with the Database of Genomic Variants (DGV; 71.7% overlap), and quantitative PCR (qPCR). The identified variations, which encompassed functional genes, were significantly enriched in the cellular part, in the membrane-bound organelle, in the development process, in cell communication, in signal transduction, and in biological regulation.

Conclusion/Significance

We expect that the methods and findings in this study will contribute in particular to genome studies of Asian populations.     

Molecular and clinical study of 18 families with ADCA type II: evidence for genetic heterogeneity and de novo mutation.

The SCA7 mutation has been found in 54 patients and 7 at-risk subjects from 17 families who have autosomal dominant cerebellar ataxia (ADCA) II with progressive pigmentary maculopathy. In one isolated case, haplotype reconstruction through three generations confirmed a de novo mutation owing to paternal meiotic instability. Different disease-associated haplotypes segregated among the SCA7-positive kindreds, which indicated a multiple origin of the mutation. One family with the clinical phenotype of ADCA type II did not have the CAG expansion that indicated locus heterogeneity. The distribution of the repeat size in 944 independent normal chromosomes from controls, unaffected at-risk subjects, and one affected individual fell into two ranges. The majority of the alleles were in the first range of 7-19 CAG repeats. A second range could be identified with 28-35 repeats, and we provide evidence that these repeats represent intermediate alleles that are prone to further expansion. The repeat size of the pathological allele, the widest reported for all CAG-repeat disorders, ranged from 37 to approximately 220. The repeat size showed significant negative correlation with both age at onset and age at death. Analysis of the clinical features in the patients with SCA7 confirmed that the most frequently associated features are pigmentary maculopathy, pyramidal tract involvement, and slow saccades. The subjects with <49 repeats tended to have a less complicated neurological phenotype and a longer disease duration, whereas the converse applied to subjects with >/=49 repeats. The degree of instability during meiotic transmission was greater than in all other CAG-repeat disorders and was particularly striking in paternal transmission, in which a median increase in repeat size of 6 and an interquartile range of 12 were observed, versus a median increase of 3 and interquartile range of 3.5 in maternal transmission.     

Fragile X syndrome and the (CGG)n mutation: two families with discordant MZ twins.

The fragile X phenotype has been found, in the majority of cases, to be due to the expansion of a CGG repeat in the 5'-UTR region of the FMR-1 gene, accompanied by methylation of the adjacent CpG island and inactivation of the FMR-1 gene. Although several important aspects of the genetics of fragile X have been resolved, it remains to be elucidated at which stage in development the transition from the premutation to the full mutation occurs. We present two families in which discordance between two sets of MZ twins illustrates two important genetic points. In one family, two affected MZ brothers differed in the number of CGG repeats, demonstrating in vivo mitotic instability of this CGG repeat and suggesting that the transition to the full mutation occurred postzygotically. In the second family, two MZ sisters had the same number of repeats, but only one was mentally retarded. When the methylation status of the FMR-1 CpG island was studied, we found that the majority of normal chromosomes had been inactivated in the affected twin, thus leading to the expression of the fragile X phenotype.     

Owl pellets as a source of DNA for genetic studies of small mammals

Owl pellets contain a good skeletal record of the small mammals consumed, and correspond to the undigested portions of prey which are regurgitated. These pellets are easy to find at the roosting site of owls. As it has been demonstrated that amplifiable DNA can be isolated from ancient bone remains, the possibility of using owl pellets as a source of DNA for small mammal genetics studies via the polymerase chain reaction has been investigated. The main uncertainties when isolating DNA from such a material are firstly the possibility that the extracted DNA would be too degraded during the digestion in the stomach of the owl, and secondly that extensive cross-contaminations could occur among the different prey consumed. The results obtained clearly demonstrate that cross-contamination does not occur, and that mitochondrial and nuclear DNA can be amplified using skulls of small mammals found in owl pellets as a source of DNA. The relative efficiency of two methods of DNA extraction is estimated and discussed. Thus, owl pellets represent a non-invasive sampling technique which provides a valuable source of DNA for studying population genetics of small mammals.     

Comparative-high resolution melting: a novel method of simultaneous screening for small mutations and copy number variations

Efficient and cost-effective screening for DNA sequence changes, both small mutations and copy number variations (CNVs), is a crucial aspect for routine genetic diagnostics as well as for basic research. In this study we present a development and evaluation of comparative-high resolution melting (C-HRM), a new approach for the simultaneous screening of small DNA changes and gene CNVs. In contrast to other methods, relative quantification in C-HRM is based on the results obtained during the melting process and calculations of the melting peak height ratio in the multiplex reaction. Validation of the method was conducted on DNA samples from 50 individuals from Duchenne muscular dystrophy (DMD) families, 50 probands diagnosed with familial adenomatous polyposis and a control group of 36 women and 36 men. The results of analyses conducted on fragments of the DMD and APC genes correspond completely (100 %) with the results of previous studies. C-HRM sensitivity in CNV detection was assessed through the analysis of mixed DNA samples with different proportions of a deletion carrier and wild type control. The results are presented as a linear regression with R ² of 0.9974 and imply the capability of the method to detect mosaics. C-HRM is an attractive and powerful alternative to other methods of point mutations and CNV detection with 100 % accuracy in our studied group.     

Model-building studies of Inovirus: genetic variations on a geometric theme

Inovirus (filamentous bacteriophage) is a simple system for studying the rules by which protein primary structure (amino acid sequence) controls secondary and higher order structure, and thereby function. The virus occurs naturally as a number of different strains with similar secondary and higher order structure, but the protein subunit that assembles to form the virion coat has quite different primary structures in different virus strains. Despite these differences in primary structure, the subunits of all strains have much the same size, about 50 residues, which are distributed by type in much the same way into three domains of primary structure: a collection of acidic residues in the N-terminal region, a hydrophobic domain of about 19 residues near the middle, and a collection of basic residues near the C-terminus. Each subunit can be closely approximated by an alpha-helix with its long axis roughly parallel to the fibre axis, sloping from large to small radius in the virion and interleaving between subunits in the next turn or level. The acidic residues near the N-terminus of the subunit face outwards on the virion surface, and explain the low isoelectric point of the virion; the basic residues near the C-terminus face inwards, where they neutralize the charge on the DNA at the core of the virion; and the hydrophobic central domain is involved in interactions which bind neighbouring subunits. Detailed X-ray fibre diffraction analysis of one strain gives the subunit structure. Comparative model-building studies of different strains illustrate the common structural principles.     

The use of ligand-based de novo design for scaffold hopping and sidechain optimization: two case studies

This paper describes the application of de novo design utilizing exclusively ligand information. In the current approach, ligand design criteria, including pharmacophores, similarity and desired properties are applied as part of a fitness function driving the design process, instead of using them as filters after the process. This allows relevant parts of chemical space to be explored more efficiently. Two case studies of successful ligand design are also presented, one aimed at scaffold hopping, the other for exploring substitution patterns around a novel scaffold.     

Genome-wide association study of copy number variations associated with pulmonary function measures in Korea Associated Resource (KARE) cohorts

Copy number variation (CNV) is an attractive emerging approach to study the association with various diseases. We performed a CNV-based genome-wide association study of pulmonary function measures (FEV(1), FVC, and FEV(1)/FVC) in KARE cohorts. Affymetrix Genome-wide Human SNP Array 5.0 was used to measure genome-wide variation and CNV segmentation was performed using Golden Helix SVS 7.0. Single and multivariate regressions were used for the association study using the R statistical package and the Dabatase for Annotation, Visualization and Integrated (DAVID v6.7b) tool for the functional annotation. We identified significantly associated 1260 CNVs with pulmonary function measures of FEV(1) and FVC. Functional gene classification and annotation analysis found 5 highly enriched clusters, the BPI/LBP/Plunc superfamily, myosin, serpin peptidase inhibitor, protein tyrosine phosphatase, and olfactory receptors. According to the functional annotation, gene-based CNVs are likely to be involved in the pathogenesis and inflammatory responsiveness of pulmonary diseases.     

cn.MOPS: mixture of Poissons for discovering copy number variations in next-generation sequencing data with a low false discovery rate

Quantitative analyses of next-generation sequencing (NGS) data, such as the detection of copy number variations (CNVs), remain challenging. Current methods detect CNVs as changes in the depth of coverage along chromosomes. Technological or genomic variations in the depth of coverage thus lead to a high false discovery rate (FDR), even upon correction for GC content. In the context of association studies between CNVs and disease, a high FDR means many false CNVs, thereby decreasing the discovery power of the study after correction for multiple testing. We propose 'Copy Number estimation by a Mixture Of PoissonS' (cn.MOPS), a data processing pipeline for CNV detection in NGS data. In contrast to previous approaches, cn.MOPS incorporates modeling of depths of coverage across samples at each genomic position. Therefore, cn.MOPS is not affected by read count variations along chromosomes. Using a Bayesian approach, cn.MOPS decomposes variations in the depth of coverage across samples into integer copy numbers and noise by means of its mixture components and Poisson distributions, respectively. The noise estimate allows for reducing the FDR by filtering out detections having high noise that are likely to be false detections. We compared cn.MOPS with the five most popular methods for CNV detection in NGS data using four benchmark datasets: (i) simulated data, (ii) NGS data from a male HapMap individual with implanted CNVs from the X chromosome, (iii) data from HapMap individuals with known CNVs, (iv) high coverage data from the 1000 Genomes Project. cn.MOPS outperformed its five competitors in terms of precision (1-FDR) and recall for both gains and losses in all benchmark data sets. The software cn.MOPS is publicly available as an R package at http://www.bioinf.jku.at/software/cnmops/ and at Bioconductor.     

Murine DNA cytosine C(5)-methyltransferase: in vitro studies of de novo methylation spreading

The preference of murine DNA (cytosine-5)-methyltransferase (Dnmt1) for single stranded DNA substrates is increased up to 50-fold by the presence of a proximal 5-methyl cytosine (5(me)C). This modulation is distance-dependent and is due to an enhanced binding affinity and minor changes in catalytic efficiency. No modulation was observed with double stranded DNA. Modulation requires that the 5(me)C moiety be attached to the DNA strand containing the CpG methylation target. Our results support a model in which 5(me)C binding by the enzyme occurs to at least one site outside the region involved in CpG recognition. No modulation in response to 5(me)C is observed with the bacterial enzyme M.SssI, which lacks the large N-terminal regulatory domain found in Dnmt1. We suggest that this allosteric modulation involves the N-terminal domain of Dnmt1.     

Loggerhead turtle eggshells as a source of maternal nuclear genomic DNA for population genetic studies

Tagging studies on nesting beaches are commonly used to estimate nesting frequency, remigration interval and nesting population size for marine turtle rookeries. Estimates of these demographic parameters from tagging projects may be biased because of the small scale of tagging efforts relative to female nest site fidelity and the logistical difficulty of intercepting all nesting females. Therefore, alternative and supplemental means of individual identification of nesting females are required. We demonstrate that maternal nuclear microsatellite DNA can be isolated from unincubated eggshells of the loggerhead sea turtle (Caretta caretta) through comparison of DNA extracted from 59 eggs collected within 15 h of oviposition and DNA derived from skin samples from respective nesting females. Scorable microsatellite genotypes were produced in 897 of 994 (90.2%) single-locus egg amplifications attempted. Among eggs from known females, 730 of 748 (97.6%) single-locus, egg-derived genotypes matched the respective skin-derived genotypes. Allelic dropout was the most common type of error, followed by the presence of nonmaternal, presumably paternal, alleles. Genotypes derived from unincubated eggshells permit individual assignment of nests and therefore demographic parameter estimates for loggerhead turtle nesting populations, despite genotyping errors that require further optimization. Although sampling unincubated eggs is destructive, this technique is noninvasive to nesting females and is applicable in marine turtle population genetics studies when individual resolution is required but direct interception of nesting females is undesirable or logistically infeasible.     

Amplified Fragment Length Polymorphism (AFLP) as a source of genetic markers for red algae

A recent PCR-based fingerprinting technique, amplified fragment length polymorphism (AFLP), was successfully applied to the red alga Chondrus crispus Stackh. This is apparently the first account to describe the application of AFLP methodology to an alga. Six isolates of C. crispus were analyzed by AFLP. A total of twenty-five primer pairs were screened and six primer pairs were selected for further investigation. Both conservative and variable markers were identified within and between populations; some markers were unique to individuals. As such, AFLP should prove useful as a source of genetic markers in algae for applications as diverse as genome mapping to population genetic investigations. This revised version was published online in June 2006 with corrections to the Cover Date.     

Origin pairing ('handcuffing') as a mode of negative control of P1 plasmid copy number.

In one family of bacterial plasmids, multiple initiator binding sites, called iterons, are used for initiation of plasmid replication as well as for the control of plasmid copy number. Iterons can also pair in vitro via the bound initiators. This pairing, called handcuffing, has been suggested to cause steric hindrance to initiation and thereby control the copy number. To test this hypothesis, we have compared copy numbers of isogenic miniP1 plasmid monomer and dimer. The dimer copy number was only one-quarter that of the monomer, suggesting that the higher local concentration of origins in the dimer facilitated their pairing. Physical evidence consistent with iteron-mediated pairing of origins preferentially in the dimer was obtained in vivo. Thus, origin handcuffing can be a mechanism to control P1 plasmid replication.     

Sensitive and specific real-time polymerase chain reaction assays to accurately determine copy number variations (CNVs) of human complement C4A, C4B, C4-long, C4-short, and RCCX modules: elucidation of C4 CNVs in 50 consanguineous subjects with defined HLA genotypes

Recent comparative genome hybridization studies revealed that hundreds to thousands of human genomic loci can have interindividual copy number variations (CNVs). One of such CNV loci in the HLA codes for the immune effector protein complement component C4. Sensitive, specific, and accurate assays to interrogate the C4 CNV and its associated polymorphisms by using submicrogram quantities of genomic DNA are needed for high throughput epidemiologic studies of C4 CNVs in autoimmune, infectious, and neurological diseases. Quantitative real-time PCR (qPCR) assays were developed using TaqMan chemistry and based on sequences specific for C4A and C4B genes, structural characteristics corresponding to the long and short forms of C4 genes, and the breakpoint region of RP-C4-CYP21-TNX (RCCX) modular duplication. Assignments for gene copy numbers were achieved by relative standard curve methods using cloned C4 genomic DNA covering 6 logs of DNA concentrations for calibrations. The accuracies of test results were cross-confirmed internally in each sample, as the sum of C4A plus C4B equals to the sum of C4L plus C4S or the total copy number of RCCX modules. These qPCR assays were applied to determine C4 CNVs from samples of 50 consanguineous subjects who were mostly homozygous in HLA genotypes. The results revealed eight HLA haplotypes with single C4 genes in monomodular RCCX that are associated with multiple autoimmune and infectious diseases and 32 bimodular, 4 trimodular, and one quadrimodular RCCX. These C4 qPCR assays are proven to be robust, sensitive, and reliable, as they have contributed to the elucidation of C4 CNVs in >1000 human samples with autoimmune and neurological diseases.