Search

Asthenoteratozoospermia characterized by multiple morphological abnormalities of the flagella (MMAF) has been identified as a sub-type of male infertility. Recent progress has identified several MMAF-associated genes with an autosomal recessive inheritance in human affected individuals, but the etiology in approximately 40% of affected individuals remains unknown. Here, we conducted whole-exome sequencing (WES) and identified hemizygous missense variants in the X-linked CFAP47 in three unrelated Chinese individuals with MMAF. These three CFAP47 variants were absent in human control population genome databases and were predicted to be deleterious by multiple bioinformatic tools. CFAP47 encodes a cilia- and flagella-associated protein that is highly expressed in testis. Immunoblotting and immunofluorescence assays revealed obviously reduced levels of CFAP47 in spermatozoa from all three men harboring deleterious missense variants of CFAP47. Furthermore, WES data from an additional cohort of severe asthenoteratozoospermic men originating from Australia permitted the identification of a hemizygous Xp21.1 deletion removing the entire CFAP47 gene. All men harboring hemizygous CFAP47 variants displayed typical MMAF phenotypes. We also generated a Cfap47-mutated mouse model, the adult males of which were sterile and presented with reduced sperm motility and abnormal flagellar morphology and movement. However, fertility could be rescued by the use of intra-cytoplasmic sperm injections (ICSIs). Altogether, our experimental observations in humans and mice demonstrate that hemizygous mutations in CFAP47 can induce X-linked MMAF and asthenoteratozoospermia, for which good ICSI prognosis is suggested. These findings will provide important guidance for genetic counseling and assisted reproduction treatments.

IntroductionDe novo mutations (DNMs) are known to play a prominent role in sporadic disorders with reduced fitness1. We hypothesize that DNMs play an important role in male infertility and explain a significant fraction of the genetic causes of this understudied disorder. To test this hypothesis, we performed trio-based exome-sequencing in a unique cohort of 185 infertile males and their unaffected parents. Following a systematic analysis, 29 of 145 rare protein altering DNMs were classified as possibly causative of the male infertility phenotype. We observed a significant enrichment of Loss-of-Function (LoF) DNMs in LoF-intolerant genes (p-value=1.00×10-5) as well as predicted pathogenic missense DNMs in missense-intolerant genes (p-value=5.01×10-4). One DNM gene identified, RBM5, is an essential regulator of male germ cell pre-mRNA splicing2. In a follow-up study, 5 rare pathogenic missense mutations affecting this gene were observed in a cohort of 2,279 infertile patients, with no such mutations found in a cohort of 5,784 fertile men (p-value=0.009). Our results provide the first evidence for the role of DNMs in severe male infertility and point to many new candidate genes affecting fertility.

STUDY QUESTION What are the causative genetic variants in patients with male infertility due to severe sperm motility disorders? SUMMARY ANSWER We identified high confidence disease-causing variants in multiple genes previously associated with severe sperm motility disorders in 10 out of 21 patients (48%) and variants in novel candidate genes in seven additional patients (33%). WHAT IS KNOWN ALREADY Severe sperm motility disorders are a form of male infertility characterised by immotile sperm often in combination with a spectrum of structural abnormalities of the sperm flagellum that do not affect viability. Currently, depending on the clinical sub-categorisation, up to 50% of causality in patients with severe sperm motility disorders can be explained by pathogenic variants in at least 22 genes. STUDY DESIGN, SIZE, DURATION We performed exome sequencing in 21 patients with severe sperm motility disorders from two different clinics. PARTICIPANTS/MATERIALS, SETTING, METHOD Two groups of infertile men, one from Argentina (n = 9) and one from Australia (n = 12), with clinically defined severe sperm motility disorders (motility MAIN RESULTS AND ROLE OF CHANCE In 10 of 21 patients (48%), we identified pathogenic variants in known sperm assembly genes: CFAP43 (3 patients); CFAP44 (2 patients), CFAP58 (1 patient), QRICH2 (2 patients), DNAH1 (1 patient) and DNAH6 (1 patient). The diagnostic rate did not differ markedly between the Argentinian and the Australian cohort (55% and 42%, respectively). Furthermore, we identified patients with variants in the novel human candidate sperm motility genes: DNAH12, DRC1, MDC1, PACRG, SSPL2C and TPTE2. One patient presented with variants in four candidate genes and it remains unclear which variants were responsible for the severe sperm motility defect in this patient. LARGE SCALE DATA N/A LIMITATIONS, REASONS FOR CAUTION In this study, we described patients with either a homozygous or two heterozygous candidate pathogenic variants in genes linked to sperm motility disorders. Due to unavailability of parental DNA, we have not assessed the frequency of de novo or maternally inherited dominant variants and could not determine the parental origin of the mutations to establish in all cases that the mutations are present on both alleles. WIDER IMPLICATIONS OF THE FINDINGS Our results confirm the likely causal role of variants in six known genes for sperm motility and we demonstrate that exome sequencing is an effective method to diagnose patients with severe sperm motility disorders (10/21 diagnosed; 48%). Furthermore, our analysis revealed six novel candidate genes for severe sperm motility disorders. Genome-wide sequencing of additional patient cohorts and re-analysis of exome data of currently unsolved cases may reveal additional variants in these novel candidate genes. STUDY FUNDING/COMPETING INTEREST(S) This project was supported in part by funding from the Australian National Health and Medical Research Council (APP1120356) to M.K.O.B., J.A.V. and R.I.M.L., The Netherlands Organisation for Scientific Research (918-15-667) to J.A.V., the Royal Society and Wolfson Foundation (WM160091) to J.A.V., as well as an Investigator Award in Science from the Wellcome Trust (209451) to J.A.V. and Grants from the National Research Council of Argentina (PIP 0900 and 4584) and ANPCyT (PICT 9591) to H.E.C. and a UUKi Rutherford Fund Fellowship awarded to B.J.H.

Background Circular RNAs (circRNAs) are predominantly derived from protein coding genes, and some can act as microRNA sponges or transcriptional regulators. Changes in circRNA levels have been identified during human development which may be functionally important, but lineage-specific analyses are currently lacking. To address this, we performed RNAseq analysis of human embryonic stem (ES) cells differentiated for 90 days towards 3D laminated retina. Results A transcriptome-wide increase in circRNA expression, size, and exon count was observed, with circRNA levels reaching a plateau by day 45. Parallel statistical analyses, controlling for sample and locus specific effects, identified 239 circRNAs with expression changes distinct from the transcriptome-wide pattern, but these all also increased in abundance over time. Surprisingly, circRNAs derived from long non-coding RNAs (lncRNAs) were found to account for a significantly larger proportion of transcripts from their loci of origin than circRNAs from coding genes. The most abundant, circRMST:E12-E6, showed a > 100X increase during differentiation accompanied by an isoform switch, and accounts for > 99% of RMST transcripts in many adult tissues. The second most abundant, circFIRRE:E10-E5, accounts for > 98% of FIRRE transcripts in differentiating human ES cells, and is one of 39 FIRRE circRNAs, many of which include multiple unannotated exons. Conclusions Our results suggest that during human ES cell differentiation, changes in circRNA levels are primarily globally controlled. They also suggest that RMST and FIRRE, genes with established roles in neurogenesis and topological organisation of chromosomal domains respectively, are processed as circular lncRNAs with only minor linear species. Electronic supplementary material The online version of this article (10.1186/s12864-018-4660-7) contains supplementary material, which is available to authorized users.