The Mouse Genome Database (MGD): facilitating mouse as a model for human biology and disease

Developmental coordination disorder: What can we learn from RI mice using motor learning tasks and QTL analysis

Developmental Coordination Disorder (DCD) is a neurodevelopmental disorder of unknown etiology that affects one in 20 children. There is an indication that DCD has an underlying genetic component due to its high heritability. Therefore, we explored the use of a recombinant inbred family of mice known as the BXD panel to understand the genetic basis of complex traits (i.e., motor learning) through identification of quantitative trait loci (QTLs). The overall aim of this study was to utilize the QTL approach to evaluate the genome-to-phenome correlation in BXD strains of mice in order to better understand the human presentation of DCD. Results of this current study confirm differences in motor learning in selected BXD strains and strains with altered cerebellar volume. Five strains - BXD15, BXD27, BXD28, BXD75, and BXD86 - exhibited the most DCD-like phenotype when compared with other BXD strains of interest. Results indicate that BXD15 and BXD75 struggled primarily with gross motor skills, BXD28 primarily had difficulties with fine motor skills, and BXD27 and BXD86 strains struggled with both fine and gross motor skills. The functional roles of genes within significant QTLs were assessed in relation to DCD-like behavior. Only Rab3a (Ras-related protein Rab-3A) emerged as a high likelihood candidate gene for the horizontal ladder rung task. This gene is associated with brain and skeletal muscle development, but lacked nonsynonymous polymorphisms. This study along with Gill et al. (same issue) is the first studies to specifically examine the genetic linkage of DCD using BXD strains of mice.

Sequence-based machine learning reveals 3D genome differences between bonobos and chimpanzees

Phenotypic divergence between closely related species, including bonobos and chimpanzees (genus Pan), is largely driven by variation in gene regulation. The 3D structure of the genome mediates gene expression; however, genome folding differences in Pan are not well understood. Here, we apply machine learning to predict genome-wide 3D genome contact maps from DNA sequence for 56 bonobos and chimpanzees, encompassing all five extant lineages. We use a pairwise approach to estimate 3D divergence between individuals from the resulting contact maps in 4,420 1 Mb genomic windows. While most pairs were similar, ∼17% were predicted to be substantially divergent in genome folding. The most dissimilar maps were largely driven by single individuals with rare variants that produce unique 3D genome folding in a region. We also identified 89 genomic windows where bonobo and chimpanzee contact maps substantially diverged, including several windows harboring genes associated with traits implicated in Pan phenotypic divergence. We used in silico mutagenesis to identify 51 3D-modifying variants in these bonobo-chimpanzee divergent windows, finding that 34 or 66.67% induce genome folding changes via CTCF binding motif disruption. Our results reveal 3D genome variation at the population-level and identify genomic regions where changes in 3D folding may contribute to phenotypic differences in our closest living relatives.

A pathogenic variant in the uncharacterized RNF212B gene results in severe aneuploidy male infertility and repeated IVF failure

Quantitative and qualitative spermatogenic impairments are major causes of men's infertility. Although in vitro fertilization (IVF) is effective, some couples persistently fail to conceive. To identify causal variants in patients with severe male infertility factor and repeated IVF failures, we sequenced the exome of two consanguineous family members who underwent several failed IVF cycles and were diagnosed with low sperm count and motility. We identified a rare homozygous nonsense mutation in a previously uncharacterized gene, RNF212B, as the causative variant. Recurrence was identified in another unrelated, infertile patient who also faced repeated failed IVF treatments. scRNA-seq demonstrated meiosis-specific expression of RNF212B. Sequence analysis located a protein domain known to be associated with aneuploidy, which can explain multiple IVF failures. Accordingly, FISH analysis revealed a high aneuploidy rate in the patients' sperm cells and their IVF embryos. Finally, inactivation of the Drosophila orthologs significantly reduced male fertility. Given that members of the evolutionary conserved RNF212 gene family are involved in meiotic recombination and crossover maturation, our findings indicate a critical role of RNF212B in meiosis, genome stability, and in human fertility. Since recombination is completely absent in Drosophila males, our findings may indicate an additional unrelated role for the RNF212-like paralogs in spermatogenesis.

Genome-wide association analyses define pathogenic signaling pathways and prioritize drug targets for IgA nephropathy

IgA nephropathy (IgAN) is a progressive form of kidney disease defined by glomerular deposition of IgA. Here we performed a genome-wide association study of 10,146 kidney-biopsy-diagnosed IgAN cases and 28,751 controls across 17 international cohorts. We defined 30 genome-wide significant risk loci explaining 11% of disease risk. A total of 16 loci were new, including TNFSF4/TNFSF18, REL, CD28, PF4V1, LY86, LYN, ANXA3, TNFSF8/TNFSF15, REEP3, ZMIZ1, OVOL1/RELA, ETS1, IGH, IRF8, TNFRSF13B and FCAR. The risk loci were enriched in gene orthologs causing abnormal IgA levels when genetically manipulated in mice. We also observed a positive genetic correlation between IgAN and serum IgA levels. High polygenic score for IgAN was associated with earlier onset of kidney failure. In a comprehensive functional annotation analysis of candidate causal genes, we observed convergence of biological candidates on a common set of inflammatory signaling pathways and cytokine ligand-receptor pairs, prioritizing potential new drug targets.

Biallelic variants in KCTD19 associated with male factor infertility and oligoasthenoteratozoospermia

STUDY QUESTION

Can whole-exome sequencing (WES) reveal new genetic factors responsible for male infertility characterized by oligozoospermia?

SUMMARY ANSWER

We identified biallelic missense variants in the Potassium Channel Tetramerization Domain Containing 19 gene (KCTD19) and confirmed it to be a novel pathogenic gene for male infertility.

WHAT IS KNOWN ALREADY

KCTD19 is a key transcriptional regulator that plays an indispensable role in male fertility by regulating meiotic progression. Kctd19 gene-disrupted male mice exhibit infertility due to meiotic arrest.

STUDY DESIGN, SIZE, DURATION

We recruited a cohort of 536 individuals with idiopathic oligozoospermia from 2014 to 2022 and focused on five infertile males from three unrelated families. Semen analysis data and ICSI outcomes were collected. WES and homozygosity mapping were performed to identify potential pathogenic variants. The pathogenicity of the identified variants was investigated in silico and in vitro.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Male patients diagnosed with primary infertility were recruited from the Reproductive and Genetic Hospital of CITIC-Xiangya. Genomic DNA extracted from affected individuals was used for WES and Sanger sequencing. Sperm phenotype, sperm nuclear maturity, chromosome aneuploidy, and sperm ultrastructure were assessed using hematoxylin and eosin staining and toluidine blue staining, FISH and transmission electron microscopy. The functional effects of the identified variants in HEK293T cells were investigated via western blotting and immunofluorescence.

MAIN RESULTS AND THE ROLE OF CHANCE

We identified three homozygous missense variants (NM_001100915, c.G628A:p.E210K, c.C893T:p.P298L, and c.G2309A:p.G770D) in KCTD19 in five infertile males from three unrelated families. Abnormal morphology of the sperm heads with immature nuclei and/or nuclear aneuploidy were frequently observed in individuals with biallelic KCTD19 variants, and ICSI was unable to rescue these deficiencies. These variants reduced the abundance of KCTD19 due to increased ubiquitination and impaired its nuclear colocalization with its functional partner, zinc finger protein 541 (ZFP541), in HEK293T cells.

LIMITATIONS, REASONS FOR CAUTION

The exact pathogenic mechanism remains unclear, and warrants further studies using knock-in mice that mimic the missense mutations found in individuals with biallelic KCTD19 variants.

WIDER IMPLICATIONS OF THE FINDINGS

Our study is the first to report a likely causal relationship between KCTD19 deficiency and male infertility, confirming the critical role of KCTD19 in human reproduction. Additionally, this study provided evidence for the poor ICSI clinical outcomes in individuals with biallelic KCTD19 variants, which may guide clinical treatment strategies.

STUDY FUNDING/COMPETING INTEREST(S)

This work was supported by the National Key Research and Developmental Program of China (2022YFC2702604 to Y.-Q.T.), the National Natural Science Foundation of China (81971447 and 82171608 to Y.-Q.T., 82101961 to C.T.), a key grant from the Prevention and Treatment of Birth Defects from Hunan Province (2019SK1012 to Y.-Q.T.), a Hunan Provincial Grant for Innovative Province Construction (2019SK4012), and the China Postdoctoral Science Foundation (2022M721124 to W.W.). The authors declare no conflicts of interest.

