Allelic recombination between distinct genomic locations generates copy number diversity in human beta-defensins

HiCanu: accurate assembly of segmental duplications, satellites, and allelic variants from high-fidelity long reads

Complete and accurate genome assemblies form the basis of most downstream genomic analyses and are of critical importance. Recent genome assembly projects have relied on a combination of noisy long-read sequencing and accurate short-read sequencing, with the former offering greater assembly continuity and the latter providing higher consensus accuracy. The recently introduced Pacific Biosciences (PacBio) HiFi sequencing technology bridges this divide by delivering long reads (>10 kbp) with high per-base accuracy (>99.9%). Here we present HiCanu, a modification of the Canu assembler designed to leverage the full potential of HiFi reads via homopolymer compression, overlap-based error correction, and aggressive false overlap filtering. We benchmark HiCanu with a focus on the recovery of haplotype diversity, major histocompatibility complex (MHC) variants, satellite DNAs, and segmental duplications. For diploid human genomes sequenced to 30× HiFi coverage, HiCanu achieved superior accuracy and allele recovery compared to the current state of the art. On the effectively haploid CHM13 human cell line, HiCanu achieved an NG50 contig size of 77 Mbp with a per-base consensus accuracy of 99.999% (QV50), surpassing recent assemblies of high-coverage, ultralong Oxford Nanopore Technologies (ONT) reads in terms of both accuracy and continuity. This HiCanu assembly correctly resolves 337 out of 341 validation BACs sampled from known segmental duplications and provides the first preliminary assemblies of nine complete human centromeric regions. Although gaps and errors still remain within the most challenging regions of the genome, these results represent a significant advance toward the complete assembly of human genomes.

The Dichotomous Responses Driven by β-Defensins

Defensins are short, rapidly evolving, cationic antimicrobial host defence peptides with a repertoire of functions, still incompletely realised, that extends beyond direct microbial killing. They are released or secreted at epithelial surfaces, and in some cases, from immune cells in response to infection and inflammation. Defensins have been described as endogenous alarmins, alerting the body to danger and responding to inflammatory signals by promoting both local innate and adaptive systemic immune responses. However, there is now increasing evidence that they exert variable control on the response to danger; creating a dichotomous response that can suppress inflammation in some circumstances but exacerbate the response to danger and damage in others and, at higher levels, lead to a cytotoxic effect. Focussing in this review on human β-defensins, we discuss the evidence for their functions as proinflammatory, immune activators amplifying the response to infection or damage signals and/or as mediators of resolution of damage, contributing to a return to homeostasis. Finally, we consider their involvement in the development of autoimmune diseases.

Characterization and evolutionary dynamics of complex regions in eukaryotic genomes

Complex regions in eukaryotic genomes are typically characterized by duplications of chromosomal stretches that often include one or more genes repeated in a tandem array or in relatively close proximity. Nevertheless, the repetitive nature of these regions, together with the often high sequence identity among repeats, have made complex regions particularly recalcitrant to proper molecular characterization, often being misassembled or completely absent in genome assemblies. This limitation has prevented accurate functional and evolutionary analyses of these regions. This is becoming increasingly relevant as evidence continues to support a central role for complex genomic regions in explaining human disease, developmental innovations, and ecological adaptations across phyla. With the advent of long-read sequencing technologies and suitable assemblers, the development of algorithms that can accommodate sample heterozygosity, and the adoption of a pangenomic-like view of these regions, accurate reconstructions of complex regions are now within reach. These reconstructions will finally allow for accurate functional and evolutionary studies of complex genomic regions, underlying the generation of genotype-phenotype maps of unprecedented resolution.

β-Defensins Coordinate In Vivo to Inhibit Bacterial Infections of the Trachea

β-defensins are predicted to play an important role in innate immunity against bacterial infections in the airway. We previously observed that a type III-secretion product of Bordetella bronchiseptica inhibits the NF-κB-mediated induction of a β-defensin in airway epithelial cells in vitro. To confirm this in vivo and to examine the relative roles of other β-defensins in the airway, we infected wild-type C57BL/6 mice and mice with a deletion of the mBD-1 gene with B. bronchiseptica wild-type strain, RB50 and its mutant strain lacking the type III-secretion system, WD3. The bacteria were quantified in the trachea and the nasal tissue and mRNA levels of mouse β-defensin-3 (mBD-3) were assessed after 24 h. Infection with the wild-type bacterial strain resulted in lower mBD-3 mRNA levels in the trachea than in mice infected with the type III-deficient strain. Furthermore, we observed an increase in bacterial numbers of RB50 only in the tracheas of mBD-1-deficient mice. Neutrophils were also more abundant on the trachea in RB50 infected WT mice but not in the bronchiolar lavage fluid (BAL), compared with WD3 infected WT and mBD-1^−/− mice, indicating that the coordination of β-defensin chemotactic effects may be confined to tracheal epithelial cells (TEC). RB50 decreased the ability of mice to mount an early specific antibody response, seven days after infection in both WT and mBD-1^−/− mice but there were no differences in titers between RB50-infected WT and mBD-1^−/− mice or between WD3-infected WT and mBD-1^−/− mice, indicating mBD-1 was not involved in induction of the humoral immune response to the B. bronchiseptica. Challenge of primary mouse TEC in vitro with RB50 and WD3, along with IL-1β, further corroborated the in vivo studies. The results demonstrate that at least two β-defensins can coordinate early in an infection to limit the growth of bacteria in the trachea.

Human beta defensin (HBD) gene copy number affects HBD2 protein levels: impact on cervical bactericidal immunity in pregnancy

Human beta defensin 2 (HBD2) is an endogenous mucosal antimicrobial peptide (AMP) upregulated during infection and inflammation. HBD2 is encoded by the DEFB4 gene, which exhibits extensive copy number variation. Previous studies have demonstrated a relationship between HBD copy number and serum HBD2 protein levels; however, our current understanding of the influence of copy number on mucosal AMP function remains limited. This study explores the relationship between HBD copy number, cervicovaginal HBD2 protein levels and antimicrobial activity in 203 women with risk factors for preterm birth. We provide evidence that suggests HBD copy number modulates cervical antimicrobial immunity.

Analysis of Copy Number Variation Using the Paralogue Ratio Test (PRT)

Copy number variation (CNV), where a segment of DNA differs in copy number between different individuals, is an extensive and often underappreciated source of genetic variation within species. However, reliably determining copy number of a particular DNA sequence for a large number of samples can be challenging. Here, I describe and review the paralogue ratio test (PRT) in detail. PRT was developed to robustly type the CNV of the beta-defensin locus using small amounts of genomic DNA in a high-throughput manner, and has been applied successfully at many other loci. I discuss the strategies for designing successful PRT assays using both manual and bioinformatics methods, how to optimize experimental conditions, and approaches for analyzing the data. I discuss strengths and weaknesses of the approach, and how to troubleshoot results, as well as the range of problems to which PRT can be a potential solution.

Interchromosomal core duplicons drive both evolutionary instability and disease susceptibility of the Chromosome 8p23.1 region

