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

GLP-1 receptor agonists show no detrimental effect on sperm quality in mouse models and cell lines

PURPOSE

Glucagon-like peptide-1 receptor (GLP-1R) agonists exert multiple beneficial effects. However, their effects on reproduction system are controversial. Here, we aimed to investigate their effects on male reproduction and provide safety evidence for future clinical use.

METHODS

Male diabetic mice and aged mice were treated with liraglutide or vehicle, and sperm concentration and motility were assessed. The expression and location of GLP-1R in testicular tissues and in four testicular cell lines (spermatogonia, spermatocytes, Leydig cells, and Sertoli cells) were detected. Cauda epididymis and testicular cells were treated with liraglutide, semaglutide or vehicle, and sperm motility and cell proliferation were detected to determine the direct effect of GLP-1R agonists. Global Glp1r knockout mice were constructed, and testicular morphology, sperm concentration and motility were detected to confirm the effects of GLP-1R signaling on male reproduction.

RESULTS

Liraglutide significantly reduced blood glucose levels, but did not improve sperm parameters in diabetic mice. No significant differences were observed between liraglutide and control group in aged mice. GLP-1R was expressed in testicular tissues and all four cell lines, with the highest expression in Leydig cells. Liraglutide or semaglutide had no impacts on sperm count and motility in vitro, and had no effects on cell proliferation in four cell lines. The Glp1r knockout mice exhibited higher blood glucose levels and preserved normal testicular morphology, but their sperm concentration was higher than that in wildtype mice.

CONCLUSION

GLP-1R agonists have no detrimental effect on sperm concentration and motility in vivo and in vitro, while GLP-1R absence increase the sperm concentration.

Novel mutations in ZMYND15 are associated with male infertility with oligozoospermia/azoospermia

PURPOSE

This study aimed to identify the genetic causes of male infertility associated with oligozoospermia/azoospermia in two unrelated Chinese families.

METHODS

Whole-exome sequencing (WES) and Sanger sequencing were performed on peripheral blood samples from three infertile individuals with reduced sperm counts. Semen analysis data were collected, and sperm morphology was evaluated using hematoxylin and eosin staining, along with transmission electron microscopy. Acidic aniline staining and fluorescence in situ hybridization (FISH) were employed to assess sperm nuclear maturity and chromosome aneuploidy. In vitro analyses were performed to determine the effect of the identified variants.

RESULTS

We identified two novel homozygous variants in ZMYND15: a frameshift variant (NM_001136046.3:c.828-2_833dupAGAGAGCT) in family 1 and a missense variant (c.2051 T > A:p.Met684Lys) in family 2. Both variants were absent in public databases, and the missense variant was predicted to be deleterious. In vitro analyses confirmed that the frameshift variants likely impact protein function. Abnormal sperm head morphologies, characterized by reduced chromatin condensation and nuclear aneuploidy, were frequently observed in ZMYND15 mutant individual.

CONCLUSION

Our findings reveal two novel ZMYND15 variants in three infertile patients with oligozoospermia/azoospermia, expanding the mutational spectrum of ZMYND15 and providing valuable insights for genetic counseling and the diagnosis in cases of male infertility.

Investigating Motor Coordination Using BXD Recombinant Inbred Mice to Model the Genetic Underpinnings of Developmental Coordination Disorder

The fundamental skills for motor coordination and motor control emerge through development. Neurodevelopmental disorders such as developmental coordination disorder (DCD) lead to impaired acquisition of motor skills. This study investigated motor behaviors that reflect the core symptoms of human DCD through the use of BXD recombinant inbred strains of mice that are known to have divergent phenotypes in many behavioral traits, including motor activity. We sought to correlate behavior in basic motor control tasks with the known genotypes of these reference populations of mice using quantitative trait locus (QTL) mapping. We used 12 BXD strains with an average of 16 mice per group to assess the onset of reflexes during the early neonatal stage of life and differences in motor coordination using the tests for open field, rotarod, and gait behaviors during the adolescent/young adulthood period. Results indicated significant variability between strains in when neonatal reflexes appeared and significant strain differences for all measures of motor coordination. Five strains (BXD15, BXD27, BXD28, BXD75, BXD86) struggled with sensorimotor coordination as seen in gait analysis, rotarod, and open field, similar to human presentation of DCD. We identified three significant quantitative trait loci for gait on proximal Chr 3, Chr 4, and distal Chr 6. Based on expression, function, and polymorphism within the mapped QTL intervals, seven candidate genes (Gpr63, Spata5, Trpc3, Cntn6, Chl1, Grm7, Ogg1) emerged. This study offers new insights into mouse motor behavior, which promises to be a first murine model to explore the genetics and neural correlates of DCD.

