Identification of stock units of Mullus barbatus (Linnaeus, 1758) in the Italian seas: development of genomic 2b-RAD markers and biological variation analysis.

This thesis is settled within the STOCKMAPPING project, which represents one of the studies that were developed in the framework of RITMARE Flagship project. The main goals of STOCKMAPPING were the creation of a genomic mapping for stocks of demersal target species and the assembling of a database of population genomic, in order to identify stocks and stocks boundaries. The thesis focuses on three main objectives representing the core for the initial assessment of the methodologies and structure that would be applied to the entire STOCKMAPPING project: individuation of an analytical design to identify and locate stocks and stocks boundaries of Mullus barbatus, application of a multidisciplinary approach to validate biological methods and an initial assessment and improvement for the genotyping by sequencing technique utilized (2b-RAD). The first step is the individuation of an analytical design that has to take in to account the biological characteristics of red mullet and being representative for STOCKMAPPING commitments. In this framework a reduction and selection steps was needed due to budget reduction. Sampling areas were ranked according the individuation of four priorities. To guarantee a multidisciplinary approach the biological data associated to the collected samples were used to investigate differences between sampling areas and GSAs. Genomic techniques were applied to red mullet for the first time so an initial assessment of molecular protocols for DNA extraction and 2b-RAD processing were needed. At the end 192 good quality DNAs have been extracted and eight samples have been processed with 2b-RAD. Utilizing the software Stacks for sequences analyses a great number of SNPs markers among the eight samples have been identified. Several tests have been performed changing the main parameter of the Stacks pipeline in order to identify the most explicative and functional sets of parameters.

Stock boundaries of Parapenaeus longirostris (Lucas, 1846) populations in the Italian seas: Development of genomic 2b-RAD markers and biological variation analysis.

This thesis is developed in the contest of Ritmare project WP1, which main objective is the development of a sustainable fishery through the identification of populations boundaries in commercially important species in Italian Seas. Three main objectives are discussed in order to help reach the main purpose of identification of stock boundaries in Parapenaeus longirostris: 1 -Development of a representative sampling design for Italian seas; 2 -Evaluation of 2b-RAD protocol; 3 -Investigation of populations through biological data analysis. First of all we defined and accomplished a sampling design which properly represents all Italian seas. Then we used information and data about nursery areas distribution, abundance of populations and importance of P. longirostris in local fishery, to develop an experimental design that prioritize the most important areas to maximize the results with actual project funds. We introduced for the first time the use of 2b-RAD on this species, a genotyping method based on sequencing the uniform fragments produced by type IIB restriction endonucleases. Thanks to this method we were able to move from genetics to the more complex genomics. In order to proceed with 2b-RAD we performed several tests to identify the best DNA extraction kit and protocol and finally we were able to extract 192 high quality DNA extracts ready to be processed. We tested 2b-RAD with five samples and after high-throughput sequencing of libraries we used the software “Stacks” to analyze the sequences. We obtained positive results identifying a great number of SNP markers among the five samples. To guarantee a multidisciplinary approach we used the biological data associated to the collected samples to investigate differences between geographical samples. Such approach assures continuity with other project, for instance STOCKMED, which utilize a combination of molecular and biological analysis as well.

Gespeicherte Ionen in der Paulfalle: Voruntersuchungen zur Lebensdauermessung des metastabilen 6D-3/2-Zustands des 226 Ra +, Speichereigenschaften von 138 Ba + und 28 N 2 + bei Variation des Puffergasdrucks