Recurrent rearrangements of Chromosome 8p23.1 are associated with congenital heart defects and developmental delay. The complexity of this region has led to inconsistencies in the current reference assembly, confounding studies of genetic variation. Using comparative sequence-based approaches, we generated a high-quality 6.3-Mbp alternate reference assembly of an inverted Chromosome 8p23.1 haplotype. Comparison with nonhuman primates reveals a 746-kbp duplicative transposition and two separate inversion events that arose in the last million years of human evolution. The breakpoints associated with these rearrangements map to an ape-specific interchromosomal core duplicon that clusters at sites of evolutionary inversion (P = 7.8 × 10⁻⁵). Refinement of microdeletion breakpoints identifies a subgroup of patients that map to the same interchromosomal core involved in the evolutionary formation of the duplication blocks. Our results define a higher-order genomic instability element that has shaped the structure of specific chromosomes during primate evolution contributing to rearrangements associated with inversion and disease.

Low α-defensin gene copy number increases the risk for IgA nephropathy and renal dysfunction

Although a major source of genetic variation, copy number variations (CNVs) and their involvement in disease development have not been well studied. Immunoglobulin A nephropathy (IgAN) is the most common primary glomerulonephritis worldwide. We performed association analysis of the DEFA1A3 CNV locus in two independent IgAN cohorts of southern Chinese Han (total of 1189 cases and 1187 controls). We discovered three independent copy number associations within the locus: DEFA1A3 [P = 3.99 × 10−9; odds ratio (OR), 0.88], DEFA3 (P = 6.55 × 10−5; OR, 0.82), and a noncoding deletion variant (211bp) (P = 3.50 × 10−16; OR, 0.75) (OR per copy, fixed-effects meta-analysis). While showing strong association with an increased risk for IgAN (P = 9.56 × 10−20), low total copy numbers of the three variants also showed significant association with renal dysfunction in patients with IgAN (P = 0.03; hazards ratio, 3.69; after controlling for the effects of known prognostic factors) and also with increased serum IgA1 (P = 0.02) and galactose-deficient IgA1 (P = 0.03). For replication, we confirmed the associations of DEFA1A3 (P = 4.42 × 10−4; OR, 0.82) and DEFA3 copy numbers (P = 4.30 × 10−3; OR, 0.74) with IgAN in a Caucasian cohort (531 cases and 198 controls) and found the 211bp variant to be much rarer in Caucasians. We also observed an association of the 211bp copy number with membranous nephropathy (P = 1.11 × 10−7; OR, 0.74; in 493 Chinese cases and 500 matched controls), but not with diabetic kidney disease (in 806 Chinese cases and 786 matched controls). By explaining 4.96% of disease risk and influencing renal dysfunction in patients with IgAN, the DEFA1A3 CNV locus may be a potential therapeutic target for developing treatments for this disease.

Novel phenotype of mouse spermatozoa following deletion of nine β-defensin genes

β-defensin peptides are a large family of antimicrobial peptides. Although they kill microbes in vitro and interact with immune cells, the precise role of these genes in vivo remains uncertain. Despite their inducible presence at mucosal surfaces, their main site of expression is the epididymis. Recent evidence suggests that a major function of these peptides is in sperm maturation. In addition to previous work suggesting this, work at the MRC Human Genetics Unit, Edinburgh, has shown that homozygous deletion of a cluster of nine β-defensin genes in the mouse results in profound male sterility. The spermatozoa derived from the mutants had reduced motility and increased fragility. Epididymal spermatozoa isolated from the cauda region of the homozygous mutants demonstrated precocious capacitation and increased spontaneous acrosome reactions compared with those from wild-types. Despite this, these mutant spermatozoa had reduced ability to bind to the zona pellucida of oocytes. Ultrastructural examination revealed a disintegration of the microtubule structure of mutant-derived spermatozoa isolated from the epididymal cauda region, but not from the caput. Consistent with premature acrosome reaction and hyperactivation, spermatozoa from mutant animals had significantly increased intracellular calcium content. This work demonstrates that in vivo β-defensins are essential for successful sperm maturation, and that their disruption alters intracellular calcium levels, which most likely leads to premature activation and spontaneous acrosome reactions that result in hyperactivation and loss of microtubule structure of the axoneme. Determining which of the nine genes are responsible for the phenotype and the relevance to human sperm function is important for future work on male infertility.

Determining multiallelic complex copy number and sequence variation from high coverage exome sequencing data

BACKGROUND

Copy number variation (CNV) is a major component of genomic variation, yet methods to accurately type genomic CNV lag behind methods that type single nucleotide variation. High-throughput sequencing can contribute to these methods by using sequence read depth, which takes the number of reads that map to a given part of the reference genome as a proxy for copy number of that region, and compares across samples. Furthermore, high-throughput sequencing also provides information on the sequence differences between copies within and between individuals.

METHODS

In this study we use high-coverage phase 3 exome sequences of the 1000 Genomes project to infer diploid copy number of the beta-defensin genomic region, a well-studied CNV that carries several beta-defensin genes involved in the antimicrobial response, signalling, and fertility. We also use these data to call sequence variants, a particular challenge given the multicopy nature of the region.

RESULTS

We confidently call copy number and sequence variation of the beta-defensin genes on 1285 samples from 26 global populations, validate copy number using Nanostring nCounter and triplex paralogue ratio test data. We use the copy number calls to verify the genomic extent of the CNV and validate sequence calls using analysis of cloned PCR products. We identify novel variation, mostly individually rare, predicted to alter amino-acid sequence in the beta-defensin genes. Such novel variants may alter antimicrobial properties or have off-target receptor interactions, and may contribute to individuality in immunological response and fertility.

CONCLUSIONS

Given that 81% of identified sequence variants were not previously in dbSNP, we show that sequence variation in multiallelic CNVs represent an unappreciated source of genomic diversity.

Adaptive potential of genomic structural variation in human and mammalian evolution

Because phenotypic innovations must be genetically heritable for biological evolution to proceed, it is natural to consider new mutation events as well as standing genetic variation as sources for their birth. Previous research has identified a number of single-nucleotide polymorphisms that underlie a subset of adaptive traits in organisms. However, another well-known class of variation, genomic structural variation, could have even greater potential to produce adaptive phenotypes, due to the variety of possible types of alterations (deletions, insertions, duplications, among others) at different genomic positions and with variable lengths. It is from these dramatic genomic alterations, and selection on their phenotypic consequences, that adaptations leading to biological diversification could be derived. In this review, using studies in humans and other mammals, we highlight examples of how phenotypic variation from structural variants might become adaptive in populations and potentially enable biological diversification. Phenotypic change arising from structural variants will be described according to their immediate effect on organismal metabolic processes, immunological response and physical features. Study of population dynamics of segregating structural variation can therefore provide a window into understanding current and historical biological diversification.

An evolutionary history of defensins: a role for copy number variation in maximizing host innate and adaptive immune responses

Defensins represent an evolutionary ancient family of antimicrobial peptides that play diverse roles in human health and disease. Defensins are cationic cysteine-containing multifunctional peptides predominantly expressed by epithelial cells or neutrophils. Defensins play a key role in host innate immune responses to infection and, in addition to their classically described role as antimicrobial peptides, have also been implicated in immune modulation, fertility, development, and wound healing. Aberrant expression of defensins is important in a number of inflammatory diseases as well as modulating host immune responses to bacteria, unicellular pathogens, and viruses. In parallel with their role in immunity, in other species, defensins have evolved alternative functions, including the control of coat color in dogs. Defensin genes reside in complex genomic regions that are prone to structural variations and some defensin family members exhibit copy number variation (CNV). Structural variations have mediated, and continue to influence, the diversification and expression of defensin family members. This review highlights the work currently being done to better understand the genomic architecture of the β-defensin locus. It evaluates current evidence linking defensin CNV to autoimmune disease (i.e., Crohn’s disease and psoriasis) as well as the contribution CNV has in influencing immune responses to HIV infection.

