High-throughput discovery of novel developmental phenotypes

Ciliary and non-ciliary functions of Rab34 during craniofacial bone development

The primary cilium is a hair-like projection that controls cell development and tissue homeostasis. Although accumulated studies identify the molecular link between cilia and cilia-related diseases, the underlying etiology of ciliopathies has not been fully understood. In this paper, we determine the function of Rab34, a small GTPase, as a key regulator for controlling ciliogenesis and type I collagen trafficking in craniofacial development. Mechanistically, Rab34 is required to form cilia that control osteogenic proliferation, survival, and differentiation via cilia-mediated Hedgehog signaling. In addition, Rab34 is indispensable for regulating type I collagen trafficking from the ER to the Golgi. These results demonstrate that Rab34 has both ciliary and non-ciliary functions to regulate osteogenesis. Our study highlights the critical function of Rab34, which may contribute to understanding the novel etiology of ciliopathies that are associated with the dysfunction of RAB34 in humans.

Homeostatic regulation of NAD(H) and NADP(H) in cells

Nicotinamide adenine dinucleotide (NAD+)/reduced NAD+ (NADH) and nicotinamide adenine dinucleotide phosphate (NADP+)/reduced NADP+ (NADPH) are essential metabolites involved in multiple metabolic pathways and cellular processes. NAD+ and NADH redox couple plays a vital role in catabolic redox reactions, while NADPH is crucial for cellular anabolism and antioxidant responses. Maintaining NAD(H) and NADP(H) homeostasis is crucial for normal physiological activity and is tightly regulated through various mechanisms, such as biosynthesis, consumption, recycling, and conversion between NAD(H) and NADP(H). The conversions between NAD(H) and NADP(H) are controlled by NAD kinases (NADKs) and NADP(H) phosphatases [specifically, metazoan SpoT homolog-1 (MESH1) and nocturnin (NOCT)]. NADKs facilitate the synthesis of NADP+ from NAD+, while MESH1 and NOCT convert NADP(H) into NAD(H). In this review, we summarize the physiological roles of NAD(H) and NADP(H) and discuss the regulatory mechanisms governing NAD(H) and NADP(H) homeostasis in three key aspects: the transcriptional and posttranslational regulation of NADKs, the role of MESH1 and NOCT in maintaining NAD(H) and NADP(H) homeostasis, and the influence of the circadian clock on NAD(H) and NADP(H) homeostasis. In conclusion, NADKs, MESH1, and NOCT are integral to various cellular processes, regulating NAD(H) and NADP(H) homeostasis. Dysregulation of these enzymes results in various human diseases, such as cancers and metabolic disorders. Hence, strategies aiming to restore NAD(H) and NADP(H) homeostasis hold promise as novel therapeutic approaches for these diseases.

ZSCAN10 deficiency causes a neurodevelopmental disorder with characteristic oto-facial malformations

Neurodevelopmental disorders are major indications for genetic referral and have been linked to more than 1500 loci including genes encoding transcriptional regulators. The dysfunction of transcription factors often results in characteristic syndromic presentations; however, at least half of these patients lack a genetic diagnosis. The implementation of machine learning approaches has the potential to aid in the identification of new disease genes and delineate associated phenotypes. Next generation sequencing was performed in seven affected individuals with neurodevelopmental delay and dysmorphic features. Clinical characterization included reanalysis of available neuroimaging datasets and 2D portrait image analysis with GestaltMatcher. The functional consequences of ZSCAN10 loss were modelled in mouse embryonic stem cells (mESCs), including a knockout and a representative ZSCAN10 protein truncating variant. These models were characterized by gene expression and western blot analyses, chromatin immunoprecipitation and quantitative PCR (ChIP-qPCR) and immunofluorescence staining. Zscan10 knockout mouse embryos were generated and phenotyped. We prioritized bi-allelic ZSCAN10 loss-of-function variants in seven affected individuals from five unrelated families as the underlying molecular cause. RNA-sequencing analyses in Zscan10-/- mESCs indicated dysregulation of genes related to stem cell pluripotency. In addition, we established in mESCs the loss-of-function mechanism for a representative human ZSCAN10 protein truncating variant by showing alteration of its expression levels and subcellular localization, interfering with its binding to DNA enhancer targets. Deep phenotyping revealed global developmental delay, facial asymmetry and malformations of the outer ear as consistent clinical features. Cerebral MRI showed dysplasia of the semicircular canals as an anatomical correlate of sensorineural hearing loss. Facial asymmetry was confirmed as a clinical feature by GestaltMatcher and was recapitulated in the Zscan10 mouse model along with inner and outer ear malformations. Our findings provide evidence of a novel syndromic neurodevelopmental disorder caused by bi-allelic loss-of-function variants in ZSCAN10.

Congenital syndromic Chiari-like malformation (CSCM) in Holstein cattle: towards unravelling of possible genetic causes

BACKGROUND

Chiari malformation type II (CMII) was originally reported in humans as a rare disorder characterized by the downward herniation of the hindbrain and towering cerebellum. The congenital brain malformation is usually accompanied by spina bifida, a congenital spinal anomaly resulting from incomplete closure of the dorsal aspect of the spinal neural tube, and occasionally by other lesions. A similar disorder has been reported in several animal species, including cattle, particularly as a congenital syndrome. A cause of congenital syndromic Chiari-like malformation (CSCM) in cattle has not been reported to date. We collected a series of 14 CSCM-affected Holstein calves (13 purebred, one Red Danish Dairy F1 cross) and performed whole-genome sequencing (WGS). WGS was performed on 33 cattle, including eight cases with parents (trio-based; group 1), three cases with one parent (group 2), and three single cases (solo-based; group 3).

RESULTS

Sequencing-based genome-wide association study of the 13 Holstein calves with CSCM and 166 controls revealed no significantly associated genome region. Assuming a single Holstein breed-specific recessive allele, no region of shared homozygosity was detected suggesting heterogeneity. Subsequent filtering for protein-changing variants that were only homozygous in the genomes of the individual cases allowed the identification of two missense variants affecting different genes, SHC4 in case 4 in group 1 and WDR45B in case 13 in group 3. Furthermore, these two variants were only observed in Holstein cattle when querying WGS data of > 5,100 animals. Alternatively, potential de novo mutational events were assessed in each case. Filtering for heterozygous private protein-changing variants identified one DYNC1H1 frameshift variant as a candidate causal dominant acting allele in case 12 in group 3. Finally, the presence of larger structural DNA variants and chromosomal abnormalities was investigated in all cases. Depth of coverage analysis revealed two different partial monosomies of chromosome 2 segments in cases 1 and 7 in group 1 and a trisomy of chromosome 12 in the WDR45B homozygous case 13 in group 3.

CONCLUSIONS

This study presents for the first time a detailed genomic evaluation of CSCM in Holstein cattle and suggests an unexpected genetic and allelic heterogeneity considering the mode of inheritance, as well as the type of variant. For the first time, we propose candidate causal variants that may explain bovine CSCM in a certain proportion of affected calves. We present cattle as a large animal model for human CMII and propose new genes and genomic variants as possible causes for related diseases in both animals and humans.

Morphological simulation tests the limits on phenotype discovery in 3D image analysis

In the past few decades, advances in 3D imaging have created new opportunities for reverse genetic screens. Rapidly growing datasets of 3D images of genetic knockouts require high-throughput, automated computational approaches for identifying and characterizing new phenotypes. However, exploratory, discovery-oriented image analysis pipelines used to discover these phenotypes can be difficult to validate because, by their nature, the expected outcome is not known a priori . Introducing known morphological variation through simulation can help distinguish between real phenotypic differences and random variation; elucidate the effects of sample size; and test the sensitivity and reproducibility of morphometric analyses. Here we present a novel approach for 3D morphological simulation that uses open-source, open-access tools available in 3D Slicer, SlicerMorph, and Advanced Normalization Tools in R (ANTsR). While we focus on diffusible-iodine contrast-enhanced micro-CT (diceCT) images, this approach can be used on any volumetric image. We then use our simulated datasets to test whether tensor-based morphometry (TBM) can recover our introduced differences; to test how effect size and sample size affect detectability; and to determine the reproducibility of our results. In our approach to morphological simulation, we first generate a simulated deformation based on a reference image and then propagate this deformation to subjects using inverse transforms obtained from the registration of subjects to the reference. This produces a new dataset with a shifted population mean while retaining individual variability because each sample deforms more or less based on how different or similar it is from the reference. TBM is a widely-used technique that statistically compares local volume differences associated with local deformations. Our results showed that TBM recovered our introduced morphological differences, but that detectability was dependent on the effect size, the sample size, and the region of interest (ROI) included in the analysis. Detectability of subtle phenotypes can be improved both by increasing the sample size and by limiting analyses to specific body regions. However, it is not always feasible to increase sample sizes in screens of essential genes. Therefore, methodical use of ROIs is a promising way to increase the power of TBM to detect subtle phenotypes. Generating known morphological variation through simulation has broad applicability in developmental, evolutionary, and biomedical morphometrics and is a useful way to distinguish between a failure to detect morphological difference and a true lack of morphological difference. Morphological simulation can also be applied to AI-based supervised learning to augment datasets and overcome dataset limitations.