TRIAL REGISTRATION NUMBER

N/A.

Defined Pig Microbiota with a Potential Protective Effect against Infection with Salmonella Typhimurium

A balanced microbiota is a main prerequisite for the host's health. The aim of the present work was to develop defined pig microbiota (DPM) with the potential ability to protect piglets against infection with Salmonella Typhimurium, which causes enterocolitis. A total of 284 bacterial strains were isolated from the colon and fecal samples of wild and domestic pigs or piglets using selective and nonselective cultivation media. Isolates belonging to 47 species from 11 different genera were identified by MALDI-TOF mass spectrometry (MALDI-TOF MS). The bacterial strains for the DPM were selected for anti-Salmonella activity, ability to aggregate, adherence to epithelial cells, and to be bile and acid tolerant. The selected combination of 9 strains was identified by sequencing of the 16S rRNA gene as Bacillus sp., Bifidobacterium animalis subsp. lactis, B. porcinum, Clostridium sporogenes, Lactobacillus amylovorus, L. paracasei subsp. tolerans, Limosilactobacillus reuteri subsp. suis, and Limosilactobacillus reuteri (two strains) did not show mutual inhibition, and the mixture was stable under freezing for at least 6 months. Moreover, strains were classified as safe without pathogenic phenotype and resistance to antibiotics. Future experiments with Salmonella-infected piglets are needed to test the protective effect of the developed DPM.

Inheritance of perturbed methylation and metabolism caused by uterine malnutrition via oocytes

BACKGROUND

Undernourishment in utero has deleterious effects on the metabolism of offspring, but the mechanism of the transgenerational transmission of metabolic disorders is not well known. In the present study, we found that undernourishment in utero resulted in metabolic disorders of female F1 and F2 in mouse model.

RESULTS

Undernutrition in utero induced metabolic disorders of F1 females, which was transmitted to F2 females. The global methylation in oocytes of F1 exposed to undernutrition in utero was decreased compared with the control. KEGG analysis showed that genes with differential methylation regions (DMRs) in promoters were significantly enriched in metabolic pathways. The altered methylation of some DMRs in F1 oocytes located at the promoters of metabolic-related genes were partially observed in F2 tissues, and the expressions of these genes were also changed. Meanwhile, the abnormal DNA methylation of the validated DMRs in F1 oocytes was also observed in F2 oocytes.

CONCLUSIONS

These results indicate that DNA methylation may mediate the transgenerational inheritance of metabolic disorders induced by undernourishment in utero via female germline.

Exosomes induce neurogenesis of pluripotent P19 cells

Exosomes play a role in tissue/organ development and differentiation. Retinoic acid induces differentiation of P19 cells (UD-P19) to P19 neurons (P19N) that behave like cortical neurons and express characteristic neuronal genes such as NMDA receptor subunits. Here we report P19N exosome-mediated differentiation of UD-P19 to P19N. Both UD-P19 and P19N released exosomes with characteristic exosome morphology, size, and common protein markers. P19N internalized significantly higher number of Dil-P19N exosomes as compared to UD-P19 with accumulation in the perinuclear region. Continuous exposure of UD-P19 to P19N exosomes for six days induced formation of small-sized embryoid bodies that differentiated into MAP2-/GluN2B-positive neurons recapitulating RA-induction of neurogenesis. Incubation with UD-P19 exosomes for six days did not affect UD-P19. Small RNA-seq identified enrichment of P19N exosomes with pro-neurogenic non-coding RNAs (ncRNAs) such as miR-9, let-7, MALAT1 and depleted with ncRNAs involved in maintenance of stem cell characteristics. UD-P19 exosomes were rich with ncRNAs required for maintenance of stemness. P19N exosomes provide an alternative method to genetic modifications for cellular differentiation of neurons. Our novel findings on exosomes-mediated differentiation of UD-P19 to P19 neurons provide tools to study pathways directing neuron development/differentiation and develop novel therapeutic strategies in neuroscience.

Identification of ocular regulatory functions of core histone variant H3.2

The posttranscriptional modifications (PTM) of the Histone H3 family play an important role in ocular system differentiation. However, there has been no study on the nature of specific Histone H3 subtype carrying these modifications. Fortuitously, we had previously identified a dominant small-eye mutant Aey69 mouse with a mutation in the H3.2 encoding Hist2h3c1 gene (Vetrivel et al., 2019). In continuation, in the present study, the role of Histone H3.2 with relation to the microphtalmic Aey69 has been elaborated. Foremost, a transgenic mouse line expressing the fusion protein H3.2-GFP was generated using Crispr/Cas9. The approach was intended to confer a unique tag to the Hist2h3c1 gene which is similar in sequence and encoded protein structure to other histones. The GFP tag was then used for ChIP Seq analysis of the genes regulated by H3.2. The approach revealed ocular specific H3.2 targets including Ephrin family genes. Altered enrichment of H3.2 was found in the mutant Aey69 mouse, specifically around the ligand Efna5 and the receptor Ephb2. The effect of this altered enrichment on Ephrin signaling was further analysed by QPCR and immunohistochemistry. This study identifies Hist2h3c1 encoded H3.2 as an important epigenetic player in ocular development. By binding to specific regions of ocular developmental factors Histone H3.2 facilitates the function of these genes for successful early ocular development.

Hidden diversity of the genus Trinomys (Rodentia: Echimyidae): phylogenetic and populational structure analyses uncover putative new lineages

Trinomys, one of the most species-rich spiny rat genera in Brazil, is widely distributed in Caatinga, Cerrado and Atlantic Forest biomes, and currently includes ten recognized species, three of which are polytypic. Although some studies employing molecular data have been conducted to better characterize phylogenetic relationships among species, 19 nominal taxa have been suggested, implying considerable incongruence regarding species boundaries. We addressed this incongruence by intensively sampling all species across the geographic distribution of the genus. In addition to publicly available data, we generated 182 mt-Cytb gene sequences, and employed phylogenetic and computational species delimitation methods to obtain a clearer picture of the genus diversity. Moreover, we evaluated populational diversity within each accepted species, considering their geographical distribution and a timescale for the evolution of the genus. Beyond confirming the general patterns described for the evolution of the group, this new analysis suggests that Trinomys is comprised of at least 16 evolutionary lineages, 13 of them recognized as species or subspecies, and three never before characterized. This study highlights the importance of increased sample sizes and computational species delimitation methods in uncovering hidden diversity in Trinomys.

Cross-Disorder Analysis of Genic and Regulatory Copy Number Variations in Bipolar Disorder, Schizophrenia, and Autism Spectrum Disorder

BACKGROUND

We aimed to determine the similarities and differences in the roles of genic and regulatory copy number variations (CNVs) in bipolar disorder (BD), schizophrenia (SCZ), and autism spectrum disorder (ASD).

METHODS

Based on high-resolution CNV data from 8708 Japanese samples, we performed to our knowledge the largest cross-disorder analysis of genic and regulatory CNVs in BD, SCZ, and ASD.

RESULTS

In genic CNVs, we found an increased burden of smaller (<100 kb) exonic deletions in BD, which contrasted with the highest burden of larger (>500 kb) exonic CNVs in SCZ/ASD. Pathogenic CNVs linked to neurodevelopmental disorders were significantly associated with the risk for each disorder, but BD and SCZ/ASD differed in terms of the effect size (smaller in BD) and subtype distribution of CNVs linked to neurodevelopmental disorders. We identified 3 synaptic genes (DLG2, PCDH15, and ASTN2) as risk factors for BD. Whereas gene set analysis showed that BD-associated pathways were restricted to chromatin biology, SCZ and ASD involved more extensive and similar pathways. Nevertheless, a correlation analysis of gene set results indicated weak but significant pathway similarities between BD and SCZ or ASD (r = 0.25-0.31). In SCZ and ASD, but not BD, CNVs were significantly enriched in enhancers and promoters in brain tissue.

CONCLUSIONS

BD and SCZ/ASD differ in terms of CNV burden, characteristics of CNVs linked to neurodevelopmental disorders, and regulatory CNVs. On the other hand, they have shared molecular mechanisms, including chromatin biology. The BD risk genes identified here could provide insight into the pathogenesis of BD.

Comprehensive Assessment of Indian Variations in the Druggable Kinome Landscape Highlights Distinct Insights at the Sequence, Structure and Pharmacogenomic Stratum