Large multiallelic copy number variations in humans

Thousands of genomic segments appear to be present in widely varying copy numbers in different human genomes. We developed ways to use increasingly abundant whole-genome sequence data to identify the copy numbers, alleles and haplotypes present at most large multiallelic CNVs (mCNVs). We analyzed 849 genomes sequenced by the 1000 Genomes Project to identify most large (>5-kb) mCNVs, including 3,878 duplications, of which 1,356 appear to have 3 or more segregating alleles. We find that mCNVs give rise to most human variation in gene dosage-seven times the combined contribution of deletions and biallelic duplications-and that this variation in gene dosage generates abundant variation in gene expression. We describe 'runaway duplication haplotypes' in which genes, including HPR and ORM1, have mutated to high copy number on specific haplotypes. We also describe partially successful initial strategies for analyzing mCNVs via imputation and provide an initial data resource to support such analyses.

Evidence of convergent evolution in humans and macaques supports an adaptive role for copy number variation of the β-defensin-2 gene

β-defensins are a family of important peptides of innate immunity, involved in host defense, immunomodulation, reproduction, and pigmentation. Genes encoding β-defensins show evidence of birth-and-death evolution, adaptation by amino acid sequence changes, and extensive copy number variation (CNV) within humans and other species. The role of CNV in the adaptation of β-defensins to new functions remains unclear, as does the adaptive role of CNV in general. Here, we fine-map CNV of a cluster of β-defensins in humans and rhesus macaques. Remarkably, we found that the structure of the CNV is different between primates, with distinct mutational origins and CNV boundaries defined by retroviral long terminal repeat elements. Although the human β-defensin CNV region is 322 kb and encompasses several genes, including β-defensins, a long noncoding RNA gene, and testes-specific zinc-finger transcription factors, the orthologous region in the rhesus macaque shows CNV of a 20-kb region, containing only a single gene, the ortholog of the human β-defensin-2 gene. Despite its independent origins, the range of gene copy numbers in the rhesus macaque is similar to humans. In addition, the rhesus macaque gene has been subject to divergent positive selection at the amino acid level following its initial duplication event between 3 and 9.5 Ma, suggesting adaptation of this gene as the macaque successfully colonized novel environments outside Africa. Therefore, the molecular phenotype of β-defensin-2 CNV has undergone convergent evolution, and this gene shows evidence of adaptation at the amino acid level in rhesus macaques.

Human gene copy number variation and infectious disease

Variability in the susceptibility to infectious disease and its clinical manifestation can be determined by variation in the environment and by genetic variation in the pathogen and the host. Despite several successes based on candidate gene studies, defining the host variation affecting infectious disease has not been as successful as for other multifactorial diseases. Both single nucleotide variation and copy number variation (CNV) of the host contribute to the host's susceptibility to infectious disease. In this review we focus on CNV, particularly on complex multiallelic CNV that is often not well characterised either directly by hybridisation methods or indirectly by analysis of genotypes and flanking single nucleotide variants. We summarise the well-known examples, such as α-globin deletion and susceptibility to severe malaria, as well as more recent controversies, such as the extensive CNV of the chemokine gene CCL3L1 and HIV infection. We discuss the potential biological mechanisms that could underly any genetic association and reflect on the extensive complexity and functional variation generated by a combination of CNV and sequence variation, as illustrated by the Fc gamma receptor genes FCGR3A, FCGR3B and FCGR2C. We also highlight some understudied areas that might prove fruitful areas for further research.

Copy number variation of the beta defensin gene cluster on chromosome 8p influences the bacterial microbiota within the nasopharynx of otitis-prone children