Lethal phenotypes in Mendelian disorders

PURPOSE

Existing resources that characterize the essentiality status of genes are based on either proliferation assessment in human cell lines, viability evaluation in mouse knockouts, or constraint metrics derived from human population sequencing studies. Several repositories document phenotypic annotations for rare disorders; however, there is a lack of comprehensive reporting on lethal phenotypes.

METHODS

We queried Online Mendelian Inheritance in Man for terms related to lethality and classified all Mendelian genes according to the earliest age of death recorded for the associated disorders, from prenatal death to no reports of premature death. We characterized the genes across these lethality categories, examined the evidence on viability from mouse models and explored how this information could be used for novel gene discovery.

RESULTS

We developed the Lethal Phenotypes Portal to showcase this curated catalog of human essential genes. Differences in the mode of inheritance, physiological systems affected, and disease class were found for genes in different lethality categories, as well as discrepancies between the lethal phenotypes observed in mouse and human.

CONCLUSION

We anticipate that this resource will aid clinicians in the diagnosis of early lethal conditions and assist researchers in investigating the properties that make these genes essential for human development.

TATA-binding associated factors have distinct roles during early mammalian development

Early embryonic development is a finely orchestrated process that requires precise regulation of gene expression coordinated with morphogenetic events. TATA-box binding protein-associated factors (TAFs), integral components of transcription initiation coactivators like TFIID and SAGA, play a crucial role in this intricate process. Here we show that disruptions in TAF5, TAF12 and TAF13 individually lead to embryonic lethality in the mouse, resulting in overlapping yet distinct phenotypes. Taf5 and Taf12 mutant embryos exhibited a failure to implant post-blastocyst formation, and Taf5 mutants have aberrant lineage specification within the inner cell mass. In contrast, Taf13 mutant embryos successfully implant and form egg-cylinder stages but fail to initiate gastrulation. Strikingly, we observed a depletion of pluripotency factors in TAF13-deficient embryos, including OCT4, NANOG and SOX2, highlighting an indispensable role of TAF13 in maintaining pluripotency. Transcriptomic analysis revealed distinct gene targets affected by the loss of TAF5, TAF12 and TAF13. Thus, we propose that TAF5, TAF12 and TAF13 convey locus specificity to the TFIID complex throughout the mouse genome.

The 36th International Mammalian Genome Conference: A scientific gathering under the cherry blossoms in Tsukuba

The 36th International Mammalian Genome Conference (IMGC) was held in a hybrid format at the Tsukuba International Congress Center in Tsukuba, Ibaraki, Japan, for 4 days from March 28 to 31, 2023. This international conference on functional genomics of mouse, human, and other mammalian species attracted 246 participants in total, of which 129 were from outside Japan, including Europe, the United States and Asia, and 117 participants were from Japan. The conference included three technical workshops, keynote lectures by domestic researchers, commemorative lectures for the conference awards, 57 oral presentations, and 97 poster presentations. The event was a great success. Topics included the establishment and analysis of disease models using genetically engineered or spontaneous mutant mice, systems genetic analysis using mouse strains such as wild-derived mice and recombinant inbred mouse strains, infectious diseases, immunology, and epigenetics. In addition, as a joint program, a two-day RIKEN Symposium was held, and active discussions continued over the four-day period. Also, there was a trainee symposium, in which young researchers were encouraged to participate, and excellent papers were selected as oral presentations in the main session.

Our current understanding of the biological impact of endometrial cancer mtDNA genome mutations and their potential use as a biomarker

Endometrial cancer (EC) is a devastating and common disease affecting women's health. The NCI Surveillance, Epidemiology, and End Results Program predicted that there would be >66,000 new cases in the United States and >13,000 deaths from EC in 2023, and EC is the sixth most common cancer among women worldwide. Regulation of mitochondrial metabolism plays a role in tumorigenesis. In proliferating cancer cells, mitochondria provide the necessary building blocks for biosynthesis of amino acids, lipids, nucleotides, and glucose. One mechanism causing altered mitochondrial activity is mitochondrial DNA (mtDNA) mutation. The polyploid human mtDNA genome is a circular double-stranded molecule essential to vertebrate life that harbors genes critical for oxidative phosphorylation plus mitochondrial-derived peptide genes. Cancer cells display aerobic glycolysis, known as the Warburg effect, which arises from the needs of fast-dividing cells and is characterized by increased glucose uptake and conversion of glucose to lactate. Solid tumors often contain at least one mtDNA substitution. Furthermore, it is common for cancer cells to harbor mixtures of wild-type and mutant mtDNA genotypes, known as heteroplasmy. Considering the increase in cancer cell energy demand, the presence of functionally relevant carcinogenesis-inducing or environment-adapting mtDNA mutations in cancer seems plausible. We review 279 EC tumor-specific mtDNA single nucleotide variants from 111 individuals from different studies. Many transition mutations indicative of error-prone DNA polymerase γ replication and C to U deamination events were present. We examine the spectrum of mutations and their heteroplasmy and discuss the potential biological impact of recurrent, non-synonymous, insertion, and deletion mutations. Lastly, we explore current EC treatments, exploiting cancer cell mitochondria for therapy and the prospect of using mtDNA variants as an EC biomarker.

Systems Genetics Analyses Reveals Key Genes Related to Behavioral Traits in the Striatum of CFW Mice

The striatum plays a central role in directing many complex behaviors ranging from motor control to action choice and reward learning. In our study, we used 55 male CFW mice with rapid decay linkage disequilibrium to systematically mine the striatum-related behavioral functional genes by analyzing their striatal transcriptomes and 79 measured behavioral phenotypic data. By constructing a gene coexpression network, we clustered the genes into 13 modules, with most of them being positively correlated with motor traits. Based on functional annotations as well as Fisher's exact and hypergeometric distribution tests, brown and magenta modules were identified as core modules. They were significantly enriched for striatal-related functional genes. Subsequent Mendelian randomization analysis verified the causal relationship between the core modules and dyskinesia. Through the intramodular gene connectivity analysis, Adcy5 and Kcnma1 were identified as brown and magenta module hub genes, respectively. Knock outs of both Adcy5 and Kcnma1 lead to motor dysfunction in mice, and KCNMA1 acts as a risk gene for schizophrenia and smoking addiction in humans. We also evaluated the cellular composition of each module and identified oligodendrocytes in the striatum to have a positive role in motor regulation.

PP1β opposes classic PP1 function, inhibiting spine maturation and promoting LTP

Protein phosphatase 1 (PP1) regulates synaptic plasticity and has been described as a molecular constraint on learning and memory. There are three neuronal isoforms, PP1α, PP1β, and PP1γ, but little is known about their individual functions. PP1α and PP1γ are assumed to mediate the effects of PP1 on learning and memory based on their enrichment at dendritic spines and their preferential binding to neurabin and spinophilin, major PP1 synaptic scaffolding proteins. However, it was recently discovered that human de novo PP1β mutations cause intellectual disability, suggesting an important but ill-defined role for PP1β. In this study, we investigated the functions of each PP1 isoform in hippocampal synaptic physiology using conditional CA1-specific knockout mice. In stark contrast to classic PP1 function, we found that PP1β promotes synaptic plasticity as well as spatial memory. These changes in synaptic plasticity and memory are accompanied by changes in GluA1 phosphorylation, GluN2A levels, and dendritic spine density and morphology, including silent synapse number. These functions of PP1β reveal a previously unidentified signaling pathway regulating spine maturation and plasticity, broadening our understanding of the complex role of PP1 in synaptic physiology.

