Assortative mating and parental genetic relatedness contribute to the pathogenicity of variably expressive variants.

We examined more than 97,000 families from four neurodevelopmental disease cohorts and the UK Biobank to identify phenotypic and genetic patterns in parents contributing to neurodevelopmental disease risk in children. We identified within- and cross-disorder correlations between six phenotypes in parents and children, such as obsessive-compulsive disorder (R = 0.32-0.38, p < 10-126). We also found that measures of sub-clinical autism features in parents are associated with several autism severity measures in children, including biparental mean Social Responsiveness Scale scores and proband Repetitive Behaviors Scale scores (regression coefficient = 0.14, p = 3.38 × 10-4). We further describe patterns of phenotypic similarity between spouses, where spouses show correlations for six neurological and psychiatric phenotypes, including a within-disorder correlation for depression (R = 0.24-0.68, p < 0.001) and a cross-disorder correlation between anxiety and bipolar disorder (R = 0.09-0.22, p < 10-92). Using a simulated population, we also found that assortative mating can lead to increases in disease liability over generations and the appearance of "genetic anticipation" in families carrying rare variants. We identified several families in a neurodevelopmental disease cohort where the proband inherited multiple rare variants in disease-associated genes from each of their affected parents. We further identified parental relatedness as a risk factor for neurodevelopmental disorders through its inverse relationship with variant pathogenicity and propose that parental relatedness modulates disease risk by increasing genome-wide homozygosity in children (R = 0.05-0.26, p < 0.05). Our results highlight the utility of assessing parent phenotypes and genotypes toward predicting features in children who carry rare variably expressive variants and implicate assortative mating as a risk factor for increased disease severity in these families.

Understanding the new BRD4-related syndrome: Clinical and genomic delineation with an international cohort study.

BRD4 is part of a multiprotein complex involved in loading the cohesin complex onto DNA, a fundamental process required for cohesin-mediated loop extrusion and formation of Topologically Associating Domains. Pathogenic variations in this complex have been associated with a growing number of syndromes, collectively known as cohesinopathies, the most classic being Cornelia de Lange syndrome. However, no cohort study has been conducted to delineate the clinical and molecular spectrum of BRD4-related disorder. We formed an international collaborative study, and collected 14 new patients, including two fetuses. We performed phenotype and genotype analysis, integrated prenatal findings from fetopathological examinations, phenotypes of pediatric patients and adults. We report the first cohort of patients with BRD4-related disorder and delineate the dysmorphic features at different ages. This work extends the phenotypic spectrum of cohesinopathies and characterize a new clinically relevant and recognizable pattern, distinguishable from the other cohesinopathies.

Novel Loss of Function Variant in BCKDK Causes a Treatable Developmental and Epileptic Encephalopathy.

Branched-chain amino acids (BCAA) are essential amino acids playing crucial roles in protein synthesis and brain neurotransmission. Branched-chain ketoacid dehydrogenase (BCKDH), the flux-generating step of BCAA catabolism, is tightly regulated by reversible phosphorylation of its E1α-subunit. BCKDK is the kinase responsible for the phosphorylation-mediated inactivation of BCKDH. In three siblings with severe developmental delays, microcephaly, autism spectrum disorder and epileptic encephalopathy, we identified a new homozygous in-frame deletion (c.999_1001delCAC; p.Thr334del) of BCKDK. Plasma and cerebrospinal fluid concentrations of BCAA were markedly reduced. Hyperactivity of BCKDH and over-consumption of BCAA were demonstrated by functional tests in cells transfected with the mutant BCKDK. Treatment with pharmacological doses of BCAA allowed the restoring of BCAA concentrations and greatly improved seizure control. Behavioral and developmental skills of the patients improved to a lesser extent. Importantly, a retrospective review of the newborn screening results allowed the identification of a strong decrease in BCAA concentrations on dried blood spots, suggesting that BCKDK is a new treatable metabolic disorder probably amenable to newborn screening programs.

16p13.11 microduplication in 45 new patients: refined clinical significance and genotype-phenotype correlations.