As there is increasing evidence that aberrant defensin expression is related to susceptibility for infectious disease and inflammatory disorders, we sought to determine if copy number of the beta-defensin gene cluster located on chromosome 8p23.1 (DEFB107, 106, 105, 104, 103, DEFB4 and SPAG11), that shows copy number variation as a block, was associated with susceptibility to otitis media (OM). The gene DEFB103 within this complex encodes human beta defensin-3 (hBD-3), an antimicrobial peptide (AP) expressed by epithelial cells that line the mammalian airway, important for defense of mucosal surfaces and previously shown to have bactericidal activity in vitro against multiple human pathogens, including the three that predominate in OM. To this end, we conducted a retrospective case-control study of 113 OM prone children and 267 controls aged five to sixty months. We identified the copy number of the above defined beta-defensin gene cluster (DEFB-CN) in each study subject by paralogue ratio assays. The mean DEFB-CN was indistinguishable between subjects classified as OM prone based on a recent history of multiple episodes of OM and control subjects who had no history of OM (4.4 ± 0.96 versus 4.4 ± 1.08, respectively: Odds Ratio [OR]: 1.16 (95% CI: 0.61, 2.20). Despite a lack of direct association, we observed a statistically significant correlation between DEFB-CN and nasopharyngeal bacterial colonization patterns. Collectively, our findings suggested that susceptibility to OM might be mediated by genetic variation among individuals, wherein a DEFB-CN less than 4 exerts a marked influence on the microbiota of the nasopharynx, specifically with regard to colonization by the three predominant bacterial pathogens of OM.

Narrowing down the distal border of the copy number variable beta-defensin gene cluster on human 8p23

Background

Copy number variation (CNV) in the range from 2 to 12 per diploid genome is an outstanding feature of the beta-defensin gene (DEFB) cluster on human chromosome 8p23.1 numerously demonstrated by different methods. So far, CNV was proven for a 115 kb region between DEFB4 and 21 kb proximal of DEFB107 but the borders for the entire CNV repeat unit are still unknown. Our study aimed to narrow down the distal border of the DEFB cluster.

Results

We established tests for length polymorphisms based on amplification and capillary electrophoresis with laser-induced fluorescence (CE-LIF) analysis of seven insertion/deletion (indel) containing regions spread over the entire cluster. The tests were carried out with 25 genomic DNAs with different previously determined cluster copy numbers. CNV was demonstrated for six indels between ~1 kb distal of DEFB108P and 10 kb proximal of DEFB107. In contrast, the most distal indel is not affected by CNV.

Conclusion

Our analysis fixes the minimal length of proven CNV to 157 kb including DEFB108P but excluding DEFB109P. The distal border between CNV and non-CNV part of the DEF cluster is located in the 59 kb interval chr8:7,171,082-7,230,128.

The challenges and importance of structural variation detection in livestock

Recent studies in humans and other model organisms have demonstrated that structural variants (SVs) comprise a substantial proportion of variation among individuals of each species. Many of these variants have been linked to debilitating diseases in humans, thereby cementing the importance of refining methods for their detection. Despite progress in the field, reliable detection of SVs still remains a problem even for human subjects. Many of the underlying problems that make SVs difficult to detect in humans are amplified in livestock species, whose lower quality genome assemblies and incomplete gene annotation can often give rise to false positive SV discoveries. Regardless of the challenges, SV detection is just as important for livestock researchers as it is for human researchers, given that several productive traits and diseases have been linked to copy number variations (CNVs) in cattle, sheep, and pig. Already, there is evidence that many beneficial SVs have been artificially selected in livestock such as a duplication of the agouti signaling protein gene that causes white coat color in sheep. In this review, we will list current SV and CNV discoveries in livestock and discuss the problems that hinder routine discovery and tracking of these polymorphisms. We will also discuss the impacts of selective breeding on CNV and SV frequencies and mention how SV genotyping could be used in the future to improve genetic selection.

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.

Copy number variation of the beta-defensin genes in europeans: no supporting evidence for association with lung function, chronic obstructive pulmonary disease or asthma

Lung function measures are heritable, predict mortality and are relevant in diagnosis of chronic obstructive pulmonary disease (COPD). COPD and asthma are diseases of the airways with major public health impacts and each have a heritable component. Genome-wide association studies of SNPs have revealed novel genetic associations with both diseases but only account for a small proportion of the heritability. Complex copy number variation may account for some of the missing heritability. A well-characterised genomic region of complex copy number variation contains beta-defensin genes (DEFB103, DEFB104 and DEFB4), which have a role in the innate immune response. Previous studies have implicated these and related genes as being associated with asthma or COPD. We hypothesised that copy number variation of these genes may play a role in lung function in the general population and in COPD and asthma risk. We undertook copy number typing of this locus in 1149 adult and 689 children using a paralogue ratio test and investigated association with COPD, asthma and lung function. Replication of findings was assessed in a larger independent sample of COPD cases and smoking controls. We found evidence for an association of beta-defensin copy number with COPD in the adult cohort (OR = 1.4, 95%CI:1.02-1.92, P = 0.039) but this finding, and findings from a previous study, were not replicated in a larger follow-up sample(OR = 0.89, 95%CI:0.72-1.07, P = 0.217). No robust evidence of association with asthma in children was observed. We found no evidence for association between beta-defensin copy number and lung function in the general populations. Our findings suggest that previous reports of association of beta-defensin copy number with COPD should be viewed with caution. Suboptimal measurement of copy number can lead to spurious associations. Further beta-defensin copy number measurement in larger sample sizes of COPD cases and children with asthma are needed.

Partial deletion of chromosome 8 β-defensin cluster confers sperm dysfunction and infertility in male mice

β-defensin peptides are a family of antimicrobial peptides present at mucosal surfaces, with the main site of expression under normal conditions in the male reproductive tract. Although they kill microbes in vitro and interact with immune cells, the precise role of these genes in vivo remains uncertain. We show here that homozygous deletion of a cluster of nine β-defensin genes (DefbΔ9) in the mouse results in male sterility. The sperm derived from the mutants have reduced motility and increased fragility. Epididymal sperm isolated from the cauda should require capacitation to induce the acrosome reaction but sperm from the mutants demonstrate precocious capacitation and increased spontaneous acrosome reaction compared to wild-types but have reduced ability to bind the zona pellucida of oocytes. Ultrastructural examination reveals a defect in microtubule structure of the axoneme with increased disintegration in mutant derived sperm present in the epididymis cauda region, but not in caput region or testes. Consistent with premature acrosome reaction, sperm from mutant animals have significantly increased intracellular calcium content. Thus we demonstrate in vivo that β-defensins are essential for successful sperm maturation, and their disruption leads to alteration in intracellular calcium, inappropriate spontaneous acrosome reaction and profound male infertility.

β-defensins are small molecules, considered primarily to be antimicrobials and important in the first defence response to invading organisms. They are predominantly produced at surfaces in contact with the outside environment and these include skin, airway and reproductive tract. We show here that when we delete from the mouse a subset of nine β-defensin genes, surprisingly the main consequence is that the male mice are completely infertile. When normal sperm leave the male and enter the female reproductive tract they are triggered to undergo a reaction that alters the membrane properties of the sperm and allows fertilisation. We show here that sperm isolated from the male mice, that no longer make these β-defensins, are prematurely ready to fertilise an egg. It is far too early for this to happen and as a consequence the sperm are severely reduced in their ability to move and have a major defect in the structure of their tail. We provide evidence that the reason this has happened is due to a dysregulation of calcium transport. This work is important for understanding defensin gene function in a living organism and may enable the design of novel contraceptives with additional antibiotic ability.

Accurate measurement of gene copy number for human alpha-defensin DEFA1A3

Background

Multi-allelic copy number variants include examples of extensive variation between individuals in the copy number of important genes, most notably genes involved in immune function. The definition of this variation, and analysis of its impact on function, has been hampered by the technical difficulty of large-scale but accurate typing of genomic copy number. The copy-variable alpha-defensin locus DEFA1A3 on human chromosome 8 commonly varies between 4 and 10 copies per diploid genome, and presents considerable challenges for accurate high-throughput typing.

Results

In this study, we developed two paralogue ratio tests and three allelic ratio measurements that, in combination, provide an accurate and scalable method for measurement of DEFA1A3 gene number. We combined information from different measurements in a maximum-likelihood framework which suggests that most samples can be assigned to an integer copy number with high confidence, and applied it to typing 589 unrelated European DNA samples. Typing the members of three-generation pedigrees provided further reassurance that correct integer copy numbers had been assigned. Our results have allowed us to discover that the SNP rs4300027 is strongly associated with DEFA1A3 gene copy number in European samples.

Conclusions

We have developed an accurate and robust method for measurement of DEFA1A3 copy number. Interrogation of rs4300027 and associated SNPs in Genome-Wide Association Study SNP data provides no evidence that alpha-defensin copy number is a strong risk factor for phenotypes such as Crohn’s disease, type I diabetes, HIV progression and multiple sclerosis.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-719) contains supplementary material, which is available to authorized users.

A potential relationship among beta-defensins haplotype, SOX7 duplication and cardiac defects

OBJECTIVE

To determine the pathogenesis of a patient born with congenital heart defects, who had appeared normal in prenatal screening.

METHODS

In routine prenatal screening, G-banding was performed to analyse the karyotypes of the family and fluorescence in situ hybridization was used to investigate the 22q11.2 deletion in the fetus. After birth, the child was found to be suffering from heart defects by transthoracic echocardiography. In the following study, sequencing was used to search for potential mutations in pivotal genes. SNP-array was employed for fine mapping of the aberrant region and quantitative real-time PCR was used to confirm the results. Furthermore, other patients with a similar phenotype were screened for the same genetic variations. To compare with a control, these variations were also assessed in the general population.

RESULTS

The child and his mother each had a region that was deleted in the beta-defensin repeats, which are usually duplicated in the general population. Besides, the child carried a SOX7-gene duplication. While this duplication was not detected in his mother, it was found in two other patients with cardiac defects who also had the similar deletion in the beta-defensin repeats.

CONCLUSION

The congenital heart defects of the child were probably caused by a SOX7-gene duplication, which may be a consequence of the partial haplotype of beta-defensin regions at 8p23.1. To our knowledge, this is the first congenital heart defect case found to have the haplotype of beta-defensin and the duplication of SOX7.

Population genetic tools for dissecting innate immunity in humans

Innate immunity involves direct interactions between the host and microorganisms, both pathogenic and symbiotic, so natural selection is expected to strongly influence genes involved in these processes. Population genetics investigates the impact of past natural selection events on the genome of present-day human populations, and it complements immunological as well as clinical and epidemiological genetic studies. Recent data show that the impact of selection on the different families of innate immune receptors and their downstream signalling molecules varies considerably. This Review discusses these findings and highlights how they help to delineate the relative functional importance of innate immune pathways, which can range from being essential to being redundant.