Epigenetic signature of very low birth weight in young adult life

BACKGROUND

Globally, one in ten babies is born preterm (<37 weeks), and 1-2% preterm at very low birth weight (VLBW, <1500 g). As adults, they are at increased risk for a plethora of health conditions, e.g., cardiometabolic disease, which may partly be mediated by epigenetic regulation. We compared blood DNA methylation between young adults born at VLBW and controls.

METHODS

157 subjects born at VLBW and 161 controls born at term, from the Helsinki Study of Very Low Birth Weight Adults, were assessed for peripheral venous blood DNA methylation levels at mean age of 22 years. Significant CpG-sites (5'-C-phosphate-G-3') were meta-analyzed against continuous birth weight in four independent cohorts (pooled n = 2235) with cohort mean ages varying from 0 to 31 years.

RESULTS

In the discovery cohort, 66 CpG-sites were differentially methylated between VLBW adults and controls. Top hits were located in HIF3A, EBF4, and an intergenic region nearest to GLI2 (distance 57,533 bp). Five CpG-sites, all in proximity to GLI2, were hypermethylated in VLBW and associated with lower birth weight in the meta-analysis.

CONCLUSION

We identified differentially methylated CpG-sites suggesting an epigenetic signature of preterm birth at VLBW present in adult life.

IMPACT

Being born preterm at very low birth weight has major implications for later health and chronic disease risk factors. The mechanism linking preterm birth to later outcomes remains unknown. Our cohort study of 157 very low birth weight adults and 161 controls found 66 differentially methylated sites at mean age of 22 years. Our findings suggest an epigenetic mark of preterm birth present in adulthood, which opens up opportunities for mechanistic studies.

Influence of extracellular ATP on mammalian sperm physiology

In addition to its central role in cellular metabolism, adenosine 5'-triphosphate (ATP) is an important extracellular signalling molecule involved in various physiological processes. In reproduction, extracellular ATP participates in both autocrine and paracrine paths regulating gametogenesis, gamete maturation and fertilisation. This review focusses on how extracellular ATP modulates sperm physiology with emphasis on the mammalian acrosome reaction. The presence of extracellular ATP in the reproductive tract is primarily determined by the ion channels and transporters that influence its movement within the cells comprising the tract. The main targets of extracellular ATP in spermatozoa are its own transporters, particularly species-specific sperm purinergic receptors. We also discuss notable phenotypes from knock-out mouse models and human Mendelian inheritance related to ATP release mechanisms, along with immunological, proteomic, and functional observations regarding sperm purinergic receptors and their involvement in sperm signalling.

Computational identification of disease models through cross-species phenotype comparison

The use of standardised phenotyping screens to identify abnormal phenotypes in mouse knockouts, together with the use of ontologies to describe such phenotypic features, allows the implementation of an automated and unbiased pipeline to identify new models of disease by performing phenotype comparisons across species. Using data from the International Mouse Phenotyping Consortium (IMPC), approximately half of mouse mutants are able to mimic, at least partially, the human ortholog disease phenotypes as computed by the PhenoDigm algorithm. We found the number of phenotypic abnormalities in the mouse and the corresponding Mendelian disorder, the pleiotropy and severity of the disease, and the viability and zygosity status of the mouse knockout to be associated with the ability of mouse models to recapitulate the human disorder. An analysis of the IMPC impact on disease gene discovery through a publication-tracking system revealed that the resource has been implicated in at least 109 validated rare disease-gene associations over the last decade.

Trefoil factor 1 suppresses stemness and enhances chemosensitivity of pancreatic cancer

BACKGROUND AND AIMS

Pancreatic cancer is one of the most lethal malignancies, partly due to resistance to conventional chemotherapy. The chemoresistance of malignant tumors is associated with epithelial-mesenchymal transition (EMT) and the stemness of cancer cells. The aim of this study is to investigate the availability and functional mechanisms of trefoil factor family 1 (TFF1), a tumor-suppressive protein in pancreatic carcinogenesis, to treat pancreatic cancer.

METHODS

To investigate the role of endogenous TFF1 in human and mice, specimens of human pancreatic cancer and genetically engineered mouse model of pancreatic cancer (KPC/TFF1KO; Pdx1-Cre/LSL-KRASG12D/LSL-p53R172H/TFF1-/-) were analyzed by immunohistochemistry (IHC). To explore the efficacy of extracellular administration of TFF1, recombinant and chemically synthesized TFF1 were administered to pancreatic cancer cell lines, a xenograft mouse model and a transgenic mouse model.

RESULTS

The deficiency of TFF1 was associated with increased EMT of cancer cells in mouse models of pancreatic cancer, KPC. The expression of TFF1 in cancer cells was associated with better survival rate of the patients who underwent chemotherapy, and loss of TFF1 deteriorated the benefit of gemcitabine in KPC mice. Extracellular administration of TFF1 inhibited gemcitabine-induced EMT, Wnt pathway activation and cancer stemness, eventually increased apoptosis of pancreatic cancer cells in vitro. In vivo, combined treatment of gemcitabine and subcutaneous administration of TFF1 arrested tumor growth in xenograft mouse model and resulted in the better survival of KPC mice by inhibiting EMT and cancer stemness.

CONCLUSION

These results indicate that TFF1 can contribute to establishing a novel strategy to treat pancreatic cancer patients by enhancing chemosensitivity.

Loss of the mitochondrial carrier, SLC25A1, during embryogenesis induces a unique senescence program controlled by p53

Germline inactivating mutations of the SLC25A1 gene contribute to various human developmental disorders, including combined D/L-2-hydroxyglutaric aciduria (D/L-2HGA), a severe systemic syndrome characterized by the accumulation of both enantiomers of 2-hydroxyglutaric acid (2HG). The mechanisms by which SLC25A1 deficiency leads to this disease and the role of 2HG are unclear and no therapies exist. We now show that mice lacking both Slc25a1 alleles display a spectrum of alterations that resemble human D/L-2HGA. Mechanistically, SLC25A1 loss results in a proliferation defect and activates two distinct senescence pathways, oncogene-induced senescence (OIS) and mitochondrial dysfunction-induced senescence (MiDAS), both involving the p53 tumor suppressor and driven by two discernible signals: the accumulation of 2HG, inducing OIS, and mitochondrial dysfunction, triggering MiDAS. Inhibiting these senescence programs or blocking p53 activity reverses the growth defect caused by SLC25A1 dysfunction and restores proliferation. These findings reveal novel pathogenic roles of senescence in human disorders and suggest potential strategies to correct the molecular alterations caused by SLC25A1 loss.

Causal network perturbation analysis identifies known and novel type-2 diabetes driver genes

The molecular pathogenesis of diabetes is multifactorial, involving genetic predisposition and environmental factors that are not yet fully understood. However, pancreatic β-cell failure remains among the primary reasons underlying the progression of type-2 diabetes (T2D) making targeting β-cell dysfunction an attractive pathway for diabetes treatment. To identify genetic contributors to β-cell dysfunction, we investigated single-cell gene expression changes in β-cells from healthy (C57BL/6J) and diabetic (NZO/HlLtJ) mice fed with normal or high-fat, high-sugar diet (HFHS). Our study presents an innovative integration of the causal network perturbation assessment (ssNPA) framework with meta-cell transcriptome analysis to explore the genetic underpinnings of type-2 diabetes (T2D). By generating a reference causal network and in silico perturbation, we identified novel genes implicated in T2D and validated our candidates using the Knockout Mouse Phenotyping (KOMP) Project database.

Linking regulatory variants to target genes by integrating single-cell multiome methods and genomic distance