India confines more than 17% of the world’s population and has a diverse genetic makeup with several clinically relevant rare mutations belonging to many sub-group which are undervalued in global sequencing datasets like the 1000 Genome data (1KG) containing limited samples for Indian ethnicity. Such databases are critical for the pharmaceutical and drug development industry where diversity plays a crucial role in identifying genetic disposition towards adverse drug reactions. A qualitative and comparative sequence and structural study utilizing variant information present in the recently published, largest curated Indian genome database (IndiGen) and the 1000 Genome data was performed for variants belonging to the kinase coding genes, the second most targeted group of drug targets. The sequence-level analysis identified similarities and differences among different populations based on the nsSNVs and amino acid exchange frequencies whereas a comparative structural analysis of IndiGen variants was performed with pathogenic variants reported in UniProtKB Humsavar data. The influence of these variations on structural features of the protein, such as structural stability, solvent accessibility, hydrophobicity, and the hydrogen-bond network was investigated. In-silico screening of the known drugs to these Indian variation-containing proteins reveals critical differences imparted in the strength of binding due to the variations present in the Indian population. In conclusion, this study constitutes a comprehensive investigation into the understanding of common variations present in the second largest population in the world and investigating its implications in the sequence, structural and pharmacogenomic landscape. The preliminary investigation reported in this paper, supporting the screening and detection of ADRs specific to the Indian population could aid in the development of techniques for pre-clinical and post-market screening of drug-related adverse events in the Indian population.

Creation and evaluation of full-text literature-derived, feature-weighted disease models of genetically determined developmental disorders

There are >2500 different genetically determined developmental disorders (DD), which, as a group, show very high levels of both locus and allelic heterogeneity. This has led to the wide-spread use of evidence-based filtering of genome-wide sequence data as a diagnostic tool in DD. Determining whether the association of a filtered variant at a specific locus is a plausible explanation of the phenotype in the proband is crucial and commonly requires extensive manual literature review by both clinical scientists and clinicians. Access to a database of weighted clinical features extracted from rigorously curated literature would increase the efficiency of this process and facilitate the development of robust phenotypic similarity metrics. However, given the large and rapidly increasing volume of published information, conventional biocuration approaches are becoming impractical. Here, we present a scalable, automated method for the extraction of categorical phenotypic descriptors from the full-text literature. Papers identified through literature review were downloaded and parsed using the Cadmus custom retrieval package. Human Phenotype Ontology terms were extracted using MetaMap, with 76-84% precision and 65-73% recall. Mean terms per paper increased from 9 in title + abstract, to 68 using full text. We demonstrate that these literature-derived disease models plausibly reflect true disease expressivity more accurately than widely used manually curated models, through comparison with prospectively gathered data from the Deciphering Developmental Disorders study. The area under the curve for receiver operating characteristic (ROC) curves increased by 5-10% through the use of literature-derived models. This work shows that scalable automated literature curation increases performance and adds weight to the need for this strategy to be integrated into informatic variant analysis pipelines. Database URL: https://doi.org/10.1093/database/baac038.

Bi-allelic variants in SHOC1 cause non-obstructive azoospermia with meiosis arrest in humans and mice

Meiosis is pivotal to gametogenesis and fertility. Meiotic recombination is a mandatory process that ensures faithful chromosome segregation and generates genetic diversity in gametes. Non-obstructive azoospermia (NOA) caused by meiotic arrest is a common cause of male infertility and has many genetic origins, including chromosome abnormalities, Y chromosome microdeletion and monogenic mutations. However, the genetic causes of the majority of NOA cases remain to be elucidated. Here, we report our findings of three Shortage in chiasmata 1 (SHOC1) bi-allelic variants in three NOA patients, of which two are homozygous for the same loss-of-function variant (c.231_232del: p. L78Sfs*9), and one is heterozygous for two different missense variants (c.1978G>A: p.A660T; c.4274G>A: p.R1425H). Testicular biopsy of one patient revealed impairment of spermatocyte maturation. Both germ-cell-specific and general Shoc1-knockout mice exhibited similar male infertility phenotypes. Subsequent analysis revealed comprehensive defects in homologous pairing and synapsis along with abnormal expression of DMC1, RAD51 and RPA2 in Shoc1-defective spermatocyte spreads. These findings imply that SHOC1 may have a presynaptic function during meiotic recombination apart from its previously identified role in crossover formation. Overall, our results provide strong evidence for the clinical relevance of SHOC1 mutations in patients with NOA and contribute to a deeper mechanistic understanding of the role of SHOC1 during meiotic recombination.

Modelling the genetic aetiology of complex disease: human-mouse conservation of noncoding features and disease-associated loci

Understanding the genetic aetiology of loci associated with a disease is crucial for developing preventative measures and effective treatments. Mouse models are used extensively to understand human pathobiology and mechanistic functions of disease-associated loci. However, the utility of mouse models is limited in part by evolutionary divergence in transcription regulation for pathways of interest. Here, we summarize the alignment of genomic (exonic and multi-cell regulatory) annotations alongside Mendelian and complex disease-associated variant sites between humans and mice. Our results highlight the importance of understanding evolutionary divergence in transcription regulation when interpreting functional studies using mice as models for human disease variants.

Invertebrate Model Organisms as a Platform to Investigate Rare Human Neurological Diseases

Patients suffering from rare human diseases often go through a painful journey for finding a definite molecular diagnosis prerequisite of appropriate cures. With a novel variant isolated from a single patient, determination of its pathogenicity to end such "diagnostic odyssey" requires multi-step processes involving experts in diverse areas of interest, including clinicians, bioinformaticians and research scientists. Recent efforts in building large-scale genomic databases and in silico prediction platforms have facilitated identification of potentially pathogenic variants causative of rare human diseases of a Mendelian basis. However, the functional significance of individual variants remains elusive in many cases, thus requiring incorporation of versatile and rapid model organism (MO)-based platforms for functional analyses. In this review, the current scope of rare disease research is briefly discussed. In addition, an overview of invertebrate MOs for their key features relevant to rare neurological diseases is provided, with the characteristics of two representative invertebrate MOs, Drosophila melanogaster and Caenorhabditis elegans, as well as the challenges against them. Finally, recently developed research networks integrating these MOs in collaborative research are portraited with an array of bioinformatical analyses embedded. A comprehensive survey of MO-based research activities provided in this review will help us to design a wellstructured analysis of candidate genes or potentially pathogenic variants for their roles in rare neurological diseases in future.

miRNA-Profiling in Ejaculated and Epididymal Pig Spermatozoa and Their Relation to Fertility after Artificial Insemination

Simple Summary

The present study searched for the presence and abundance of porcine spermatozoa small RNA sequences (microRNAs) that have the potential to alter gene expression patterns. Four different sperm sources were compared: spermatozoa from three different sections of the ejaculate and from the caudal epididymis, also classed as spermatozoa from higher (HF) or lower (LF) fertility boars. Sperm miRNAs were compared using high-output small RNA sequencing. We identified five sperm miRNAs not previously reported in pigs. Differences in abundance of four miRNAs known to affect the expression of genes with key roles in fertility were related to boar fertility. These miRNAs could be used as fertility markers in artificial insemination programs.

Abstract

MicroRNAs (miRNAs) are short non-coding RNAs (20–25 nucleotides in length) capable of regulating gene expression by binding -fully or partially- to the 3’-UTR of target messenger RNA (mRNA). To date, several studies have investigated the role of sperm miRNAs in spermatogenesis and their remaining presence toward fertilization and early embryo development. However, little is known about the miRNA cargo in the different sperm sources and their possible implications in boar fertility. Here, we characterized the differential abundance of miRNAs in spermatozoa from the terminal segment of the epididymis and three different fractions of the pig ejaculate (sperm-peak, sperm-rich, and post-sperm rich) comparing breeding boars with higher (HF) and lower (LF) fertility after artificial insemination (AI) using high-output small RNA sequencing. We identified five sperm miRNAs that, to our knowledge, have not been previously reported in pigs (mir-10386, mir-10390, mir-6516, mir-9788-1, and mir-9788-2). Additionally, four miRNAs (mir-1285, mir-92a, mir-34c, mir-30), were differentially expressed among spermatozoa sourced from ejaculate fractions and the cauda epididymis, and also different abundance was found between HF and LF groups in mir-182, mir-1285, mir-191, and mir-96. These miRNAs target genes with key roles in fertility, sperm survival, immune tolerance, or cell cycle regulation, among others. Linking the current findings with the expression of specific sperm proteins would help predict fertility in future AI-sires.

DNA methylation aging and transcriptomic studies in horses

