Consensus Recommendations to Optimize the Detection and Reporting of NTRK Gene Fusions by RNA-Based Next-Generation Sequencing

The detection of gene fusions by RNA-based next-generation sequencing (NGS) is an emerging method in clinical genetic laboratories for oncology biomarker testing to direct targeted therapy selections. A recent Canadian study (CANTRK study) comparing the detection of NTRK gene fusions on different NGS assays to determine subjects’ eligibility for tyrosine kinase TRK inhibitor therapy identified the need for recommendations for best practices for laboratory testing to optimize RNA-based NGS gene fusion detection. To develop consensus recommendations, representatives from 17 Canadian genetic laboratories participated in working group discussions and the completion of survey questions about RNA-based NGS. Consensus recommendations are presented for pre-analytic, analytic and reporting aspects of gene fusion detection by RNA-based NGS.

CANTRK: A Canadian Ring Study to Optimize Detection of NTRK Gene Fusions by Next-Generation RNA Sequencing

The Canadian NTRK (CANTRK) study is an interlaboratory comparison ring study to optimize testing for neurotrophic receptor tyrosine kinase (NTRK) fusions in Canadian laboratories. Sixteen diagnostic laboratories used next-generation sequencing (NGS) for NTRK1, NTRK2, or NTRK3 fusions. Each laboratory received 12 formalin-fixed, paraffin-embedded tumor samples with unique NTRK fusions and two control non-NTRK fusion samples (one ALK and one ROS1). Laboratories used validated protocols for NGS fusion detection. Panels included Oncomine Comprehensive Assay v3, Oncomine Focus Assay, Oncomine Precision Assay, AmpliSeq for Illumina Focus, TruSight RNA Pan-Cancer Panel, FusionPlex Lung, and QIAseq Multimodal Lung. One sample was withdrawn from analysis because of sample quality issues. Of the remaining 13 samples, 6 of 11 NTRK fusions and both control fusions were detected by all laboratories. Two fusions, WNK2::NTRK2 and STRN3::NTRK2, were not detected by 10 laboratories using the Oncomine Comprehensive or Focus panels, due to absence of WNK2 and STRN3 in panel designs. Two fusions, TPM3::NTRK1 and LMNA::NTRK1, were challenging to detect on the AmpliSeq for Illumina Focus panel because of bioinformatics issues. One ETV6::NTRK3 fusion at low levels was not detected by two laboratories using the TruSight Pan-Cancer Panel. Panels detecting all fusions included FusionPlex Lung, Oncomine Precision, and QIAseq Multimodal Lung. The CANTRK study showed competency in detection of NTRK fusions by NGS across different panels in 16 Canadian laboratories and identified key test issues as targets for improvements.

Evaluation of DNA Methylation Array for Glioma Tumor Profiling and Description of a Novel Epi-Signature to Distinguish IDH1/IDH2 Mutant and Wild-Type Tumors

Molecular biomarkers, such as IDH1/IDH2 mutations and 1p19q co-deletion, are included in the histopathological and clinical criteria currently used to diagnose and classify gliomas. IDH1/IDH2 mutation is a common feature of gliomas and is associated with a glioma-CpG island methylator phenotype (CIMP). Aberrant genomic methylation patterns can also be used to extrapolate information about copy number variation in a tumor. This project's goal was to assess the feasibility of DNA methylation array for the simultaneous detection of glioma biomarkers as a more effective testing strategy compared to existing single analyte tests.

METHODS

Whole-genome methylation array (WGMA) testing was performed using 48 glioma DNA samples to detect methylation aberrations and chromosomal gains and losses. The analyzed samples include 39 tumors in the discovery cohort and 9 tumors in the replication cohort. Methylation profiles for each sample were correlated with IDH1 p.R132G mutation, immunohistochemistry (IHC), and previous 1p19q clinical testing to assess the sensitivity and specificity of the WGMA assay for the detection of these variants.

RESULTS

We developed a DNA methylation signature to specifically distinguish a IDH1/IDH2 mutant tumor from normal samples. This signature is composed of 11 CpG sites that were significantly hypermethylated in the IDH1/IDH2 mutant group. Copy number analysis using WGMA data was able to identify five of five positive samples for 1p19q co-deletion and was concordant for all negative samples.

CONCLUSIONS