Methods that analyze single-cell paired RNA-seq and ATAC-seq multiome data have shown great promise in linking regulatory elements to genes. However, existing methods differ in their modeling assumptions and approaches to account for biological and technical noise-leading to low concordance in their linking scores-and do not capture the effects of genomic distance. We propose pgBoost, an integrative modeling framework that trains a non-linear combination of existing linking strategies (including genomic distance) on fine-mapped eQTL data to assign a probabilistic score to each candidate SNP-gene link. We applied pgBoost to single-cell multiome data from 85k cells representing 6 major immune/blood cell types. pgBoost attained higher enrichment for fine-mapped eSNP-eGene pairs (e.g. 21x at distance >10kb) than existing methods (1.2-10x; p-value for difference = 5e-13 vs. distance-based method and < 4e-35 for each other method), with larger improvements at larger distances (e.g. 35x vs. 0.89-6.6x at distance >100kb; p-value for difference < 0.002 vs. each other method). pgBoost also outperformed existing methods in enrichment for CRISPR-validated links (e.g. 4.8x vs. 1.6-4.1x at distance >10kb; p-value for difference = 0.25 vs. distance-based method and < 2e-5 for each other method), with larger improvements at larger distances (e.g. 15x vs. 1.6-2.5x at distance >100kb; p-value for difference < 0.009 for each other method). Similar improvements in enrichment were observed for links derived from Activity-By-Contact (ABC) scores and GWAS data. We further determined that restricting pgBoost to features from a focal cell type improved the identification of SNP-gene links relevant to that cell type. We highlight several examples where pgBoost linked fine-mapped GWAS variants to experimentally validated or biologically plausible target genes that were not implicated by other methods. In conclusion, a non-linear combination of linking strategies, including genomic distance, improves power to identify target genes underlying GWAS associations.

Integrating Multiple Database Resources to Elucidate the Gene Flow in Southeast Asian Pig Populations

The domestic pig (Sus scrofa) and its subfamilies have experienced long-term and extensive gene flow, particularly in Southeast Asia. Here, we analyzed 236 pigs, focusing on Yunnan indigenous, European commercial, East Asian, and Southeast Asian breeds, using the Pig Genomics Reference Panel (PGRP v1) of Pig Genotype-Tissue Expression (PigGTEx) to investigate gene flow and associated complex traits by integrating multiple database resources. In this study, we discovered evidence of admixtures from European pigs into the genome of Yunnan indigenous pigs. Additionally, we hypothesized that a potential conceptual gene flow route that may have contributed to the genetic composition of the Diannan small-ear pig is a gene exchange from the Vietnamese pig. Based on the most stringent gene introgression scan using the fd statistic, we identified three specific loci on chromosome 8, ranging from 51.65 to 52.45 Mb, which exhibited strong signatures of selection and harbored the NAF1, NPY1R, and NPY5R genes. These genes are associated with complex traits, such as fat mass, immunity, and litter weight, in pigs, as supported by multiple bio-functionalization databases. We utilized multiple databases to explore the potential dynamics of genetic exchange in Southeast Asian pig populations and elucidated specific gene functionalities.

Mitochondrial citrate carrier SLC25A1 is a dosage-dependent regulator of metabolic reprogramming and morphogenesis in the developing heart

The developing mammalian heart undergoes an important metabolic shift from glycolysis toward mitochondrial oxidation, such that oxidative phosphorylation defects may present with cardiac abnormalities. Here, we describe a new mechanistic link between mitochondria and cardiac morphogenesis, uncovered by studying mice with systemic loss of the mitochondrial citrate carrier SLC25A1. Slc25a1 null embryos displayed impaired growth, cardiac malformations, and aberrant mitochondrial function. Importantly, Slc25a1 heterozygous embryos, which are overtly indistinguishable from wild type, exhibited an increased frequency of these defects, suggesting Slc25a1 haploinsuffiency and dose-dependent effects. Supporting clinical relevance, we found a near-significant association between ultrarare human pathogenic SLC25A1 variants and pediatric congenital heart disease. Mechanistically, SLC25A1 may link mitochondria to transcriptional regulation of metabolism through epigenetic control of gene expression to promote metabolic remodeling in the developing heart. Collectively, this work positions SLC25A1 as a novel mitochondrial regulator of ventricular morphogenesis and cardiac metabolic maturation and suggests a role in congenital heart disease.

A comprehensive analysis of spermatozoal RNA elements in idiopathic infertile males undergoing fertility treatment

Current approaches to diagnosing male infertility inadequately assess the complexity of the male gamete. Beyond the paternal haploid genome, spermatozoa also deliver coding and non-coding RNAs to the oocyte. While sperm-borne RNAs have demonstrated potential involvement in embryo development, the underlying mechanisms remain unclear. In this study, 47 sperm samples from normozoospermic males undergoing fertility treatment using donor oocytes were sequenced and analyzed to evaluate associations between sperm RNA elements (exon-sized sequences) and blastocyst progression. A total of 366 RNA elements (REs) were significantly associated with blastocyst rate (padj < 0.05), some of which were linked to genes related to critical developmental processes, including mitotic spindle formation and both ectoderm and mesoderm specification. Of note, 27 RE-associated RNAs are predicted targets of our previously reported list of developmentally significant miRNAs. Inverse RE-miRNA expression patterns were consistent with miRNA-mediated down-regulation. This study provides a comprehensive set of REs which differ by the patient's ability to produce blastocysts. This knowledge can be leveraged to improve clinical screening of male infertility and ultimately reduce time to pregnancy.

Best practices for managing and disseminating resources and outreach and evaluating the impact of the IDG Consortium

The Illuminating the Druggable Genome (IDG) consortium generated reagents, biological model systems, data, informatic databases, and computational tools. The Resource Dissemination and Outreach Center (RDOC) played a central administrative role, organized internal meetings, fostered collaboration, and coordinated consortium-wide efforts. The RDOC developed and deployed a Resource Management System (RMS) to enable efficient workflows for collecting, accessing, validating, registering, and publishing resource metadata. IDG policies for repositories and standardized representations of resources were established, adopting the FAIR (findable, accessible, interoperable, reusable) principles. The RDOC also developed metrics of IDG impact. Outreach initiatives included digital content, the Protein Illumination Timeline (representing milestones in generating data and reagents), the Target Watch publication series, the e-IDG Symposium series, and leveraging social media platforms.

Loss of KANSL3 leads to defective inner cell mass and early embryonic lethality

e-Lysine acetylation is a prominent histone mark found at transcriptionally active loci. Among many lysine acetyl transferases, nonspecific lethal complex (NSL) members are known to mediate the modification of histone H4. In addition to histone modifications, the KAT8 regulatory complex subunit 3 gene (Kansl3), a core member of NSL complex, has been shown to be involved in several other cellular processes such as mitosis and mitochondrial activity. Although functional studies have been performed on NSL complex members, none of the four core proteins, including Kansl3, have been studied during early mouse development. Here we show that homozygous knockout Kansl3 embryos are lethal at peri-implantation stages, failing to hatch out of the zona pellucida. When the zona pellucida is removed in vitro, Kansl3 null embryos form an abnormal outgrowth with significantly disrupted inner cell mass (ICM) morphology. We document lineage-specific defects at the blastocyst stage with significantly reduced ICM cell number but no difference in trophectoderm cell numbers. Both epiblast and primitive endoderm lineages are altered with reduced cell numbers in null mutants. These results show that Kansl3 is indispensable during early mouse embryonic development and with defects in both ICM and trophectoderm lineages.

Checkpoint kinase 2 controls insulin secretion and glucose homeostasis

After the discovery of insulin, a century ago, extensive work has been done to unravel the molecular network regulating insulin secretion. Here we performed a chemical screen and identified AZD7762, a compound that potentiates glucose-stimulated insulin secretion (GSIS) of a human β cell line, healthy and type 2 diabetic (T2D) human islets and primary cynomolgus macaque islets. In vivo studies in diabetic mouse models and cynomolgus macaques demonstrated that AZD7762 enhances GSIS and improves glucose tolerance. Furthermore, genetic manipulation confirmed that ablation of CHEK2 in human β cells results in increased insulin secretion. Consistently, high-fat-diet-fed Chk2-/- mice show elevated insulin secretion and improved glucose clearance. Finally, untargeted metabolic profiling demonstrated the key role of the CHEK2-PP2A-PLK1-G6PD-PPP pathway in insulin secretion. This study successfully identifies a previously unknown insulin secretion regulating pathway that is conserved across rodents, cynomolgus macaques and human β cells in both healthy and T2D conditions.

Critical assessment of variant prioritization methods for rare disease diagnosis within the rare genomes project

BACKGROUND

