Further delineation of the 15q13 microdeletion and duplication syndromes: a clinical spectrum varying from non-pathogenic to a severe outcome

Background: Recurrent 15q13.3 microdeletions were recently identified with identical proximal (BP4) and distal (BP5) breakpoints and associated with mild to moderate mental retardation and epilepsy. Methods: To assess further the clinical implications of this novel 15q13.3 microdeletion syndrome, 18 new probands with a deletion were molecularly and clinically characterised. In addition, we evaluated the characteristics of a family with a more proximal deletion between BP3 and BP4. Finally, four patients with a duplication in the BP3-BP4-BP5 region were included in this study to ascertain the clinical significance of duplications in this region. Results: The 15q13.3 microdeletion in our series was associated with a highly variable intra-and inter-familial phenotype. At least 11 of the 18 deletions identified were inherited. Moreover, 7 of 10 siblings from four different families also had this deletion: one had a mild developmental delay, four had only learning problems during childhood, but functioned well in daily life as adults, whereas the other two had no learning problems at all. In contrast to previous findings, seizures were not a common feature in our series (only 2 of 17 living probands). Three patients with deletions had cardiac defects and deletion of the KLF13 gene, located in the critical region, may contribute to these abnormalities. The limited data from the single family with the more proximal BP3-BP4 deletion suggest this deletion may have little clinical significance. Patients with duplications of the BP3-BP4-BP5 region did not share a recognisable phenotype, but psychiatric disease was noted in 2 of 4 patients. Conclusions: Overall, our findings broaden the phenotypic spectrum associated with 15q13.3 deletions and suggest that, in some individuals, deletion of 15q13.3 is not sufficient to cause disease. The existence of microdeletion syndromes, associated with an unpredictable and variable phenotypic outcome, will pose the clinician with diagnostic difficulties and challenge the commonly used paradigm in the diagnostic setting that aberrations inherited from a phenotypically normal parent are usually without clinical consequences.

High frequency of submicroscopic chromosomal imbalances in patients with syndromic craniosynostosis detected by a combined approach of microsatellite segregation analysis, multiplex ligation-dependent probe amplification and array-based comparative genome hybridisation

We present the first comprehensive study, to our knowledge, on genomic chromosomal analysis in syndromic craniosynostosis. In total, 45 patients with craniosynostotic disorders were screened with a variety of methods including conventional karyotype, microsatellite segregation analysis, subtelomeric multiplex ligation-dependent probe amplification) and whole-genome array-based comparative genome hybridisation. Causative abnormalities were present in 42.2% (19/45) of the samples, and 27.8% (10/36) of the patients with normal conventional karyotype carried submicroscopic imbalances. Our results include a wide variety of imbalances and point to novel chromosomal regions associated with craniosynostosis. The high incidence of pure duplications or trisomies suggests that these are important mechanisms in craniosynostosis, particularly in cases involving the metopic suture.

A study of adrenocortical tumors in children and adolescents by a comparative genomic hybridization technique

Adrenocortical tumors (ACT) are rare neoplasms of the adrenal glands accounting for 0.2% of all pediatric cancers. However, the incidence of ACT in South Brazilian children is 10 to 15 times greater than the worldwide incidence. Comparative genomic hybridization studies have revealed the presence of a high degree of chromosomal instability in ACT. We evaluated 16 ACT, 8 of them carcinomas and 8 adenomas. The presence of changes in DNA copy numbers was determined by comparative genomic hybridization, and the findings were validated by real-time polymerase chain reaction on the basis of IGF-II gene expression. The adenomas showed a mean of 19.7 imbalances per case, with the most frequent gain and loss being 4p15.1-p15.3 and 20p11.2-p13.2, respectively. The carcinomas presented with a mean of 35.5 imbalances per case, with the more frequent gain being 2q14.1-q24.3 and the more frequent losses being 3q21-q26.2, 20q12-qter, and 22q11.2-q13.3. The most frequent imbalances in both adenomas and carcinomas were gains of 1p21-p31.2, 2p12-p21 and loss of 20p11.2-p12. The expression of IGF-II mRNA (11p15.5) was higher in samples that presented with a gain of this region. It has been established that great genomic instability exists in pediatric ACT.

Influence of dentin on pH of 2% chlorhexidine gel and calcium hydroxide alone or in combination

