The Nrf family of transcription factors plays a critical role in mediating adaptive responses to cellular stress and defends against neurodegeneration, aging, and cancer. Here, we report a novel role for the Caenorhabditis elegans Nrf homolog SKN-1 in regulating synaptic transmission at neuromuscular junctions (NMJs). Activation of SKN-1, either by acute pharmacological treatment with the mitochondrial toxin sodium arsenite or by mutations that cause constitutive SKN-1 activation, results in defects in neuromuscular function. Additionally, elimination of the conserved WD40 repeat protein WDR-23, a principal negative regulator of SKN-1, results in impaired locomotion and synaptic vesicle and neuropeptide release from cholinergic motor axons. Mutations that abolish skn-1 activity restore normal neuromuscular function to wdr-23 mutants and animals treated with toxin. We show that negative regulation of SKN-1 by WDR-23 in the intestine, but not at neuromuscular junctions, is necessary and sufficient for proper neuromuscular function. WDR-23 isoforms differentially localize to the outer membranes of mitochondria and to nuclei, and the effects of WDR-23 on neuromuscular function are dependent on its interaction with cullin E3 ubiquitin ligase. Finally, whole-transcriptome RNA sequencing of wdr-23 mutants reveals an increase in the expression of known SKN-1/Nrf2-regulated stress-response genes, as well as neurotransmission genes not previously implicated in SKN-1/Nrf2 responses. Together, our results indicate that SKN-1/Nrf2 activation may be a mechanism through which cellular stress, detected in one tissue, affects cellular function of a distal tissue through endocrine signaling. These results provide insight into how SKN-1/Nrf2 might protect the nervous system from damage in response to oxidative stress.

The large yellow croaker Larimichthys crocea (L. crocea) is one of the most economically important marine fish in China and East Asian countries. It also exhibits peculiar behavioral and physiological characteristics, especially sensitive to various environmental stresses, such as hypoxia and air exposure. These traits may render L. crocea a good model for investigating the response mechanisms to environmental stress. To understand the molecular and genetic mechanisms underlying the adaptation and response of L. crocea to environmental stress, we sequenced and assembled the genome of L. crocea using a bacterial artificial chromosome and whole-genome shotgun hierarchical strategy. The final genome assembly was 679 Mb, with a contig N50 of 63.11 kb and a scaffold N50 of 1.03 Mb, containing 25,401 protein-coding genes. Gene families underlying adaptive behaviours, such as vision-related crystallins, olfactory receptors, and auditory sense-related genes, were significantly expanded in the genome of L. crocea relative to those of other vertebrates. Transcriptome analyses of the hypoxia-exposed L. crocea brain revealed new aspects of neuro-endocrine-immune/metabolism regulatory networks that may help the fish to avoid cerebral inflammatory injury and maintain energy balance under hypoxia. Proteomics data demonstrate that skin mucus of the air-exposed L. crocea had a complex composition, with an unexpectedly high number of proteins (3,209), suggesting its multiple protective mechanisms involved in antioxidant functions, oxygen transport, immune defence, and osmotic and ionic regulation. Our results reveal the molecular and genetic basis of fish adaptation and response to hypoxia and air exposure. The data generated by this study will provide valuable resources for the genetic improvement of stress resistance and yield potential in L. crocea.