A major obstacle faced by families with rare diseases is obtaining a genetic diagnosis. The average "diagnostic odyssey" lasts over five years and causal variants are identified in under 50%, even when capturing variants genome-wide. To aid in the interpretation and prioritization of the vast number of variants detected, computational methods are proliferating. Knowing which tools are most effective remains unclear. To evaluate the performance of computational methods, and to encourage innovation in method development, we designed a Critical Assessment of Genome Interpretation (CAGI) community challenge to place variant prioritization models head-to-head in a real-life clinical diagnostic setting.

METHODS

We utilized genome sequencing (GS) data from families sequenced in the Rare Genomes Project (RGP), a direct-to-participant research study on the utility of GS for rare disease diagnosis and gene discovery. Challenge predictors were provided with a dataset of variant calls and phenotype terms from 175 RGP individuals (65 families), including 35 solved training set families with causal variants specified, and 30 unlabeled test set families (14 solved, 16 unsolved). We tasked teams to identify causal variants in as many families as possible. Predictors submitted variant predictions with estimated probability of causal relationship (EPCR) values. Model performance was determined by two metrics, a weighted score based on the rank position of causal variants, and the maximum F-measure, based on precision and recall of causal variants across all EPCR values.

RESULTS

Sixteen teams submitted predictions from 52 models, some with manual review incorporated. Top performers recalled causal variants in up to 13 of 14 solved families within the top 5 ranked variants. Newly discovered diagnostic variants were returned to two previously unsolved families following confirmatory RNA sequencing, and two novel disease gene candidates were entered into Matchmaker Exchange. In one example, RNA sequencing demonstrated aberrant splicing due to a deep intronic indel in ASNS, identified in trans with a frameshift variant in an unsolved proband with phenotypes consistent with asparagine synthetase deficiency.

CONCLUSIONS

Model methodology and performance was highly variable. Models weighing call quality, allele frequency, predicted deleteriousness, segregation, and phenotype were effective in identifying causal variants, and models open to phenotype expansion and non-coding variants were able to capture more difficult diagnoses and discover new diagnoses. Overall, computational models can significantly aid variant prioritization. For use in diagnostics, detailed review and conservative assessment of prioritized variants against established criteria is needed.

Morphomics via next-generation electron microscopy

The living body is composed of innumerable fine and complex structures. Although these structures have been studied in the past, a vast amount of information pertaining to them still remains unknown. When attempting to observe these ultra-structures, the use of electron microscopy (EM) has become indispensable. However, conventional EM settings are limited to a narrow tissue area, which can bias observations. Recently, new trends in EM research have emerged, enabling coverage of far broader, nano-scale fields of view for two-dimensional wide areas and three-dimensional large volumes. Moreover, cutting-edge bioimage informatics conducted via deep learning has accelerated the quantification of complex morphological bioimages. Taken together, these technological and analytical advances have led to the comprehensive acquisition and quantification of cellular morphology, which now arises as a new omics science termed 'morphomics'.

Network-based analysis predicts interacting genetic modifiers from a meta-mapping study of spike-wave discharge in mice

Absence seizures are characterized by brief lapses in awareness accompanied by a hallmark spike-and-wave discharge (SWD) electroencephalographic pattern and are common to genetic generalized epilepsies (GGEs). While numerous genes have been associated with increased risk, including some Mendelian forms with a single causal allele, most cases of GGE are idiopathic and there are many unknown genetic modifiers of GGE influencing risk and severity. In a previous meta-mapping study, crosses between transgenic C57BL/6 and C3HeB/FeJ strains, each carrying one of three SWD-causing mutations (Gabrg2tm1Spet(R43Q) , Scn8a8j or Gria4spkw1 ), demonstrated an antagonistic epistatic interaction between loci on mouse chromosomes 2 and 7 influencing SWD. These results implicate universal modifiers in the B6 background that mitigate SWD severity through a common pathway, independent of the causal mutation. In this study, we prioritized candidate modifiers in these interacting loci. Our approach integrated human genome-wide association results with gene interaction networks and mouse brain gene expression to prioritize candidate genes and pathways driving variation in SWD outcomes. We considered candidate genes that are functionally associated with human GGE risk genes and genes with evidence for coding or non-coding allele effects between the B6 and C3H backgrounds. Our analyses output a summary ranking of gene pairs, one gene from each locus, as candidates for explaining the epistatic interaction. Our top-ranking gene pairs implicate microtubule function, cytoskeletal stability and cell cycle regulation as novel hypotheses about the source of SWD variation across strain backgrounds, which could clarify underlying mechanisms driving differences in GGE severity in humans.

Systems genetics of influenza A virus-infected mice identifies TRIM21 as a critical regulator of pulmonary innate immune response

Tripartite motif 21 (TRIM21) is a cytosolic Fc receptor that targets antibody-bound, internalized pathogens for destruction. Apart from this intrinsic defense role, TRIM21 is implicated in autoimmune diseases, inflammation, and autophagy. Whether TRIM21 participates in host interactions with influenza A virus (IAV), however, is unknown. By computational modeling of body weight and lung transcriptome data from the BXD parents (C57BL/6 J (B6) and DBA/2 J (D2)) and 41 BXD mouse strains challenged by IAV, we reveal that a Trim21-associated gene network modulates the early host responses to IAV infection. Trim21 transcripts were significantly upregulated in infected mice of both B6 and D2 backgrounds. Its expression was significantly higher in infected D2 than in infected B6 early after infection and significantly correlated with body weight loss. We identified significant trans-eQTL on chromosome 14 that regulates Trim21 expression. Nr1d2 and Il3ra were among the strongest candidate genes. Pathway analysis found Trim21 to be involved in inflammation and immunity related pathways, such as inflammation signaling pathways (TNF, IL-17, and NF-κB), viral detection signaling pathways (NOD-like and RIG-I-like), influenza, and other respiratory viral infections. Knockdown of TRIM21 in human lung epithelial A549 cells significantly augmented IAV-induced expression of IFNB1, IFNL1, CCL5, CXCL10, and IFN-stimulated genes including DDX58 and IFIH1, among others. Our data suggest that a TRIM21-associated gene network is involved in several aspects of inflammation and viral detection mechanisms during IAV infection. We identify and validate TRIM21 as a critical regulator of innate immune responses to IAV in human lung epithelial cells.

Chemical Inhibition of Sterol Biosynthesis

Cholesterol is an essential molecule of life, and its synthesis can be inhibited by both genetic and nongenetic mechanisms. Hundreds of chemicals that we are exposed to in our daily lives can alter sterol biosynthesis. These also encompass various classes of FDA-approved medications, including (but not limited to) commonly used antipsychotic, antidepressant, antifungal, and cardiovascular medications. These medications can interfere with various enzymes of the post-lanosterol biosynthetic pathway, giving rise to complex biochemical changes throughout the body. The consequences of these short- and long-term homeostatic disruptions are mostly unknown. We performed a comprehensive review of the literature and built a catalogue of chemical agents capable of inhibiting post-lanosterol biosynthesis. This process identified significant gaps in existing knowledge, which fall into two main areas: mechanisms by which sterol biosynthesis is altered and consequences that arise from the inhibitions of the different steps in the sterol biosynthesis pathway. The outcome of our review also reinforced that sterol inhibition is an often-overlooked mechanism that can result in adverse consequences and that there is a need to develop new safety guidelines for the use of (novel and already approved) medications with sterol biosynthesis inhibiting side effects, especially during pregnancy.

A novel mouse model for N-terminal truncated Aβ2-x generation through meprin β overexpression in astrocytes

Neurotoxic amyloid-β (Aβ) peptides cause neurodegeneration in Alzheimer's disease (AD) patients' brains. They are released upon proteolytic processing of the amyloid precursor protein (APP) extracellularly at the β-secretase site and intramembranously at the γ-secretase site. Several AD mouse models were developed to conduct respective research in vivo. Most of these classical models overexpress human APP with mutations driving AD-associated pathogenic APP processing. However, the resulting pattern of Aβ species in the mouse brains differs from those observed in AD patients' brains. Particularly mutations proximal to the β-secretase cleavage site (e.g., the so-called Swedish APP (APPswe) fostering Aβ1-x formation) lead to artificial Aβ production, as N-terminally truncated Aβ peptides are hardly present in these mouse brains. Meprin β is an alternative β-secretase upregulated in brains of AD patients and capable of generating N-terminally truncated Aβ2-x peptides. Therefore, we aimed to generate a mouse model for the production of so far underestimated Aβ2-x peptides by conditionally overexpressing meprin β in astrocytes. We chose astrocytes as meprin β was detected in this cell type in close proximity to Aβ plaques in AD patients' brains. The meprin β-overexpressing mice showed elevated amyloidogenic APP processing detected with a newly generated neo-epitope-specific antibody. Furthermore, we observed elevated Aβ production from endogenous APP as well as AD-related behavior changes (hyperlocomotion and deficits in spatial memory). The novel mouse model as well as the established tools and methods will be helpful to further characterize APP cleavage and the impact of different Aβ species in future studies.