Variation in human β-defensin genes: new insights from a multi-population study

Human β-defensin 2 (hBD-2) and hBD-3, encoded by DEFB4 and DEFB103A, respectively, have shown anti-HIV activity, and both genes exhibit copy number variation (CNV). Although the role of hBD-1, encoded by DEFB1, in HIV-1 infection is less clear, single nucleotide polymorphisms (SNPs) in DEFB1 may influence viral loads and disease progression. We examined the distribution of DEFB1 SNPs and DEFB4/103A CNV, and the relationship between DEFB1 SNPs and DEFB4/103A CNV using samples from two HIV/AIDS cohorts from the United States (n = 150) and five diverse populations from the Coriell Cell Repositories (n = 46). We determined the frequencies of 10 SNPs in DEFB1 using a post-PCR, oligonucleotide ligation detection reaction-fluorescent microsphere assay, and CNV in DEFB4/103A by real-time quantitative PCR. There were noticeable differences in the frequencies of DEFB1 SNP alleles and haplotypes among various racial/ethnic groups. The DEFB4/103A copy numbers varied from 2 to 8 (median, 4), and there was a significant difference between the copy numbers of self-identified whites and blacks in the US cohorts (Mann-Whitney U-test P = 0.04). A significant difference was observed in the distribution of DEFB4/103A CNV among DEFB1 -52G/A and -390T/A genotypes (Kruskal-Wallis P = 0.017 and 0.026, respectively), while not in the distribution of DEFB4/103A CNV among -52G/A_-44C/G_-20G/A diplotypes. These observations provide additional insights for further investigating the complex interplay between β-defensin genetic polymorphisms and susceptibility to, or the progression or severity of, HIV infection/disease.

Evolutionary history of copy-number-variable locus for the low-affinity Fcγ receptor: mutation rate, autoimmune disease, and the legacy of helminth infection

Both sequence variation and copy-number variation (CNV) of the genes encoding receptors for immunoglobulin G (Fcγ receptors) have been genetically and functionally associated with a number of autoimmune diseases. However, the molecular nature and evolutionary context of this variation is unknown. Here, we describe the structure of the CNV, estimate its mutation rate and diversity, and place it in the context of the known functional alloantigen variation of these genes. Deletion of Fcγ receptor IIIB, associated with systemic lupus erythematosus, is a result of independent nonallelic homologous recombination events with a frequency of approximately 0.1%. We also show that pathogen diversity, in particular helminth diversity, has played a critical role in shaping the functional variation at these genes both between mammalian species and between human populations. Positively selected amino acids are involved in the interaction with IgG and include some amino acids that are known polymorphic alloantigens in humans. This supports a genetic contribution to the hygiene hypothesis, which states that past evolution in the context of helminth diversity has left humans with an array of susceptibility alleles for autoimmune disease in the context of a helminth-free environment. This approach shows the link between pathogens and autoimmune disease at the genetic level and provides a strategy for interrogating the genetic variation underlying autoimmune-disease risk and infectious-disease susceptibility.

Exploring the role of copy number variants in human adaptation

Over the past decade, the ubiquity of copy number variants (CNVs, the gain or loss of genomic material) in the genomes of healthy humans has become apparent. Although some of these variants are associated with disorders, a handful of studies documented an adaptive advantage conferred by CNVs. In this review, we propose that CNVs are substrates for human evolution and adaptation. We discuss the possible mechanisms and evolutionary processes in which CNVs are selected, outline the current challenges in identifying these loci, and highlight that copy number variable regions allow for the creation of novel genes that may diversify the repertoire of such genes in response to rapidly changing environments. We expect that many more adaptive CNVs will be discovered in the coming years, and we believe that these new findings will contribute to our understanding of human-specific phenotypes.

The challenges of studying complex and dynamic regions of the human genome

Recent work has emphasised that the human genome is not simple and static, but complex and dynamic. This review focuses on the regions that are particularly hard to dissect and analyse, yet hold clues to how the genome changes during evolution and disease. I begin by summarising recent key advances in the understanding of the variable structure of our genome, and then I discuss a medley of methods that may allow us to analyse this structure in fine detail. In the final part, I describe potential future developments in this field, and make an argument that, just as we routinely genotype single-nucleotide polymorphisms now and will routinely re-sequence genomes in the near future, we should be aiming to physically re-map the individual human genome for each individual we study.

β-Defensins: multifunctional modulators of infection, inflammation and more?

Defensins comprise one of the largest groups of host defence peptides, present throughout evolution, in fungi and flowering plants as well as in invertebrates and vertebrates. These cysteine-rich, cationic peptides have a common ability to kill a broad range of microorganisms including bacteria, yeast and viruses. As such, they are a strong component of the arsenal that is an organism's innate immunity. It is becoming increasingly clear, however, that antimicrobial action is only one of the numerous roles of these multifunctional peptides. In recent years, the functions of defensins in immunomodulation have been widely investigated, and their involvement in other processes (such as fertility) is becoming evident. This review addresses recent advances in the immunomodulatory activity of β-defensins as well as the involvement of β-defensins in fertility, development, wound healing and cancer.

Big defensins, a diverse family of antimicrobial peptides that follows different patterns of expression in hemocytes of the oyster Crassostrea gigas

Background

Big defensin is an antimicrobial peptide composed of a highly hydrophobic N-terminal region and a cationic C-terminal region containing six cysteine residues involved in three internal disulfide bridges. While big defensin sequences have been reported in various mollusk species, few studies have been devoted to their sequence diversity, gene organization and their expression in response to microbial infections.

Findings

Using the high-throughput Digital Gene Expression approach, we have identified in Crassostrea gigas oysters several sequences coding for big defensins induced in response to a Vibrio infection. We showed that the oyster big defensin family is composed of three members (named Cg-BigDef1, Cg-BigDef2 and Cg-BigDef3) that are encoded by distinct genomic sequences. All Cg-BigDefs contain a hydrophobic N-terminal domain and a cationic C-terminal domain that resembles vertebrate β-defensins. Both domains are encoded by separate exons. We found that big defensins form a group predominantly present in mollusks and closer to vertebrate defensins than to invertebrate and fungi CSαβ-containing defensins. Moreover, we showed that Cg-BigDefs are expressed in oyster hemocytes only and follow different patterns of gene expression. While Cg-BigDef3 is non-regulated, both Cg-BigDef1 and Cg-BigDef2 transcripts are strongly induced in response to bacterial challenge. Induction was dependent on pathogen associated molecular patterns but not damage-dependent. The inducibility of Cg-BigDef1 was confirmed by HPLC and mass spectrometry, since ions with a molecular mass compatible with mature Cg-BigDef1 (10.7 kDa) were present in immune-challenged oysters only. From our biochemical data, native Cg-BigDef1 would result from the elimination of a prepropeptide sequence and the cyclization of the resulting N-terminal glutamine residue into a pyroglutamic acid.

Conclusions

We provide here the first report showing that big defensins form a family of antimicrobial peptides diverse not only in terms of sequences but also in terms of genomic organization and regulation of gene expression.

A worldwide analysis of beta-defensin copy number variation suggests recent selection of a high-expressing DEFB103 gene copy in East Asia