The DNA methylation signature presented here has the potential to refine the utility of WGMA to predict IDH1/IDH2 mutation status of gliomas, thus improving diagnostic yield and efficiency of laboratory testing compared to single analyte IDH1/IDH2 or 1p19q tests.

Practice guidelines for BRCA1/2 tumour testing in ovarian cancer

The purpose of this document is to provide pre-analytical, analytical and post-analytical considerations and recommendations to Canadian clinical laboratories developing, validating and offering next-generation sequencing (NGS)-based BRCA1 and BRCA2 (BRCA1/2) tumour testing in ovarian cancers. This document was drafted by the members of the Canadian College of Medical Geneticists (CCMG) somatic BRCA Ad Hoc Working Group, and representatives from the Canadian Association of Pathologists. The document was circulated to the CCMG members for comment. Following incorporation of feedback, this document has been approved by the CCMG board of directors. The CCMG is a Canadian organisation responsible for certifying medical geneticists and clinical laboratory geneticists, and for establishing professional and ethical standards for clinical genetics services in Canada. The current CCMG Practice Guidelines were developed as a resource for clinical laboratories in Canada; however, they are not inclusive of all information laboratories should consider in the validation and use of NGS for BRCA1/2 tumour testing in ovarian cancers.

Multisite verification of the accuracy of a multi-gene next generation sequencing panel for detection of mutations and copy number alterations in solid tumours

Molecular variants including single nucleotide variants (SNVs), copy number variants (CNVs) and fusions can be detected in the clinical setting using deep targeted sequencing. These assays support low limits of detection using little genomic input material. They are gaining in popularity in clinical laboratories, where sample volumes are limited, and low variant allele fractions may be present. However, data on reproducibility between laboratories is limited. Using a ring study, we evaluated the performance of 7 Ontario laboratories using targeted sequencing panels. All laboratories analysed a series of control and clinical samples for SNVs/CNVs and gene fusions. High concordance was observed across laboratories for measured CNVs and SNVs. Over 97% of SNV calls in clinical samples were detected by all laboratories. Whilst only a single CNV was detected in the clinical samples tested, all laboratories were able to reproducibly report both the variant and copy number. Concordance for information derived from RNA was lower than observed for DNA, due largely to decreased quality metrics associated with the RNA components of the assay, suggesting that the RNA portions of comprehensive NGS assays may be more vulnerable to variations in approach and workflow. Overall the results of this study support the use of the OFA for targeted sequencing for testing of clinical samples and suggest specific internal quality metrics that can be reliable indicators of assay failure. While we believe this evidence can be interpreted to support deep targeted sequencing in general, additional studies should be performed to confirm this.

A Pan-Canadian Validation Study for the Detection of EGFR T790M Mutation Using Circulating Tumor DNA From Peripheral Blood

INTRODUCTION

Genotyping circulating tumor DNA (ctDNA) is a promising noninvasive clinical tool to identify the EGFR T790M resistance mutation in patients with advanced NSCLC with resistance to EGFR inhibitors. To facilitate standardization and clinical adoption of ctDNA testing across Canada, we developed a 2-phase multicenter study to standardize T790M mutation detection using plasma ctDNA testing.

METHODS

In phase 1, commercial reference standards were distributed to participating clinical laboratories, to use their existing platforms for mutation detection. Baseline performance characteristics were established using known and blinded engineered plasma samples spiked with predetermined concentrations of T790M, L858R, and exon 19 deletion variants. In phase II, peripheral blood collected from local patients with known EGFR activating mutations and progressing on treatment were assayed for the presence of EGFR variants and concordance with a clinically validated test at the reference laboratory.

RESULTS

All laboratories in phase 1 detected the variants at 0.5 % and 5.0 % allele frequencies, with no false positives. In phase 2, the concordance with the reference laboratory for detection of both the primary and resistance mutation was high, with next-generation sequencing and droplet digital polymerase chain reaction exhibiting the best overall concordance. Data also suggested that the ability to detect mutations at clinically relevant limits of detection is generally not platform-specific, but rather impacted by laboratory-specific practices.

CONCLUSIONS

Discrepancies among sending laboratories using the same assay suggest that laboratory-specific practices may impact performance. In addition, a negative or inconclusive ctDNA test should be followed by tumor testing when possible.

DNA Methylation of the Oxytocin Receptor Across Neurodevelopmental Disorders