During ischemic stroke, occlusion of the cerebrovasculature causes neuronal cell death (infarction), but naturally occurring genetic factors modulating infarction have been difficult to identify in human populations. In a surgically induced mouse model of ischemic stroke, we have previously mapped Civq1 to distal chromosome 7 as a quantitative trait locus determining infarct volume. In this study, genome-wide association mapping using 32 inbred mouse strains and an additional linkage scan for infarct volume confirmed that the size of the infarct is determined by ancestral alleles of the causative gene(s). The genetically isolated Civq1 locus in reciprocal recombinant congenic mice refined the critical interval and demonstrated that infarct size is determined by both vascular (collateral vessel anatomy) and non-vascular (neuroprotection) effects. Through the use of interval-specific SNP haplotype analysis, we further refined the Civq1 locus and identified integrin alpha L (Itgal) as one of the causative genes for Civq1. Itgal is the only gene that exhibits both strain-specific amino acid substitutions and expression differences. Coding SNPs, a 5-bp insertion in exon 30b, and increased mRNA and protein expression of a splice variant of the gene (Itgal-003, ENSMUST00000120857), all segregate with infarct volume. Mice lacking Itgal show increased neuronal cell death in both ex vivo brain slice and in vivo focal cerebral ischemia. Our data demonstrate that sequence variation in Itgal modulates ischemic brain injury, and that infarct volume is determined by both vascular and non-vascular mechanisms.

Significant advances have been made in the discovery of genes affecting bone mineral density (BMD); however, our understanding of its genetic basis remains incomplete. In the current study, genome-wide association (GWA) and co-expression network analysis were used in the recently described Hybrid Mouse Diversity Panel (HMDP) to identify and functionally characterize novel BMD genes. In the HMDP, a GWA of total body, spinal, and femoral BMD revealed four significant associations (-log10P>5.39) affecting at least one BMD trait on chromosomes (Chrs.) 7, 11, 12, and 17. The associations implicated a total of 163 genes with each association harboring between 14 and 112 genes. This list was reduced to 26 functional candidates by identifying those genes that were regulated by local eQTL in bone or harbored potentially functional non-synonymous (NS) SNPs. This analysis revealed that the most significant BMD SNP on Chr. 12 was a NS SNP in the additional sex combs like-2 (Asxl2) gene that was predicted to be functional. The involvement of Asxl2 in the regulation of bone mass was confirmed by the observation that Asxl2 knockout mice had reduced BMD. To begin to unravel the mechanism through which Asxl2 influenced BMD, a gene co-expression network was created using cortical bone gene expression microarray data from the HMDP strains. Asxl2 was identified as a member of a co-expression module enriched for genes involved in the differentiation of myeloid cells. In bone, osteoclasts are bone-resorbing cells of myeloid origin, suggesting that Asxl2 may play a role in osteoclast differentiation. In agreement, the knockdown of Asxl2 in bone marrow macrophages impaired their ability to form osteoclasts. This study identifies a new regulator of BMD and osteoclastogenesis and highlights the power of GWA and systems genetics in the mouse for dissecting complex genetic traits.

Females with polycystic ovary syndrome (PCOS), the most common endocrine disorder in women, have an increased risk of developing cardiometabolic disorders such as insulin resistance, obesity, and type 2 diabetes (T2D). While only diagnosable in females, males with a family history of PCOS can also exhibit a poor cardiometabolic profile. Therefore, we aimed to elucidate the role of sex in the cardiometabolic comorbidities observed in PCOS by conducting bidirectional genetic risk score analyses in both sexes. We first conducted a phenome-wide association study (PheWAS) using PCOS polygenic risk scores (PCOSPRS) to identify potential pleiotropic effects of PCOSPRS across 1,380 medical conditions recorded in the Vanderbilt University Medical Center electronic health record (EHR) database, in females and males. After adjusting for age and genetic ancestry, we found that European (EUR)-ancestry males with higher PCOSPRS were significantly more likely to develop hypertensive diseases than females at the same level of genetic risk. We performed the same analysis in an African (AFR)-ancestry population, but observed no significant associations, likely due to poor trans-ancestry performance of the PRS. Based on observed significant associations in the EUR-ancestry population, we then tested whether the PRS for comorbid conditions (e.g., T2D, body mass index (BMI), hypertension, etc.) also increased the odds of a PCOS diagnosis. Only BMIPRS and T2DPRS were significantly associated with a PCOS diagnosis in EUR-ancestry females. We then further adjusted the T2DPRS for measured BMI and BMIresidual (regressed on the BMIPRS and enriched for the environmental contribution to BMI). Results demonstrated that genetically regulated BMI primarily accounted for the relationship between T2DPRS and PCOS. Overall, our findings show that the genetic architecture of PCOS has distinct sex differences in associations with genetically correlated cardiometabolic traits. It is possible that the cardiometabolic comorbidities observed in PCOS are primarily explained by their shared genetic risk factors, which can be further influenced by biological variables including sex and BMI.

