Human defensin gene copy number polymorphisms: comprehensive analysis of independent variation in alpha- and beta-defensin regions at 8p22-p23

To investigate defensin gene copy number polymorphisms, a quantitative real-time PCR assay was developed and used to study DNA from 27 unrelated individuals of diverse ethnic and racial backgrounds. The DEFB4 and DEFB103A genes varied in tandem, with copy numbers 2 to 8, with a mode of 6 per diploid genome (PDG). The combined copy numbers of the DEFA1 and DEFA3 genes ranged from 5 to 14, with a mode of 10 copies PDG. The copy numbers of the DEFA1/3 genes varied independently of those of the DEFB4 and DEFB103A genes. The amount of HNP-1 and HNP-3 peptides expressed in neutrophils was found to be proportional to the combined copy number of DEFA1 and DEFA3. The DEFA3 allele was absent in 7/27 subjects. The highly copy-number-variable DEFA1 and DEFA3 genes are flanked by other defensin genes present uniformly at 2 copies PDG. The remarkable variability in defensin gene copy numbers could contribute to differences in individual resistance to infections.     

Exome sequencing and arrayCGH detection of gene sequence and copy number variation between ILS and ISS mouse strains

It has been well documented that genetic factors can influence predisposition to develop alcoholism. While the underlying genomic changes may be of several types, two of the most common and disease associated are copy number variations (CNVs) and sequence alterations of protein coding regions. The goal of this study was to identify CNVs and single-nucleotide polymorphisms that occur in gene coding regions that may play a role in influencing the risk of an individual developing alcoholism. Toward this end, two mouse strains were used that have been selectively bred based on their differential sensitivity to alcohol: the Inbred long sleep (ILS) and Inbred short sleep (ISS) mouse strains. Differences in initial response to alcohol have been linked to risk for alcoholism, and the ILS/ISS strains are used to investigate the genetics of initial sensitivity to alcohol. Array comparative genomic hybridization (arrayCGH) and exome sequencing were conducted to identify CNVs and gene coding sequence differences, respectively, between ILS and ISS mice. Mouse arrayCGH was performed using catalog Agilent 1 × 244 k mouse arrays. Subsequently, exome sequencing was carried out using an Illumina HiSeq 2000 instrument. ArrayCGH detected 74 CNVs that were strain-specific (38 ILS/36 ISS), including several ISS-specific deletions that contained genes implicated in brain function and neurotransmitter release. Among several interesting coding variations detected by exome sequencing was the gain of a premature stop codon in the alpha-amylase 2B (AMY2B) gene specifically in the ILS strain. In total, exome sequencing detected 2,597 and 1,768 strain-specific exonic gene variants in the ILS and ISS mice, respectively. This study represents the most comprehensive and detailed genomic comparison of ILS and ISS mouse strains to date. The two complementary genome-wide approaches identified strain-specific CNVs and gene coding sequence variations that should provide strong candidates to contribute to the alcohol-related phenotypic differences associated with these strains.     

A 450-kb contig of defensin genes on human chromosome 8p23.

Defensins are a large family of host defense peptides expressed in leukocytes and epithelia. Using P1 and BAC clones, we have determined the organization of the human alpha-defensin genes and the beta-defensin gene HDEFB1 on chromosome 8p23. From the telomere, the order of the genes (with encoded peptides in parentheses) is HDEFA5 (HD-5), HDEFA1/1A (HNP-1/3), HDEFA4 (HNP-4), HDEFA6 (HD-6), and HDEFB1 (HBD-1). These genes span a region of approximately 450kb. Genes encoding intestinal Paneth cell defensins (HDEFA5 and HDEFA6) flank the myeloid defensin gene cluster (HDEFA1, HDEFA1A, HDEFA4). Based on our previous studies, the remaining known defensin gene, HDEFB2 (HBD-2), is about 400kb centromeric to HDEFB1. This map supports the hypothesis, originally proposed because of sequence similarities, that myeloid alpha-defensin genes evolved by reduplication and divergence from Paneth cell defensin genes, and identifies regions and clones, which should be useful in the search for new defensin genes.     

A distance-type measure approach to the analysis of copy number variation in DNA sequencing data

Background

The next generation sequencing technology allows us to obtain a large amount of short DNA sequence (DNA-seq) reads at a genome-wide level. DNA-seq data have been increasingly collected during the recent years. Count-type data analysis is a widely used approach for DNA-seq data. However, the related data pre-processing is based on the moving window method, in which a window size need to be defined in order to obtain count-type data. Furthermore, useful information can be reduced after data pre-processing for count-type data.

Results