Auf Paulfallen basierende Experimente spielen eine wichtige Rolle in verschiedenen Bereichen der Physik, z.B. der Atomphysik zum Test theoretischer Modelle und der Massenspektroskopie. Die vorliegende Arbeit widmet sich beiden Themengebieten und gliedert sich entsprechend in zwei Teilbereiche: 1) Erdalkali-Ionen sind aufgrund ihrer Energieniveaus optimale Kandidaten für Laserspektroskopie-Experimente mit Ionenfallen und bestens geeignet, um mittels der spektroskopischen Daten die theoretischen Modelle zu testen. Lediglich für Ra+ fehlen bislang als einzigem Erdalkali-Ion diese Daten wie z.B. die Lebensdauern der metastabilen Niveaus. Diese wären auch von Interesse für bereits geplante Radium-Experimente zur Paritätsverletzung. Im ersten Teil dieser Arbeit wird der Aufbau eines Laser-Paulfallenexperiments zur Messung der Lebensdauer des 6D3/2 Zustands von 226Ra+ dokumentiert und es werden Testmessungen mit 138Ba+ vorgestellt. 2) Für die Verwendung der Paulfalle in der Massenspektroskopie und zur Analyse von Reaktionsprodukten ist die Kenntnis der Lage der im Speicherbereich auftretenden nichtlinearen Resonanzen wesentlich, ebenso wie deren Veränderung durch Dämpfung und Raumladung. Im zweiten Teil dieser Arbeit werden detaillierte Untersuchungen der Speicherung großer puffergasgekühlter Ionenwolken an zwei unterschiedlichen Paulfallen-Experimenten vorgestellt. Am ersten wurden 138Ba+-Ionenwolken kontinuierlich durch Laserspektroskopie bzw. über einen elektronischen Nachweis beobachtet, während das zweite N2+-Molekülionen automatisiert destruktiv nachwies. Am N2+-Experiment wurden zwei hochaufgelöste Messungen des ersten Speicherbereichs durchgeführt, die erstmals eine direkte Überprüfung der theoretisch berechneten Verläufe der Resonanzen mit experimentellen Daten erlauben. Die Nachweiseichung ermöglichte dabei zum ersten Mal die Angabe absoluter Ionenzahlen. Im Gegensatz zu vergleichbaren früheren Messungen wurden hierbei die sich überlagernden Speicherbereiche von 4 simultan gespeicherten Ionensorten beobachtet und zur Analyse der Resonanzen herangezogen. Die nichtlinearen Resonanzen wurden untersucht bei Variation von Puffergasdruck und Ionenzahl, wobei kollektive Resonanzen ohne zusätzliche externe Anregung beobachtet wurden. Die gemessenen Raumladungsverschiebungen wurden mit theoretischen Modellen verglichen. Bei Variation des Puffergasdrucks wurde mit Bariumionen die räumliche Ausdehnung der Ionenwolke gemessen und mit Stickstoffionen die Verschiebung des Punktes optimaler Speicherung bestimmt. Dabei wurde festgestellt, daß der zum Ioneneinfang optimale Puffergasdruck kleiner ist als der, bei dem die längsten Speicherdauern erzielt werden. Mit gespeicherten N2+-Ionen wurde die Position extern angeregter kollektiver und individueller Resonanzen im Frequenzspektrum bei Änderung der Parameter Ionenzahl, Puffergasdruck und Dauer der Anregung untersucht, ebenso wie die Resonanzform, die mit theoretischen Linienformen verglichen wurde. Bei Änderung der Fallenparameter wurden verstärkende Effekte zwischen nahen kollektiven Resonanzen festgestellt. Die Beobachtung, welche der im Frequenzspektrum vorher identifizierten Bewegungs-Resonanzen sich bei Variation der Fallenparameter a bzw. q überlagern, ermöglicht eine bislang nicht beschriebene einfache Methode der Bestimmung von nichtlinearen Resonanzen im Stabilitätsdiagramm.

Analysis of Copy Number Variants identifies new candidate genes for Autism Spectrum Disorder and Intellectual Disability

Autism spectrum disorder (ASD) and Intellectual Disability (ID) are complex neuropsychiatric disorders characterized by extensive clinical and genetic heterogeneity and with overlapping risk factors. The aim of my project was to further investigate the role of Copy Numbers Variants (CNVs), identified through genome-wide studies performed by the Autism Geome Project (AGP) and the CHERISH consortium in large cohorts of ASD and ID cases, respectively. Specifically, I focused on four rare genic CNVs, selected on the basis of their impact on interesting ASD/ID candidate genes: a) a compound heterozygous deletion involving CTNNA3, predicted to cause the lack of functional protein; b) a 15q13.3 duplication containing CHRNA7; c) a 2q31.1 microdeletion encompassing KLHL23, SSB and METTL5; d) Lastly, I investigated the putative imprinting regulation of the CADPS2 gene, disrupted by a maternal deletion in two siblings with ASD and ID. This study provides further evidence for the role of CTNNA3, CHRNA7, KLHL23 and CADPS2 as ASD and/or ID susceptibility genes, and highlights that rare genetic variation contributes to disease risk in different ways: some rare mutations, such as those impacting CTNNA3, act in a recessive mode of inheritance, while other CNVs, such as those occurring in the 15q13.3 region, are implicated in multiple developmental and/or neurological disorders possibly interacting with other susceptibility variants elsewhere in the genome. On the other hand, the discovery of a tissue-specific monoallelic expression for the CADPS2 gene, implicates the involvement of epigenetic regulatory mechanisms as risk factors conferring susceptibility to ASD/ID.