BACKGROUND: The clinical significance of 16p13.11 duplications remains controversial while frequently detected in patients with developmental delay (DD), intellectual deficiency (ID) or autism spectrum disorder (ASD). Previously reported patients were not or poorly characterised. The absence of consensual recommendations leads to interpretation discrepancy and makes genetic counselling challenging. This study aims to decipher the genotype-phenotype correlations to improve genetic counselling and patients' medical care. METHODS: We retrospectively analysed data from 16 013 patients referred to 12 genetic centers for DD, ID or ASD, and who had a chromosomal microarray analysis. The referring geneticists of patients for whom a 16p13.11 duplication was detected were asked to complete a questionnaire for detailed clinical and genetic data for the patients and their parents. RESULTS: Clinical features are mainly speech delay and learning disabilities followed by ASD. A significant risk of cardiovascular disease was noted. About 90% of the patients inherited the duplication from a parent. At least one out of four parents carrying the duplication displayed a similar phenotype to the propositus. Genotype-phenotype correlations show no impact of the size of the duplicated segment on the severity of the phenotype. However, NDE1 and miR-484 seem to have an essential role in the neurocognitive phenotype. CONCLUSION: Our study shows that 16p13.11 microduplications are likely pathogenic when detected in the context of DD/ID/ASD and supports an essential role of NDE1 and miR-484 in the neurocognitive phenotype. Moreover, it suggests the need for cardiac evaluation and follow-up and a large study to evaluate the aortic disease risk.

Estimating the effect size of the 15Q11.2 BP1-BP2 deletion and its contribution to neurodevelopmental symptoms: recommendations for practice.

BACKGROUND: The 15q11.2 deletion is frequently identified in the neurodevelopmental clinic. Case-control studies have associated the 15q11.2 deletion with neurodevelopmental disorders, and clinical case series have attempted to delineate a microdeletion syndrome with considerable phenotypic variability. The literature on this deletion is extensive and confusing, which is a challenge for genetic counselling. The aim of this study was to estimate the effect size of the 15q11.2 deletion and quantify its contribution to neurodevelopmental disorders. METHODS: We performed meta-analyses on new and previously published case-control studies and used statistical models trained in unselected populations with cognitive assessments. We used new (n=241) and previously published (n=150) data from a clinically referred group of deletion carriers. 15q11.2 duplications (new n=179 and previously published n=35) were used as a neutral control variant. RESULTS: The deletion decreases IQ by 4.3 points. The estimated ORs and respective frequencies in deletion carriers for intellectual disabilities, schizophrenia and epilepsy are 1.7 (3.4%), 1.5 (2%) and 3.1 (2.1%), respectively. There is no increased risk for heart malformations and autism. In the clinically referred group, the frequency and nature of symptoms in deletions are not different from those observed in carriers of the 15q11.2 duplication suggesting that most of the reported symptoms are due to ascertainment bias. CONCLUSIONS: We recommend that the deletion should be classified as 'pathogenic of mild effect size'. Since it explains only a small proportion of the phenotypic variance in carriers, it is not worth discussing in the developmental clinic or in a prenatal setting.

Rare variants in the genetic background modulate cognitive and developmental phenotypes in individuals carrying disease-associated variants.

PURPOSE: To assess the contribution of rare variants in the genetic background toward variability of neurodevelopmental phenotypes in individuals with rare copy-number variants (CNVs) and gene-disruptive variants. METHODS: We analyzed quantitative clinical information, exome sequencing, and microarray data from 757 probands and 233 parents and siblings who carry disease-associated variants. RESULTS: The number of rare likely deleterious variants in functionally intolerant genes ("other hits") correlated with expression of neurodevelopmental phenotypes in probands with 16p12.1 deletion (n=23, p=0.004) and in autism probands carrying gene-disruptive variants (n=184, p=0.03) compared with their carrier family members. Probands with 16p12.1 deletion and a strong family history presented more severe clinical features (p=0.04) and higher burden of other hits compared with those with mild/no family history (p=0.001). The number of other hits also correlated with severity of cognitive impairment in probands carrying pathogenic CNVs (n=53) or de novo pathogenic variants in disease genes (n=290), and negatively correlated with head size among 80 probands with 16p11.2 deletion. These co-occurring hits involved known disease-associated genes such as SETD5, AUTS2, and NRXN1, and were enriched for cellular and developmental processes. CONCLUSION: Accurate genetic diagnosis of complex disorders will require complete evaluation of the genetic background even after a candidate disease-associated variant is identified.

Exome copy number variation detection: Use of a pool of unrelated healthy tissue as reference sample.

