Search

SETD2 belongs to the family of histone methyltransferase proteins and has been associated with three nosologically distinct entities with different clinical and molecular features: Luscan–Lumish syndrome (LLS), intellectual developmental disorder, autosomal dominant 70 (MRD70), and Rabin–Pappas syndrome (RAPAS). LLS [MIM #616831] is an overgrowth disorder with multisystem involvement including intellectual disability, speech delay, autism spectrum disorder (ASD), macrocephaly, tall stature, and motor delay. RAPAS [MIM #6201551] is a recently reported multisystemic disorder characterized by severely impaired global and intellectual development, hypotonia, feeding difficulties with failure to thrive, microcephaly, and dysmorphic facial features. Other neurologic findings may include seizures, hearing loss, ophthalmologic defects, and brain imaging abnormalities. There is variable involvement of other organ systems, including skeletal, genitourinary, cardiac, and potentially endocrine. Three patients who carried the missense variant p.Arg1740Gln in SETD2 were reported with a moderately impaired intellectual disability, speech difficulties, and behavioral abnormalities. More variable findings included hypotonia and dysmorphic features. Due to the differences with the two previous phenotypes, this association was then named intellectual developmental disorder, autosomal dominant 70 [MIM 620157]. These three disorders seem to be allelic and are caused either by loss-of-function, gain-of-function, or missense variants in the SETD2 gene. Here we describe 18 new patients with variants in SETD2, most of them with the LLS phenotype, and reviewed 33 additional patients with variants in SETD2 that have been previously reported in the scientific literature. This article offers an expansion of the number of reported individuals with LLS and highlights the clinical features and the similarities and differences among the three phenotypes associated with SETD2.

PURPOSE CTR9 is a subunit of the PAF1 complex (PAF1C) that plays a crucial role in transcription regulation by binding CTR9 to RNA polymerase II. It is involved in transcription-coupled histone modification through promoting H3K4 and H3K36 methylation. We describe the clinical and molecular studies in 13 probands, harboring likely pathogenic CTR9 missense variants, collected through GeneMatcher. METHODS Exome sequencing was performed in all individuals. CTR9 variants were assessed through 3-dimensional modeling of the activated human transcription complex Pol II-DSIF-PAF-SPT6 and the PAF1/CTR9 complex. H3K4/H3K36 methylation analysis, mitophagy assessment based on tetramethylrhodamine ethyl ester perchlorate immunofluorescence, and RNA-sequencing in skin fibroblasts from 4 patients was performed. RESULTS Common clinical findings were variable degrees of intellectual disability, hypotonia, joint hyperlaxity, speech delay, coordination problems, tremor, and autism spectrum disorder. Mild dysmorphism and cardiac anomalies were less frequent. For 11 CTR9 variants, de novo occurrence was shown. Three-dimensional modeling predicted a likely disruptive effect of the variants on local CTR9 structure and protein interaction. Additional studies in fibroblasts did not unveil the downstream functional consequences of the identified variants. CONCLUSION We describe a neurodevelopmental disorder caused by (mainly) de novo variants in CTR9, likely affecting PAF1C function.

The study aimed at widening the clinical and genetic spectrum of ASXL3-related syndrome, a neurodevelopmental disorder, caused by truncating variants in the ASXL3 gene. In this international collaborative study, we have undertaken a detailed clinical and molecular analysis of 45 previously unpublished individuals with ASXL3-related syndrome, as well as a review of all previously published individuals. We have reviewed the rather limited functional characterization of pathogenic variants in ASXL3 and discuss current understanding of the consequences of the different ASXL3 variants. In this comprehensive analysis of ASXL3-related syndrome, we define its natural history and clinical evolution occurring with age. We report familial ASXL3 pathogenic variants, characterize the phenotype in mildly affected individuals and discuss nonpenetrance. We also discuss the role of missense variants in ASXL3. We delineate a variable but consistent phenotype. The most characteristic features are neurodevelopmental delay with consistently limited speech, significant neuro-behavioral issues, hypotonia, and feeding difficulties. Distinctive features include downslanting palpebral fissures, hypertelorism, tubular nose with a prominent nasal bridge, and low-hanging columella. The presented data will inform clinical management of individuals with ASXL3-related syndrome and improve interpretation of new ASXL3 sequence variants.

Purpose To expand the recent description of a new neurodevelopmental syndrome related to alterations in CDK19. Methods Individuals were identified through international collaboration. Functional studies included autophosphorylation assays for CDK19 Gly28Arg and Tyr32His variants and in vivo zebrafish assays of the CDK19(G28R) and CDK19(Y32H). Results We describe 11 unrelated individuals (age range: 9 months to 14 years) with de novo missense variants mapped to the kinase domain of CDK19, including two recurrent changes at residues Tyr32 and Gly28. In vitro autophosphorylation and substrate phosphorylation assays revealed that kinase activity of protein was lower for p.Gly28Arg and higher for p.Tyr32His substitutions compared with that of the wild-type protein. Injection of CDK19 messenger RNA (mRNA) with either the Tyr32His or the Gly28Arg variants using in vivo zebrafish model significantly increased fraction of embryos with morphological abnormalities. Overall, the phenotype of the now 14 individuals with CDK19-related disorder includes universal developmental delay and facial dysmorphism, hypotonia (79%), seizures (64%), ophthalmologic anomalies (64%), and autism/autistic traits (56%). Conclusion CDK19 de novo missense variants are responsible for a novel neurodevelopmental disorder. Both kinase assay and zebrafish experiments showed that the pathogenetic mechanism may be more diverse than previously thought.

THG1L-associated autosomal recessive ataxia belongs to a group of disorders that occur due to abnormal mitochondrial tRNA modification. The product of THG1L is the tRNA-histidine guanylyltransferase 1-like enzyme that catalyzes the 3'-5"addition of guanine to the 5"-end of tRNA-histidine in the mitochondrion. To date, five individuals with homozygosity for p.(Val55Ala) in THG1L have been reported and presented with mild delays or normal development and cerebellar dysfunction. We present seven individuals with biallelic variants in THG1L. Three individuals were compound heterozygous for the p.(Cys51Trp) and p.(Val55Ala) variants and presented with profound developmental delays, microcephaly, intractable epilepsy, and cerebellar hypoplasia. Four siblings were homozygous for the p.(Val55Ala) variant and presented with cerebellar ataxia with cerebellar vermis hypoplasia, dysarthria, mild developmental delays, and normal/near-normal cognition. All seven patients were of Ashkenazi Jewish descent. Carrier rates for the two variants were calculated in a cohort of 26,731 Ashkenazi Jewish individuals tested by the Dor Yeshorim screening program. The p.(Cys51Trp) variant is novel and was found in 40 of the Ashkenazi Jewish individuals tested, with a carrier rate of 1 in 668 (0.15%). The p.(Val55Ala) variant was found in 229 of the Ashkenazi Jewish individuals tested, with a carrier rate of 1 in 117 (0.85%). The individuals with compound heterozygosity of the p.(Val55Ala) and p.(Cys51Trp) variants expand the phenotypic spectrum of THG1L-related disorders to include severe epileptic encephalopathy. The individuals with homozygosity of the p.(V55A) variant further establish the associated mild and slowly progressive or nonprogressive neurodevelopmental phenotype.

Nonsyndromic hearing loss is genetically heterogeneous. Despite comprehensive genetic testing, many cases remain unsolved because the clinical significance of identified variants is uncertain or because biallelic pathogenic variants are not identified for presumed autosomal recessive cases. Common synonymous variants are often disregarded. Determining the pathogenicity of synonymous variants may improve genetic diagnosis. We report a synonymous variant c.9861 C > T/p.(Gly3287=) in MYO15A in homozygosity or compound heterozygosity with another pathogenic or likely pathogenic MYO15A variant in 10 unrelated families with nonsyndromic sensorineural hearing loss. Biallelic variants in MYO15A were identified in 21 affected and were absent in 22 unaffected siblings. A mini-gene assay confirms that the synonymous variant leads to abnormal splicing. The variant is enriched in the Ashkenazi Jewish population. Individuals carrying biallelic variants involving c.9861 C > T often exhibit progressive post-lingual hearing loss distinct from the congenital profound deafness typically associated with biallelic loss-of-function MYO15A variants. This study establishes the pathogenicity of the c.9861 C > T variant in MYO15A and expands the phenotypic spectrum of MYO15A-related hearing loss. Our work also highlights the importance of multicenter collaboration and data sharing to establish the pathogenicity of a relatively common synonymous variant for improved diagnosis and management of hearing loss.

ZMYND8 encodes a multidomain protein that serves as a central interactive hub for coordinating critical roles in transcription regulation, chromatin remodeling, regulation of super-enhancers, DNA damage response and tumor suppression. We delineate a novel neurocognitive disorder caused by variants in the ZMYND8 gene.An international collaboration, exome sequencing, molecular modeling, yeast two-hybrid assays, analysis of available transcriptomic data and a knockdown Drosophila model were used to characterize the ZMYND8 variants.ZMYND8 variants were identified in 11 unrelated individuals; 10 occurred de novo and one suspected de novo; 2 were truncating, 9 were missense, of which one was recurrent. The disorder is characterized by intellectual disability with variable cardiovascular, ophthalmologic and minor skeletal anomalies. Missense variants in the PWWP domain of ZMYND8 abolish the interaction with Drebrin and missense variants in the MYND domain disrupt the interaction with GATAD2A. ZMYND8 is broadly expressed across cell types in all brain regions and shows highest expression in the early stages of brain development. Neuronal knockdown of the DrosophilaZMYND8 ortholog results in decreased habituation learning, consistent with a role in cognitive function.We present genomic and functional evidence for disruption of ZMYND8 as a novel etiology of syndromic intellectual disability.

The SET domain containing 2, histone lysine methyltransferase encoded by SETD2 is a dual-function methyltransferase for histones and microtubules and plays an important role for transcriptional regulation, genomic stability, and cytoskeletal functions. Specifically, SETD2 is associated with trimethylation of histone H3 at lysine 36 (H3K36me3) and methylation of α-tubulin at lysine 40. Heterozygous loss of function and missense variants have previously been described with Luscan-Lumish syndrome (LLS), which is characterized by overgrowth, neurodevelopmental features, and absence of overt congenital anomalies. We have identified 15 individuals with de novo variants in codon 1740 of SETD2 whose features differ from those with LLS. Group 1 consists of 12 individuals with heterozygous variant c.5218C>T p.(Arg1740Trp) and Group 2 consists of 3 individuals with heterozygous variant c.5219G>A p.(Arg1740Gln). The phenotype of Group 1 includes microcephaly, profound intellectual disability, congenital anomalies affecting several organ systems, and similar facial features. Individuals in Group 2 had moderate to severe intellectual disability, low normal head circumference, and absence of additional major congenital anomalies. While LLS is likely due to loss of function of SETD2, the clinical features seen in individuals with variants affecting codon 1740 are more severe suggesting an alternative mechanism, such as gain of function, effects on epigenetic regulation, or posttranslational modification of the cytoskeleton. Our report is a prime example of different mutations in the same gene causing diverging phenotypes and the features observed in Group 1 suggest a new clinically recognizable syndrome uniquely associated with the heterozygous variant c.5218C>T p.(Arg1740Trp) in SETD2.

Significance Approximately 400 United States children 1 y of age and older die suddenly from unexplained causes annually. We studied whole-exome sequence data from 124 “trios” (decedent child and living parents) to identify genetic risk factors. Nonsynonymous mutations, mostly de novo (present in child but absent in both biological parents), were highly enriched in genes associated with cardiac and seizure disorders relative to controls, and contributed to 9% of deaths. We found significant overtransmission of loss-of-function or pathogenic missense variants in cardiac and seizure disorder genes. Most pathogenic variants were de novo in origin, highlighting the importance of trio studies. Many of these pathogenic de novo mutations altered a protein network regulating calcium-related excitability at submembrane junctions in cardiomyocytes and neurons., Sudden unexplained death in childhood (SUDC) is an understudied problem. Whole-exome sequence data from 124 “trios” (decedent child, living parents) was used to test for excessive de novo mutations (DNMs) in genes involved in cardiac arrhythmias, epilepsy, and other disorders. Among decedents, nonsynonymous DNMs were enriched in genes associated with cardiac and seizure disorders relative to controls (odds ratio = 9.76, P = 2.15 × 10−4). We also found evidence for overtransmission of loss-of-function (LoF) or previously reported pathogenic variants in these same genes from heterozygous carrier parents (11 of 14 transmitted, P = 0.03). We identified a total of 11 SUDC proband genotypes (7 de novo, 1 transmitted parental mosaic, 2 transmitted parental heterozygous, and 1 compound heterozygous) as pathogenic and likely contributory to death, a genetic finding in 8.9% of our cohort. Two genes had recurrent missense DNMs, RYR2 and CACNA1C. Both RYR2 mutations are pathogenic (P = 1.7 × 10−7) and were previously studied in mouse models. Both CACNA1C mutations lie within a 104-nt exon (P = 1.0 × 10−7) and result in slowed L-type calcium channel inactivation and lower current density. In total, six pathogenic DNMs can alter calcium-related regulation of cardiomyocyte and neuronal excitability at a submembrane junction, suggesting a pathway conferring susceptibility to sudden death. There was a trend for excess LoF mutations in LoF intolerant genes, where ≥1 nonhealthy sample in denovo-db has a similar variant (odds ratio = 6.73, P = 0.02); additional uncharacterized genetic causes of sudden death in children might be discovered with larger cohorts.

Warsaw breakage syndrome (WABS), caused by bi-allelic variants in the DDX11 gene, is a rare cohesinopathy characterized by pre- and postnatal growth retardation, microcephaly, intellectual disability, facial dysmorphia, and sensorineural hearing loss due to cochlear hypoplasia. The DDX11 gene codes for an iron-sulfur DNA helicase in the Superfamily 2 helicases and plays an important role in genomic stability and maintenance. Fourteen individuals with WABS have been previously reported in the medical literature. Affected individuals have been of various ethnic backgrounds with different pathogenic variants. We report two unrelated individuals of Ashkenazi Jewish descent affected with WABS, who are homozygous for the c.1763-1G>C variant in the DDX11 gene. Their phenotype is consistent with previously reported individuals. RNA studies showed that this variant causes an alternative splice acceptor site leading to a frameshift in the open reading frame. Carrier screening of the c.1763-1G>C variant in the Jewish population revealed a high carrier frequency of 1 in 68 in the Ashkenazi Jewish population. Due to the high carrier frequency and the low number of affected individuals, we hypothesize a high rate of miscarriage of homozygous fetuses and/or subfertility for carrier couples. If the carrier frequency is reproducible in additional Ashkenazi Jewish populations, we suggest including DDX11 to Ashkenazi Jewish carrier screening panels.

PPP3CA encodes the catalytic subunit of calcineurin, a calcium-calmodulin-regulated serine-threonine phosphatase. Loss-of-function (LoF) variants in the catalytic domain have been associated with epilepsy, while gain-of-function (GoF) variants in the auto-inhibitory domain cause multiple congenital abnormalities. We herein report five new patients with de novo PPP3CA variants. Interestingly, the two frameshift variants in this study and the six truncating variants reported previously are all located within a 26-amino acid region in the regulatory domain (RD). Patients with a truncating variant had more severe earlier onset seizures compared to patients with a LoF missense variant, while autism spectrum disorder was a more frequent feature in the latter. Expression studies of a truncating variant showed apparent RNA expression from the mutant allele, but no detectable mutant protein. Our data suggest that PPP3CA truncating variants clustered in the RD, causing more severe early-onset refractory epilepsy and representing a type of variants distinct from LoF or GoF missense variants.

ARGONAUTE-2 and associated miRNAs form the RNA-induced silencing complex (RISC), which targets mRNAs for translational silencing and degradation as part of the RNA interference pathway. Despite the essential nature of this process for cellular function, there is little information on the role of RISC components in human development and organ function. We identify 13 heterozygous mutations in AGO2 in 21 patients affected by disturbances in neurological development. Each of the identified single amino acid mutations result in impaired shRNA-mediated silencing. We observe either impaired RISC formation or increased binding of AGO2 to mRNA targets as mutation specific functional consequences. The latter is supported by decreased phosphorylation of a C-terminal serine cluster involved in mRNA target release, increased formation of dendritic P-bodies in neurons and global transcriptome alterations in patient-derived primary fibroblasts. Our data emphasize the importance of gene expression regulation through the dynamic AGO2-RNA association for human neuronal development., AGO2 binds to miRNAs to repress expression of cognate target mRNAs. Here the authors report that heterozygous AGO2 mutations result in defects in neurological development and impair RNA interference.

The GATA zing finger domain-containing protein 2B (GATAD2B) gene encodes the p66beta protein, a subunit of the MeCP1-Mi2/ nucleosome remodeling and deacetylase complex, which is involved in transcription regulation and chromatin remodeling. Pathogenic variants in the GATAD2B gene have recently been associated with a recognizable neurodevelopmental syndrome, characterized by intellectual disability, limited speech, childhood hypotonia, and dysmorphic features. The majority of reported patients resulted from de novo loss of function (LOF) variants. We report a patient identified through whole exome sequencing analysis where a healthy mother was found to be low level mosaic for the pathogenic LOF variant found in her child, who is affected with GATAD2B-associated neurodevelopmental disorder (GAND). This variant was only found with the use of next generation sequencing technology in the mother and confirmed by digital droplet PCR. We summarize additional patients with GATAD2B LOF variants from a literature review and with our patient we contribute to delineate the phenotypic spectrum of GAND. We highlight the importance of detailed genetic testing, testing method, and counseling for cases of somatic mosaicism in an unaffected parent of children with GAND. This inheritance may be underreported and has a direct impact on reproductive planning and prenatal diagnosis.

Background Sudden death in children is a tragic event that often remains unexplained after comprehensive investigation. We report two asymptomatic siblings who died unexpectedly at approximately 1 year of age found to have biallelic (compound heterozygous) variants in PPA2. Methods The index case, parents, and sister were enrolled in the Sudden Unexplained Death in Childhood Registry and Research Collaborative, which included next‐generation genetic screening. Prior published cases of PPA2 variants, along with the known biology of PPA2, were also summarized. Results Whole exome sequencing in both siblings revealed biallelic rare missense variants in PPA2: c.182C > T (p.Ser61Phe) and c.380G > T (p.Arg127Leu). PPA2 encodes a mitochondrially located inorganic pyrophosphatase implicated in progressive and lethal cardiomyopathies. As a regulator and supplier of inorganic phosphate, PPA2 is central to phosphate metabolism. Biological roles include the following: mtDNA maintenance; oxidative phosphorylation and generation of ATP; reactive oxygen species homeostasis; mitochondrial membrane potential regulation; and possibly, retrograde signaling between mitochondria and nucleus. Conclusions Two healthy and asymptomatic sisters died unexpectedly at ages 12 and 10 months, and were diagnosed by molecular autopsy to carry biallelic variants in PPA2. Our cases add additional details to those reported thus far, and broaden the spectrum of clinical and molecular features of PPA2 variants., Two healthy and asymptomatic sisters died unexpectedly at ages 12 and 10 months, and were diagnosed by molecular autopsy to carry biallelic variants in PPA2. This gene encodes a mitochondrially located inorganic pyrophosphatase implicated in progressive and lethal cardiomyopathies.

We describe unrelated individuals with ichthyosis, failure to thrive, thrombocytopenia, photophobia, and progressive hearing loss. Each have bi-allelic mutations in AP1B1, the gene encoding the β subunit of heterotetrameric adaptor protein 1 (AP-1) complexes, which mediate endomembrane polarization, sorting, and transport. In affected keratinocytes the AP-1 β subunit is lost, and the γ subunit is greatly reduced, demonstrating destabilization of the AP-1 complex. Affected cells and tissue contain an abundance of abnormal vesicles and show hyperproliferation, abnormal epidermal differentiation, and derangement of intercellular junction proteins. Transduction of affected cells with wild-type AP1B1 rescues the vesicular phenotype, conclusively establishing that loss of AP1B1 function causes this disorder.

Background: Sudden unexplained death in childhood (SUDC) is an understudied public health issue with multiple potential risk factors, including undiagnosed cardiac and seizure disorders, for which genetic risk factors can contribute. Methods: We studied whole exome sequence data from 124 consecutively ascertained trios (decedent child and unaffected parents) to test for excessive de novo mutations in cardiac, epilepsy and other genes. Findings: Among SUDC decedents, nonsynonymous mutations were enriched in genes associated with cardiac and seizure disorders relative to non-sudden death controls (odds ratio=9.76, p=2.15x10-4). Beyond these genes there was a trend for excess loss-of-function mutations in genes that are loss-of-function intolerant (odds ratio=2.34, p=0.08). In our cohort, two genes had recurrent missense mutations, RYR2 and CACNA1C. Both RYR2 mutations are known as pathogenic (p=1.7x10-7). Both CACNA1C mutations are within a short (104 nucleotide) exon (p=1.0x10-7). We found evidence for over-transmission of loss-of-function or pathogenic variants within cardiac and seizure disorder genes (11/14 transmitted, p=0.03). Interpretation: Eleven of 124 (8.9%) decedents carried likely pathogenic or pathogenic variants associated with sudden death via a lethal cardiac or seizure disorder. Of these, 7 were de novo in origin, 1 was a transmitted parental mosaic mutation, and 3 were transmitted from heterozygous parents. Funding Statement: This study was supported by funding from the SUDC Foundation and Finding A Cure for Epilepsy and Seizures (NYU). RWT was supported by NIH grants R01 DA040484 and R01 MH71739. MJA and DJT were supported by the Mayo Clinic Windland Smith Rice Comprehensive Sudden Cardiac Death Program. Declaration of Interests: MJA is a consultant for Audentes Therapeutics, Boston Scientific, Gilead Sciences, Invitae, Medtronic, Myokardia, and St. Jude Medical. MJA and Mayo Clinic have an equity/royalty relationship with AliveCor, Blue Ox Health, and StemoniX. However, none of these entities are related to study. DG reports he is cofounder and holds equity in Q State – Pairnomix and Praxis Therapeutics, and equity in Apostle Inc. He receives personal fees from AstraZeneca, outside the submitted work. RT reports his participation on the Scientific Advisory boards of Howard Hughes Medical Institute, the Institute of Neuroscience (Shanghai), and the Institute of Science and Technology (Vienna). RT also reports his participation in the Actelion Scientific Advisory Board (Basel) with no compensation during 3 year span. OD reports equity interest in Empatica and receive funding from NINDS on SUDEP research. MJA, OD and PL are members of the SUDC Foundation’s scientific advisory board. LGC is President of the SUDC Foundation. MJA, OD, PL and LGC do not participate in grant funding decisions or negotiations by the SUDC Foundation to NYU Langone Health. PL, DMH, RR, XW, JP have no declarations of interests. Ethics Approval Statement: Parental consent was obtained in this NYU Institutional Review Board approved study (# S14-01061).

Purpose: To define the clinical characteristics of patients with variants in TCF20, we describe 27 patients, 26 of whom were identified via exome sequencing. We compare detailed clinical data with 17 previously reported patients. Methods: Patients were ascertained through molecular testing laboratories performing exome sequencing (and other testing) with orthogonal confirmation; collaborating referring clinicians provided detailed clinical information. Results: The cohort of 27 patients all had novel variants, and ranged in age from two to 68 years. All had developmental delay/intellectual disability. Autism spectrum disorders/autistic features were reported in 69%, attention disorders or hyperactivity in 67%, craniofacial features (no recognizable facial gestalt) in 67%, structural brain anomalies in 24%, and seizures in 12%. Additional features affecting various organ systems were described in 93%. In a majority of patients, we did not observe previously reported findings of postnatal overgrowth or craniosynostosis, in comparison to earlier reports. Conclusion: We provide valuable data regarding the prognosis and clinical manifestations of patients with variants in TCF20.

The ELISpot, a heterogeneous immunoassay, is widely used for detection of low abundant analytes. It is a reliable and robust assay to monitor responses of the immune system at the single-cell level by capturing secreted molecules of interest with specific, membrane-bound antibodies. Those molecules are then made visible by a cascade of ELISA-related development steps. The final results are distinct spots on the membrane as an imprint of the cell secreting the captured molecules, not only allowing their quantification but also providing insight on the kinetics and strength of secretion. This chapter describes the optimized protocol steps of the ELISpot technique, important improvements and tools available for the community, and the current expansion of the technique into polyfunctional cell analysis.

The ELISpot, a heterogeneous immunoassay, is widely used for detection of low abundant analytes. It is a reliable and robust assay to monitor responses of the immune system at the single-cell level by capturing secreted molecules of interest with specific, membrane-bound antibodies. Those molecules are then made visible by a cascade of ELISA-related development steps. The final results are distinct spots on the membrane as an imprint of the cell secreting the captured molecules, not only allowing their quantification but also providing insight on the kinetics and strength of secretion. This chapter describes the optimized protocol steps of the ELISpot technique, important improvements and tools available for the community, and the current expansion of the technique into polyfunctional cell analysis.