Temporal variation of genetic diversity and differentiation of albacore, Thunnus alalunga, in the Mediterranean Sea

This study is on albacore (Thunnus alalunga, Bonnaterre 1788), an epi- and mesopelagic oceanic tuna species cosmopolitan in the tropical and temperate waters of all oceans including the Mediterranean Sea, extending in a broad band between 40°N and 40°S. What it’s known about albacore population structure is based on different studies that used fisheries data, RFLP, mtDNA control region and nuDNA markers, blood lectins analysis, individual tags and microsatellite. At the moment, for T. alalunga six management units are recognized: the North Pacific, South Pacific, Indian, North Atlantic, South Atlantic and Mediterranean stocks. In this study I have done a temporal and spatial comparison of genetic variability between different Mediterranean populations of Thunnus alalunga matching an historical dataset ca. from 1920s composed of 43 individuals divided in 3 populations (NADR, SPAIN and CMED) with a modern dataset composed of 254 individuals and 7 populations (BAL, CYP, LIG, TYR, TUR, ADR, ALB). The investigation was possible using a panel of 94 nuclear SNPs, built specifically for the target species at the University of Basque Country UPV/EHU. First analysis done was the Hardy-Weinberg, then the number of clusters (K) was determined using STRUCTURE and to assess the genetic variability, allele frequencies, the average number of alleles per locus, expected (He) and observed (Ho) heterozygosis, and the index of polymorphism (P) was used the software Genetix. Historical and modern samples gives different results, showing a clear loss of genetic diversity over time leading to a single cluster in modern albacore instead of the two found in historical samples. What this study reveals is very important for conservation concerns, and additional research endeavours are needed.

Are there genetic breaks between Atlantic and Pacific Yellowfin tuna (Thunnus albacares) populations? A preliminary study based on microsatellites gene variation