Cytosine methylation patterns have not yet been thoroughly studied in horses. Here, we profile n = 333 samples from 42 horse tissue types at loci that are highly conserved between mammalian species using a custom array (HorvathMammalMethylChip40). Using the blood and liver tissues from horses, we develop five epigenetic aging clocks: a multi-tissue clock, a blood clock, a liver clock and two dual-species clocks that apply to both horses and humans. In addition, using blood methylation data from three additional equid species (plains zebra, Grevy's zebras and Somali asses), we develop another clock that applies across all equid species. Castration does not significantly impact the epigenetic aging rate of blood or liver samples from horses. Methylation and RNA data from the same tissues define the relationship between methylation and RNA expression across horse tissues. We expect that the multi-tissue atlas will become a valuable resource.

Genetically Modified Mouse Models of Congenital Diaphragmatic Hernia: Opportunities and Limitations for Studying Altered Lung Development

Congenital diaphragmatic hernia (CDH) is a relatively common and life-threatening birth defect, characterized by an abnormal opening in the primordial diaphragm that interferes with normal lung development. As a result, CDH is accompanied by immature and hypoplastic lungs, being the leading cause of morbidity and mortality in patients with this condition. In recent decades, various animal models have contributed novel insights into the pathogenic mechanisms underlying CDH and associated pulmonary hypoplasia. In particular, the generation of genetically modified mouse models, which show both diaphragm and lung abnormalities, has resulted in the discovery of multiple genes and signaling pathways involved in the pathogenesis of CDH. This article aims to offer an up-to-date overview on CDH-implicated transcription factors, molecules regulating cell migration and signal transduction as well as components contributing to the formation of extracellular matrix, whilst also discussing the significance of these genetic models for studying altered lung development with regard to the human situation.

Complementary evolution of coding and noncoding sequence underlies mammalian hairlessness

Body hair is a defining mammalian characteristic, but several mammals, such as whales, naked mole-rats, and humans, have notably less hair. To find the genetic basis of reduced hair quantity, we used our evolutionary-rates-based method, RERconverge, to identify coding and noncoding sequences that evolve at significantly different rates in so-called hairless mammals compared to hairy mammals. Using RERconverge, we performed a genome-wide scan over 62 mammal species using 19,149 genes and 343,598 conserved noncoding regions. In addition to detecting known and potential novel hair-related genes, we also discovered hundreds of putative hair-related regulatory elements. Computational investigation revealed that genes and their associated noncoding regions show different evolutionary patterns and influence different aspects of hair growth and development. Many genes under accelerated evolution are associated with the structure of the hair shaft itself, while evolutionary rate shifts in noncoding regions also included the dermal papilla and matrix regions of the hair follicle that contribute to hair growth and cycling. Genes that were top ranked for coding sequence acceleration included known hair and skin genes KRT2, KRT35, PKP1, and PTPRM that surprisingly showed no signals of evolutionary rate shifts in nearby noncoding regions. Conversely, accelerated noncoding regions are most strongly enriched near regulatory hair-related genes and microRNAs, such as mir205, ELF3, and FOXC1, that themselves do not show rate shifts in their protein-coding sequences. Such dichotomy highlights the interplay between the evolution of protein sequence and regulatory sequence to contribute to the emergence of a convergent phenotype.

Pathogenic variations in Germ Cell Nuclear Acidic Peptidase (GCNA) are associated with human male infertility

Infertility affects one in six couples, half of which are caused by a male factor. Male infertility can be caused by both, qualitative and quantitative defects, leading to Oligo- astheno-terato-zoospermia (OAT; impairment in ejaculate sperm cell concentration, motility and morphology). Azoospermia defined as complete absence of sperm cells in the ejaculation. While hundreds of genes are involved in spermatogenesis the genetic etiology of men's infertility remains incomplete.We identified a hemizygous stop gain pathogenic variation (PV) in the X-linked Germ Cell Nuclear Acidic Peptidase (GCNA), in an Azoospermic patient by exome sequencing. Assessment of the prevalence of pathogenic variations in this gene in infertile males by exome sequence data of 11 additional unrelated patients identified a probable hemizygous causative missense PV in GCNA in a severe OAT patient. Expression of GCNA in the patients' testes biopsies and the stage of spermatogonial developmental arrest were determined by immunofluorescence and immunohistochemistry. The Azoospermic patient presented spermatogenic maturation arrest with an almost complete absence of early and late primary spermatocytes and thus the complete absence of sperm. GCNA is critical for genome integrity and its loss results in genomic instability and infertility in Drosophila, C. elegans, zebrafish, and mouse. PVs in GCNA appear to be incompatible with male fertility in humans as well: A stop-gain PV caused Azoospermia and a missense PV caused severe OAT with very low fertilization rates and no pregnancy in numerous IVF treatments.

Homozygous duplication identified by whole genome sequencing causes LRBA deficiency

In more than one-third of primary immunodeficiency (PID) patients, extensive genetic analysis including whole-exome sequencing (WES) fails to identify the genetic defect. Whole-genome sequencing (WGS) is able to detect variants missed by other genomics platforms, enabling the molecular diagnosis of otherwise unresolved cases. Here, we report two siblings, offspring of consanguineous parents, who experienced similar severe events encompassing early onset of colitis, lymphoproliferation, and hypogammaglobulinemia, typical of lipopolysaccharide-responsive and beige-like anchor (LRBA) or cytotoxic T lymphocyte antigen 4 (CTLA4) deficiencies. Gene-panel sequencing, comparative genomic hybridization (CGH) array, and WES failed to reveal a genetic aberration in relevant genes. WGS of these patients detected a 12.3 kb homozygous tandem duplication that was absent in control cohorts and is predicted to disrupt the reading frame of the LRBA gene. The variant was validated by PCR and Sanger sequencing, demonstrating the presence of the junction between the reference and the tandem-duplicated sequence. Droplet digital PCR (ddPCR) further confirmed the copy number in the unaffected parents (CN = 3, heterozygous) and affected siblings (CN = 4, homozygous), confirming the expected segregation pattern. In cases of suspected inherited immunodeficiency, WGS may reveal a mutation when other methods such as microarray and WES analysis failed to detect an aberration.

Mouse Genome Informatics (MGI): latest news from MGD and GXD

The Mouse Genome Informatics (MGI) database system combines multiple expertly curated community data resources into a shared knowledge management ecosystem united by common metadata annotation standards. MGI's mission is to facilitate the use of the mouse as an experimental model for understanding the genetic and genomic basis of human health and disease. MGI is the authoritative source for mouse gene, allele, and strain nomenclature and is the primary source of mouse phenotype annotations, functional annotations, developmental gene expression information, and annotations of mouse models with human diseases. MGI maintains mouse anatomy and phenotype ontologies and contributes to the development of the Gene Ontology and Disease Ontology and uses these ontologies as standard terminologies for annotation. The Mouse Genome Database (MGD) and the Gene Expression Database (GXD) are MGI's two major knowledgebases. Here, we highlight some of the recent changes and enhancements to MGD and GXD that have been implemented in response to changing needs of the biomedical research community and to improve the efficiency of expert curation. MGI can be accessed freely at http://www.informatics.jax.org .

ELMO1 signaling is a promoter of osteoclast function and bone loss

Osteoporosis affects millions worldwide and is often caused by osteoclast induced bone loss. Here, we identify the cytoplasmic protein ELMO1 as an important ‘signaling node’ in osteoclasts. We note that ELMO1 SNPs associate with bone abnormalities in humans, and that ELMO1 deletion in mice reduces bone loss in four in vivo models: osteoprotegerin deficiency, ovariectomy, and two types of inflammatory arthritis. Our transcriptomic analyses coupled with CRISPR/Cas9 genetic deletion identify Elmo1 associated regulators of osteoclast function, including cathepsin G and myeloperoxidase. Further, we define the ‘ELMO1 interactome’ in osteoclasts via proteomics and reveal proteins required for bone degradation. ELMO1 also contributes to osteoclast sealing zone on bone-like surfaces and distribution of osteoclast-specific proteases. Finally, a 3D structure-based ELMO1 inhibitory peptide reduces bone resorption in wild type osteoclasts. Collectively, we identify ELMO1 as a signaling hub that regulates osteoclast function and bone loss, with relevance to osteoporosis and arthritis.

Osteoporosis and bone fractures affect millions of patients worldwide and are often due to increased bone resorption. Here the authors identify the cytoplasmic protein ELMO1 as an important ‘signaling node’ promoting the bone resorption function of osteoclasts.

Dynamic changes in RNA-protein interactions and RNA secondary structure in mammalian erythropoiesis

Two features of eukaryotic RNA molecules that regulate their post-transcriptional fates are RNA secondary structure and RNA-binding protein (RBP) interaction sites. However, a comprehensive global overview of the dynamic nature of these sequence features during erythropoiesis has never been obtained. Here, we use our ribonuclease-mediated structure and RBP-binding site mapping approach to reveal the global landscape of RNA secondary structure and RBP-RNA interaction sites and the dynamics of these features during this important developmental process. We identify dynamic patterns of RNA secondary structure and RBP binding throughout the process and determine a set of corresponding protein-bound sequence motifs along with their dynamic structural and RBP-binding contexts. Finally, using these dynamically bound sequences, we identify a number of RBPs that have known and putative key functions in post-transcriptional regulation during mammalian erythropoiesis. In total, this global analysis reveals new post-transcriptional regulators of mammalian blood cell development.

Protein Interaction Network-based Deep Learning Framework for Identifying Disease-Associated Human Proteins

Infectious diseases in humans appear to be one of the most primary public health issues. Identification of novel disease-associated proteins will furnish an efficient recognition of the novel therapeutic targets. Here, we develop a Graph Convolutional Network (GCN)-based model called PINDeL to identify the disease-associated host proteins by integrating the human Protein Locality Graph and its corresponding topological features. Because of the amalgamation of GCN with the protein interaction network, PINDeL achieves the highest accuracy of 83.45% while AUROC and AUPRC values are 0.90 and 0.88, respectively. With high accuracy, recall, F1-score, specificity, AUROC, and AUPRC, PINDeL outperforms other existing machine-learning and deep-learning techniques for disease gene/protein identification in humans. Application of PINDeL on an independent dataset of 24320 proteins, which are not used for training, validation, or testing purposes, predicts 6448 new disease-protein associations of which we verify 3196 disease-proteins through experimental evidence like disease ontology, Gene Ontology, and KEGG pathway enrichment analyses. Our investigation informs that experimentally-verified 748 proteins are indeed responsible for pathogen-host protein interactions of which 22 disease-proteins share their association with multiple diseases such as cancer, aging, chem-dependency, pharmacogenomics, normal variation, infection, and immune-related diseases. This unique Graph Convolution Network-based prediction model is of utmost use in large-scale disease-protein association prediction and hence, will provide crucial insights on disease pathogenesis and will further aid in developing novel therapeutics.

Evolutionary conservation and divergence of the human brain transcriptome

Background

Mouse models have allowed for the direct interrogation of genetic effects on molecular, physiological, and behavioral brain phenotypes. However, it is unknown to what extent neurological or psychiatric traits may be human- or primate-specific and therefore which components can be faithfully recapitulated in mouse models.

Results

We compare conservation of co-expression in 116 independent data sets derived from human, mouse, and non-human primate representing more than 15,000 total samples. We observe greater changes occurring on the human lineage than mouse, and substantial regional variation that highlights cerebral cortex as the most diverged region. Glia, notably microglia, astrocytes, and oligodendrocytes are the most divergent cell type, three times more on average than neurons. We show that cis-regulatory sequence divergence explains a significant fraction of co-expression divergence. Moreover, protein coding sequence constraint parallels co-expression conservation, such that genes with loss of function intolerance are enriched in neuronal, rather than glial modules. We identify dozens of human neuropsychiatric and neurodegenerative disease risk genes, such as COMT, PSEN-1, LRRK2, SHANK3, and SNCA, with highly divergent co-expression between mouse and human and show that 3D human brain organoids recapitulate in vivo co-expression modules representing several human cell types.

Conclusions

We identify robust co-expression modules reflecting whole-brain and regional patterns of gene expression. Compared with those that represent basic metabolic processes, cell-type-specific modules, most prominently glial modules, are the most divergent between species. These data and analyses serve as a foundational resource to guide human disease modeling and its interpretation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-020-02257-z.

Mutation in PHACTR1 associated with multifocal epilepsy with infantile spasms and hypsarrhythmia

A young boy with multifocal epilepsy with infantile spasms and hypsarrhythmia with minimal organic lesions of brain structures underwent DNA diagnosis using whole‐exome sequencing. A heterozygous amino‐acid substitution p.L519R in a PHACTR1 gene was identified. PHACTR1 belongs to a protein family of G‐actin binding protein phosphatase 1 (PP1) cofactors and was not previously associated with a human disease. The missense single nucleotide variant in the proband was shown to occur de novo in the paternal allele. The mutation was shown in vitro to reduce the affinity of PHACTR1 for G‐actin, and to increase its propensity to form complexes with the catalytic subunit of PP1. These properties are associated with altered subcellular localization of PHACTR1 and increased ability to induce cytoskeletal rearrangements. Although the molecular role of the PHACTR1 in neuronal excitability and differentiation remains to be defined, PHACTR1 has been previously shown to be involved in Slack channelopathy pathogenesis, consistent with our findings. We conclude that this activating mutation in PHACTR1 causes a severe type of sporadic multifocal epilepsy in the patient.

Comparative analysis on the anti-inflammatory/immune effect of mesenchymal stem cell therapy for the treatment of pulmonary arterial hypertension

Despite the advancement of targeted therapy for pulmonary arterial hypertension (PAH), poor prognosis remains a reality. Mesenchymal stem cells (MSCs) are one of the most clinically feasible alternative treatment options. We compared the treatment effects of adipose tissue (AD)-, bone marrow (BD)-, and umbilical cord blood (UCB)-derived MSCs in the rat monocrotaline-induced pulmonary hypertension (PH) model. The greatest improvement in the right ventricular function was observed in the UCB-MSCs treated group. The UCB-MSCs treated group also exhibited the greatest improvement in terms of the largest decrease in the medial wall thickness, perivascular fibrosis, and vascular cell proliferation, as well as the lowest levels of recruitment of innate and adaptive immune cells and associated inflammatory cytokines. Gene expression profiling of lung tissue confirmed that the UCB-MSCs treated group had the most notably attenuated immune and inflammatory profiles. Network analysis further revealed that the UCB-MSCs group had the greatest therapeutic effect in terms of the normalization of all three classical PAH pathways. The intravenous injection of the UCB-MSCs, compared with those of other MSCs, showed superior therapeutic effects in the PH model for the (1) right ventricular function, (2) vascular remodeling, (3) immune/inflammatory profiles, and (4) classical PAH pathways.

Bi-allelic BRWD1 variants cause male infertility with asthenoteratozoospermia and likely primary ciliary dyskinesia

Genetics-associated asthenoteratozoospermia is often seen in patients with multiple morphological abnormalities of the sperm flagella (MMAF). Although 24 causative genes have been identified, these explain only approximately half of patients with MMAF. Since sperm flagella and motile cilia (especially respiratory cilia) have similar axonemal structures, many patients with MMAF also exhibit respiratory symptoms, such as recurrent airway infection, chronic sinusitis, and bronchiectasis, which are frequently associated with primary ciliary dyskinesia (PCD), another recessive disorder. Here, exome sequencing was conducted to evaluate the genetic cause in 53 patients with MMAF and classic PCD/PCD-like symptoms. Two homozygous missense variants and a compound-heterozygous variant in the BRWD1 gene were identified in three unrelated individuals. BRWD1 staining was detected in the whole flagella and respiratory cilia of normal controls but was absent in BRWD1-mutated individuals. Transmission electron microscopy and immunostaining demonstrated that BRWD1 deficiency in human affected respiratory cilia and sperm flagella differently, as the absence of outer and inner dynein arms in sperm flagellum and respiratory cilia, while with a decreased number and outer doublet microtubule defects of respiratory cilia. To our knowledge, this is the first report of a BRWD1-variant-related disease in humans, manifesting as an autosomal recessive form of MMAF and PCD/PCD-like symptoms. Our data provide a basis for further exploring the molecular mechanism of BRWD1 gene during spermatogenesis and ciliogenesis.

Subacute cannabidiol alters genome-wide DNA methylation in adult mouse hippocampus

Use of cannabidiol (CBD), the most abundant non-psychoactive compound found in cannabis (Cannabis sativa), has recently increased as a result of widespread availability of CBD-containing products. CBD is FDA-approved for the treatment of epilepsy and exhibits anxiolytic, antipsychotic, prosocial, and other behavioral effects in animal studies and clinical trials, however, the underlying mechanisms governing these phenotypes are still being elucidated. The epigenome, particularly DNA methylation, is responsive to environmental input and can govern persistent patterns of gene regulation affecting phenotype across the life course. In order to understand the epigenomic activity of cannabidiol exposure in the adult brain, 12-week-old male wild-type a/a Agouti viable yellow (A^vy ) mice were exposed to either 20 mg/kg CBD or vehicle daily by oral administration for 14 days. Hippocampal tissue was collected and reduced-representation bisulfite sequencing (RRBS) was performed. Analyses revealed 3,323 differentially methylated loci (DMLs) in CBD-exposed animals with a small skew toward global hypomethylation. Genes for cell adhesion and migration, dendritic spine development, and excitatory postsynaptic potential were found to be enriched in a gene ontology term analysis of DML-containing genes, and disease ontology enrichment revealed an overrepresentation of DMLs in gene sets associated with autism spectrum disorder, schizophrenia, and other phenotypes. These results suggest that the epigenome may be a key substrate for CBD's behavioral effects and provides a wealth of gene regulatory information for further study.

Anatomical deadspace during resuscitation of infants with congenital diaphragmatic hernia

BACKGROUND

Congenital diaphragmatic hernia (CDH) has a high mortality and morbidity related to pulmonary hypoplasia.

AIMS

To test the hypothesis that CDH infants who survived would have a greater anatomical deadspace reflecting less severe pulmonary hypoplasia. Furthermore, infants with CDH who had undergone feto-tracheal occlusion (FETO) would have a greater anatomical deadspace.

STUDY DESIGN

Infants were studied during resuscitation in the delivery suite. They were all intubated immediately at delivery, given a neuromuscular blocking agent and underwent respiratory monitoring. The anatomical deadspace was calculated from volumetric capnography measurements.

SUBJECTS

Thirty infants born at 32 weeks of gestation or greater and diagnosed antenatally with a CDH were studied. Eleven had undergone FETO and overall five died.

OUTCOME MEASURES

Anatomical deadspace (Vd_ANA) and survival to discharge.

RESULTS

The median (IQR) gestational age of the infants was 38.1 (35.2-39.3) weeks and birthweight 2.8 (2.3-3.3) kg. The anatomical deadspace was higher in those infants who survived (2.9 (2.8-3.3) mls/kg) compared to those who died (2.2 (2.1-2.7) mls/kg; p = 0.003) and was higher in those who had undergone FETO (3.0 (2.8-3.8) mls/kg) compared to those who had not (2.8 (2.4-3.0) mls/kg; p = 0.032). In predicting survival to discharge, the anatomical deadspace had an AUC of 0.90 (p = 0.006).

CONCLUSIONS

CDH infants who survived had a larger anatomical deadspace than those who died suggesting they had less lung hypoplasia. In addition, infants who had undergone FETO had greater anatomical deadspace possibly reflecting distension of the conducting airways.

Transgenic animal models of congenital diaphragmatic hernia: a comprehensive overview of candidate genes and signaling pathways

Congenital diaphragmatic hernia (CDH) is a relatively common and life-threatening birth defect, characterized by incomplete formation of the diaphragm. Because CDH herniation occurs at the same time as preacinar airway branching, normal lung development becomes severely disrupted, resulting almost invariably in pulmonary hypoplasia. Despite various research efforts over the past decades, the pathogenesis of CDH and associated lung hypoplasia remains poorly understood. With the advent of molecular techniques, transgenic animal models of CDH have generated a large number of candidate genes, thus providing a novel basis for future research and treatment. This review article offers a comprehensive overview of genes and signaling pathways implicated in CDH etiology, whilst also discussing strengths and limitations of transgenic animal models in relation to the human condition.

miR155 regulation of behavior, neuropathology, and cortical transcriptomics in Alzheimer's disease

MicroRNAs are recognized as important regulators of many facets of physiological brain function while also being implicated in the pathogenesis of several neurological disorders. Dysregulation of miR155 is widely reported across a variety of neurodegenerative conditions, including Alzheimer's disease (AD), Parkinson's disease, amyotrophic lateral sclerosis, and traumatic brain injury. In previous work, we observed that experimentally validated miR155 gene targets were consistently enriched among genes identified as differentially expressed across multiple brain tissue and disease contexts. In particular, we found that human herpesvirus-6A (HHV-6A) suppressed miR155, recapitulating reports of miR155 inhibition by HHV-6A in infected T-cells, thyrocytes, and natural killer cells. In earlier studies, we also reported the effects of constitutive deletion of miR155 on accelerating the accumulation of Aβ deposits in 4-month-old APP/PSEN1 mice. Herein, we complete the cumulative characterization of transcriptomic, electrophysiological, neuropathological, and learning behavior profiles from 4-, 8- and 10-month-old WT and APP/PSEN1 mice in the absence or presence of miR155. We also integrated human post-mortem brain RNA-sequences from four independent AD consortium studies, together comprising 928 samples collected from six brain regions. We report that gene expression perturbations associated with miR155 deletion in mouse cortex are in aggregate observed to be concordant with AD-associated changes across these independent human late-onset AD (LOAD) data sets, supporting the relevance of our findings to human disease. LOAD has recently been formulated as the clinicopathological manifestation of a multiplex of genetic underpinnings and pathophysiological mechanisms. Our accumulated data are consistent with such a formulation, indicating that miR155 may be uniquely positioned at the intersection of at least four components of this LOAD "multiplex": (1) innate immune response pathways; (2) viral response gene networks; (3) synaptic pathology; and (4) proamyloidogenic pathways involving the amyloid β peptide (Aβ).

Computational resources associating diseases with genotypes, phenotypes and exposures

The causes of a disease and its therapies are not only related to genotypes, but also associated with other factors, including phenotypes, environmental exposures, drugs and chemical molecules. Distinguishing disease-related factors from many neutral factors is critical as well as difficult. Over the past two decades, bioinformaticians have developed many computational resources to integrate the omics data and discover associations among these factors. However, researchers and clinicians are experiencing difficulties in choosing appropriate resources from hundreds of relevant databases and software tools. Here, in order to assist the researchers and clinicians, we systematically review the public computational resources of human diseases related to genotypes, phenotypes, environment factors, drugs and chemical exposures. We briefly describe the development history of these computational resources, followed by the details of the relevant databases and software tools. We finally conclude with a discussion of current challenges and future opportunities as well as prospects on this topic.

A new MEIOB mutation is a recurrent cause for azoospermia and testicular meiotic arrest

STUDY QUESTION

Are there genetic variants that can be used for the clinical evaluation of azoospermic men?

SUMMARY ANSWER

A novel homozygous frame-shift mutation in the MEIOB gene was identified in three azoospermic patients from two different families.

WHAT IS KNOWN ALREADY

Up to 1% of all men have complete absence of sperm in the semen, a condition known as azoospermia. There are very few tools for determining the etiology of azoospermia and the likelihood of sperm cells in the testis. The MEIOB gene codes for a single-strand DNA binding protein required for DNA double-strand breaks repair during meiosis. MEIOB appears to be exclusively expressed in human and mouse testis, and MeioB knockout mice are azoospermic due to meiotic arrest.

STUDY DESIGN, SIZE, DURATION

Two brothers with non-obstructive azoospermia (NOA) underwent whole-exome sequencing followed by comprehensive bioinformatics analyses. Candidate variations were further screened in infertile and fertile men, as well as in public and local reference databases.

PARTICIPANTS/MATERIALS, SETTING, METHODS

This study included 159 infertile and 77 fertile men. The exomes of two Arab men were completely sequenced. In addition, 213 other men of the same Arab ethnicity (136 infertile and 77 fertile men) underwent restriction fragment length polymorphism (RFLP) screening, as did 21 NOA men, of other ethnicities, with testicular impairment of spermatocyte arrest. All of the infertile men underwent Y-chromosome microdeletion and CFTR gene mutation assessments. Comprehensive bioinformatics analyses were designed to uncover candidate mutations associated with azoospermia.

MAIN RESULTS AND THE ROLE OF CHANCE

A novel homozygous frame-shift mutation in the MEIOB gene was identified in two brothers of Arab ethnicity. This frame-shift is predicted to result in a truncated MEIOB protein, which lacks the conserved C-terminal DNA binding domain. RFLP screening of the mutation in 157 infertile men, including 112 NOA patients of Arab ethnicity, identified an additional unrelated NOA patient with the same homozygous mutation and a similar testicular impairment. This mutation was not found in available public databases (n > 160 000), nor in the 77 proven fertile men, nor in our database of local Israeli population variations derived from exome and genome sequencing data (n = 500).

LIMITATIONS, REASONS FOR CAUTION

We have thus far screened for only two specific MEIOB probable pathogenic mutations in a relatively small local cohort. Therefore, the relative incidence of MEIOB mutations in azoospermia should be further assessed in larger and diverse cohorts in order to determine the efficiency of MEIOB sequence screening for clinical evaluations.

WIDER IMPLICATIONS OF THE FINDINGS

The relatively high incidence of likely NOA-causing mutations in MEIOB that was found in our cohort supports the idea that a complete screening of this gene might be beneficial for clinical evaluation of NOA patients.

STUDY FUNDING/COMPETING INTEREST(S)

This research was supported in part by a grant to EA from the European Research Council under the European Union's Seventh Framework Programme (FP/2007-2013)/ERC grant agreement (616088). There are no competing interests.

TRIAL REGISTRATION NUMBER

N/A.

Dynamic rewiring of the human interactome by interferon signaling

BACKGROUND

The type I interferon (IFN) response is an ancient pathway that protects cells against viral pathogens by inducing the transcription of hundreds of IFN-stimulated genes. Comprehensive catalogs of IFN-stimulated genes have been established across species and cell types by transcriptomic and biochemical approaches, but their antiviral mechanisms remain incompletely characterized. Here, we apply a combination of quantitative proteomic approaches to describe the effects of IFN signaling on the human proteome, and apply protein correlation profiling to map IFN-induced rearrangements in the human protein-protein interaction network.

RESULTS

We identify > 26,000 protein interactions in IFN-stimulated and unstimulated cells, many of which involve proteins associated with human disease and are observed exclusively within the IFN-stimulated network. Differential network analysis reveals interaction rewiring across a surprisingly broad spectrum of cellular pathways in the antiviral response. We identify IFN-dependent protein-protein interactions mediating novel regulatory mechanisms at the transcriptional and translational levels, with one such interaction modulating the transcriptional activity of STAT1. Moreover, we reveal IFN-dependent changes in ribosomal composition that act to buffer IFN-stimulated gene protein synthesis.

CONCLUSIONS

Our map of the IFN interactome provides a global view of the complex cellular networks activated during the antiviral response, placing IFN-stimulated genes in a functional context, and serves as a framework to understand how these networks are dysregulated in autoimmune or inflammatory disease.

Mouse Genome Database (MGD) 2019

The Mouse Genome Database (MGD; http://www.informatics.jax.org) is the community model organism genetic and genome resource for the laboratory mouse. MGD is the authoritative source for biological reference data sets related to mouse genes, gene functions, phenotypes, and mouse models of human disease. MGD is the primary outlet for official gene, allele and mouse strain nomenclature based on the guidelines set by the International Committee on Standardized Nomenclature for Mice. In this report we describe significant enhancements to MGD, including two new graphical user interfaces: (i) the Multi Genome Viewer for exploring the genomes of multiple mouse strains and (ii) the Phenotype-Gene Expression matrix which was developed in collaboration with the Gene Expression Database (GXD) and allows researchers to compare gene expression and phenotype annotations for mouse genes. Other recent improvements include enhanced efficiency of our literature curation processes and the incorporation of Transcriptional Start Site (TSS) annotations from RIKEN's FANTOM 5 initiative.

The mutational constraint spectrum quantified from variation in 141,456 humans

Genetic variants that inactivate protein-coding genes are a powerful source of information about the phenotypic consequences of gene disruption: genes that are crucial for the function of an organism will be depleted of such variants in natural populations, whereas non-essential genes will tolerate their accumulation. However, predicted loss-of-function variants are enriched for annotation errors, and tend to be found at extremely low frequencies, so their analysis requires careful variant annotation and very large sample sizes1. Here we describe the aggregation of 125,748 exomes and 15,708 genomes from human sequencing studies into the Genome Aggregation Database (gnomAD). We identify 443,769 high-confidence predicted loss-of-function variants in this cohort after filtering for artefacts caused by sequencing and annotation errors. Using an improved model of human mutation rates, we classify human protein-coding genes along a spectrum that represents tolerance to inactivation, validate this classification using data from model organisms and engineered human cells, and show that it can be used to improve the power of gene discovery for both common and rare diseases.

Biomolecular Databases and Subnetwork Identification Approaches of Interest to Big Data Community: An Expert Review

Next-generation sequencing approaches and genome-wide studies have become essential for characterizing the mechanisms of human diseases. Consequently, many researchers have applied these approaches to discover the genetic/genomic causes of common complex and rare human diseases, generating multiomics big data that span the continuum of genomics, proteomics, metabolomics, and many other system science fields. Therefore, there is a significant and unmet need for biological databases and tools that enable and empower the researchers to analyze, integrate, and make sense of big data. There are currently large number of databases that offer different types of biological information. In particular, the integration of gene expression profiles and protein-protein interaction networks provides a deeper understanding of the complex multilayered molecular architecture of human diseases. Therefore, there has been a growing interest in developing methodologies that integrate and contextualize big data from molecular interaction networks to identify biomarkers of human diseases at a subnetwork resolution as well. In this expert review, we provide a comprehensive summary of most popular biomolecular databases for molecular interactions (e.g., Biological General Repository for Interaction Datasets, Kyoto Encyclopedia of Genes and Genomes and Search Tool for The Retrieval of Interacting Genes/Proteins), gene-disease associations (e.g., Online Mendelian Inheritance in Man, Disease-Gene Network, MalaCards), and population-specific databases (e.g., Human Genetic Variation Database), and describe some examples of their usage and potential applications. We also present the most recent subnetwork identification approaches and discuss their main advantages and limitations. As the field of data science continues to emerge, the present analysis offers a deeper and contextualized understanding of the available databases in molecular biomedicine.

RNA-seq analysis of chondrocyte transcriptome reveals genetic heterogeneity in LG/J and SM/J murine strains

OBJECTIVE

To investigate the transcriptomic differences in chondrocytes obtained from LG/J (large, healer) and SM/J (small, non-healer) murine strains in an attempt to discern the molecular pathways implicated in cartilage regeneration and susceptibility to osteoarthritis (OA).

DESIGN

We performed RNA-sequencing on chondrocytes derived from LG/J (n = 16) and SM/J (n = 16) mice. We validated the expression of candidate genes and compared single nucleotide polymorphisms (SNPs) between the two mouse strains. We also examined gene expression of positional candidates for ear pinna regeneration and long bone length quantitative trait loci (QTLs) that display differences in cartilaginous expression.

RESULTS

We observed a distinct genetic heterogeneity between cells derived from LG/J and SM/J mouse strains. We found that gene ontologies representing cell development, cartilage condensation, and regulation of cell differentiation were enriched in LG/J chondrocytes. In contrast, gene ontologies enriched in the SM/J chondrocytes were mainly related to inflammation and degeneration. Moreover, SNP analysis revealed that multiple validated genes vary in sequence between LG/J and SM/J in coding and highly conserved noncoding regions. Finally, we showed that most QTLs have 20-30% of their positional candidates displaying differential expression between the two mouse strains.

CONCLUSIONS

While the enrichment of pathways related to cell differentiation, cartilage development and cartilage condensation infers superior healing potential of LG/J strain, the enrichment of pathways related to cytokine production, immune cell activation and inflammation entails greater susceptibility of SM/J strain to OA. These data provide novel insights into chondrocyte transcriptome and aid in identification of the quantitative trait genes and molecular differences underlying the phenotypic differences associated with individual QTLs.

Novel DNAAF6 variants identified by whole-exome sequencing cause male infertility and primary ciliary dyskinesia

PURPOSE

To identify the genetic cause of patients with primary ciliary dyskinesia (PCD) and male infertility from two unrelated Han Chinese families.

METHODS

We conducted whole-exome sequencing of three individuals with PCD and male infertility from two unrelated Chinese families, and performed a targeted look-up for DNAAF6 variants in our previously reported cohort of 442 individuals (219 with isolated oligoasthenospermia and 223 fertile controls). Ultrastructural and immunostaining analyses of patients' spermatozoa were performed. The pathogenicity of the variants was validated using patient's spermatozoa and HEK293T cells. Intracytoplasmic sperm injection (ICSI) treatment was conducted in two patients.

RESULTS

We identified one novel hemizygous frameshift variant (NM_173494, c.319_329del: p.R107fs) of DNAAF6 gene (previously named PIH1D3) in family 1 and one novel hemizygous missense variant (c.290G>T: p.G97V) in family 2. No hemizygous deleterious variants in DNAAF6 were detected in the control cohort of 442 individuals. Ultrastructural and immunostaining analyses of patients' spermatozoa showed the absence of outer and inner dynein arms in sperm flagella. Both variants were proven to lead to DNAAF6 protein degradation in HEK293T cells. Both patients carrying DNAAF6 variants underwent one ICSI cycle and delivered one healthy child each.

CONCLUSION

We identified novel DNAAF6 variants causing male infertility and PCD in Han Chinese patients. This finding extended the spectrum of variants in DNAAF6 and revealed new light on the impact of DNAAF6 variants in sperm flagella.

An atlas of evidence-based phenotypic associations across the mouse phenome

To date, reliable relationships between mammalian phenotypes, based on diagnostic test measurements, have not been reported on a large scale. The purpose of this study was to present a large mouse phenotype-phenotype relationships dataset as a reference resource, alongside detailed evaluation of the resource. We used bias-minimized comprehensive mouse phenotype data and applied association rule mining to a dataset consisting of only binary (normal and abnormal phenotypes) data to determine relationships among phenotypes. We present 3,686 evidence-based significant associations, comprising 345 phenotypes covering 60 biological systems (functions), and evaluate their characteristics in detail. To evaluate the relationships, we defined a set of phenotype-phenotype association pairs (PPAPs) as a module of phenotypic expression for each of the 345 phenotypes. By analyzing each PPAP, we identified phenotype sub-networks consisting of the largest numbers of phenotypes and distinct biological systems. Furthermore, using hierarchical clustering based on phenotype similarities among the 345 PPAPs, we identified seven community types within a putative phenome-wide association network. Moreover, to promote leverage of these data, we developed and published web-application tools. These mouse phenome-wide phenotype-phenotype association data reveal general principles of relationships among mammalian phenotypes and provide a reference resource for biomedical analyses.

Pan-mammalian analysis of molecular constraints underlying extended lifespan

Although lifespan in mammals varies over 100-fold, the precise evolutionary mechanisms underlying variation in longevity remain unknown. Species-specific genetic changes have been observed in long-lived species including the naked mole-rat, bats, and the bowhead whale, but these adaptations do not generalize to other mammals. We present a novel method to identify associations between rates of protein evolution and continuous phenotypes across the entire mammalian phylogeny. Unlike previous analyses that focused on individual species, we treat absolute and relative longevity as quantitative traits and demonstrate that these lifespan traits affect the evolutionary constraint on hundreds of genes. Specifically, we find that genes related to cell cycle, DNA repair, cell death, the IGF1 pathway, and immunity are under increased evolutionary constraint in large and long-lived mammals. For mammals exceptionally long-lived for their body size, we find increased constraint in inflammation, DNA repair, and NFKB-related pathways. Strikingly, these pathways have considerable overlap with those that have been previously reported to have potentially adaptive changes in single-species studies, and thus would be expected to show decreased constraint in our analysis. This unexpected finding of increased constraint in many longevity-associated pathways underscores the power of our quantitative approach to detect patterns that generalize across the mammalian phylogeny.

Variation in mouse pelvic morphology maps to locations enriched in Sox9 Class II and Pitx1 regulatory features

Variation in pelvic morphology has a complex genetic basis and its patterning and specification is governed by conserved developmental pathways. Whether the mechanisms underlying the differentiation and specification of the pelvis also produce the morphological covariation on which natural selection may act, is still an open question in evolutionary developmental biology. We use high-resolution quantitative trait locus (QTL) mapping in the F₃₄ generation of an advanced intercross experiment (LG,SM-G₃₄ ) to characterize the genetic architecture of the mouse pelvis. We test the prediction that genomic features linked to developmental patterning and differentiation of the hind limb and pelvis and the regulation of chondrogenesis are overrepresented in QTL. We find 31 single QTL trait associations at the genome- or chromosome-wise significance level coalescing to 27 pleiotropic loci. We recover further QTL at a more relaxed significance threshold replicating locations found in a previous experiment in an earlier generation of the same population. QTL were more likely than chance to harbor Pitx1 and Sox9 Class II chromatin immunoprecipitation-seq features active during development of skeletal features. There was weak or no support for the enrichment of seven more categories of developmental features drawn from the literature. Our results suggest that genotypic variation is channeled through a subset of developmental processes involved in the generation of phenotypic variation in the pelvis. This finding indicates that the evolvability of complex traits may be subject to biases not evident from patterns of covariance among morphological features or developmental patterning when either is considered in isolation.

Matching whole genomes to rare genetic disorders: Identification of potential causative variants using phenotype-weighted knowledge in the CAGI SickKids5 clinical genomes challenge

Precise identification of causative variants from whole-genome sequencing data, including both coding and noncoding variants, is challenging. The Critical Assessment of Genome Interpretation 5 SickKids clinical genome challenge provided an opportunity to assess our ability to extract such information. Participants in the challenge were required to match each of the 24 whole-genome sequences to the correct phenotypic profile and to identify the disease class of each genome. These are all rare disease cases that have resisted genetic diagnosis in a state-of-the-art pipeline. The patients have a range of eye, neurological, and connective-tissue disorders. We used a gene-centric approach to address this problem, assigning each gene a multiphenotype-matching score. Mutations in the top-scoring genes for each phenotype profile were ranked on a 6-point scale of pathogenicity probability, resulting in an approximately equal number of top-ranked coding and noncoding candidate variants overall. We were able to assign the correct disease class for 12 cases and the correct genome to a clinical profile for five cases. The challenge assessor found genes in three of these five cases as likely appropriate. In the postsubmission phase, after careful screening of the genes in the correct genome, we identified additional potential diagnostic variants, a high proportion of which are noncoding.

Reproducible Colonization of Germ-Free Mice With the Oligo-Mouse-Microbiota in Different Animal Facilities

The Oligo-Mouse-Microbiota (OMM¹²) is a recently developed synthetic bacterial community for functional microbiome research in mouse models (Brugiroux et al., 2016). To date, the OMM¹² model has been established in several germ-free mouse facilities world-wide and is employed to address a growing variety of research questions related to infection biology, mucosal immunology, microbial ecology and host-microbiome metabolic cross-talk. The OMM¹² consists of 12 sequenced and publically available strains isolated from mice, representing five bacterial phyla that are naturally abundant in the murine gastrointestinal tract (Lagkouvardos et al., 2016). Under germ-free conditions, the OMM¹² colonizes mice stably over multiple generations. Here, we investigated whether stably colonized OMM¹² mouse lines could be reproducibly established in different animal facilities. Germ-free C57Bl/6J mice were inoculated with a frozen mixture of the OMM¹² strains. Within 2 weeks after application, the OMM¹² community reached the same stable composition in all facilities, as determined by fecal microbiome analysis. We show that a second application of the OMM¹² strains after 72 h leads to a more stable community composition than a single application. The availability of such protocols for reliable de novo generation of gnotobiotic rodents will certainly contribute to increasing experimental reproducibility in biomedical research.

LeukmiR: a database for miRNAs and their targets in acute lymphoblastic leukemia

Acute lymphoblastic leukemia (ALL) is one of the most common hematological malignancies in children. Recent studies suggest the involvement of multiple microRNAs in the tumorigenesis of various leukemias. However, until now, no comprehensive database exists for miRNAs and their cognate target genes involved specifically in ALL. Therefore, we developed 'LeukmiR' a dynamic database comprising in silico predicted microRNAs, and experimentally validated miRNAs along with the target genes they regulate in mouse and human. LeukmiR is a user-friendly platform with search strings for ALL-associated microRNAs, their sequences, description of target genes, their location on the chromosomes and the corresponding deregulated signaling pathways. For the user query, different search modules exist where either quick search can be carried out using any fuzzy term or by providing exact terms in specific modules. All entries for both human and mouse genomes can be retrieved through multiple options such as miRNA ID, their accession number, sequence, target genes, Ensemble-ID or Entrez-ID. User can also access miRNA: mRNA interaction networks in different signaling pathways, the genomic location of the targeted regions such as 3'UTR, 5'UTR and exons with their gene ontology and disease ontology information in both human and mouse systems. Herein, we also report 51 novel microRNAs which are not described earlier for ALL. Thus, LeukmiR database will be a valuable source of information for researchers to understand and investigate miRNAs and their targets with diagnostic and therapeutic potential in ALL. Database URL: http://tdb.ccmb.res.in/LeukmiR/.