The liver and pancreas originate from overlapping embryonic regions, and single-cell lineage tracing in zebrafish has shown that Bone morphogenetic protein 2b (Bmp2b) signaling is essential for determining the fate of bipotential hepatopancreatic progenitors towards the liver or pancreas. Despite its pivotal role, the gene regulatory networks functioning downstream of Bmp2b signaling in this process are poorly understood. We have identified four and a half LIM domains 1b (fhl1b), which is primarily expressed in the prospective liver anlage, as a novel target of Bmp2b signaling. fhl1b depletion compromised liver specification and enhanced induction of pancreatic cells from endodermal progenitors. Conversely, overexpression of fhl1b favored liver specification and inhibited induction of pancreatic cells. By single-cell lineage tracing, we showed that fhl1b depletion led lateral endodermal cells, destined to become liver cells, to become pancreatic cells. Reversely, when fhl1b was overexpressed, medially located endodermal cells, fated to differentiate into pancreatic and intestinal cells, contributed to the liver by directly or indirectly modulating the discrete levels of pdx1 expression in endodermal progenitors. Moreover, loss of fhl1b increased the regenerative capacity of β-cells by increasing pdx1 and neurod expression in the hepatopancreatic ductal system. Altogether, these data reveal novel and critical functions of Fhl1b in the hepatic versus pancreatic fate decision and in β-cell regeneration.

Danforth's short tail mutant (Sd) mouse, first described in 1930, is a classic spontaneous mutant exhibiting defects of the axial skeleton, hindgut, and urogenital system. We used meiotic mapping in 1,497 segregants to localize the mutation to a 42.8-kb intergenic segment on chromosome 2. Resequencing of this region identified an 8.5-kb early retrotransposon (ETn) insertion within the highly conserved regulatory sequences upstream of Pancreas Specific Transcription Factor, 1a (Ptf1a). This mutation resulted in up to tenfold increased expression of Ptf1a as compared to wild-type embryos at E9.5 but no detectable changes in the expression levels of other neighboring genes. At E9.5, Sd mutants exhibit ectopic Ptf1a expression in embryonic progenitors of every organ that will manifest a developmental defect: the notochord, the hindgut, and the mesonephric ducts. Moreover, at E 8.5, Sd mutant mice exhibit ectopic Ptf1a expression in the lateral plate mesoderm, tail bud mesenchyme, and in the notochord, preceding the onset of visible defects such as notochord degeneration. The Sd heterozygote phenotype was not ameliorated by Ptf1a haploinsufficiency, further suggesting that the developmental defects result from ectopic expression of Ptf1a. These data identify disruption of the spatio-temporal pattern of Ptf1a expression as the unifying mechanism underlying the multiple congenital defects in Danforth's short tail mouse. This striking example of an enhancer mutation resulting in profound developmental defects suggests that disruption of conserved regulatory elements may also contribute to human malformation syndromes.

Experimental autoimmune orchitis (EAO), the principal model of non-infectious testicular inflammatory disease, can be induced in susceptible mouse strains by immunization with autologous testicular homogenate and appropriate adjuvants. As previously established, the genome of DBA/2J mice encodes genes that are capable of conferring dominant resistance to EAO, while the genome of BALB/cByJ mice does not and they are therefore susceptible to EAO. In a genome scan, we previously identified Orch3 as the major quantitative trait locus controlling dominant resistance to EAO and mapped it to chromosome 11. Here, by utilizing a forward genetic approach, we identified kinesin family member 1C (Kif1c) as a positional candidate for Orch3 and, using a transgenic approach, demonstrated that Kif1c is Orch3. Mechanistically, we showed that the resistant Kif1c(D2) allele leads to a reduced antigen-specific T cell proliferative response as a consequence of decreased MHC class II expression by antigen presenting cells, and that the L(578) → P(578) and S(1027) → P(1027) polymorphisms distinguishing the BALB/cByJ and DBA/2J alleles, respectively, can play a role in transcriptional regulation. These findings may provide mechanistic insight into how polymorphism in other kinesins such as KIF21B and KIF5A influence susceptibility and resistance to human autoimmune diseases.