Yellowfin tuna (Thunnus albacares, YFT, Bonnaterre 1788) is one of the most important market tuna species in the world. The high mortality of juveniles is in part caused by their bycatch. Indeed, if unregulated, it could permanently destabilize stocks health. For this reason investigating and better knowing the stock boundaries represent a crucial concern. Aim of this thesis was to preliminary investigate the YFT population structure within and between Atlantic and Pacific Oceans through the analysis of genetic variation at eight microsatellite loci and assess the occurrence of barriers to the gene flow between Oceans. For this propouse we collected 4 geographical samples coming from Atlantic and Pacific Ocean and selected a panel of 8 microsatellites loci developped by Antoni et al., (2014). Samples 71-2-Y and 77-2-Y, came from rispectively west central pacific ocean (WCPO) and east central pacific ocean (ECPO), instead samples 41-1-Y and 34-2-Y derive from west central atlantic ocean (WCAO) and east central atlantic ocean (ECAO). Total 160 specimens were analyzed (40 per sample) and were carried out several genetic information as allele frequencies, allele number, allelic richness, HWE (using He and Ho) and pairwise Fst genetic distance. Results obtained, may support the panmictic theory of this species, only one of pairwise Fst obtained is statistically significant (Fst= 0.00927; pV= 0.00218) between 41-1-Y and 71-2-Y samples. Results suggest low genetic differentiation and consequent high level of gene flow between Atlantic and Pacific populations. Furthermore, we performed an analysis of molecular taxonomy through the use of ATCO (the flaking region between ATPse6 and cytochrome oxidase subunit III genes mt DNA, to discriminate within the gener Thunnus two of the related species (Yellofin and bigeye tuna) according with their difficult recognition at certain size (<40 cm). ATCO analysis in this thesis, has provided strong discriminate evidence between the target species proving to be one of the most reliable genetic tools capable to indagate within the genus Thunnus. Thus, our study has provided useful information for possible use of this protocol for conservation plans and management of this fish stocks.

First assessment on genetic structure and phylogeography of Mediterranean blue shark (prionace glauca, L. 1758) population using mitochondrial gene variation: a comparison with the Atlantic.

The blue shark, Prionace glauca, is one of the most vagile shark species worldwide distributed. The particular body shape allows blue sharks make transoceanic movements, leading to a circumglobal distribution. Due to its reproductive cycle, an extraordinarily high number of specimens is globally registered but, even if it is still a major bycatch of longline fishery rather than a commercial target, it is characterized by a high vulnerability. In this perspective it is important to increase the amount of informations regarding its population extent in the different worldwide areas, evaluating the possible phylogeographic patterns between different locations. This study, included in the "MedBlueSGen" European project, aims exactly at filling a gap in knowledges regarding the genetic population structure of the Mediterranean blue sharks, which has never been investigated before, with a comparison with the North-Eastern Atlantic blue shark population. To reach this objective, we used a dataset of samples from different Mediterranean areas implementing it with some samples from North-Eastern Atlantic. Analyzing the variability of the two mitochondrial markers control region and cytochrome b, with the design of new species-specific primer pairs, we assessed the mitochondrial genetic structure of Mediterranean and North-Eastern Atlantic samples, focusing on the analysis of their possible connectivity, and we tried to reconstruct their demographic history and population size. Data analyses highlighted the absence of a genetic structuring within the Mediterranean and among it and North-Eastern Atlantic, suggesting that the Strait of Gibraltar doesn't represent a phylogeographic barrier. These results are coherent to what has been found in similar investigations on other worldwide blue shark populations. Analysis of the historical demographic trend revealed a general stable pattern for the cytochrome-b and a slightly population expansion for the control region marker.

Hybridization between distant lineages increases adaptive variation during a biological invasion: stickleback in Switzerland

The three-spined stickleback is a widespread Holarctic species complex that radiated from the sea into freshwaters after the retreat of the Pleistocene ice sheets. In Switzerland, sticklebacks were absent with the exception of the far northwest, but different introduced populations have expanded to occupy a wide range of habitats since the late 19th century. A well-studied adaptive phenotypic trait in sticklebacks is the number of lateral plates. With few exceptions, freshwater and marine populations in Europe are fixed for either the low plated phenotype or the fully plated phenotype, respectively. Switzerland, in contrast, harbours in close proximity the full range of phenotypic variation known from across the continent. We addressed the phylogeographic origins of Swiss sticklebacks using mitochondrial partial cytochrome b and control region sequences. We found only five different haplotypes but these originated from three distinct European regions, fixed for different plate phenotypes. These lineages occur largely in isolation at opposite ends of Switzerland, but co-occur in a large central part. Across the country, we found a strong correlation between a microsatellite linked to the high plate ectodysplasin allele and the mitochondrial haplotype from a region where the fully plated phenotype is fixed. Phylogenomic and population genomic analysis of 481 polymorphic amplified fragment length polymorphism loci indicate genetic admixture in the central part of the country. The same part of the country also carries elevated within-population phenotypic variation. We conclude that during the recent invasive range expansion of sticklebacks in Switzerland, adaptive and neutral between-population genetic variation was converted into within-population variation, raising the possibility that hybridization between colonizing lineages contributed to the ecological success of sticklebacks in Switzerland.

Variation in partner notification outcomes for chlamydia in UK genitourinary medicine clinics: multilevel study

Objectives To compare different ways of measuring partner notification (PN) outcomes with published audit standards, examine variability between clinics and examine factors contributing to variation in PN outcomes in genitourinary medicine (GUM) clinics in the UK. Methods Reanalysis of the 2007 BASHH national chlamydia audit. The primary outcome was the number of partners per index case tested for chlamydia, as verified by a healthcare worker or, if missing, reported by the patient. Control charts were used to examine variation between clinics considering missing values as zero or excluding missing values. Hierarchical logistic regression was used to investigate factors contributing to variation in outcomes. Results Data from 4616 individuals in 169 genitourinary medicine clinics were analysed. There was no information about the primary outcome in 41% of records. The mean number of partners tested for chlamydia ranged from 0 to 1.5 per index case per clinic. The median across all clinics was 0.47 when missing values were assumed to be zero and 0.92 per index case when missing values were excluded. Men who have sex with men were less likely than heterosexual men and patients with symptoms (4-week look-back period) were less likely than asymptomatic patients (6-month look-back) to report having one or more partners tested for chlamydia. There was no association between the primary outcome and the type of the health professional giving the PN advice. Conclusions The completeness of PN outcomes recorded in clinical notes needs to improve. Further research is needed to identify auditable measures that are associated with successful PN that prevents repeated chlamydia in index cases.

Estimating the net contribution of interleukin-28B variation to spontaneous hepatitis C virus clearance

The identification of associations between interleukin-28B (IL-28B) variants and the spontaneous clearance of hepatitis C virus (HCV) raises the issues of causality and the net contribution of host genetics to the trait. To estimate more precisely the net effect of IL-28B genetic variation on HCV clearance, we optimized genotyping and compared the host contributions in multiple- and single-source cohorts to control for viral and demographic effects. The analysis included individuals with chronic or spontaneously cleared HCV infections from a multiple-source cohort (n = 389) and a single-source cohort (n = 71). We performed detailed genotyping in the coding region of IL-28B and searched for copy number variations to identify the genetic variant or haplotype carrying the strongest association with viral clearance. This analysis was used to compare the effects of IL-28B variation in the two cohorts. Haplotypes characterized by carriage of the major alleles at IL-28B single-nucleotide polymorphisms (SNPs) were highly overrepresented in individuals with spontaneous clearance versus those with chronic HCV infections (66.1% versus 38.6%, P = 6 × 10(-9) ). The odds ratios for clearance were 2.1 [95% confidence interval (CI) = 1.6-3.0] and 3.9 (95% CI = 1.5-10.2) in the multiple- and single-source cohorts, respectively. Protective haplotypes were in perfect linkage (r(2) = 1.0) with a nonsynonymous coding variant (rs8103142). Copy number variants were not detected. CONCLUSION: We identified IL-28B haplotypes highly predictive of spontaneous HCV clearance. The high linkage disequilibrium between IL-28B SNPs indicates that association studies need to be complemented by functional experiments to identify single causal variants. The point estimate for the genetic effect was higher in the single-source cohort, which was used to effectively control for viral diversity, sex, and coinfections and, therefore, offered a precise estimate of the net host genetic contribution.

Genetic maps of the proximal half of chromosome arm 2L of Drosophila melanogaster

Drosophila mutants have played an important role in elucidating the physiologic function of genes. Large-scale projects have succeeded in producing mutations in a large proportion of Drosophila genes. Many mutant fly lines have also been produced through the efforts of individual laboratories over the past century. In an effort to make some of these mutants more useful to the research community, we systematically mapped a large number of mutations affecting genes in the proximal half of chromosome arm 2L to more precisely defined regions, defined by deficiency intervals, and, when possible, by individual complementation groups. To further analyze regions 36 and 39-40, we produced 11 new deficiencies with gamma irradiation, and we constructed 6 new deficiencies in region 30-33, using the DrosDel system. trans-heterozygous combinations of deficiencies revealed 5 additional functions, essential for viability or fertility.

Losses of chromosome arms 4q, 8p, 13q and gain of 8q are correlated with increasing chromosomal instability in hepatocellular carcinoma

OBJECTIVE: Chromosomal instability is a key feature in hepatocellular carcinoma (HCC). Array comparative genomic hybridization (aCGH) revealed recurring structural aberrations, whereas fluorescence in situ hybridization (FISH) indicated an increasing number of numerical aberrations in dedifferentiating HCC. Therefore, we examined whether there was a correlation between structural and numerical aberrations of chromosomal instability in HCC. METHODS AND RESULTS: 27 HCC (5 well, 10 moderately, 12 lower differentiated) already cytogenetically characterized by aCGH were analyzed. FISH analysis using probes for chromosomes 1, 3, 7, 8 and 17 revealed 1.46-4.24 signals/nucleus, which correlated with the histological grade (well vs. moderately,p < 0.0003; moderately vs. lower, p < 0.004). The number of chromosomes to each other was stable with exceptions only seen for chromosome 8. Loss of 4q and 13q, respectively, were correlated with the number of aberrations detected by aCGH (p < 0.001, p < 0.005; Mann-Whitney test). Loss of 4q and gain of 8q were correlated with an increasing number of numerical aberrations detected by FISH (p < 0.020, p < 0.031). Loss of 8p was correlated with the number of structural imbalances seen in aCGH (p < 0.048), but not with the number of numerical changes seen in FISH. CONCLUSION: We found that losses of 4q, 8p and 13q were closely correlated with an increasing number of aberrations detected by aCGH, whereas a loss of 4q and a gain of 8q were also observed in the context of polyploidization, the cytogenetic correlate of morphological dedifferentiation.

Space trajectories optimization using variable-chromosome-length genetic algorithms

The problem of optimal design of a multi-gravity-assist space trajectories, with free number of deep space maneuvers (MGADSM) poses multi-modal cost functions. In the general form of the problem, the number of design variables is solution dependent. To handle global optimization problems where the number of design variables varies from one solution to another, two novel genetic-based techniques are introduced: hidden genes genetic algorithm (HGGA) and dynamic-size multiple population genetic algorithm (DSMPGA). In HGGA, a fixed length for the design variables is assigned for all solutions. Independent variables of each solution are divided into effective and ineffective (hidden) genes. Hidden genes are excluded in cost function evaluations. Full-length solutions undergo standard genetic operations. In DSMPGA, sub-populations of fixed size design spaces are randomly initialized. Standard genetic operations are carried out for a stage of generations. A new population is then created by reproduction from all members based on their relative fitness. The resulting sub-populations have different sizes from their initial sizes. The process repeats, leading to increasing the size of sub-populations of more fit solutions. Both techniques are applied to several MGADSM problems. They have the capability to determine the number of swing-bys, the planets to swing by, launch and arrival dates, and the number of deep space maneuvers as well as their locations, magnitudes, and directions in an optimal sense. The results show that solutions obtained using the developed tools match known solutions for complex case studies. The HGGA is also used to obtain the asteroids sequence and the mission structure in the global trajectory optimization competition (GTOC) problem. As an application of GA optimization to Earth orbits, the problem of visiting a set of ground sites within a constrained time frame is solved. The J2 perturbation and zonal coverage are considered to design repeated Sun-synchronous orbits. Finally, a new set of orbits, the repeated shadow track orbits (RSTO), is introduced. The orbit parameters are optimized such that the shadow of a spacecraft on the Earth visits the same locations periodically every desired number of days.

Pollinator responses to variation in floral display and flower size in dioecious Sagittaria latifolia (Alismataceae)

In animal-pollinated plants with unisexual flowers, sexual dimorphism in floral traits may be the consequence of pollinator-mediated selection. Experimental investigations of the effects of variation in flower size and floral display on pollinator visitation can provide insights into the evolution of floral dimorphism in dioecious plants. Here, we investigated pollinator responses to experimental arrays of dioecious Sagittaria latifolia in which we manipulated floral display and flower size. We also examined whether there were changes in pollinator visitation with increasing dimorphism in flower size. In S. latifolia, males have larger flowers and smaller floral displays than females. Visitation by pollinators, mainly flies and bees, was more frequent for male than for female inflorescences and increased with increasing flower size, regardless of sex. The number of insect visits per flower decreased with increasing floral display in males but remained constant in females. Greater sexual dimorphism in flower size increased visits to male inflorescences but had no influence on the number of visits to female inflorescences. These results suggest that larger flower sizes would be advantageous to both females and males, and no evidence was found that females suffer from increased flower-size dimorphism. Small daily floral displays may benefit males by allowing extended flowering periods and greater opportunities for effective pollen dispersal.

Comparison of multiplex ligation-dependent probe amplification and real-time PCR accuracy for gene copy number quantification using the beta-defensin locus

The reliable quantification of gene copy number variations is a precondition for future investigations regarding their functional relevance. To date, there is no generally accepted gold standard method for copy number quantification, and methods in current use have given inconsistent results in selected cohorts. In this study, we compare two methods for copy number quantification. beta-defensin gene copy numbers were determined in parallel in 80 genomic DNA samples by real-time PCR and multiplex ligation-dependent probe amplification (MLPA). The pyrosequencing-based paralog ratio test (PPRT) was used as a standard of comparison in 79 out of 80 samples. Realtime PCR and MPLA results confirmed concordant DEFB4, DEFB103A, and DEFB104A copy numbers within samples. These two methods showed identical results in 32 out of 80 samples; 29 of these 32 samples comprised four or fewer copies. The coefficient of variation of MLPA is lower compared with PCR. In addition, the consistency between MLPA and PPRT is higher than either PCR/MLPA or PCR/PPRT consistency. In summary, these results suggest that MLPA is superior to real-time PCR in beta-defensin copy number quantification.