An increasing number of bioinformatic tools designed to detect CNVs (copy number variants) in tumor samples based on paired exome data where a matched healthy tissue constitutes the reference have been published in the recent years. The idea of using a pool of unrelated healthy DNA as reference has previously been formulated but not thoroughly validated. As of today, the gold standard for CNV calling is still aCGH but there is an increasing interest in detecting CNVs by exome sequencing. We propose to design a metric allowing the comparison of two CNV profiles, independently of the technique used and assessed the validity of using a pool of unrelated healthy DNA instead of a matched healthy tissue as reference in exome-based CNV detection. We compared the CNV profiles obtained with three different approaches (aCGH, exome sequencing with a matched healthy tissue as reference, exome sequencing with a pool of eight unrelated healthy tissue as reference) on three multiple myeloma samples. We show that the usual analyses performed to compare CNV profiles (deletion/amplification ratios and CNV size distribution) lack in precision when confronted with low LRR values, as they only consider the binary status of each CNV. We show that the metric-based distance constitutes a more accurate comparison of two CNV profiles. Based on these analyses, we conclude that a reliable picture of CNV alterations in multiple myeloma samples can be obtained from whole-exome sequencing in the absence of a matched healthy sample.

Genomic studies of multiple myeloma reveal an association between X chromosome alterations and genomic profile complexity.

The genomic profile of multiple myeloma (MM) has prognostic value by dividing patients into a good prognosis hyperdiploid group and a bad prognosis nonhyperdiploid group with a higher incidence of IGH translocations. This classification, however, is inadequate and many other parameters like mutations, epigenetic modifications, and genomic heterogeneity may influence the prognosis. We performed a genomic study by array-based comparative genomic hybridization on a cohort of 162 patients to evaluate the frequency of genomic gains and losses. We identified a high frequency of X chromosome alterations leading to partial Xq duplication, often associated with inactive X (Xi) deletion in female patients. This partial X duplication could be a cytogenetic marker of aneuploidy as it is correlated with a high number of chromosomal breakages. Patient with high level of chromosomal breakage had reduced survival regardless the region implicated. A higher transcriptional level was shown for genes with potential implication in cancer and located in this altered region. Among these genes, IKBKG and IRAK1 are members of the NFKB pathway which plays an important role in MM and is a target for specific treatments. © 2016 Wiley Periodicals, Inc.

Gigantism and acromegaly due to Xq26 microduplications and GPR101 mutation.

BACKGROUND: Increased secretion of growth hormone leads to gigantism in children and acromegaly in adults; the genetic causes of gigantism and acromegaly are poorly understood. METHODS: We performed clinical and genetic studies of samples obtained from 43 patients with gigantism and then sequenced an implicated gene in samples from 248 patients with acromegaly. RESULTS: We observed microduplication on chromosome Xq26.3 in samples from 13 patients with gigantism; of these samples, 4 were obtained from members of two unrelated kindreds, and 9 were from patients with sporadic cases. All the patients had disease onset during early childhood. Of the patients with gigantism who did not carry an Xq26.3 microduplication, none presented before the age of 5 years. Genomic characterization of the Xq26.3 region suggests that the microduplications are generated during chromosome replication and that they contain four protein-coding genes. Only one of these genes, GPR101, which encodes a G-protein-coupled receptor, was overexpressed in patients' pituitary lesions. We identified a recurrent GPR101 mutation (p.E308D) in 11 of 248 patients with acromegaly, with the mutation found mostly in tumors. When the mutation was transfected into rat GH3 cells, it led to increased release of growth hormone and proliferation of growth hormone-producing cells. CONCLUSIONS: We describe a pediatric disorder (which we have termed X-linked acrogigantism [X-LAG]) that is caused by an Xq26.3 genomic duplication and is characterized by early-onset gigantism resulting from an excess of growth hormone. Duplication of GPR101 probably causes X-LAG. We also found a recurrent mutation in GPR101 in some adults with acromegaly. (Funded by the Eunice Kennedy Shriver National Institute of Child Health and Human Development and others.).

Combined treatment with octreotide LAR and pegvisomant in patients with pituitary gigantism: clinical evaluation and genetic screening.