Many neurodevelopmental disorders (NDDs) share common learning and behavioural impairments, as well as features such as dysregulation of the oxytocin hormone. Here, we examined DNA methylation (DNAm) in the 1st intron of the oxytocin receptor gene, OXTR, in patients with autism spectrum (ASD), attention deficit and hyperactivity (ADHD) and obsessive compulsive (OCD) disorders. DNAm of OXTR was assessed for cohorts of ASD (blood), ADHD (saliva), OCD (saliva), which uncovered sex-specific DNAm differences compared to neurotypical, tissue-matched controls. Individuals with ASD or ADHD exhibiting extreme DNAm values had lower IQ and more social problems, respectively, than those with DNAm within normative ranges. This suggests that OXTR DNAm patterns are altered across NDDs and may be correlated with common clinical outcomes.

ATM whole gene deletion in an Italian family with hereditary pancreatic cancer: Challenges to cancer risk prediction associated with an 11q22.3 microdeletion

Hereditary pancreatic cancer has been attributed to variants of several cancer predisposition genes including ATM. While heterozygous pathogenic variants in the ATM gene are implicated as a cause of familial breast and pancreatic cancers to our knowledge ATM whole gene deletions have not been previously reported. We describe a contiguous gene deletion of the ATM locus in a multi-generation family of Italian descent with a strong family history of pancreatic cancer. A deletion of one copy of the entire ATM gene was identified by routine panel testing and further characterized by chromosomal microarray analysis. An 11q22.3 microdeletion of approximately 960 kb was identified that is predicted to result in loss of 10 genes including ATM. The deletion was identified in two additional family members including a presymptomatic daughter and an affected sibling. A normal disomic complement of the 11q22.3 region was detected in a third family member with a history of prostate and pancreatic cancer. Additional family members were not available for testing. Given available evidence that ATM haploinsufficiency can increase cancer risk, we predict that the observed copy number loss has likely contributed to hereditary cancer in this family. However, absence of the familial microdeletion in at least one affected family member suggests that ATM deletions are unlikely the sole contributing factor influencing tumor development in affected individuals. This case highlights 11q22.3 microdeletions of the ATM gene region as a possible risk factor for hereditary cancer, including pancreatic cancer. The same case provides a further cautionary tale for over interpretation of cancer risk associated tumor suppressor microdeletions and suggests that the variant may not be sufficient for tumor development or may modify the cancer risks associated with other, yet unidentified hereditary cancer genes.

Genome-Wide DNA Methylation Analysis of Chinese Patients with Systemic Lupus Erythematosus Identified Hypomethylation in Genes Related to the Type I Interferon Pathway

Background

Epigenetic variants have been shown in recent studies to be important contributors to the pathogenesis of systemic lupus erythematosus (SLE). Here, we report a 2-step study of discovery followed by replication to identify DNA methylation alterations associated with SLE in a Chinese population. Using a genome-wide DNA methylation microarray, the Illumina Infinium HumanMethylation450 BeadChip, we compared the methylation levels of CpG sites in DNA extracted from white blood cells from 12 female Chinese SLE patients and 10 healthy female controls.

Results

We identified 36 CpG sites with differential loss of DNA methylation and 8 CpG sites with differential gain of DNA methylation, representing 25 genes and 7 genes, respectively. Surprisingly, 42% of the hypomethylated CpG sites were located in CpG shores, which indicated the functional importance of the loss of DNA methylation. Microarray results were replicated in another cohort of 100 SLE patients and 100 healthy controls by performing bisulfite pyrosequencing of four hypomethylated genes, MX1, IFI44L, NLRC5 and PLSCR1. In addition, loss of DNA methylation in these genes was associated with an increase in mRNA expression. Gene ontology analysis revealed that the hypomethylated genes identified in the microarray study were overrepresented in the type I interferon pathway, which has long been implicated in the pathogenesis of SLE.

Conclusion

Our epigenetic findings further support the importance of the type I interferon pathway in SLE pathogenesis. Moreover, we showed that the DNA methylation signatures of SLE can be defined in unfractionated white blood cells.

An Update on Molecular Diagnostic Testing of Human Imprinting Disorders