Beta-defensins are a family of multifunctional genes with roles in defense against pathogens, reproduction, and pigmentation. In humans, six beta-defensin genes are clustered in a repeated region which is copy-number variable (CNV) as a block, with a diploid copy number between 1 and 12. The role in host defense makes the evolutionary history of this CNV particularly interesting, because morbidity due to infectious disease is likely to have been an important selective force in human evolution, and to have varied between geographical locations. Here, we show CNV of the beta-defensin region in chimpanzees, and identify a beta-defensin block in the human lineage that contains rapidly evolving noncoding regulatory sequences. We also show that variation at one of these rapidly evolving sequences affects expression levels and cytokine responsiveness of DEFB103, a key inhibitor of influenza virus fusion at the cell surface. A worldwide analysis of beta-defensin CNV in 67 populations shows an unusually high frequency of high-DEFB103-expressing copies in East Asia, the geographical origin of historical and modern influenza epidemics, possibly as a result of selection for increased resistance to influenza in this region. Hum Mutat 32:743–750, 2011. © 2011 Wiley-Liss, Inc.

Determination of beta-defensin genomic copy number in different populations: a comparison of three methods

Background

There have been conflicting reports in the literature on association of gene copy number with disease, including CCL3L1 and HIV susceptibility, and β-defensins and Crohn's disease. Quantification of precise gene copy numbers is important in order to define any association of gene copy number with disease. At present, real-time quantitative PCR (QPCR) is the most commonly used method to determine gene copy number, however the Paralogue Ratio Test (PRT) is being used in more and more laboratories.

Findings

In this study we compare a Pyrosequencing-based Paralogue Ratio Test (PPRT) for determining beta-defensin gene copy number with two currently used methods for gene copy number determination, QPCR and triplex PRT by typing five different cohorts (UK, Danish, Portuguese, Ghanaian and Czech) of DNA from a total of 576 healthy individuals. We found a systematic measurement bias between DNA cohorts revealed by QPCR, but not by the PRT-based methods. Using PRT, copy number ranged from 2 to 9 copies, with a modal copy number of 4 in all populations.

Conclusions

QPCR is very sensitive to quality of the template DNA, generating systematic biases that could produce false-positive or negative disease associations. Both triplex PRT and PPRT do not show this systematic bias, and type copy number within the correct range, although triplex PRT appears to be a more precise and accurate method to type beta-defensin copy number.

Genome-wide association analysis identifies three psoriasis susceptibility loci

We carried out a meta-analysis of two recent psoriasis genome-wide association studies with a combined discovery sample of 1,831 affected individuals (cases) and 2,546 controls. One hundred and two loci selected based on P value rankings were followed up in a three-stage replication study including 4,064 cases and 4,685 controls from Michigan, Toronto, Newfoundland and Germany. In the combined meta-analysis, we identified three new susceptibility loci, including one at NOS2 (rs4795067, combined P = 4 × 10⁻¹¹), one at FBXL19 (rs10782001, combined P = 9 × 10⁻¹⁰) and one near PSMA6-NFKBIA (rs12586317, combined P = 2 × 10⁻⁸). All three loci were also associated with psoriatic arthritis (rs4795067, combined P = 1 × 10⁻⁵; rs10782001, combined P = 4 × 10⁻⁸; and rs12586317, combined P = 6 × 1⁻⁵) and purely cutaneous psoriasis (rs4795067, combined P = 1 × 10⁻⁸; rs10782001, combined P = 2 × 10⁻⁶; and rs12586317, combined P = 1 × 10⁻⁶). We also replicated a recently identified association signal near RNF114 (rs495337, combined P = 2 × 10⁻⁷).

Identification of copy number variation hotspots in human populations

Copy number variants (CNVs) in the human genome contribute to both Mendelian and complex traits as well as to genomic plasticity in evolution. The investigation of mutational rates of CNVs is critical to understanding genomic instability and the etiology of the copy number variation (CNV)-related traits. However, the evaluation of the CNV mutation rate at the genome level poses an insurmountable practical challenge that requires large samples and accurate typing. In this study, we show that an approximate estimation of the CNV mutation rate could be achieved by using the phylogeny information of flanking SNPs. This allows a genome-wide comparison of mutation rates between CNVs with the use of vast, readily available data of SNP genotyping. A total of 4187 CNV regions (CNVRs) previously identified in HapMap populations were investigated in this study. We showed that the mutation rates for the majority of these CNVRs are at the order of 10⁻⁵ per generation, consistent with experimental observations at individual loci. Notably, the mutation rates of 104 (2.5%) CNVRs were estimated at the order of 10⁻³ per generation; therefore, they were identified as potential hotspots. Additional analyses revealed that genome architecture at CNV loci has a potential role in inciting mutational hotspots in the human genome. Interestingly, 49 (47%) CNV hotspots include human genes, some of which are known to be functional CNV loci (e.g., CNVs of C4 and β-defensin causing autoimmune diseases and CNVs of HYDIN with implication in control of cerebral cortex size), implicating the important role of CNV in human health and evolution, especially in common and complex diseases.

Measurement methods and accuracy in copy number variation: failure to replicate associations of beta-defensin copy number with Crohn's disease

The copy number variation in beta-defensin genes on human chromosome 8 has been proposed to underlie susceptibility to inflammatory disorders, but presents considerable challenges for accurate typing on the scale required for adequately powered case-control studies. In this work, we have used accurate methods of copy number typing based on the paralogue ratio test (PRT) to assess beta-defensin copy number in more than 1500 UK DNA samples including more than 1000 cases of Crohn's disease. A subset of 625 samples was typed using both PRT-based methods and standard real-time PCR methods, from which direct comparisons highlight potentially serious shortcomings of a real-time PCR assay for typing this variant. Comparing our PRT-based results with two previous studies based only on real-time PCR, we find no evidence to support the reported association of Crohn's disease with either low or high beta-defensin copy number; furthermore, it is noteworthy that there are disagreements between different studies on the observed frequency distribution of copy number states among European controls. We suggest safeguards to be adopted in assessing and reporting the accuracy of copy number measurement, with particular emphasis on integer clustering of results, to avoid reporting of spurious associations in future case-control studies.

Copy number variation shapes genome diversity in Arabidopsis over immediate family generational scales

Arabidopsis thaliana is the model plant and is grown worldwide by plant biologists seeking to dissect the molecular underpinning of plant growth and development. Gene copy number variation (CNV) is a common form of genome natural diversity that is currently poorly studied in plants and may have broad implications for model organism research, evolutionary biology, and crop science. Herein, comparative genomic hybridization (CGH) was used to identify and interrogate regions of gene CNV across the A. thaliana genome. A common temperature condition used for growth of A. thaliana in our laboratory and many around the globe is 22 degrees C. The current study sought to test whether A. thaliana, grown under different temperature (16 and 28 degrees C) and stress regimes (salicylic acid spray) for five generations, selecting for fecundity at each generation, displayed any differences in CNV relative to a plant lineage growing under normal conditions. Three siblings from each alternative temperature or stress lineage were also compared with the reference genome (22 degrees C) by CGH to determine repetitive and nonrepetitive CNVs. Findings document exceptional rates of CNV in the genome of A. thaliana over immediate family generational scales. A propensity for duplication and nonrepetitive CNVs was documented in 28 degrees C CGH, which was correlated with the greatest plant stress and infers a potential CNV-environmental interaction. A broad diversity of gene species were observed within CNVs, but transposable elements and biotic stress response genes were notably overrepresented as a proportion of total genes and genes initiating CNVs. Results support a model whereby segmental CNV and the genes encoded within these regions contribute to adaptive capacity of plants through natural genome variation.