INTRODUCTION: Pituitary gigantism is a rare condition caused by growth hormone secreting hypersecretion, usually by a pituitary tumor. Acromegaly and gigantism cases that have a genetic cause are challenging to treat, due to large tumor size and poor responses to some medical therapies (e.g. AIP mutation affected cases and those with X-linked acrogigantism syndrome). MATERIALS AND METHODS: We performed a retrospective study to identify gigantism cases among 160 somatotropinoma patients treated between 1985 and 2015 at the University Hospital of Caracas, Venezuela. We studied clinical details at diagnosis, hormonal responses to therapy and undertook targeted genetic testing. Among the 160 cases, eight patients (six males; 75 %) were diagnosed with pituitary gigantism and underwent genetic analysis that included array comparative genome hybridization for Xq26.3 duplications. RESULTS: All patients had GH secreting pituitary macroadenomas that were difficult to control with conventional treatment options, such as surgery or primary somatostatin receptor ligand (SRL) therapy. Combined therapy (long-acting SRL and pegvisomant) as primary treatment or after pituitary surgery and radiotherapy permitted the normalization of IGF-1 levels and clinical improvement. Novel AIP mutations were the found in three patients. None of the patients had Xq26.3 microduplications. CONCLUSIONS: Treatment of pituitary gigantism is frequently challenging; delayed control increases the harmful effects of GH excess, such as, excessive stature and symptom burden, so early diagnosis and effective treatment are particularly important in these cases.

Cytogenetic Studies of Rwandan Pediatric Patients Presenting with Global Developmental Delay, Intellectual Disability and/or Multiple Congenital Anomalies.

Global developmental delay (GDD) is defined as a significant delay in two or more developmental domains: gross or fine motor, speech/language, cognitive, social/personal and activities of daily living. Many of these children will go on to be diagnosed with intellectual disability (ID), which is most commonly defined as having an IQ <75 in addition to impairment in adaptive functioning. Cytogenetic studies have been performed in 664 Rwandan pediatric patients presenting GDD/ID and/or multiple congenital abnormalities (MCA). Karyotype analysis was performed in all patients and revealed 260 chromosomal abnormalities. The most frequent chromosomal abnormality was Down syndrome and then Edward syndrome and Patau syndrome. Other identified chromosomal abnormalities included 47,XX,+del(9)(q11), 46,XY,del(13)(q34) and 46,XX,der(22)t(10;22)(p10;p10)mat. In conclusion, our results highlight the high frequency of cytogenetically detectable abnormalities in this series, with implications for the burden on the healthcare. This study demonstrates the importance of cytogenetic analysis in patients with GDD/ID and MCA.

Mutation of the iron-sulfur cluster assembly gene IBA57 causes fatal infantile leukodystrophy.

Leukodystrophies are a heterogeneous group of severe genetic neurodegenerative disorders. A multiple mitochondrial dysfunctions syndrome was found in an infant presenting with a progressive leukoencephalopathy. Homozygosity mapping, whole exome sequencing, and functional studies were used to define the underlying molecular defect. Respiratory chain studies in skeletal muscle isolated from the proband revealed a combined deficiency of complexes I and II. In addition, western blotting indicated lack of protein lipoylation. The combination of these findings was suggestive for a defect in the iron-sulfur (Fe/S) protein assembly pathway. SNP array identified loss of heterozygosity in large chromosomal regions, covering the NFU1 and BOLA3, and the IBA57 and ABCB10 candidate genes, in 2p15-p11.2 and 1q31.1-q42.13, respectively. A homozygous c.436C > T (p.Arg146Trp) variant was detected in IBA57 using whole exome sequencing. Complementation studies in a HeLa cell line depleted for IBA57 showed that the mutant protein with the semi-conservative amino acid exchange was unable to restore the biochemical phenotype indicating a loss-of-function mutation of IBA57. In conclusion, defects in the Fe/S protein assembly gene IBA57 can cause autosomal recessive neurodegeneration associated with progressive leukodystrophy and fatal outcome at young age. In the affected patient, the biochemical phenotype was characterized by a defect in the respiratory chain complexes I and II and a decrease in mitochondrial protein lipoylation, both resulting from impaired assembly of Fe/S clusters.

Array-CGH analysis in Rwandan patients presenting development delay/intellectual disability with multiple congenital anomalies.

