Fine mapping of copy number variations on two cattle genome assemblies using high density SNP array

ABSTRACT: BACKGROUND: Btau_4.0 and UMD3.1 are two distinct cattle reference genome assemblies. In our previous study using the low density BovineSNP50 array, we reported a copy number variation (CNV) analysis on Btau_4.0 with 521 animals of 21 cattle breeds, yielding 682 CNV regions with a total length of 139.8 megabases. RESULTS: In this study using the high density BovineHD SNP array, we performed high resolution CNV analyses on both Btau_4.0 and UMD3.1 with 674 animals of 27 cattle breeds. We first compared CNV results derived from these two different SNP array platforms on Btau_4.0. With two thirds of the animals shared between studies, on Btau_4.0 we identified 3,346 candidate CNV regions representing 142.7 megabases (~4.70%) of the genome. With a similar total length but 5 times more event counts, the average CNVR length of current Btau_4.0 dataset is significantly shorter than the previous one (42.7kb vs. 205 kb). Although subsets of these two results overlapped, 64% (91.6 megabases) of current dataset was not present in the previous study. We also performed similar analyses on UMD3.1 using these BovineHD SNP array results. Approximately 50% more and 20% longer CNVs were called on UMD3.1 as compared to those on Btau_4.0. However, a comparable result of CNVRs (3,438 regions with a total length 146.9 megabases) was obtained. We suspect that these results are due to that UMD3.1's efforts of placing unplaced contigs and removing unmerged alleles. Selected CNVs were further experimentally validated, achieving a 73% PCR validation rate, which is considerably higher than the previous validation rate. About 20-45% of CNV regions overlapped with cattle RefSeq genes and Ensembl genes. Panther and IPA analyses indicated that these genes provide a wide spectrum of biological processes involving immune system, lipid metabolism, cell, organism and system development. CONCLUSION: We present a comprehensive result of cattle CNVs at a higher resolution and sensitivity. We identified over 3,000 candidate CNV regions on both Btau_4.0 and UMD3.1, further compared current datasets with previous results, and examined the impacts of genome assemblies on CNV calling.     

Dependence of the Raman signature of genomic B-DNA on nucleotide base sequence.

The vibrational spectra of four genomic and two synthetic DNAs, encompassing a wide range in base composition [poly(dA-dT). poly(dA-dT), 0% G + C; Clostridium perfringens DNA, 27% G + C; calf thymus DNA, 42% G + C; Escherichia coli DNA, 50% G + C; Micrococcus luteus DNA, 72% G + C; poly(dG-dC).poly(dG-dC), 100% G + C] (dA: deoxyadenosine; dG: deoxyguanosine; dC: deoxycytidine; dT: thymidine), have been analyzed using Raman difference methods of high sensitivity. The results show that the Raman signature of B DNA depends in detail upon both genomic base composition and sequence. Raman bands assigned to vibrational modes of the deoxyribose-phosphate backbone are among the most sensitive to base sequence, indicating that within the B family of conformations major differences occur in the backbone geometry of AT- and GC-rich domains. Raman bands assigned to in-plane vibrations of the purine and pyrimidine bases-particularly of A and T-exhibit large deviations from the patterns expected for random base distributions, establishing that Raman hypochromic effects in genomic DNA are also highly sequence dependent. The present study provides a basis for future use of Raman spectroscopy to analyze sequence-specific DNA-ligand interactions. The demonstration of sequence dependency in the Raman spectrum of genomic B DNA also implies the capability to distinguish genomic DNAs by means of their characteristic Raman signatures.     

cDNA sequence of five mouse guanine deaminase (Gda) alleles and mapping to mouse chromosome 19.