Imprinted genes are expressed in a parent of origin manner. Dysregulation of imprinted genes expression causes various disorders associated with abnormalities of growth, neurodevelopment, and metabolism. Molecular mechanisms leading to imprinting disorders and strategies for their diagnosis are discussed in this review article.

NSD1 mutations generate a genome-wide DNA methylation signature

Sotos syndrome (SS) represents an important human model system for the study of epigenetic regulation; it is an overgrowth/intellectual disability syndrome caused by mutations in a histone methyltransferase, NSD1. As layered epigenetic modifications are often interdependent, we propose that pathogenic NSD1 mutations have a genome-wide impact on the most stable epigenetic mark, DNA methylation (DNAm). By interrogating DNAm in SS patients, we identify a genome-wide, highly significant NSD1^+/−-specific signature that differentiates pathogenic NSD1 mutations from controls, benign NSD1 variants and the clinically overlapping Weaver syndrome. Validation studies of independent cohorts of SS and controls assigned 100% of these samples correctly. This highly specific and sensitive NSD1^+/− signature encompasses genes that function in cellular morphogenesis and neuronal differentiation, reflecting cardinal features of the SS phenotype. The identification of SS-specific genome-wide DNAm alterations will facilitate both the elucidation of the molecular pathophysiology of SS and the development of improved diagnostic testing.

Sotos syndrome is an growth syndrome characterized by advanced growth in childhood, characteristic facial appearance and intellectual disability. Here the authors identify a genome-wide DNA methylation signature that accurately diagnoses Sotos Syndrome and distinguishes it from similar conditions.

A CGG-repeat expansion mutation in ZNF713 causes FRA7A: association with autistic spectrum disorder in two families

We report de novo occurrence of the 7p11.2 folate-sensitive fragile site FRA7A in a male with an autistic spectrum disorder (ASD) due to a CGG-repeat expansion mutation (∼450 repeats) in a 5' intron of ZNF713. This expanded allele showed hypermethylation of the adjacent CpG island with reduced ZNF713 expression observed in a proband-derived lymphoblastoid cell line (LCL). His unaffected mother carried an unmethylated premutation (85 repeats). This CGG-repeat showed length polymorphism in control samples (five to 22 repeats). In a second unrelated family, three siblings with ASD and their unaffected father were found to carry FRA7A premutations, which were partially or mosaically methylated. In one of the affected siblings, mitotic instability of the premutation was observed. ZNF713 expression in LCLs in this family was increased in three of these four premutation carriers. A firm link cannot yet be established between ASD and the repeat expansion mutation but plausible pathogenic mechanisms are discussed.

Multilocus loss of DNA methylation in individuals with mutations in the histone H3 lysine 4 demethylase KDM5C

BACKGROUND

A number of neurodevelopmental syndromes are caused by mutations in genes encoding proteins that normally function in epigenetic regulation. Identification of epigenetic alterations occurring in these disorders could shed light on molecular pathways relevant to neurodevelopment.

RESULTS

Using a genome-wide approach, we identified genes with significant loss of DNA methylation in blood of males with intellectual disability and mutations in the X-linked KDM5C gene, encoding a histone H3 lysine 4 demethylase, in comparison to age/sex matched controls. Loss of DNA methylation in such individuals is consistent with known interactions between DNA methylation and H3 lysine 4 methylation. Further, loss of DNA methylation at the promoters of the three top candidate genes FBXL5, SCMH1, CACYBP was not observed in more than 900 population controls. We also found that DNA methylation at these three genes in blood correlated with dosage of KDM5C and its Y-linked homologue KDM5D. In addition, parallel sex-specific DNA methylation profiles in brain samples from control males and females were observed at FBXL5 and CACYBP.

CONCLUSIONS

We have, for the first time, identified epigenetic alterations in patient samples carrying a mutation in a gene involved in the regulation of histone modifications. These data support the concept that DNA methylation and H3 lysine 4 methylation are functionally interdependent. The data provide new insights into the molecular pathogenesis of intellectual disability. Further, our data suggest that some DNA methylation marks identified in blood can serve as biomarkers of epigenetic status in the brain.

Discovery of cross-reactive probes and polymorphic CpGs in the Illumina Infinium HumanMethylation450 microarray