Hexosamine biosynthesis and related pathways, protein N-glycosylation and O-GlcNAcylation: their interconnection and role in plants

N-Acetylglucosamine (GlcNAc), a fundamental amino sugar moiety, is essential for protein glycosylation, glycolipid, GPI-anchor protein, and cell wall components. Uridine diphosphate-GlcNAc (UDP-GlcNAc), an active form of GlcNAc, is synthesized through the hexosamine biosynthesis pathway (HBP). Although HBP is highly conserved across organisms, the enzymes involved perform subtly distinct functions among microbes, mammals, and plants. A complete block of HBP normally causes lethality in any life form, reflecting the pivotal role of HBP in the normal growth and development of organisms. Although HBP is mainly composed of four biochemical reactions, HBP is exquisitely regulated to maintain the homeostasis of UDP-GlcNAc content. As HBP utilizes substrates including fructose-6-P, glutamine, acetyl-CoA, and UTP, endogenous nutrient/energy metabolites may be integrated to better suit internal growth and development, and external environmental stimuli. Although the genes encoding HBP enzymes are well characterized in microbes and mammals, they were less understood in higher plants in the past. As the HBP-related genes/enzymes have largely been characterized in higher plants in recent years, in this review we update the latest advances in the functions of the HBP-related genes in higher plants. In addition, HBP's salvage pathway and GlcNAc-mediated two major co- or post-translational modifications, N-glycosylation and O-GlcNAcylation, are also included in this review. Further knowledge on the function of HBP and its product conjugates, and the mechanisms underlying their response to deleterious environments might provide an alternative strategy for agricultural biofortification in the future.

Research advances of polycomb group proteins in regulating mammalian development

Polycomb group (PcG) proteins are a subset of epigenetic factors that are highly conserved throughout evolution. In mammals, PcG proteins can be classified into two muti-proteins complexes: Polycomb repressive complex 1 (PRC1) and PRC2. Increasing evidence has demonstrated that PcG complexes play critical roles in the regulation of gene expression, genomic imprinting, chromosome X-inactivation, and chromatin structure. Accordingly, the dysfunction of PcG proteins is tightly orchestrated with abnormal developmental processes. Here, we summarized and discussed the current knowledge of the biochemical and molecular functions of PcG complexes, especially the PRC1 and PRC2 in mammalian development including embryonic development and tissue development, which will shed further light on the deep understanding of the basic knowledge of PcGs and their functions for reproductive health and developmental disorders.

High Clinical Exome Sequencing Diagnostic Rates and Novel Phenotypic Expansions for Nonisolated Microphthalmia, Anophthalmia, and Coloboma

Purpose

A molecular diagnosis is only made in a subset of individuals with nonisolated microphthalmia, anophthalmia, and coloboma (MAC). This may be due to underutilization of clinical (whole) exome sequencing (cES) and an incomplete understanding of the genes that cause MAC. The purpose of this study is to determine the efficacy of cES in cases of nonisolated MAC and to identify new MAC phenotypic expansions.

Methods

We determined the efficacy of cES in 189 individuals with nonisolated MAC. We then used cES data, a validated machine learning algorithm, and previously published expression data, case reports, and animal models to determine which candidate genes were most likely to contribute to the development of MAC.

Results

We found the efficacy of cES in nonisolated MAC to be between 32.3% (61/189) and 48.1% (91/189). Most genes affected in our cohort were not among genes currently screened in clinically available ophthalmologic gene panels. A subset of the genes implicated in our cohort had not been clearly associated with MAC. Our analyses revealed sufficient evidence to support low-penetrance MAC phenotypic expansions involving nine of these human disease genes.

Conclusions

We conclude that cES is an effective means of identifying a molecular diagnosis in individuals with nonisolated MAC and may identify putatively damaging variants that would be missed if only a clinically available ophthalmologic gene panel was obtained. Our data also suggest that deleterious variants in BRCA2, BRIP1, KAT6A, KAT6B, NSF, RAC1, SMARCA4, SMC1A, and TUBA1A can contribute to the development of MAC.

Metabolic Deficiencies Underlie Plasmacytoid Dendritic Cell Exhaustion After Viral Infection

Type I Interferons (IFN-I) are central to host protection against viral infections 1 . While any cell can produce IFN-I, Plasmacytoid Dendritic Cells (pDCs) make greater quantities and more varieties of these cytokines than any other cell type 2 . However, following an initial burst of IFN- I, pDCs lose their exceptional IFN-I production capacity and become "exhausted", a phenotype that associates with enhanced susceptibility to secondary infections 3-5 . Despite this apparent cost for the host, pDC exhaustion is conserved across multiple species and viral infections, but the underlying mechanisms and the potential evolutionary advantages are not well understood. Here we characterize pDC exhaustion and demonstrate that it is associated with a reduced capacity of pDCs to engage both oxidative and glycolytic metabolism. Mechanistically, we identify lactate dehydrogenase B (LDHB) as a novel positive regulator of pDC IFN-I production in mice and humans, show that LDHB deficiency is associated with suppressed IFN-I production, pDC metabolic capacity, and viral control following a viral infection, and demonstrate that preservation of LDHB expression is sufficient to partially restore exhausted pDC function in vitro and in vivo . Furthermore, restoring LDHB in vivo in exhausted pDCs increased IFNAR dependent infection- associated pathology. Therefore, our work identifies a novel and conserved mechanism for balancing immunity and pathology during viral infections, while also providing insight into the highly preserved but previously unexplained phenomenon of pDC exhaustion.

FBXL16 promotes cell growth and drug resistance in lung adenocarcinomas with KRAS mutation by stabilizing IRS1 and upregulating IRS1/AKT signaling

Lung cancer is the leading cause of cancer-related deaths worldwide. Lung adenocarcinomas (LUADs) are a major subtype of non-small-cell lung cancers (NSCLCs). About 25% of LUADs harbor GTPase KRAS mutations associated with poor prognosis and limited treatment options. While encouraging tumor response to novel covalent inhibitors specifically targeting KRASG12C has been shown in the clinic, either intrinsic resistance exists or acquired therapeutic resistance arises upon treatment. There is an unmet need to identify new therapeutic targets for treating LUADs with activating KRAS mutations, particularly those with resistance to KRASG12C inhibitor(s). In this study, we have revealed that F-box/LRR-repeat protein 16 (FBXL16) is selectively upregulated in LUAD with KRAS mutations. It promotes LUAD cell growth and transforms lung epithelial cells. Importantly, FBXL16 depletion greatly enhances sensitivity to the KRASG12C inhibitor (sotorasib) in resistant cells by downregulating phosphatidylinositol 3-kinase (PI3K)/protein kinase B (PKB; also known as AKT) signaling. Mechanistically, FBXL16 upregulates insulin receptor substrate 1 (IRS1) protein stability, leading to an increase of IGF1/AKT signaling, thereby promoting cell growth and migration. Taken together, our study highlights the potential of FBXL16 as a therapeutic target for treating LUAD with KRAS activating mutations.

Genetic determinants of micronucleus formation in vivo