Guanine deaminase catalyses the conversion of guanine to xanthine and ammonia, thereby irreversibly removing the guanine base from the pool of guanine-containing metabolites. We have identified five alleles at the mouse guanine deaminase locus by cDNA sequencing. These alleles were defined by single-nucleotide polymorphisms at a total of 19 positions. For each allele the representative strains are as follows: Gda(a), C57BL/6J and DBA/2J; Gda(b), A/J; Gda(c), MOLF/Ei; Gda(d), CAST/Ei; and Gda(e), SPRET-1. The only codon change resulting in an amino acid substitution was found at nucleotide 523, where GAT was replaced by AAT in Mus spretus resulting in the deduced substitution of Asp-174 by Asn. The single-nucleotide difference between the a and b alleles was also typed by allele-specific oligonucleotide amplification for 17 common strains of Mus musculus susbp. musculus. By typing the AxB and BxA recombinant inbred (RI) strain sets, Gda was mapped to mouse chromosome 19, a region syntenic with human chromosome 9q11-q22.     

Physical mapping of three minisatellite sequences in the Atlantic salmon (Salmo salar) genome

We have determined the localization of three minisatellite loci SsBglIIL.6 , SsBglIIU.20 and SsPstL.26 in the Atlantic salmon ( Salmo salar ) genome by fluorescent in situ hybridization (FISH). FISH analysis of the SsBglIIL.6 and SsBglIIU.20 minisatellites revealed that they are both located in single chromosome pairs allowing their direct identification; in contrast, the SsPstIL.26 probe hybridizes to four different chromosomal pairs. The analysis of chromosomal location of minisatellite sequences could be very useful for studying structural changes that have taken place during chromosome evolution in the karyotype of Salmo salar .     

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.     

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.     

Interspecific nucleotide sequence differences in the cytochrome B gene of indriidae (Primates, Strepsirhini)

The comparison of partial sequence of the mitochondrial cytochrome b gene nucleotides for different Indriidae allowed us to confirm the species status ofAvahi laniger, A. occidentalis, and Propithecus tattersalli. The nucleotide sequence allowed also to propose a phylogenetic tree which is discussed taking into account morphological, cytogenetic and former molecular biology data.     

Complete nucleotide sequence of the nucleoprotein gene of influenza B virus.

A DNA copy of influenza B/Singapore/222/79 viral RNA segment 5, containing the gene coding for the nucleoprotein (NP), has been cloned in Escherichia coli plasmid pBR322, and its nucleotide sequence has been determined. The influenza B NP gene contains 1,839 nucleotides and codes for a protein of 560 amino acids with a molecular weight of 61,593. Comparison of the influenza B NP amino acid sequence with that of influenza A NP (A/PR/8/34) reveals 37% direct homology in the aligned regions, indicating a common ancestor. However, influenza B NP has an additional 50 amino acids at its N-terminal end. As is the case with influenza A NP, influenza B NP is a basic protein, with its charged residues relatively evenly distributed rather than clustered. The structural homology suggests functional similarity between the NP of influenza A and B viruses.     

The cell biology of disease: FSHD: copy number variations on the theme of muscular dystrophy

In humans, copy number variations (CNVs) are a common source of phenotypic diversity and disease susceptibility. Facioscapulohumeral muscular dystrophy (FSHD) is an important genetic disease caused by CNVs. It is an autosomal-dominant myopathy caused by a reduction in the copy number of the D4Z4 macrosatellite repeat located at chromosome 4q35. Interestingly, the reduction of D4Z4 copy number is not sufficient by itself to cause FSHD. A number of epigenetic events appear to affect the severity of the disease, its rate of progression, and the distribution of muscle weakness. Indeed, recent findings suggest that virtually all levels of epigenetic regulation, from DNA methylation to higher order chromosomal architecture, are altered at the disease locus, causing the de-regulation of 4q35 gene expression and ultimately FSHD.     

Physical mapping and complete nucleotide sequence of the denV gene of bacteriophage T4.