DNA methylation, an important type of epigenetic modification in humans, participates in crucial cellular processes, such as embryonic development, X-inactivation, genomic imprinting and chromosome stability. Several platforms have been developed to study genome-wide DNA methylation. Many investigators in the field have chosen the Illumina Infinium HumanMethylation microarray for its ability to reliably assess DNA methylation following sodium bisulfite conversion. Here, we analyzed methylation profiles of 489 adult males and 357 adult females generated by the Infinium HumanMethylation450 microarray. Among the autosomal CpG sites that displayed significant methylation differences between the two sexes, we observed a significant enrichment of cross-reactive probes co-hybridizing to the sex chromosomes with more than 94% sequence identity. This could lead investigators to mistakenly infer the existence of significant autosomal sex-associated methylation. Using sequence identity cutoffs derived from the sex methylation analysis, we concluded that 6% of the array probes can potentially generate spurious signals because of co-hybridization to alternate genomic sequences highly homologous to the intended targets. Additionally, we discovered probes targeting polymorphic CpGs that overlapped SNPs. The methylation levels detected by these probes are simply the reflection of underlying genetic polymorphisms but could be misinterpreted as true signals. The existence of probes that are cross-reactive or of target polymorphic CpGs in the Illumina HumanMethylation microarrays can confound data obtained from such microarrays. Therefore, investigators should exercise caution when significant biological associations are found using these array platforms. A list of all cross-reactive probes and polymorphic CpGs identified by us are annotated in this paper.

Phenotypic spectrum associated with duplication of Xp11.22-p11.23 includes Autism Spectrum Disorder

Dup(X)(p11.22-p11.23) has been shown to be associated with intellectual disability (ID, also referred to as mental retardation). Here, we characterize a 4.64 Mb de novo duplication of the same Xp11.22-p11.23 ID region in a female, but for this reference case the diagnosis was Autism Spectrum Disorder (ASD). Besides ASD, she also had very persistent trichotillomania, anxiety symptoms and some non-specific dysmorphic features. We report the detailed clinical features, as well as refine the rearrangement breakpoints of this disease-associated copy number variation region, which encompasses more than 50 genes. We propose that in addition to ID, the phenotypic spectrum associated with dup(X)(p11.22-p11.23) can include ASD, language impairment, and/or other primary psychiatric disorders.

Assessment of methylation level prediction accuracy in methyl-DNA immunoprecipitation and sodium bisulfite based microarray platforms

In this study, we verified the accuracy of two array methods--methylated DNA immunoprecipitation coupled with CpG island microarrays (MeDIP-CGI-arrays) and sodium bisulfite conversion based microarrays (BC-arrays)--in predicting regional methylation levels as measured by pyrosequencing of bisulfite converted DNA (BC-pyrosequencing). To test the accuracy of these methods we used the Agilent Human CpG island and the Illumina HumanMethylation27 microarrays respectively, and compared microarray outputs to the data from targeted BC-pyrosequencing assays from several genomic regions of corresponding samples. We observed relatively high correlation with BC-pyrosequencing data for both array platforms, R = 0.87 for BC-Array and R = 0.79 for MeDIP-CGI array. However, MeDIP-CGI array were less reliable in predicting intermediate levels of DNA methylation. Several bioinformatics strategies, to ameliorate the performance of the MeDIP-CGI-Arrays did not improve the correlation with BC-pyrosequencing data. The high scalability, low cost and simpler analysis of BC-arrays, together with the recent extended coverage may make them a more versatile methylation analysis tool.

Isolation of MECP2-null Rett Syndrome patient hiPS cells and isogenic controls through X-chromosome inactivation

Rett syndrome (RTT) is a neurodevelopmental autism spectrum disorder that affects girls due primarily to mutations in the gene encoding methyl-CpG binding protein 2 (MECP2). The majority of RTT patients carry missense and nonsense mutations leading to a hypomorphic MECP2, while null mutations leading to the complete absence of a functional protein are rare. MECP2 is an X-linked gene subject to random X-chromosome inactivation resulting in mosaic expression of mutant MECP2. The lack of human brain tissue motivates the need for alternative human cellular models to study RTT. Here we report the characterization of a MECP2 mutation in a classic female RTT patient involving rearrangements that remove exons 3 and 4 creating a functionally null mutation. To generate human neuron models of RTT, we isolated human induced pluripotent stem (hiPS) cells from RTT patient fibroblasts. RTT-hiPS cells retained the MECP2 mutation, are pluripotent and fully reprogrammed, and retained an inactive X-chromosome in a nonrandom pattern. Taking advantage of the latter characteristic, we obtained a pair of isogenic wild-type and mutant MECP2 expressing RTT-hiPS cell lines that retained this MECP2 expression pattern upon differentiation into neurons. Phenotypic analysis of mutant RTT-hiPS cell-derived neurons demonstrated a reduction in soma size compared with the isogenic control RTT-hiPS cell-derived neurons from the same RTT patient. Analysis of isogenic control and mutant hiPS cell-derived neurons represents a promising source for understanding the pathogenesis of RTT and the role of MECP2 in human neurons.