Defensin-barbed innate immunity: clinical associations in the pediatric population

Defensins and related antimicrobial peptides serve a central role in innate immunity in all species of plants and animals. In humans, defensins are widely expressed, including in neutrophils, skin, and mucosal epithelia. Most defensins are potent antibiotics, and some have chemotactic and toxin-neutralizing activities. Results of recent studies on the homeostatic and disease-fighting activities of human defensins point to a key relevance in several pediatric disorders. Inherited variation in defensin gene expression may contribute to susceptibility to several diseases, including psoriasis and Crohn disease. We review here the recent discoveries in innate immunity that shed light on the potential roles of defensins, and other antimicrobial molecules, in the pathophysiology of common pediatric diseases such as atopic dermatitis, necrotizing enterocolitis, cystic fibrosis, and otitis media.

The role of keratinocytes in defense against infection

PURPOSE OF REVIEW

The discovery of increasing numbers of epithelial antimicrobial peptides (AMPs), cytokines that specifically induce AMPs in epithelial cells, and mechanisms of its regulation point toward a central role of the keratinocyte as effector cell of the epithelial innate defense system. This review summarizes recent progress in understanding the keratinocyte's role in combating infection that can help to understand why skin is usually covered with microbes but normally not infected.

RECENT FINDINGS

The AMP LL-37 has been identified as regulator of keratinocyte apoptosis. The hypothesis of a direct defense function of keratinocytes, combating bacterial, fungal, virus and parasite infection, is strengthened. The discovery of an IL-22-producing T-lymphocyte subpopulation implicates a role in AMP induction of keratinocytes. Multiple studies are adding to our understanding of how skin keratinocytes are interacting in skin barrier defects and with the microflora. Although in atopic patients AMP production is not generally impaired, in hyperIgE syndrome a lack of Th17 cytokines causes local Staphylococcus aureus infection due to a defective keratinocyte defense system. Ultraviolet radiation induces AMPs and thus may have beneficial effects to combat skin infection.

SUMMARY

There is better understanding of the keratinocyte's role in the skin's innate defense system, and these data can help to generate therapeutics that can activate this defense system.

Haplotyping and copy number estimation of the highly polymorphic human beta-defensin locus on 8p23 by 454 amplicon sequencing

Background

The beta-defensin gene cluster (DEFB) at chromosome 8p23.1 is one of the most copy number (CN) variable regions of the human genome. Whereas individual DEFB CNs have been suggested as independent genetic risk factors for several diseases (e.g. psoriasis and Crohn's disease), the role of multisite sequence variations (MSV) is less well understood and to date has only been reported for prostate cancer. Simultaneous assessment of MSVs and CNs can be achieved by PCR, cloning and Sanger sequencing, however, these methods are labour and cost intensive as well as prone to methodological bias introduced by bacterial cloning. Here, we demonstrate that amplicon sequencing of pooled individual PCR products by the 454 technology allows in-depth determination of MSV haplotypes and estimation of DEFB CNs in parallel.

Results

Six PCR products spread over ~87 kb of DEFB and harbouring 24 known MSVs were amplified from 11 DNA samples, pooled and sequenced on a Roche 454 GS FLX sequencer. From ~142,000 reads, ~120,000 haplotype calls (HC) were inferred that identified 22 haplotypes ranging from 2 to 7 per amplicon. In addition to the 24 known MSVs, two additional sequence variations were detected. Minimal CNs were estimated from the ratio of HCs and compared to absolute CNs determined by alternative methods. Concordance in CNs was found for 7 samples, the CNs differed by one in 2 samples and the estimated minimal CN was half of the absolute in one sample. For 7 samples and 2 amplicons, the 454 haplotyping results were compared to those by cloning/Sanger sequencing. Intrinsic problems related to chimera formation during PCR and differences between haplotyping by 454 and cloning/Sanger sequencing are discussed.

Conclusion

Deep amplicon sequencing using the 454 technology yield thousands of HCs per amplicon for an affordable price and may represent an effective method for parallel haplotyping and CN estimation in small to medium-sized cohorts. The obtained haplotypes represent a valuable resource to facilitate further studies of the biomedical impact of highly CN variable loci such as the beta-defensin locus.

The spatial expression pattern of antimicrobial peptides across the healthy bovine udder

Antimicrobial peptides are key molecules in local host defense. With the aim to better understand the possible involvement of these peptides in the prevention of bovine mastitis, we determined, for the first time to our knowledge, the spatial pattern of constitutive expression of 5 bovine beta-defensins and bovine psoriasin (S100A7) across 5 localizations of the bovine mammary gland applying a quantitative real-time PCR approach. We observed 3 different expression patterns in the healthy udder: 1) constitutive expression of the lingual and tracheal antimicrobial peptides (LAP and TAP), as well that of bovine neutrophil beta-defensins 4 and 10 (BNBD4 and BNBD10), is essentially restricted to the mammary lymph node; 2) bovine beta-defensin 1 (DEFB1) is mainly expressed in the cisternal epithelium and the Rosette of Fürstenberg; 3) strong constitutive mRNA expression of the calcium-binding protein S100A7, which is also known as psoriasin and which has been reported to be highly active against Escherichia coli, was detected in the streak canal. These results indicate a crucial role of S100A7 in the early-stage prevention of coliform mastitis, and the analyzed beta-defensins might be regarded as inducible weapons against already invading pathogens.

8p23.1 duplication syndrome differentiated from copy number variation of the defensin cluster at prenatal diagnosis in four new families

BACKGROUND

The 8p23.1 duplication syndrome and copy number variation of the 8p23.1 defensin gene cluster are cytogenetically indistinguishable but distinct at the molecular level. To our knowledge, the 8p23.1 duplication syndrome has been described at prenatal diagnosis only once and we report our experience with four further apparent duplications ascertained at prenatal diagnosis.

METHODS

Additional material at band 8p23.1 was detected using conventional G-banded cytogenetics in each case. Multiplex Ligation-dependent Probe Amplification (MLPA) or Fluorescence In Situ Hybridisation (FISH) were used depending on whether only DNA (Cases 1 and 4) or cytogenetic preparations (Cases 2 and 3) were available from the laboratory of origin. The extent of the duplication in Case 1 was retrospectively determined using array Comparative Genomic Hybridisation (array CGH).

RESULTS

Three cases of 8p23.1 duplication syndrome were found (Cases 1 to 3). Two were de novo and continued to term and the third, a paternally transmitted duplication, was terminated because of a previous child with psychomotor delay and 8p23.1 duplication syndrome. Case 1 was ascertained with a hypoplastic left heart but the ventricular septal and interventricular defects, in Cases 2 and 3 respectively, were found after ascertainment for advanced maternal age. By contrast, case 4 was a maternally transmitted copy number variation of the defensin cluster with normal outcome.

CONCLUSIONS