BACKGROUND: Array-CGH is considered as the first-tier investigation used to identify copy number variations. Right now, there is no available data about the genetic etiology of patients with development delay/intellectual disability and congenital malformation in East Africa. METHODS: Array comparative genomic hybridization was performed in 50 Rwandan patients with development delay/intellectual disability and multiple congenital abnormalities, using the Agilent's 180 K microarray platform. RESULTS: Fourteen patients (28%) had a global development delay whereas 36 (72%) patients presented intellectual disability. All patients presented multiple congenital abnormalities. Clinically significant copy number variations were found in 13 patients (26%). Size of CNVs ranged from 0,9 Mb to 34 Mb. Six patients had CNVs associated with known syndromes, whereas 7 patients presented rare genomic imbalances. CONCLUSION: This study showed that CNVs are present in African population and show the importance to implement genetic testing in East-African countries.

Genetic diagnosis of Duchenne and Becker muscular dystrophy using multiplex ligation-dependent probe amplification in Rwandan patients.

Duchenne and Becker muscular dystrophies are the most common clinical forms of muscular dystrophies. They are genetically X-linked diseases caused by a mutation in the dystrophin (DMD) gene. A genetic diagnosis was carried out in six Rwandan patients presenting a phenotype of Duchenne and Becker muscular dystrophies and six asymptomatic female carrier relatives using multiplex ligation-dependent probe amplification (MLPA). Our results revealed deletion of the exons 48-51 in one patient, an inherited deletion of the exons 8-21 in two brothers and a de novo deletion of the exons 46-50 in the fourth patient. No copy number variation was found in two patients. Only one female carrier presented exon deletion in the DMD gene. This is the first cohort of genetic analysis in Rwandan patients affected by Duchenne and Becker muscular dystrophies. This report confirmed that MLPA assay can be easily implemented in low-income countries.

Assortative mating and parental genetic relatedness drive the pathogenicity of variably expressive variants.

We examined more than 38,000 spouse pairs from four neurodevelopmental disease cohorts and the UK Biobank to identify phenotypic and genetic patterns in parents associated with neurodevelopmental disease risk in children. We identified correlations between six phenotypes in parents and children, including correlations of clinical diagnoses such as obsessive-compulsive disorder (R=0.31-0.49, p<0.001), and two measures of sub-clinical autism features in parents affecting several autism severity measures in children, such as bi-parental mean Social Responsiveness Scale (SRS) scores affecting proband SRS scores (regression coefficient=0.11, p=0.003). We further describe patterns of phenotypic and genetic similarity between spouses, where spouses show both within- and cross-disorder correlations for seven neurological and psychiatric phenotypes, including a within-disorder correlation for depression (R=0.25-0.72, p<0.001) and a cross-disorder correlation between schizophrenia and personality disorder (R=0.20-0.57, p<0.001). Further, these spouses with similar phenotypes were significantly correlated for rare variant burden (R=0.07-0.57, p<0.0001). We propose that assortative mating on these features may drive the increases in genetic risk over generations and the appearance of "genetic anticipation" associated with many variably expressive variants. We further identified parental relatedness as a risk factor for neurodevelopmental disorders through its inverse correlations with burden and pathogenicity of rare variants and propose that parental relatedness drives disease risk by increasing genome-wide homozygosity in children (R=0.09-0.30, p<0.001). Our results highlight the utility of assessing parent phenotypes and genotypes in predicting features in children carrying variably expressive variants and counseling families carrying these variants.

Newborn screening of duchenne muscular dystrophy specifically targeting deletions amenable to exon-skipping therapy.

Duchenne Muscular Dystrophy (DMD) is a lethal progressive muscle-wasting disease. New treatment strategies relying on DMD gene exon-skipping therapy have recently been approved and about 30% of patients could be amenable to exon 51, 53 or 45 skipping. We evaluated the spectrum of deletions reported in DMD registries, and designed a method to screen newborns and identify DMD deletions amenable to exon 51, 53 and 45 skipping. We developed a multiplex qPCR assay identifying hemi(homo)-zygotic deletions of the flanking exons of these therapeutic targets in DMD exons (i.e. exons 44, 46, 50, 52 and 54). We conducted an evaluation of our new method in 51 male patients with a DMD phenotype, 50 female carriers of a DMD deletion and 19 controls. Studies were performed on dried blood spots with patient's consent. We analyzed qPCR amplification curves of controls, carriers, and DMD patients to discern the presence or the absence of the target exons. Analysis of the exons flanking the exon-skipping targets permitted the identification of patients that could benefit from exon-skipping. All samples were correctly genotyped, with either presence or absence of amplification of the target exon. This proof-of-concept study demonstrates that this new assay is a highly sensitive method to identify DMD patients carrying deletions that are rescuable by exon-skipping treatment. The method is easily scalable to population-based screening. This targeted screening approach could address the new management paradigm in DMD, and could help to optimize the beneficial therapeutic effect of DMD therapies by permitting pre-symptomatic care.