In this study, we propose to analyze DNA-seq data based on the related distance-type measure. Distances are measured in base pairs (bps) between two adjacent alignments of short reads mapped to a reference genome. Our experimental data based simulation study confirms the advantages of distance-type measure approach in both detection power and detection accuracy. Furthermore, we propose artificial censoring for the distance data so that distances larger than a given value are considered potential outliers. Our purpose is to simplify the pre-processing of DNA-seq data. Statistically, we consider a mixture of right censored geometric distributions to model the distance data. Additionally, to reduce the GC-content bias, we extend the mixture model to a mixture of generalized linear models (GLMs). The estimation of model can be achieved by the Newton-Raphson algorithm as well as the Expectation-Maximization (E-M) algorithm. We have conducted simulations to evaluate the performance of our approach. Based on the rank based inverse normal transformation of distance data, we can obtain the related z-values for a follow-up analysis. For an illustration, an application to the DNA-seq data from a pair of normal and tumor cell lines is presented with a change-point analysis of z-values to detect DNA copy number alterations.

Conclusion

Our distance-type measure approach is novel. It does not require either a fixed or a sliding window procedure for generating count-type data. Its advantages have been demonstrated by our simulation studies and its practical usefulness has been illustrated by an experimental data application.

     

Proliferation and copy number variation of BEL-like long terminal repeat retrotransposons within the Diabrotica virgifera virgifera genome

The proliferation of retrotransposons within a genome can contribute to increased size and affect the function of eukaryotic genes. BEL/Pao-like long-terminal repeat (LTR) retrotransposons were annotated from the highly adaptable insect species Diabrotica virgifera virgifera, the Western corn rootworm, using survey sequences from bacterial artificial chromosome (BAC) inserts and contigs derived from a low coverage next-generation genome sequence assembly. Eleven unique D. v. virgifera BEL elements were identified that contained full-length gag–pol coding sequences, whereas 88 different partial coding regions were characterized from partially assembled elements. Estimated genome copy number for full and partial BEL-like elements ranged from ~ 8 to 1582 among individual contigs using a normalized depth of coverage (DOC) among Illumina HiSeq reads (total genome copy number ~ 8821). BEL element copy number was correlated among different D. v. virgifera populations (R² = 0.9846), but individual element numbers varied ≤ 1.68-fold and the total number varied by ~ 527 copies. These data indicate that BEL element proliferation likely contributed to a large genome size, and suggest that differences in copy number are a source of genetic variability among D. v. virgifera.     

General assessment of copy number variation in normal and tumor tissues of the domestic dog (Canis lupus familiaris)

In recent years, characterization of a copy number variation (CNV) of the genomic DNA has provided evidence for the relationship of this type of genetic variation with the occurrence of a broad spectrum of diseases, including cancer lesions. Copy number variants (CNVs) also occur in the genomes of healthy individuals as a result of abnormal recombination processes in germ cells and have a hereditary character contributing to the natural genetic diversity. Recent image analysis methods and advanced computational techniques allow for identification of CNVs using SNPs genotyping microarrays based on the analysis of signal intensity observed for markers located in the specific genomic regions. In this study we used CanineHD BeadChip assay (Illumina) to identify both natural and cancer-induced CNVs in the genomes of different dog breeds and in different cancer types occurring in this species. The obtained results showed that structural aberrations are a common phenomenon arising during a tumor progression and are more complex and widespread in tumors of mesenchymal tissue origin than in epithelial tissue originating tumors. The tumor derived CNVs, in comparison to healthy samples, were characterized by larger sizes of regions, higher number of amplifications, and in some cases encompassed genes with potential effect on tumor progression.     

Analysis of copy number variation in the normal human population within a region containing complex segmental duplications on 22q11 using high-resolution array-CGH

A previously detected copy number polymorphism (Ep CNP) in patients affected with neuroectodermal tumors led us to investigate its frequency and length in the normal population. For this purpose, a program called Sequence Allocator was developed and applied for the construction of an array that consisted of unique and duplicated fragments, allowing the assessment of copy number variation within regions of segmental duplications. The average resolution of this array was 11 kb and we determined the size of the Ep CNP to be 290 kb. Analysis of normal controls identified 7.7 and 7.1% gains in peripheral blood and lymphoblastoid cell line (LCL) DNA, respectively, while deletions were found only in the LCL group (7.1%). This array platform allows the detection of DNA copy number variation within regions of pronounced genomic complexity, which constitutes an improvement over available technologies.     

Increased paternal age and the influence on burden of genomic copy number variation in the general population