A novel approach identifies new differentially methylated regions (DMRs) associated with imprinted genes

Imprinted genes are critical for normal human growth and neurodevelopment. They are characterized by differentially methylated regions (DMRs) of DNA that confer parent of origin-specific transcription. We developed a new strategy to identify imprinted gene-associated DMRs. Using genome-wide methylation profiling of sodium bisulfite modified DNA from normal human tissues of biparental origin, candidate DMRs were identified by selecting CpGs with methylation levels consistent with putative allelic differential methylation. In parallel, the methylation profiles of tissues of uniparental origin, i.e., paternally-derived androgenetic complete hydatidiform moles (AnCHMs), and maternally-derived mature cystic ovarian teratoma (MCT), were examined and then used to identify CpGs with parent of origin-specific DNA methylation. With this approach, we found known DMRs associated with imprinted genomic regions as well as new DMRs for known imprinted genes, NAP1L5 and ZNF597, and novel candidate imprinted genes. The paternally methylated DMR for one candidate, AXL, a receptor tyrosine kinase, was also validated in experiments with mouse embryos that demonstrated Axl was expressed preferentially from the maternal allele in a DNA methylation-dependent manner.

Autism spectrum disorders and epigenetics

OBJECTIVE

Current research suggests that the causes of autism spectrum disorders (ASD) are multifactorial and include both genetic and environmental factors. Several lines of evidence suggest that epigenetics also plays an important role in ASD etiology and that it might, in fact, integrate genetic and environmental influences to dysregulate neurodevelopmental processes. The objective of this review is to illustrate how epigenetic modifications that are known to alter gene expression without changing primary DNA sequence may play a role in the etiology of ASD.

METHOD

In this review, we summarize current knowledge about epigenetic modifications to genes and genomic regions possibly involved in the etiology of ASD.

RESULTS

Several genetic syndromes comorbid with ASD, which include Rett, Fragile X, Prader-Willi, Angelman, and CHARGE (Coloboma of the eye, Heart defects, Atresia of the nasal choanae, Retardation of growth and/or development, Genital and/or urinary abnormalities, and Ear abnormalities and deafness), all demonstrate dysregulation of epigenetic marks or epigenetic mechanisms. We report also on genes or genomic regions exhibiting abnormal epigenetic regulation in association with either syndromic (15q11-13 maternal duplication) or nonsyndromic forms of ASD. Finally, we discuss the state of current knowledge regarding the etiologic role of environmental factors linked to both the development of ASD and epigenetic dysregulation.

CONCLUSION

Data reviewed in this article highlight a variety of situations in which epigenetic dysregulation is associated with the development of ASD, thereby supporting a role for epigenetics in the multifactorial etiologies of ASD.

Search for the sex-determining switch in monotremes: mapping WT1, SF1, LHX1, LHX2, FGF9, WNT4, RSPO1 and GATA4 in platypus

The duck-billed platypus has five pairs of sex chromosomes, but there is no information about the primary sex-determining switch in this species. As there is no apparent SRY orthologue in platypus, another gene must acquire the function of a key regulator of the gonadal male or female fate. SOX9 was ruled out from being this key regulator as it maps to an autosome in platypus. To check whether other genes in mammalian gonadogenesis could be the primary switch in monotremes, we have mapped a number of candidates in platypus. We report here the autosomal location of WT1, SF1, LHX1, LHX9, FGF9, WNT4 and RSPO1 in platypus, thus excluding these from being key regulators of sex determination in this species. We found that GATA4 maps to sex chromosomes Y1 and X2; however, it lies in the pairing region shown by chromosome painting to be homologous, so is unlikely to be either male-specific or differentially dosed in male and female.