Genomic instability arising from defective responses to DNA damage1 or mitotic chromosomal imbalances2 can lead to the sequestration of DNA in aberrant extranuclear structures called micronuclei (MN). Although MN are a hallmark of ageing and diseases associated with genomic instability, the catalogue of genetic players that regulate the generation of MN remains to be determined. Here we analyse 997 mouse mutant lines, revealing 145 genes whose loss significantly increases (n = 71) or decreases (n = 74) MN formation, including many genes whose orthologues are linked to human disease. We found that mice null for Dscc1, which showed the most significant increase in MN, also displayed a range of phenotypes characteristic of patients with cohesinopathy disorders. After validating the DSCC1-associated MN instability phenotype in human cells, we used genome-wide CRISPR-Cas9 screening to define synthetic lethal and synthetic rescue interactors. We found that the loss of SIRT1 can rescue phenotypes associated with DSCC1 loss in a manner paralleling restoration of protein acetylation of SMC3. Our study reveals factors involved in maintaining genomic stability and shows how this information can be used to identify mechanisms that are relevant to human disease biology1.

Exosome complex components 1 and 2 are vital for early mammalian development

Exosome Complex Components 1 and 2 (EXOSC1 and 2) are two proteins in the RNA Exosome complex whose main function is 5' → 3' RNA degradation and processing. The RNA exosome complex is comprised of nine subunits that form two separate components: the S1/KH cap and the PH-core. EXOSC1 and 2 are both part of the S1/KH cap and are involved in binding nascent RNA. As part of a systemic characterization of early lethal alleles produced by the Knockout Mouse Project, we have examined Exosc1 and Exosc2 homozygous null (mutant) embryos to determine developmental and molecular phenotypes of embryos lacking their functions. Our studies reveal that Exosc1 null embryos implant and form an egg cylinder but are developmentally delayed and fail to initiate gastrulation by embryonic day 7.5. In contrast, Exosc2 null embryos are lethal during peri-implantation stages, and while they do form a morphologically normal blastocyst at E3.5, they cannot be recovered at post-implantation stages. We show the absence of stage-specific developmental and altered lineage-specification in both Exosc1 and Exosc2 mutant embryos and conclude that these genes are essential for the successful progression through early mammalian development.

Functional genomics in inborn errors of immunity

Inborn errors of immunity (IEI) comprise a diverse spectrum of 485 disorders as recognized by the International Union of Immunological Societies Committee on Inborn Error of Immunity in 2022. While IEI are monogenic by definition, they illuminate various pathways involved in the pathogenesis of polygenic immune dysregulation as in autoimmune or autoinflammatory syndromes, or in more common infectious diseases that may not have a significant genetic basis. Rapid improvement in genomic technologies has been the main driver of the accelerated rate of discovery of IEI and has led to the development of innovative treatment strategies. In this review, we will explore various facets of IEI, delving into the distinctions between PIDD and PIRD. We will examine how Mendelian inheritance patterns contribute to these disorders and discuss advancements in functional genomics that aid in characterizing new IEI. Additionally, we will explore how emerging genomic tools help to characterize new IEI as well as how they are paving the way for innovative treatment approaches for managing and potentially curing these complex immune conditions.

Female behavior drives the formation of distinct social structures in C57BL/6J versus wild-derived outbred mice in field enclosures

BACKGROUND

Social behavior and social organization have major influences on individual health and fitness. Yet, biomedical research focuses on studying a few genotypes under impoverished social conditions. Understanding how lab conditions have modified social organizations of model organisms, such as lab mice, relative to natural populations is a missing link between socioecology and biomedical science.

RESULTS

Using a common garden design, we describe the formation of social structure in the well-studied laboratory mouse strain, C57BL/6J, in replicated mixed-sex populations over 10-day trials compared to control trials with wild-derived outbred house mice in outdoor field enclosures. We focus on three key features of mouse social systems: (i) territory establishment in males, (ii) female social relationships, and (iii) the social networks formed by the populations. Male territorial behaviors were similar but muted in C57 compared to wild-derived mice. Female C57 sharply differed from wild-derived females, showing little social bias toward cage mates and exploring substantially more of the enclosures compared to all other groups. Female behavior consistently generated denser social networks in C57 than in wild-derived mice.

CONCLUSIONS

C57 and wild-derived mice individually vary in their social and spatial behaviors which scale to shape overall social organization. The repeatable societies formed under field conditions highlights opportunities to experimentally study the interplay between society and individual biology using model organisms.

Compromised transcription-mRNA export factor THOC2 causes R-loop accumulation, DNA damage and adverse neurodevelopment

We implicated the X-chromosome THOC2 gene, which encodes the largest subunit of the highly-conserved TREX (Transcription-Export) complex, in a clinically complex neurodevelopmental disorder with intellectual disability as the core phenotype. To study the molecular pathology of this essential eukaryotic gene, we generated a mouse model based on a hypomorphic Thoc2 exon 37-38 deletion variant of a patient with ID, speech delay, hypotonia, and microcephaly. The Thoc2 exon 37-38 deletion male (Thoc2Δ/Y) mice recapitulate the core phenotypes of THOC2 syndrome including smaller size and weight, and significant deficits in spatial learning, working memory and sensorimotor functions. The Thoc2Δ/Y mouse brain development is significantly impacted by compromised THOC2/TREX function resulting in R-loop accumulation, DNA damage and consequent cell death. Overall, we suggest that perturbed R-loop homeostasis, in stem cells and/or differentiated cells in mice and the patient, and DNA damage-associated functional alterations are at the root of THOC2 syndrome.

Unraveling the genetic architecture of congenital vertebral malformation with reference to the developing spine

Congenital vertebral malformation, affecting 0.13-0.50 per 1000 live births, has an immense locus heterogeneity and complex genetic architecture. In this study, we analyze exome/genome sequencing data from 873 probands with congenital vertebral malformation and 3794 control individuals. Clinical interpretation identifies Mendelian etiologies in 12.0% of the probands and reveals a muscle-related disease mechanism. Gene-based burden test of ultra-rare variants identifies risk genes with large effect sizes (ITPR2, TBX6, TPO, H6PD, and SEC24B). To further investigate the biological relevance of the genetic association signals, we perform single-nucleus RNAseq on human embryonic spines. The burden test signals are enriched in the notochord at early developmental stages and myoblast/myocytes at late stages, highlighting their critical roles in the developing spine. Our work provides insights into the developmental biology of the human spine and the pathogenesis of spine malformation.

The Vps13-like protein BLTP2 is pro-survival and regulates phosphatidylethanolamine levels in the plasma membrane to maintain its fluidity and function

Lipid transport proteins (LTPs) facilitate nonvesicular lipid exchange between cellular compartments and have critical roles in lipid homeostasis1. A new family of bridge-like LTPs (BLTPs) is thought to form lipid-transporting conduits between organelles2. One, BLTP2, is conserved across species but its function is not known. Here, we show that BLTP2 and its homolog directly regulate plasma membrane (PM) fluidity by increasing the phosphatidylethanolamine (PE) level in the PM. BLTP2 localizes to endoplasmic reticulum (ER)-PM contact sites34, 5, suggesting it transports PE from the ER to the PM. We find BLTP2 works in parallel with another pathway that regulates intracellular PE distribution and PM fluidity6, 7. BLTP2 expression correlates with breast cancer aggressiveness8-10. We found BLTP2 facilitates growth of a human cancer cell line and sustains its aggressiveness in an in vivo model of metastasis, suggesting maintenance of PM fluidity by BLTP2 may be critical for tumorigenesis in humans.

Miscarriage risk assessment: a bioinformatic approach to identifying candidate lethal genes and variants

PURPOSE

Miscarriage, often resulting from a variety of genetic factors, is a common pregnancy outcome. Preconception genetic carrier screening (PGCS) identifies at-risk partners for newborn genetic disorders; however, PGCS panels currently lack miscarriage-related genes. In this study, we evaluated the potential impact of both known and candidate genes on prenatal lethality and the effectiveness of PGCS in diverse populations.

METHODS

We analyzed 125,748 human exome sequences and mouse and human gene function databases. Our goals were to identify genes crucial for human fetal survival (lethal genes), to find variants not present in a homozygous state in healthy humans, and to estimate carrier rates of known and candidate lethal genes in various populations and ethnic groups.

RESULTS

This study identified 138 genes in which heterozygous lethal variants are present in the general population with a frequency of 0.5% or greater. Screening for these 138 genes could identify 4.6% (in the Finnish population) to 39.8% (in the East Asian population) of couples at risk of miscarriage. This explains the cause of pregnancy loss in approximately 1.1-10% of cases affected by biallelic lethal variants.

CONCLUSION