Genomic copy number variations (CNVs) and increased parental age are both associated with the risk to develop a variety of clinical neuropsychiatric disorders such as autism, schizophrenia and bipolar disorder. At the same time, it has been shown that the rate of transmitted de novo single nucleotide mutations is increased with paternal age. To address whether paternal age also affects the burden of structural genomic deletions and duplications, we examined various types of CNV burden in a large population sample from the Netherlands. Healthy participants with parental age information (n = 6,773) were collected at different University Medical Centers. CNVs were called with the PennCNV algorithm using Illumina genome-wide SNP array data. We observed no evidence in support of a paternal age effect on CNV load in the offspring. Our results were negative for global measures as well as several proxies for de novo CNV events in this unique sample. While recent studies suggest de novo single nucleotide mutation rate to be dominated by the age of the father at conception, our results strongly suggest that at the level of global CNV burden there is no influence of increased paternal age. While it remains possible that local genomic effects may exist for specific phenotypes, this study indicates that global CNV burden and increased father’s age may be independent disease risk factors.     

Effect of plasmid copy number and lac operator sequence on antibiotic-free plasmid selection by operator-repressor titration in Escherichia coli

The Escherichia coli strain DH1lacdapD enables plasmid selection and maintenance that is free from antibiotics and selectable marker genes. This is achieved by using only the lac operator sequence as a selectable element. This strain is currently used to generate high copy number plasmids with no antibiotic resistance genes for use as DNA vaccines and for expression of recombinant proteins. Until now these have been limited to pUC-based plasmids containing a high copy number pMB1-derived origin of replication, and the principle lacO(1) and auxiliary lacO(3) operators. In this study we have shown that this system can also be used to select and maintain pBR322-based plasmids with the lower copy number pMB1 origin of replication, and that lacO(1) alone or a palindromic version of lacO(1) can provide a sufficient level of repressor titration for plasmid selection. This is advantageous for recombinant protein production, where low copy number plasmids are often used and plasmid maintenance is important. The degree of repressor titration due to these plasmids was measured using the natural lactose operon in E. coli DH1 as a model.     

Evidence for duplication of the human defensin gene DEFB4 in chromosomal region 8p22-23 and implications for the analysis of SNP allele distribution

Defensins constitute a primary mechanism in the innate immune system of humans and all mammals. Defensins are short, processed peptide molecules that are classified by structure into three groups: alpha-defensins, beta-defensins and theta-defensins. In humans, four beta-defensins have been described so far, corresponding to the products of the genes DEFB1 (hBD1, NM_005218), DEFB4 (hBD2, NM_004942.2), DEFB103 (hBD3, NM_018661), and DEFB104 (hBD4, NM_080389), respectively. All these genes have been mapped to chromosome 8p22-23. Much interest has been shown in genetic variation in the population at defensin loci to understand individual differences in disease susceptibility and severity. In this study, we have used an electronic search and then fluorescence in situ hybridization (FISH) on elongated chromosomes to demonstrate that the region containing the DEFB4 gene is duplicated on human chromosome 8p, making difficult the discovery of new SNPs in this gene and compromising the assessment of their allelic distribution in various ethnic populations for disease association studies.     

Mutation detection by pyrosequencing: sequencing of exons 5-8 of the p53 tumor suppressor gene

The ability to sequence a large number of DNA samples rapidly and accurately for detection of all possible mutations is a critical goal for the future application of DNA sequencing in routine medical diagnostics. Pyrosequencing() is a non-electrophoretic real-time DNA sequencing method that uses the luciferase-luciferin light release as the detection signal for nucleotide incorporation into target DNA. For pyrosequencing of the human p53 gene, a nested multiplex PCR method for amplification of exons 5-8 was prepared. In order to investigate the use of pyrosequencing in mutation detection, DNA samples from skin-cancer patients were used. Two forms of nucleotide dispensation strategy were used, cyclic and programmed. Bi-directional pyrosequencing was performed and the overlapping sequence data produced were assembled to determine the sequence of the gene. Reliable sequencing data were obtained with both dispensation strategies, but some advantages were obtained using the programmed nucleotide dispensation approach, such as longer and faster reads, and fewer out-of-phase problems. The accuracy of pyrosequencing for detection of p53 mutations and allele distribution was demonstrated.     

Quality assessment of the Affymetrix U133A&B probesets by target sequence mapping and expression data analysis