Three years pilot of spinal muscular atrophy newborn screening turned into official program in Southern Belgium.

Three new therapies for spinal muscular atrophy (SMA) have been approved by the United States Food and Drug Administration and the European Medicines Agency since 2016. Although these new therapies improve the quality of life of patients who are symptomatic at first treatment, administration before the onset of symptoms is significantly more effective. As a consequence, newborn screening programs have been initiated in several countries. In 2018, we launched a 3-year pilot program to screen newborns for SMA in the Belgian region of Liège. This program was rapidly expanding to all of Southern Belgium, a region of approximately 55,000 births annually. During the pilot program, 136,339 neonates were tested for deletion of exon 7 of SMN1, the most common cause of SMA. Nine SMA cases with homozygous deletion were identified through this screen. Another patient was identified after presenting with symptoms and was shown to be heterozygous for the SMN1 exon 7 deletion and a point mutation on the opposite allele. These ten patients were treated. The pilot program has now successfully transitioned into the official neonatal screening program in Southern Belgium. The lessons learned during implementation of this pilot program are reported.

Surface Mucin-1 does not play a role in dendritic cell migration.

Mucin-1 (MUC1) is a transmembrane glycoprotein that is upregulated upon maturation of dendritic cells (DC) in vitro or in vivo. One of the proposed functions of surface expressed MUC1 is its involvement in migration of cells. We hypothesized that MUC1 is involved in DC migration since mature DC (mDC) are highly migratory cells and MUC1 is upregulated on the surface of DC upon maturation. In this study we cultured DC using two maturation cocktails, one cocktail containing IL-4, GM-CSF, TNFalpha, PGE2, IL-1 beta and IL-6 (TP1,6-DC) and the other IL-13, GM-CSF, Ribomunyl and IFN-gamma (RI-DC). Both maturation cocktails render DC with a similar surface phenotype including CCR7 expression, but only the former induces a migratory capacity of DC to a CCL19 gradient. To analyze the role of surface-expression of MUC1 on TP1,6-DC, that are capable of migration, expression of MUC1 was prevented by adding an anti-MUC1 antibody (Ab) during the maturation process. Compared with matured DC in the absence of the Ab, no difference was observed in chemokine-induced migratory behaviour between the MUC1+ and MUC1- DC populations in a standard Transwell chemotaxis assay, nor in organotypic cultures. Our data clearly demonstrate that surface MUC1 on DC does not influence intrinsic cell-motility, nor is it involved in cell-cell and cell-matrix dependent migration.

Increased migration of Langerhans cells in response to HPV16 E6 and E7 oncogene silencing: role of CCL20.

Human papillomavirus (HPV) infection, particularly type 16, is causally associated with cancer of the uterine cervix. The persistence or progression of cervical lesions suggests that viral antigens are not adequately presented to the immune system. This hypothesis is reinforced by the observation that most squamous intraepithelial lesions (SILs) show quantitative and functional alterations of Langerhans cells (LC). The infiltration of immature LC in the squamous epithelium is mainly controlled by Macrophage Inflammatory Protein 3alpha/CCL20. After having shown that CCL20 production is altered in HPV-transformed keratinocytes (KC), the possible role of HPV16 E6 and E7 viral oncoproteins in the reduced CCL20 levels observed in SILs was investigated by silencing HPV16 E6 and E7 oncogenes by RNA interference (siRNA). This treatment not only increased CCL20 secretion but also resulted in the modulation of NF-kappaB p50, p52 and p65 precursor localization. Moreover, silencing of E6 and E7 oncogenes in HPV16-transformed KC induced a significantly higher migratory capacity of LC in a Boyden chamber assay and in an in vitro formed (pre)neoplastic epithelium reminiscent of high-grade SILs. Anti-CCL20 neutralizing antibody experiments showed that the increased migration of LC is due to the re-expression of CCL20 in E6 and E7 siRNA transfected KC. These data suggest that HPV16 E6/E7-induced down-regulation of CCL20 observed during the cervical carcinogenesis may contribute to a diminished capacity of the immune system to control HPV infection.