Phage T4 deletion mutants that are folate analog resistant (far) and contain deletions in the region of the T4 genome near denV have been isolated previously. We showed that one of these mutants (T4farP12) expressed normal denV gene activity, whereas another mutant (T4farP13) was defective in the denV gene. The rII-distal (right) physical endpoints of these deletions defined the limits of the interval in which the rII-proximal (left) endpoint of the denV gene should be located. The deletion endpoints were identified by restriction and Southern hybridization analyses of phage derivatives containing deoxycytidine instead of hydroxymethyldeoxycytidine in their DNAs. The results of these analyses localized the rII-proximal (left) end of the denV gene to a region between 62.4 and 64.3 kilobases on the T4 physical map. denV+ phage resulted from marker rescue with two of five denV- alleles tested, using plasmids containing a 1.8-kilobase fragment from this region or a 179-base-pair terminal fragment derived from it. Sequencing of the 179-base-pair fragment from wild-type DNA showed a 130-base-pair open reading frame with its termination codon at the rII-proximal end. Confirmation that this open reading frame is part of the denV coding sequence was obtained by identifying a TAG amber codon in the homologous DNA derived from a denV amber mutant strain. This mutant strain rescued the denV+ allele from plasmids containing the wild-type sequence. An adjacent overlapping restriction fragment was also cloned, permitting determination of the remaining denV gene sequence. Based on these results, the 3' end of the coding region of the denV locus was mapped to kilobase position 64.07 on the T4 physical map, and the 5' end was mapped to position 64.48.     

Echinomycin binding to the sequence CG(AT)nCG alters the structure of the central AT region.

DNA fragments containing the sequence CG(AT)nCG have been used in footprinting experiments to assess the effect of echinomycin, which binds to CG steps, on the structure of the central AT region. DNAase I normally cuts ApT much better than TpA; in the presence of the drug this preference is retained but cleavage at TpA is enhanced. Changes in cleavage by micrococcal nuclease have also been observed. Echinomycin renders alternate adenines hyperreactive to diethylpyrocarbonate. The results suggest that echinomycin induces structural changes in regions surrounding its binding site and that these can be cooperatively propagated over several turns of the DNA helix.     

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.     

Genome-wide analysis of copy number variations reveals that aging processes influence body fat distribution in Korea Associated Resource (KARE) cohorts

Many anthropometric measures, including body mass index (BMI), waist-to-hip ratio (WHR), and subcutaneous fat thickness, are used as indicators of nutritional status, fertility and predictors of future health outcomes. While BMI is currently the best available estimate of body adiposity, WHR and skinfold thickness at various sites (biceps, triceps, suprailiac, and subscapular) are used as indices of body fat distribution. Copy number variation (CNV) is an attractive emerging approach to the study of associations with various diseases. In this study, we investigated the dosage effect of genes in the CNV genome widely associated with fat distribution phenotypes in large cohorts. We used the Affymetrix genome-wide human SNP Array 5.0 data of 8,842 healthy unrelated adults in KARE cohorts and identified CNVs associated with BMI and fat distribution-related traits including WHR and subcutaneous skinfold thickness at suprailiac (SUP) and subscapular (SUB) sites. CNV segmentation of each chromosome was performed using Golden Helix SVS 7.0, and single regression analysis was used to identify CNVs associated with each phenotype. We found one CNV for BMI, 287 for WHR, 2,157 for SUP, and 2,102 for SUB at the 5?% significance level after Holm–Bonferroni correction. Genes included in the CNV were used for the analysis of functional annotations using the Database for Annotation, Visualization and Integrated Discovery (DAVID v6.7b) tool. Functional gene classification analysis identified five significant gene clusters (metallothionein, ATP-binding proteins, ribosomal proteins, kinesin family members, and zinc finger proteins) for SUP, three (keratin-associated proteins, zinc finger proteins, keratins) for SUB, and one (protamines) for WHR. BMI was excluded from this analysis because the entire structure of no gene was identified in the CNV. Based on the analysis of genes enriched in the clusters, the fat distribution traits of KARE cohorts were related to the fat redistribution associated with the aging process. In addition to structural variation, dosage effect analysis of genes based on CNV is useful to gain an understanding of the comprehensive biological phenomena underlying particular phenotypes and/or diseases.