Careful analysis of microarray probe design should be an obligatory component of MicroArray Quality Control (MACQ) project [Patterson et al., 2006; Shi et al., 2006] initiated by the FDA (USA) in order to provide quality control tools to researchers of gene expression profiles and to translate the microarray technology from bench to bedside. The identification and filtering of unreliable probesets are important preprocessing steps before analysis of microarray data. These steps may result in an essential improvement in the selection of differentially expressed genes, gene clustering and construction of co-regulatory expression networks. We revised genome localization of the Affymetrix U133A&B GeneChip initial (target) probe sequences, and evaluated the impact of erroneous and poorly annotated target sequences on the quality of gene expression data. We found about 25% of Affymetrix target sequences overlapping with interspersed repeats that could cause cross-hybridization effects. In total, discrepancies in target sequence annotation account for up to approximately 30% of 44692 Affymetrix probesets. We introduce a novel quality control algorithm based on target sequence mapping onto genome and GeneChip expression data analysis. To validate the quality of probesets we used expression data from large, clinically and genetically distinct groups of breast cancers (249 samples). For the first time, we quantitatively evaluated the effect of repeats and other sources of inadequate probe design on the specificity, reliability and discrimination ability of Affymetrix probesets. We propose that only functionally reliable Affymetrix probesets that passed our quality control algorithm (approximately 86%) for gene expression analysis should be utilized. The target sequence annotation and filtering is available upon request.     

Effects of genomic position and copy number of Acyl-ACP thioesterase transgenes on the level of the target fatty acids in Brassica napus L.

The effects of genomic position and copy number of acyl-acyl carrier protein (ACP) thioesterase (TE) transgenes on the major target fatty acid, either lauric acid (C12:0) or palmitic acid (C16:0) depending on the TE, in transgenic Brassica napus seed oil were investigated. Four transgenic parental lines, transformed individually with the bay-TE (Uc FatB1), elm-TE (Ua FatB1), nutmeg-TE (Mf FatB1) and Cuphea-TE (Ch FatB1) transgenes, were crossed with the non-transgenic recipient genotypes '212/86' or 'QO4'. Bay-TE and Cuphea-TE F₁ seeds, which carry half the number of the construct copies compared to the self-pollinated seeds of the transgenic parents, showed significantly lower levels of the target fatty acid. Doubled haploid (DH) lines were developed through microspore culture from F₁ hybrids with the elm-TE or the Cuphea-TE transgenes. DH lines carrying one to five copies of the Cuphea-TE transgene displayed a positive correlation between transgene copy number and the target fatty acid C16:0 level (r = 0.77**). DH lines with elm-TE transgene copies at four different loci showed different C16:0 levels, with one of the loci (E-II) leading to significantly higher C16:0 levels. This study supports the importance of the selection of high transgene copy number and/or the optimum genomic integration site in order to achieve maximum expression levels of the target fatty acid in transgenic oil quality modification.     

Polyadenylation can regulate ColE1 type plasmid copy number independently of any effect on RNAI decay by decreasing the interaction of antisense RNAI with its RNAII target

Replication of ColE1-type plasmids is regulated by RNAI, an antisense RNA that interacts with the replication pre-primer, RNAII. Exonucleolytic attack at the 3' end of RNAI is impeded in pcnB mutant bacteria, which lack poly(A) polymerase I-the principal RNA polyadenylase of E. coli; this leads to accumulation of an RNAI decay intermediate (RNAI(-5)) and dramatic reduction of the plasmid copy number. Here, we report that polyadenylation can also affect RNAI-mediated control of plasmid DNA replication by inhibiting interaction of RNAI(-5) with RNAII. We show that mutation of the host pcnB gene profoundly affects the plasmid copy number, even under experimental conditions that limit the effects of polyadenylation on RNAI(-5) decay. Moreover, poly(A) tails interfere with RNAI/RNAII interaction in vitro without producing any detectable alteration of RNAI secondary structure. Our results establish the existence of a previously undetected mechanism by which RNA polyadenylation can control plasmid copy number.     

Variant mapping of the Apo(B) AT rich minisatellite. Dependence on nucleotide sequence of the copy number variations. Instability of the non-canonical alleles.

Because of its variations in length, the AT rich Hyper-Variable Region (HVR) of the 3' end of the Apolipoprotein B gene is used as a polymorphic maker in genetic studies. It contains a SspI site in its repeated motif and we used this feature to precisely analyse the internal structure of the different alleles found at this locus in a Caucasian population. We performed total digestion on 194 alleles as well as Minisatellite Variant Repeat mapping (MVR mapping: partial digestion) on 54. The results show that the level of length variability (in copy number) of the 5' end of this locus is at least two times higher than that of the 3' end. This could be correlated with the difference in nucleotide sequence between the two parts of the HVR and suggests the dependence on the primary structure of the mechanism that produces length variability. A molecular model is proposed to explain this result. Moreover, the sharp analysis of the minisatellite structure by the distribution of SspI sites reveals differences between long and short alleles, indicating that in most cases, no recombination occurs between alleles of different sizes. Finally the rare alleles exhibit a non-canonical structure. These important points could explain the bimodal distribution of the frequencies of the alleles in the population.