Our data underline the need to differentiate 8p23.1 duplications from copy number variation of the defensin cluster using FISH, MLPA or array CGH. Cardiac defects were ascertained by ultrasound in only one of the three duplication 8p23.1 pregnancies but were visible in two of the three at 21 to 22 weeks gestation. Our results provide further evidence that both deletion and duplication of the GATA4 transcription factor can give rise to a variety of conotruncal heart defects with variable penetrance and expressivity.

Manual annotation and analysis of the defensin gene cluster in the C57BL/6J mouse reference genome

Background

Host defense peptides are a critical component of the innate immune system. Human alpha- and beta-defensin genes are subject to copy number variation (CNV) and historically the organization of mouse alpha-defensin genes has been poorly defined. Here we present the first full manual genomic annotation of the mouse defensin region on Chromosome 8 of the reference strain C57BL/6J, and the analysis of the orthologous regions of the human and rat genomes. Problems were identified with the reference assemblies of all three genomes. Defensins have been studied for over two decades and their naming has become a critical issue due to incorrect identification of defensin genes derived from different mouse strains and the duplicated nature of this region.

Results

The defensin gene cluster region on mouse Chromosome 8 A2 contains 98 gene loci: 53 are likely active defensin genes and 22 defensin pseudogenes. Several TATA box motifs were found for human and mouse defensin genes that likely impact gene expression. Three novel defensin genes belonging to the Cryptdin Related Sequences (CRS) family were identified. All additional mouse defensin loci on Chromosomes 1, 2 and 14 were annotated and unusual splice variants identified. Comparison of the mouse alpha-defensins in the three main mouse reference gene sets Ensembl, Mouse Genome Informatics (MGI), and NCBI RefSeq reveals significant inconsistencies in annotation and nomenclature. We are collaborating with the Mouse Genome Nomenclature Committee (MGNC) to establish a standardized naming scheme for alpha-defensins.

Conclusions

Prior to this analysis, there was no reliable reference gene set available for the mouse strain C57BL/6J defensin genes, demonstrating that manual intervention is still critical for the annotation of complex gene families and heavily duplicated regions. Accurate gene annotation is facilitated by the annotation of pseudogenes and regulatory elements. Manually curated gene models will be incorporated into the Ensembl and Consensus Coding Sequence (CCDS) reference sets. Elucidation of the genomic structure of this complex gene cluster on the mouse reference sequence, and adoption of a clear and unambiguous naming scheme, will provide a valuable tool to support studies on the evolution, regulatory mechanisms and biological functions of defensins in vivo.

Nucleotide, cytogenetic and expression impact of the human chromosome 8p23.1 inversion polymorphism

Background

The human chromosome 8p23.1 region contains a 3.8–4.5 Mb segment which can be found in different orientations (defined as genomic inversion) among individuals. The identification of single nucleotide polymorphisms (SNPs) tightly linked to the genomic orientation of a given region should be useful to indirectly evaluate the genotypes of large genomic orientations in the individuals.

Results

We have identified 16 SNPs, which are in linkage disequilibrium (LD) with the 8p23.1 inversion as detected by fluorescent in situ hybridization (FISH). The variability of the 8p23.1 orientation in 150 HapMap samples was predicted using this set of SNPs and was verified by FISH in a subset of samples. Four genes (NEIL2, MSRA, CTSB and BLK) were found differentially expressed (p<0.0005) according to the orientation of the 8p23.1 region. Finally, we have found variable levels of mosaicism for the orientation of the 8p23.1 as determined by FISH.

Conclusion

By means of dense SNP genotyping of the region, haplotype-based computational analyses and FISH experiments we could infer and verify the orientation status of alleles in the 8p23.1 region by detecting two short haplotype stretches at both ends of the inverted region, which are likely the relic of the chromosome in which the original inversion occurred. Moreover, an impact of 8p23.1 inversion on gene expression levels cannot be ruled out, since four genes from this region have statistically significant different expression levels depending on the inversion status. FISH results in lymphoblastoid cell lines suggest the presence of mosaicism regarding the 8p23.1 inversion.

Dysregulation of human beta-defensin-2 protein in inflammatory bowel disease

Background

Human β-defensin-2 (HBD2) is an antimicrobial peptide implicated in the pathogenesis of inflammatory bowel disease (IBD). Low copy number and concomitant low mRNA expression of the HBD2 gene have been implicated in susceptibility to colonic Crohn's Disease (CD). We investigated the colonic distribution of HBD2 mRNA expression, and the contributions of genetic and environmental factors on HBD2 protein production.

Methodology/Principal Findings

We examined HBD2 mRNA expression at three colonic locations by microarray analysis of biopsies from 151 patients (53 CD, 67 ulcerative colitis [UC], 31 controls). We investigated environmental and genetic influences on HBD2 protein production using ex vivo cultured sigmoid colon biopsies from 69 patients (22 CD, 26 UC, 21 controls) stimulated with lipopolysaccharide (LPS) and/or nicotine for 24 hours. HBD2 and cytokines were measured in culture supernatants. Using DNA samples from these patients, regions in the HBD2 gene promoter were sequenced for NF-κB binding-sites and HBD2 gene copy number was determined. HBD2 mRNA expression was highest in inflamed (vs. uninflamed p = 0.0122) ascending colon in CD and in inflamed (vs. uninflamed p<0.0001) sigmoid colon in UC. HBD2 protein production was increased in inflamed UC biopsies (p = 0.0078). There was no difference in HBD2 protein production from unstimulated biopsies of CD, UC and controls. LPS-induced HBD2 production was significantly increased in CD (p = 0.0375) but not UC (p = 0.2017); this LPS-induced response was augmented by nicotine in UC (p = 0.0308) but not CD (p = 0.6872). Nicotine alone did not affect HBD2 production. HBD2 production correlated with IL8 production in UC (p<0.001) and with IL10 in CD (p<0.05). Variations in the HBD2 promoter and HBD2 gene copy number did not affect HBD2 production.

Significance/Conclusions

Colonic HBD2 was dysregulated at mRNA and protein level in IBD. Inflammatory status and stimulus but not germline variations influenced these changes.

Copy number variants (CNVs) in primate species using array-based comparative genomic hybridization

A substantial amount of genomic variation is now known to exist in humans and other primate species. Single nucleotide polymorphisms (SNPs) are thought to represent the vast majority of genomic differences among individuals in a given primate species and comprise about 0.1% of the genomes of two humans. However, recent studies have now shown that structural variation msay account for as much as 0.7% of the genomic differences in humans, of which copy number variants (CNVs) are the largest component. CNVs are segments of DNA that can range in size from hundreds of bases to millions of base pairs in length and have different number of copies between individuals. Recent technological advancements in array technologies led to genome-wide identification of CNVs and consequently revealed thousands of variable loci in humans, comprising as much as 12% of the human genome [A.J. Iafrate, L. Feuk, M.N. Rivera, M.L. Listewnik, P.K. Donahoe, Y. Qi, S.W. Scherer, C. Lee, Nat. Genet. 36 (2004) 949-951, [3]]. CNVs in humans have already been associated with susceptibility to certain complex diseases, dietary adaptation, and several neurological conditions. In addition, recent studies have shown that CNVs can be successfully implemented in population genetics research, providing important insights into human genetic variation. Nevertheless, the important role of CNVs in primate evolution and genetic diversity is still largely unknown. This article aims to outline the strengths and weaknesses of current comparative genomic hybridization array technologies that have been employed to detect CNV variation and the applications of these techniques to primate genetic research.