Hybrid sterility (HS) belongs to reproductive isolation barriers that safeguard the integrity of species in statu nascendi. Although hybrid sterility occurs almost universally among animal and plant species, most of our current knowledge comes from the classical genetic studies on Drosophila interspecific crosses or introgressions. With the house mouse subspecies Mus m. musculus and Mus m. domesticus as a model, new research tools have become available for studies of the molecular mechanisms and genetic networks underlying HS. Here we used QTL analysis and intersubspecific chromosome substitution strains to identify a 4.7 Mb critical region on Chromosome X (Chr X) harboring the Hstx2 HS locus, which causes asymmetrical spermatogenic arrest in reciprocal intersubspecific F1 hybrids. Subsequently, we mapped autosomal loci on Chrs 3, 9 and 13 that can abolish this asymmetry. Combination of immunofluorescent visualization of the proteins of synaptonemal complexes with whole-chromosome DNA FISH on pachytene spreads revealed that heterosubspecific, unlike consubspecific, homologous chromosomes are predisposed to asynapsis in F1 hybrid male and female meiosis. The asynapsis is under the trans- control of Hstx2 and Hst1/Prdm9 hybrid sterility genes in pachynemas of male but not female hybrids. The finding concurred with the fertility of intersubpecific F1 hybrid females homozygous for the Hstx2(Mmm) allele and resolved the apparent conflict with the dominance theory of Haldane's rule. We propose that meiotic asynapsis in intersubspecific hybrids is a consequence of cis-acting mismatch between homologous chromosomes modulated by the trans-acting Hstx2 and Prdm9 hybrid male sterility genes.

Understanding mammalian genetic systems is predicated on the determination of the relationship between genetic variation and phenotype. Several international programmes are under way to deliver mutations in every gene in the mouse genome. The challenge for mouse geneticists is to develop approaches that will provide comprehensive phenotype datasets for these mouse mutant libraries. Several factors are critical to success in this endeavour. It will be important to catalogue assay and environment and where possible to adopt standardised procedures for phenotyping tests along with common environmental conditions to ensure comparable datasets of phenotypes. Moreover, the scale of the task underlines the need to invest in technological development improving both the speed and cost of phenotyping platforms. In addition, it will be necessary to develop new informatics standards that capture the phenotype assay as well as other factors, genetic and environmental, that impinge upon phenotype outcome.

Periventricular nodular heterotopia (PVNH) is a malformation of cortical development commonly associated with epilepsy. We exome sequenced 202 individuals with sporadic PVNH to identify novel genetic risk loci. We first performed a trio-based analysis and identified 219 de novo variants. Although no novel genes were implicated in this initial analysis, PVNH cases were found overall to have a significant excess of nonsynonymous de novo variants in intolerant genes (p = 3.27x10-7), suggesting a role for rare new alleles in genes yet to be associated with the condition. Using a gene-level collapsing analysis comparing cases and controls, we identified a genome-wide significant signal driven by four ultra-rare loss-of-function heterozygous variants in MAP1B, including one de novo variant. In at least one instance, the MAP1B variant was inherited from a parent with previously undiagnosed PVNH. The PVNH was frontally predominant and associated with perisylvian polymicrogyria. These results implicate MAP1B in PVNH. More broadly, our findings suggest that detrimental mutations likely arising in immediately preceding generations with incomplete penetrance may also be responsible for some apparently sporadic diseases.