Protocol for using Multiomics2Targets to identify targets and driver kinases for cancer cohorts profiled with multi-omics assays

The availability of multi-omics data applied to profile cancer cohorts is rapidly increasing. Here, we present a protocol for Multiomics2Targets, a computational pipeline that can identify driver cell signaling pathways, protein kinases, and cell-surface targets for immunotherapy. We describe steps for preparing the data, uploading files, and tuning parameters. We then detail procedures for running the workflow, visualizing the results, and exporting and sharing reports containing the analysis. For complete details on the use and execution of this protocol, please refer to Deng et al.1.

Petagraph: A large-scale unifying knowledge graph framework for integrating biomolecular and biomedical data

Over the past decade, there has been substantial growth in both the quantity and complexity of available biomedical data. In order to more efficiently harness this extensive data and alleviate challenges associated with integration of multi-omics data, we developed Petagraph, a biomedical knowledge graph that encompasses over 32 million nodes and 118 million relationships. Petagraph leverages more than 180 ontologies and standards in the Unified Biomedical Knowledge Graph (UBKG) to embed millions of quantitative genomics data points. Petagraph provides a cohesive data environment that enables users to efficiently analyze, annotate, and discern relationships within and across complex multi-omics datasets supported by UBKG's annotation scaffold. We demonstrate how queries on Petagraph can generate meaningful results across various research contexts and use cases.

AlzGenPred - CatBoost-based gene classifier for predicting Alzheimer's disease using high-throughput sequencing data

AD is a progressive neurodegenerative disorder characterized by memory loss. Due to the advancement in next-generation sequencing, an enormous amount of AD-associated genomics data is available. However, the information about the involvement of these genes in AD association is still a research topic. Therefore, AlzGenPred is developed to identify the AD-associated genes using machine-learning. A total of 13,504 features derived from eight sequence-encoding schemes were generated and evaluated using 16 machine learning algorithms. Network-based features significantly outperformed sequence-based features, effectively distinguishing AD-associated genes. In contrast, sequence-based features failed to classify accurately. To improve performance, we generated 24 fused features (6020 D) from sequence-based encodings, increasing accuracy by 5-7% using a two-step lightGBM-based recursive feature selection method. However, accuracy remained below 70% even after hyperparameter tuning. Therefore, network-based features were used to generate the CatBoost-based ML method AlzGenPred with 96.55% accuracy and 98.99% AUROC. The developed method is tested on the AlzGene dataset where it showed 96.43% accuracy. Then the model was validated using the transcriptomics dataset. AlzGenPred provides a reliable and user-friendly tool for identifying potential AD biomarkers, accelerating biomarker discovery, and advancing our understanding of AD. It is available at https://www.bioinfoindia.org/alzgenpred/ and https://github.com/shuklarohit815/AlzGenPred .

Sequence-Based Machine Learning Reveals 3D Genome Differences between Bonobos and Chimpanzees

