Combined Use of Whole Exome Sequencing and CRISPR/Cas9 to Study the Etiology of Non-Obstructive Azoospermia: Demonstration of the Dispensable Role of the Testis-Specific Genes C1orf185 and CCT6B

The genetic landscape of male infertility is highly complex. It is estimated that at least 4000 genes are involved in human spermatogenesis, but only few have so far been extensively studied. In this study, we investigated by whole exome sequencing two cases of idiopathic non-obstructive azoospermia (NOA) due to severe hypospermatogenesis. After variant filtering and prioritizing, we retained for each patient a homozygous loss-of-function (LoF) variant in a testis-specific gene, C1orf185 (c.250C>T; p.Gln84Ter) and CCT6B (c.615-2A>G), respectively. Both variants are rare according to the gnomAD database and absent from our local control cohort (n = 445). To verify the implication of these candidate genes in NOA, we used the CRISPR/Cas9 system to invalidate the mouse orthologs 4930522H14Rik and Cct6b and produced two knockout (KO) mouse lines. Sperm and testis parameters of homozygous KO adult male mice were analyzed and compared with those of wild-type animals. We showed that homozygous KO males were fertile and displayed normal sperm parameters and a functional spermatogenesis. Overall, these results demonstrate that not all genes highly and specifically expressed in the testes are essential for spermatogenesis, and in particular, we conclude that bi-allelic variants of C1orf185 and CCT6B are most likely not to be involved in NOA and male fertility.

Whole exome sequencing identifies rare coding variants in novel human-mouse ortholog genes in African individuals diagnosed with non-syndromic hearing impairment

Physiologically, the human and murine hearing systems are very similar, justifying the extensive use of mice in experimental models for hearing impairment (HI). About 340 murine HI genes have been reported; however, whether variants in all human-mouse ortholog genes contribute to HI has been rarely investigated. In humans, nearly 120 HI genes have been identified to date, with GJB2 and GJB6 variants accounting for half of congenital HI cases, of genetic origin, in populations of European and Asian ancestries, but not in most African populations. The contribution of variants in other known genes of HI among the populations of African ancestry is poorly studied and displays the lowest pick-up rate. We used whole exome sequencing (WES) to investigate pathogenic and likely pathogenic (PLP) variants in 34 novel human-mouse orthologs HI genes, in 40 individuals from Cameroon and South Africa diagnosed with non-syndromic hearing impairment (NSHI), and compared the data to WES data of 129 ethnically matched controls. In addition, protein modeling for selected PLP gene variants, gene enrichment, and network analyses were performed. A total of 4/38 murine genes, d6wsu163e, zfp719, grp152 and minar2, had no human orthologs. WES identified three rare PLP variants in 3/34 human-mouse orthologs genes in three unrelated Cameroonian patients, namely: OCM2, c.227G>C p.(Arg76Thr) and LRGI1, c.1657G>A p.(Gly533Arg) in a heterozygous state, and a PLP variant MCPH1, c.2311C>G p.(Pro771Ala) in a homozygous state. In silico functional analyses suggest that these human-mouse ortholog genes functionally co-expressed interactions with well-established HI genes: GJB2 and GJB6. The study found one homozygous variant in MCPH1, likely to explain HI in one patient, and suggests that human-mouse ortholog variants could contribute to the understanding of the physiology of hearing in humans. Impact statement Despite, human and murine hearing system being very similar, the contribution of variants in relevant mouse-ortholog genes to hearing impairment (HI) has not been fully investigated. The contribution of variants in the known non-syndromic hearing impairment (NSHI) genes among Africans is poorly studied, suggesting that the novel gene(s) and mutations are yet to be discovered in NSHI in the African populations. Using whole exome sequencing (WES), this study identified rare candidate pathogenic and likely pathogenic (PLP) variants in 3/34 novel human-mouse ortholog genes in 3/40 individuals, with one homozygous variant, MCPH1, c.2311C>G p.(Pro771Ala), likely to explain HI in one patient. In silico functional analyses suggest that these human-mouse ortholog genes could contribute to the understanding of the physiology of hearing in humans and thus the variants identified in those genes deserve additional investigations.

A transcription factor collective defines the HSN serotonergic neuron regulatory landscape

Cell differentiation is controlled by individual transcription factors (TFs) that together activate a selection of enhancers in specific cell types. How these combinations of TFs identify and activate their target sequences remains poorly understood. Here, we identify the cis-regulatory transcriptional code that controls the differentiation of serotonergic HSN neurons in Caenorhabditis elegans. Activation of the HSN transcriptome is directly orchestrated by a collective of six TFs. Binding site clusters for this TF collective form a regulatory signature that is sufficient for de novo identification of HSN neuron functional enhancers. Among C. elegans neurons, the HSN transcriptome most closely resembles that of mouse serotonergic neurons. Mouse orthologs of the HSN TF collective also regulate serotonergic differentiation and can functionally substitute for their worm counterparts which suggests deep homology. Our results identify rules governing the regulatory landscape of a critically important neuronal type in two species separated by over 700 million years.

Distance is not everything in imaging genomics of functional networks: reply to a commentary on Correlated gene expression supports synchronous activity in brain networks

Our 2015 paper (Richiardi et al., 2015), showed that transcriptional similarity of gene expression level is higher than expected by chance within functional brain networks (defined by functional magnetic resonance imaging), a relationship that is driven by around 140 genes. These results were replicated in vivo in adolescents, where we showed that SNPs of these genes where associated above chance with in-vivo fMRI connectivity, and in the mouse, where mouse orthologs of our genes showed above-chance association with meso-scale axonal connectivity. This paper has received a commentary on biorXiv (Pantazatos and Li, 2016), making several claims about our results and methods, mainly pointing out that Euclidean distance explains our results (“…high within-network SF is entirely attributable to proximity and is unrelated to functional brain networks…”). Here we address these claims and their weaknesses, and show that our original results stand, contrary to the claims made in the commentary.