This study has identified a set of genes and variants potentially associated with lethality across different ethnic backgrounds. The variation of these genes across ethnic groups underscores the need for a comprehensive, pan-ethnic PGCS panel that includes genes related to miscarriage.

Mapping Cell Atlases at the Single-Cell Level

Recent advancements in single-cell technologies have led to rapid developments in the construction of cell atlases. These atlases have the potential to provide detailed information about every cell type in different organisms, enabling the characterization of cellular diversity at the single-cell level. Global efforts in developing comprehensive cell atlases have profound implications for both basic research and clinical applications. This review provides a broad overview of the cellular diversity and dynamics across various biological systems. In addition, the incorporation of machine learning techniques into cell atlas analyses opens up exciting prospects for the field of integrative biology.

Integrative functional genomic analyses identify genetic variants influencing skin pigmentation in Africans

Skin color is highly variable in Africans, yet little is known about the underlying molecular mechanism. Here we applied massively parallel reporter assays to screen 1,157 candidate variants influencing skin pigmentation in Africans and identified 165 single-nucleotide polymorphisms showing differential regulatory activities between alleles. We combine Hi-C, genome editing and melanin assays to identify regulatory elements for MFSD12, HMG20B, OCA2, MITF, LEF1, TRPS1, BLOC1S6 and CYB561A3 that impact melanin levels in vitro and modulate human skin color. We found that independent mutations in an OCA2 enhancer contribute to the evolution of human skin color diversity and detect signals of local adaptation at enhancers of MITF, LEF1 and TRPS1, which may contribute to the light skin color of Khoesan-speaking populations from Southern Africa. Additionally, we identified CYB561A3 as a novel pigmentation regulator that impacts genes involved in oxidative phosphorylation and melanogenesis. These results provide insights into the mechanisms underlying human skin color diversity and adaptive evolution.

Mechanisms underlying Myosin 10's contribution to the maintenance of mitotic spindle bipolarity

Myosin 10 (Myo10) couples microtubules and integrin-based adhesions to movement along actin filaments via its microtubule-binding MyTH4 domain and integrin-binding FERM domain, respectively. Here we show that Myo10-depleted HeLa cells and mouse embryo fibroblasts (MEFs) both exhibit a pronounced increase in the frequency of multipolar spindles. Staining of unsynchronized metaphase cells showed that the primary driver of spindle multipolarity in Myo10-depleted MEFs and in Myo10-depleted HeLa cells lacking supernumerary centrosomes is pericentriolar material (PCM) fragmentation, which creates y-tubulin-positive acentriolar foci that serve as extra spindle poles. For HeLa cells possessing supernumerary centrosomes, Myo10 depletion further accentuates spindle multipolarity by impairing the clustering of the extra spindle poles. Complementation experiments show that Myo10 must interact with both microtubules and integrins to promote PCM/pole integrity. Conversely, Myo10 only needs interact with integrins to promote supernumerary centrosome clustering. Importantly, images of metaphase Halo-Myo10 knockin cells show that the myosin localizes exclusively to the spindle and the tips of adhesive retraction fibers. We conclude that Myo10 promotes PCM/pole integrity in part by interacting with spindle microtubules, and that it promotes supernumerary centrosome clustering by supporting retraction fiber-based cell adhesion, which likely serves to anchor the microtubule-based forces driving pole focusing.

Ubiquitin ligase RNF20 coordinates sequential adipose thermogenesis with brown and beige fat-specific substrates

In mammals, brown adipose tissue (BAT) and inguinal white adipose tissue (iWAT) execute sequential thermogenesis to maintain body temperature during cold stimuli. BAT rapidly generates heat through brown adipocyte activation, and further iWAT gradually stimulates beige fat cell differentiation upon prolonged cold challenges. However, fat depot-specific regulatory mechanisms for thermogenic activation of two fat depots are poorly understood. Here, we demonstrate that E3 ubiquitin ligase RNF20 orchestrates adipose thermogenesis with BAT- and iWAT-specific substrates. Upon cold stimuli, BAT RNF20 is rapidly downregulated, resulting in GABPα protein elevation by controlling protein stability, which stimulates thermogenic gene expression. Accordingly, BAT-specific Rnf20 suppression potentiates BAT thermogenic activity via GABPα upregulation. Moreover, upon prolonged cold stimuli, iWAT RNF20 is gradually upregulated to promote de novo beige adipogenesis. Mechanistically, iWAT RNF20 mediates NCoR1 protein degradation, rather than GABPα, to activate PPARγ. Together, current findings propose fat depot-specific regulatory mechanisms for temporal activation of adipose thermogenesis.

RNA-binding protein Nocte regulates Drosophila development by promoting translation reinitiation on mRNAs with long upstream open reading frames

RNA-binding proteins (RBPs) with intrinsically disordered regions (IDRs) are linked to multiple human disorders, but their mechanisms of action remain unclear. Here, we report that one such protein, Nocte, is essential for Drosophila eye development by regulating a critical gene expression cascade at translational level. Knockout of nocte in flies leads to lethality, and its eye-specific depletion impairs eye size and morphology. Nocte preferentially enhances translation of mRNAs with long upstream open reading frames (uORFs). One of the key Nocte targets, glass mRNA, encodes a transcription factor critical for differentiation of photoreceptor neurons and accessory cells, and re-expression of Glass largely rescued the eye defects caused by Nocte depletion. Mechanistically, Nocte counteracts long uORF-mediated translational suppression by promoting translation reinitiation downstream of the uORF. Nocte interacts with translation factors eIF3 and Rack1 through its BAT2 domain, and a Nocte mutant lacking this domain fails to promote translation of glass mRNA. Notably, de novo mutations of human orthologs of Nocte have been detected in schizophrenia patients. Our data suggest that Nocte family of proteins can promote translation reinitiation to overcome long uORFs-mediated translational suppression, and disruption of this function can lead to developmental defects and neurological disorders.

The Dynamics of Histone Modifications during Mammalian Zygotic Genome Activation

Mammalian fertilization initiates the reprogramming of oocytes and sperm, forming a totipotent zygote. During this intricate process, the zygotic genome undergoes a maternal-to-zygotic transition (MZT) and subsequent zygotic genome activation (ZGA), marking the initiation of transcriptional control and gene expression post-fertilization. Histone modifications are pivotal in shaping cellular identity and gene expression in many mammals. Recent advances in chromatin analysis have enabled detailed explorations of histone modifications during ZGA. This review delves into conserved and unique regulatory strategies, providing essential insights into the dynamic changes in histone modifications and their variants during ZGA in mammals. The objective is to explore recent advancements in leading mechanisms related to histone modifications governing this embryonic development phase in depth. These considerations will be useful for informing future therapeutic approaches that target epigenetic regulation in diverse biological contexts. It will also contribute to the extensive areas of evolutionary and developmental biology and possibly lay the foundation for future research and discussion on this seminal topic.

Augmenting tomato functional genomics with a genome-wide induced genetic variation resource

Induced mutations accelerate crop improvement by providing novel disease resistance and yield alleles. However, the alleles with no perceptible phenotype but have an altered function remain hidden in mutagenized plants. The whole-genome sequencing (WGS) of mutagenized individuals uncovers the complete spectrum of mutations in the genome. Genome-wide induced mutation resources can improve the targeted breeding of tomatoes and facilitate functional genomics. In this study, we sequenced 132 doubly ethyl methanesulfonate (EMS)-mutagenized lines of tomato and detected approximately 41 million novel mutations and 5.5 million short InDels not present in the parental cultivar. Approximately 97% of the genome had mutations, including the genes, promoters, UTRs, and introns. More than one-third of genes in the mutagenized population had one or more deleterious mutations predicted by Sorting Intolerant From Tolerant (SIFT). Nearly one-fourth of deleterious genes mapped on tomato metabolic pathways modulate multiple pathway steps. In addition to the reported GC>AT transition bias for EMS, our population also had a substantial number of AT>GC transitions. Comparing mutation frequency among synonymous codons revealed that the most preferred codon is the least mutagenic toward EMS. The validation of a potato leaf-like mutation, reduction in carotenoids in ζ-carotene isomerase mutant fruits, and chloroplast relocation loss in phototropin1 mutant validated the mutation discovery pipeline. Our database makes a large repertoire of mutations accessible to functional genomics studies and breeding of tomatoes.