The 3D structure of the genome is an important mediator of gene expression. As phenotypic divergence is largely driven by gene regulatory variation, comparing genome 3D contacts across species can further understanding of the molecular basis of species differences. However, while experimental data on genome 3D contacts in humans are increasingly abundant, only a handful of 3D genome contact maps exist for other species. Here, we demonstrate that human experimental data can be used to close this data gap. We apply a machine learning model that predicts 3D genome contacts from DNA sequence to the genomes from 56 bonobos and chimpanzees and identify species-specific patterns of genome folding. We estimated 3D divergence between individuals from the resulting contact maps in 4,420 1 Mb genomic windows, of which ∼17% were substantially divergent in predicted genome contacts. Bonobos and chimpanzees diverged at 89 windows, overlapping genes associated with multiple traits implicated in Pan phenotypic divergence. We discovered 51 bonobo-specific variants that individually produce the observed bonobo contact pattern in bonobo-chimpanzee divergent windows. Our results demonstrate that machine learning methods can leverage human data to fill in data gaps across species, offering the first look at population-level 3D genome variation in nonhuman primates. We also identify loci where changes in 3D folding may contribute to phenotypic differences in our closest living relatives.

Interactions between microbiota and uterine corpus endometrial cancer: A bioinformatic investigation of potential immunotherapy

Microorganisms in the gut and other niches may contribute to carcinogenesis while also altering cancer immune surveillance and therapeutic response. However, determining the impact of genetic variations and interplay with intestinal microbes' environment is difficult and unanswered. Here, we examined the frequency of thirteen mutant genes that caused aberrant gut in thirty different types of cancer using The Cancer Genomic Atlas (TCGA) database. Substantially, our findings show that all these mutated genes are quite frequent in uterine corpus endometrial cancer (UCEC). Further, these mutant genes are implicated in the infiltration of different subset of immune cells within the Tumor Microenvironment (TME) of UCEC patients. The top-ranking mutant genes that promote immune cell invasion into the TME of UCEC patients were PGLYRP2, OLFM4, and TLR5. In this regard, we used the same deconvolution of the TCGA database to analyze the microbiome that have a strong association with immune cells invasion with TME of UCEC patients. Several bacteria and viruses have been linked to the invasion of immune cells, such as B cell memory and T cell regulatory (Tregs), into the TME of UCEC patients. As a result, our findings pave the way for future research into generating novel immunizations against bacteria or viruses as immunotherapy for UCEC patients.

Evolution of ion channels in cetaceans: a natural experiment in the tree of life

Cetaceans represent a natural experiment within the tree of life in which a lineage changed from terrestrial to aquatic habitats. This shift involved phenotypic modifications, representing an opportunity to explore the genetic bases of phenotypic diversity. Among the different molecular systems that maintain cellular homeostasis, ion channels are crucial for the proper physiological functioning of all living species. This study aims to explore the evolution of ion channels during the evolutionary history of cetaceans. To do so, we created a bioinformatic pipeline to annotate the repertoire of ion channels in the genome of the species included in our sampling. Our main results show that cetaceans have, on average, fewer protein-coding genes and a higher percentage of annotated ion channels than non-cetacean mammals. Signals of positive selection were detected in ion channels related to the heart, locomotion, visual and neurological phenotypes. Interestingly, we predict that the NaV1.5 ion channel of most toothed whales (odontocetes) is sensitive to tetrodotoxin, similar to NaV1.7, given the presence of tyrosine instead of cysteine, in a specific position of the ion channel. Finally, the gene turnover rate of the cetacean crown group is more than three times faster than that of non-cetacean mammals.

SingleCellGGM enables gene expression program identification from single-cell transcriptomes and facilitates universal cell label transfer

Gene co-expression analysis of single-cell transcriptomes, aiming to define functional relationships between genes, is challenging due to excessive dropout values. Here, we developed a single-cell graphical Gaussian model (SingleCellGGM) algorithm to conduct single-cell gene co-expression network analysis. When applied to mouse single-cell datasets, SingleCellGGM constructed networks from which gene co-expression modules with highly significant functional enrichment were identified. We considered the modules as gene expression programs (GEPs). These GEPs enable direct cell-type annotation of individual cells without cell clustering, and they are enriched with genes required for the functions of the corresponding cells, sometimes at levels greater than 10-fold. The GEPs are conserved across datasets and enable universal cell-type label transfer across different studies. We also proposed a dimension-reduction method through averaging by GEPs for single-cell analysis, enhancing the interpretability of results. Thus, SingleCellGGM offers a unique GEP-based perspective to analyze single-cell transcriptomes and reveals biological insights shared by different single-cell datasets.

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.

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.

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.

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β).

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.