The aims of endodontic treatment in cases of apical periodontitis are to reduce as much as possible the number of microorganisms inside the root canal system and to inactivate toxins produced by them. Most of the times, these objectives are not achieved solely by chemomechanical preparation, and intracanal dressing may be necessary. In these cases, calcium hydroxide is used as a root canal dressing due to its well-known and recognized antimicrobial activity. Chlorhexidine has a wide spectrum of antimicrobial activity and its association with calcium hydroxide has been recommended in an attempt to amplify antimicrobial effects of calcium hydroxide. It is also known that dentin exerts a buffering effect under wide pH variations, and may be responsible for decreasing the antimicrobial activity of drugs inside the root canal. The objectives of this study were to assess the pH of 2% chlorhexidine gel and calcium hydroxide alone or in combination, as well as the influence of dentin on the pH of these compounds. Dentin powder was obtained from bovine teeth and added as 1.8% to the volume of the medications. All substances were individually stored in plastic flasks, in triplicate. A pH meter was used at five different moments to assess pH in viscous medium: immediately after preparation and after 24 h, and 7, 14, and 21 days. Results were analyzed by paired Student`s t-test. Statistically significant differences were observed in the 2% chlorhexidine gel group alone or associated with calcium hydroxide and added of dentin powder (P < 0.05). Mean pH values indicated the influence of dentin powder because of a significant increase in pH. Calcium hydroxide with propylene glycol as the vehicle always showed high pH, demonstrating that this compound was not affected by the presence of dentin.

Deletion 5q12: Delineation of a Phenotype Including Mental Retardation and Ocular Defects

Array-CGH enables the detection of submicroscopic chromosomal deletions and duplications and leads to an accurate delineation of the imbalances, raising the possibility of correlating genotype to phenotype and mapping minimal critical regions associated with particular patterns of clinical features. We report here on four patients sharing common clinical features (psychomotor retardation, coarse facies and ocular anomalies), with proximal 5q deletions identified by oligo array-CGH. The deletions range from 5.75 to 17.26-Mb in size and occurred de novo. A common 2.63-Mb region between the deletions described here can be defined in 5q12.1 (59,390,122-62,021,754 bp bp from 5pter, hg18) and includes 12 genes. Among them, KIF2A, which encodes a kinesin superfamily protein, is a particularly interesting candidate for the phenotype, as it suppresses the growth of axonal collateral branches and is involved in normal brain development. Ocular defects, albeit unspecific, seem to be common in the 5q12.1 deletion. Identification of additional cases of deletions involving the 5q12.1 region will allow more accurate genotype-phenotype correlations. (C) 2011 Wiley-Liss, Inc.

UBE2A Deficiency Syndrome: Mild to Severe Intellectual Disability Accompanied by Seizures, Absent Speech, Urogenital, and Skin Anomalies in Male Patients

We describe three patients with a comparable deletion encompassing SLC25A43, SLC25A5, CXorf56, UBE2A, NKRF, and two non-coding RNA genes, U1 and LOC100303728. Moderate to severe intellectual disability (ID), psychomotor retardation, severely impaired/absent speech, seizures, and urogenital anomalies were present in all three patients. Facial dysmorphisms include ocular hypertelorism, synophrys, and a depressed nasal bridge. These clinical features overlap with those described in two patients from a family with a similar deletion at Xq24 that also includes UBE2A, and in several patients of Brazilian and Polish families with point mutations in UBE2A. Notably, all five patients with an Xq24 deletion have ventricular septal defects that are not present inpatients with a point mutation, which might be attributed to the deletion of SLC25A5. Taken together, the UBE2A deficiency syndrome in male patients with a mutation in or a deletion of UBE2A is characterized by ID, absent speech, seizures, urogenital anomalies, frequently including a small penis, and skin abnormalities, which include generalized hirsutism, low posterior hairline, myxedematous appearance, widely spaced nipples, and hair whorls. Facial dysmorphisms include a wide face, a depressed nasal bridge, a large mouth with downturned corners, thin vermilion, and a short, broad neck. (C) 2010 Wiley-Liss, Inc.

Curiously composite structures of a retrotransposon and a complex repeat associated with chromosome ends of Rhynchosciara americana (Diptera: Sciaridae)

In Drosophila, telomere retrotransposons counterbalance the loss of telomeric DNA. The exceptional mechanism of telomere recovery characterized in Drosophila has not been found in lower dipterans (Nematocera). However, a retroelement resembling a telomere transposon and termed ""RaTART"" has been described in the nematoceran Rhynchosciara americana. In this work, DNA and protein sequence analyses, DNA cloning, and chromosomal localization of probes obtained either by PCR or by screening a genomic library were carried out in order to examine additional features of this retroelement. The analyses performed raise the possibility that RaTART represents a genomic clone composed of distinct repetitive elements, one of which is likely to be responsible for its apparent enrichment at chromosome ends. RaTART sequence in addition allowed to assess a novel subtelomeric region of R. americana chromosomes that was analyzed in this work after subcloning a DNA fragment from a phage insert. It contains a complex repeat that is located in the vicinity of simple and complex tandem repeats characterized previously. Quantification data suggest that the copy number of the repeat is significantly lower than that observed for the ribosomal DNA in the salivary gland of R. americana. A short insertion of the RaTART was identified in the cloned segment, which hybridized preferentially to subtelomeres. Like RaTART, it displays truncated sequences related to distinct retrotransposons, one of which has a conceptual translation product with significant identity with an endonuclease from a lepidopteran retrotransposon. The composite structure of this DNA stretch probably reflects mobile element activity in the subtelomeric region analyzed in this work.

Complex rearrangements in patients with duplications of MECP2 can occur by fork stalling and template switching

Duplication at the Xq28 band including the MECP2 gene is one of the most common genomic rearrangements identified in neurodevelopmentally delayed males. Such duplications are non-recurrent and can be generated by a non-homologous end joining (NHEJ) mechanism. We investigated the potential mechanisms for MECP2 duplication and examined whether genomic architectural features may play a role in their origin using a custom designed 4-Mb tiling-path oligonucleotide array CGH assay. Each of the 30 patients analyzed showed a unique duplication varying in size from similar to 250 kb to similar to 2.6 Mb. Interestingly, in 77% of these non-recurrent duplications, the distal breakpoints grouped within a 215 kb genomic interval, located 47 kb telomeric to the MECP2 gene. The genomic architecture of this region contains both direct and inverted low-copy repeat (LCR) sequences; this same region undergoes polymorphic structural variation in the general population. Array CGH revealed complex rearrangements in eight patients; in six patients the duplication contained an embedded triplicated segment, and in the other two, stretches of non-duplicated sequences occurred within the duplicated region. Breakpoint junction sequencing was achieved in four duplications and identified an inversion in one patient, demonstrating further complexity. We propose that the presence of LCRs in the vicinity of the MECP2 gene may generate an unstable DNA structure that can induce DNA strand lesions, such as a collapsed fork, and facilitate a Fork Stalling and Template Switching event producing the complex rearrangements involving MECP2.

Nonrecurrent MECP2 duplications mediated by genomic architecture-driven DNA breaks and break-induced replication repair

Recurrent submicroscopic genomic copy number changes are the result of nonallelic homologous recombination (NAHR). Nonrecurrent aberrations, however, can result from different nonexclusive recombination-repair mechanisms. We previously described small microduplications at Xq28 containing MECP2 in four male patients with a severe neurological phenotype. Here, we report on the fine-mapping and breakpoint analysis of 16 unique microduplications. The size of the overlapping copy number changes varies between 0.3 and 2.3 Mb, and FISH analysis on three patients demonstrated a tandem orientation. Although eight of the 32 breakpoint regions coincide with low-copy repeats, none of the duplications are the result of NAHR. Bioinformatics analysis of the breakpoint regions demonstrated a 2.5-fold higher frequency of Alu interspersed repeats as compared with control regions, as well as a very high GC content (53%). Unexpectedly, we obtained the junction in only one patient by long-range PCR, which revealed nonhomologous end joining as the mechanism. Breakpoint analysis in two other patients by inverse PCR and subsequent array comparative genomic hybridization analysis demonstrated the presence of a second duplicated region more telomeric at Xq28, of which one copy was inserted in between the duplicated MECP2 regions. These data suggest a two-step mechanism in which part of Xq28 is first inserted near the MECP2 locus, followed by breakage-induced replication with strand invasion of the normal sister chromatid. Our results indicate that the mechanism by which copy number changes occur in regions with a complex genomic architecture can yield complex rearrangements.

Precision of distances and ordering of microsatellite markers in consensus linkage maps of chromosomes 1, 3 and 4 from two reciprocal chicken populations using bootstrap sampling

Some factors complicate comparisons between linkage maps from different studies. This problem can be resolved if measures of precision, such as confidence intervals and frequency distributions, are associated with markers. We examined the precision of distances and ordering of microsatellite markers in the consensus linkage maps of chromosomes 1, 3 and 4 from two F 2 reciprocal Brazilian chicken populations, using bootstrap sampling. Single and consensus maps were constructed. The consensus map was compared with the International Consensus Linkage Map and with the whole genome sequence. Some loci showed segregation distortion and missing data, but this did not affect the analyses negatively. Several inversions and position shifts were detected, based on 95% confidence intervals and frequency distributions of loci. Some discrepancies in distances between loci and in ordering were due to chance, whereas others could be attributed to other effects, including reciprocal crosses, sampling error of the founder animals from the two populations, F(2) population structure, number of and distance between microsatellite markers, number of informative meioses, loci segregation patterns, and sex. In the Brazilian consensus GGA1, locus LEI1038 was in a position closer to the true genome sequence than in the International Consensus Map, whereas for GGA3 and GGA4, no such differences were found. Extending these analyses to the remaining chromosomes should facilitate comparisons and the integration of several available genetic maps, allowing meta-analyses for map construction and quantitative trait loci (QTL) mapping. The precision of the estimates of QTL positions and their effects would be increased with such information.

Karyotypic evolution trends in Rhamdia quelen (Siluriformes, Heptapteridae) with considerations about the origin and differentiation of its supernumerary chromosomes

Among catfish species of the genus Rhamdia reported for the Brazilian territory, R. quelen is the most widespread, being found in nearly all hydrographic basins of Brazil. Nowadays, R. quelen is a synonym for at least 47 other species in this genus, its taxonomic status still being controversial. The available cytogenetic reports show a wide variation in the karyotypic macrostructure, with the frequent presence of supernumerary chromosomes. The remarkable cytogenetic variability associated with taxonomic issues in this species indicates that R. quelen is actually a species complex. In order to carry out a wide comparative cytogenetic study in R. quelen from southern and southeastern Brazil and examine a species complex, we analyzed the chromosomes of 14 populations from the main hydrographic basins of these two regions. Using classic and molecular cytogenetic techniques, we found seven distinct karyotypic formulae, all bearing 2n = 58 chromosomes. Supernumerary chromosomes were present in most of the populations; their number, size and C-banding pattern allowed us to differentiate populations with similar karyotypic compositions. We examined patterns of chromosomal evolution as well as the probable mechanisms involved in the origin and morphological differentiation of their supernumerary chromosomes.

Impact of Manaus City on the Amazon Green Ocean atmosphere: ozone production, precursor sensitivity and aerosol load

As a contribution to the Large-Scale Biosphere-Atmosphere Experiment in Amazonia - Cooperative LBA Airborne Regional Experiment (LBA-CLAIRE-2001) field campaign in the heart of the Amazon Basin, we analyzed the temporal and spatial dynamics of the urban plume of Manaus City during the wet-to-dry season transition period in July 2001. During the flights, we performed vertical stacks of crosswind transects in the urban outflow downwind of Manaus City, measuring a comprehensive set of trace constituents including O(3), NO, NO(2), CO, VOC, CO(2), and H(2)O. Aerosol loads were characterized by concentrations of total aerosol number (CN) and cloud condensation nuclei (CCN), and by light scattering properties. Measurements over pristine rainforest areas during the campaign showed low levels of pollution from biomass burning or industrial emissions, representative of wet season background conditions. The urban plume of Manaus City was found to be joined by plumes from power plants south of the city, all showing evidence of very strong photochemical ozone formation. One episode is discussed in detail, where a threefold increase in ozone mixing ratios within the atmospheric boundary layer occurred within a 100 km travel distance downwind of Manaus. Observation-based estimates of the ozone production rates in the plume reached 15 ppb h(-1). Within the plume core, aerosol concentrations were strongly enhanced, with Delta CN/Delta CO ratios about one order of magnitude higher than observed in Amazon biomass burning plumes. Delta CN/Delta CO ratios tended to decrease with increasing transport time, indicative of a significant reduction in particle number by coagulation, and without substantial new particle nucleation occurring within the time/space observed. While in the background atmosphere a large fraction of the total particle number served as CCN (about 60-80% at 0.6% supersaturation), the CCN/CN ratios within the plume indicated that only a small fraction (16 +/- 12 %) of the plume particles were CCN. The fresh plume aerosols showed relatively weak light scattering efficiency. The CO-normalized CCN concentrations and light scattering coefficients increased with plume age in most cases, suggesting particle growth by condensation of soluble organic or inorganic species. We used a Single Column Chemistry and Transport Model (SCM) to infer the urban pollution emission fluxes of Manaus City, implying observed mixing ratios of CO, NO(x) and VOC. The model can reproduce the temporal/spatial distribution of ozone enhancements in the Manaus plume, both with and without accounting for the distinct (high NO(x)) contribution by the power plants; this way examining the sensitivity of ozone production to changes in the emission rates of NO(x). The VOC reactivity in the Manaus region was dominated by a high burden of biogenic isoprene from the background rainforest atmosphere, and therefore NO(x) control is assumed to be the most effective ozone abatement strategy. Both observations and models show that the agglomeration of NO(x) emission sources, like power plants, in a well-arranged area can decrease the ozone production efficiency in the near field of the urban populated cores. But on the other hand remote areas downwind of the city then bear the brunt, being exposed to increased ozone production and N-deposition. The simulated maximum stomatal ozone uptake fluxes were 4 nmol m(-2) s(-1) close to Manaus, and decreased only to about 2 nmol m(-2) s(-1) within a travel distance >1500 km downwind from Manaus, clearly exceeding the critical threshold level for broadleaf trees. Likewise, the simulated N deposition close to Manaus was similar to 70 kg N ha(-1) a(-1) decreasing only to about 30 kg N ha(-1) a(-1) after three days of simulation.

Cloud condensation nuclei in pristine tropical rainforest air of Amazonia: size-resolved measurements and modeling of atmospheric aerosol composition and CCN activity

Atmospheric aerosol particles serving as cloud condensation nuclei (CCN) are key elements of the hydrological cycle and climate. We have measured and characterized CCN at water vapor supersaturations in the range of S=0.10-0.82% in pristine tropical rainforest air during the AMAZE-08 campaign in central Amazonia. The effective hygroscopicity parameters describing the influence of chemical composition on the CCN activity of aerosol particles varied in the range of kappa approximate to 0.1-0.4 (0.16+/-0.06 arithmetic mean and standard deviation). The overall median value of kappa approximate to 0.15 was by a factor of two lower than the values typically observed for continental aerosols in other regions of the world. Aitken mode particles were less hygroscopic than accumulation mode particles (kappa approximate to 0.1 at D approximate to 50 nm; kappa approximate to 0.2 at D approximate to 200 nm), which is in agreement with earlier hygroscopicity tandem differential mobility analyzer (H-TDMA) studies. The CCN measurement results are consistent with aerosol mass spectrometry (AMS) data, showing that the organic mass fraction (f(org)) was on average as high as similar to 90% in the Aitken mode (D <= 100 nm) and decreased with increasing particle diameter in the accumulation mode (similar to 80% at D approximate to 200 nm). The kappa values exhibited a negative linear correlation with f(org) (R(2)=0.81), and extrapolation yielded the following effective hygroscopicity parameters for organic and inorganic particle components: kappa(org)approximate to 0.1 which can be regarded as the effective hygroscopicity of biogenic secondary organic aerosol (SOA) and kappa(inorg)approximate to 0.6 which is characteristic for ammonium sulfate and related salts. Both the size dependence and the temporal variability of effective particle hygroscopicity could be parameterized as a function of AMS-based organic and inorganic mass fractions (kappa(p)=kappa(org) x f(org)+kappa(inorg) x f(inorg)). The CCN number concentrations predicted with kappa(p) were in fair agreement with the measurement results (similar to 20% average deviation). The median CCN number concentrations at S=0.1-0.82% ranged from N(CCN,0.10)approximate to 35 cm(-3) to N(CCN,0.82)approximate to 160 cm(-3), the median concentration of aerosol particles larger than 30 nm was N(CN,30)approximate to 200 cm(-3), and the corresponding integral CCN efficiencies were in the range of N(CCN,0.10/NCN,30)approximate to 0.1 to N(CCN,0.82/NCN,30)approximate to 0.8. Although the number concentrations and hygroscopicity parameters were much lower in pristine rainforest air, the integral CCN efficiencies observed were similar to those in highly polluted megacity air. Moreover, model calculations of N(CCN,S) assuming an approximate global average value of kappa approximate to 0.3 for continental aerosols led to systematic overpredictions, but the average deviations exceeded similar to 50% only at low water vapor supersaturation (0.1%) and low particle number concentrations (<= 100 cm(-3)). Model calculations assuming aconstant aerosol size distribution led to higher average deviations at all investigated levels of supersaturation: similar to 60% for the campaign average distribution and similar to 1600% for a generic remote continental size distribution. These findings confirm earlier studies suggesting that aerosol particle number and size are the major predictors for the variability of the CCN concentration in continental boundary layer air, followed by particle composition and hygroscopicity as relatively minor modulators. Depending on the required and applicable level of detail, the information and parameterizations presented in this paper should enable efficient description of the CCN properties of pristine tropical rainforest aerosols of Amazonia in detailed process models as well as in large-scale atmospheric and climate models.

Knoop Microhardness and FT-Raman Spectroscopic Evaluation of a Resin-Based Dental Material Light-Cured by an Argon Ion Laser and Halogen Lamp: An in Vitro Study

Objective: The purpose of this study was to evaluate in vitro the Knoop microhardness (Knoop hardness number [KHN]) and the degree of conversion using FT-Raman spectroscopy of a light-cured microhybrid resin composite (Z350-3M-ESPE) Vita shade A3 photopolymerized with a halogen lamp or an argon ion laser. Background Data: Optimal polymerization of resin-based dental materials is important for longevity of restorations in dentistry. Materials and Methods: Thirty specimens were prepared and inserted into a disc-shaped polytetrafluoroethylene mold that was 2.0 mm thick and 3 mm in diameter. The specimens were divided into three groups (n = 10 each). Group 1 (G1) was light-cured for 20 sec with an Optilux 501 halogen light with an intensity of 1000 mW/cm(2). Group 2 (G2) was photopolymerized with an argon laser with a power of 150 mW for 10 sec, and group 3 (G3) was photopolymerized with an argon laser at 200 mW of power for 10 sec. All specimens were stored in distilled water for 24 h at 37 degrees C and kept in lightproof containers. For the KHN test five indentations were made and a depth of 100 mu m was maintained in each specimen. One hundred and fifty readings were obtained using a 25-g load for 45 sec. The degree of conversion values were measured by Raman spectroscopy. KHN and degree of conversion values were obtained on opposite sides of the irradiated surface. KHN and degree of conversion data were analyzed by one-way ANOVA and Tukey tests with statistical significance set at p < 0.05. Results: The results of KHN testing were G1 = 37.428 +/- 4.765; G2 = 23.588 +/- 6.269; and G3 = 21.652 +/- 4.393. The calculated degrees of conversion (DC%) were G1 = 48.57 +/- 2.11; G2 = 43.71 +/- 3.93; and G3 = 44.19 +/- 2.71. Conclusions: Polymerization with the halogen lamp ( G1) attained higher microhardness values than polymerization with the argon laser at power levels of 150 and 200 mW; there was no difference in hardness between the two argon laser groups. The results showed no statistically significant different degrees of conversion for the polymerization of composite samples with the two light sources tested.

The role of exon shuffling in shaping protein-protein interaction networks

Background: Physical protein-protein interaction (PPI) is a critical phenomenon for the function of most proteins in living organisms and a significant fraction of PPIs are the result of domain-domain interactions. Exon shuffling, intron-mediated recombination of exons from existing genes, is known to have been a major mechanism of domain shuffling in metazoans. Thus, we hypothesized that exon shuffling could have a significant influence in shaping the topology of PPI networks. Results: We tested our hypothesis by compiling exon shuffling and PPI data from six eukaryotic species: Homo sapiens, Mus musculus, Drosophila melanogaster, Caenorhabditis elegans, Cryptococcus neoformans and Arabidopsis thaliana. For all four metazoan species, genes enriched in exon shuffling events presented on average higher vertex degree (number of interacting partners) in PPI networks. Furthermore, we verified that a set of protein domains that are simultaneously promiscuous (known to interact to multiple types of other domains), self-interacting (able to interact with another copy of themselves) and abundant in the genomes presents a stronger signal for exon shuffling. Conclusions: Exon shuffling appears to have been a recurrent mechanism for the emergence of new PPIs along metazoan evolution. In metazoan genomes, exon shuffling also promoted the expansion of some protein domains. We speculate that their promiscuous and self-interacting properties may have been decisive for that expansion.