Meiotic DNA double-strand break (DSB) formation is pivotal for oocyte development, regulating both ovarian reserve and oocyte developmental potential. Mutations in DSB formation genes have been associated with premature ovarian insufficiency (POI) and adverse pregnancy outcomes in women. Whole exome sequencing in 1530 POI patients across two Chinese cohorts identifies loss-of-function variants in the DSB formation gene, MEI4, enriched in POI. These MEI4 variants impair DSB formation in vitro and reveal a previously unrecognized function of the MEI4 C-terminus in stabilizing the MEI4–REC114 subcomplex on the chromosome axes. Additionally, Mei4^{Arg356*/Arg356*} mice display severe defects in DSB formation, leading to massive apoptosis in oocytes triggered by the HORMAD1-dependent synapsis checkpoint in late prophase I. The few mutant oocytes surviving past the checkpoint exhibit low developmental potential, characterized by complete early embryonic arrest due to aneuploidy. Notably, heterozygous Mei4^+/Arg356* mice show intermediate follicle depletion and embryonic development arrest consistent with the phenotype of heterozygous POI and preimplantation embryonic arrest, suggesting a haploinsufficiency effect. This study defines the impacts of MEI4 mutation on oocyte quantity and quality, which can guide genetic diagnosis and intervention in patients with POI and early embryonic arrest, especially those with mutations in meiotic DSB formation genes.

Marchantia polymorpha, a model liverwort, provides a valuable system for investigating the evolution of plant sexual reproduction. To explore the cellular landscape of its reproductive structures, we generate a single-nucleus transcriptomic atlas of the antheridiophore, archegoniophore, and sporophyte. Using single-nucleus RNA sequencing (snRNA-seq), we capture over 30,000 high-quality nuclei and identify distinct cell populations. In the male organ, we characterize stages of spermatogenesis from early antheridium cells to mature sperm, revealing dynamic transcriptional programs including cell cycle regulation, chromatin remodelling, and calcium signalling. In the female organ, we define cell types including archegonial layers and secondary central cells. Sporophyte clusters are annotated as spores, elaters, capsule wall, foot, and seta cells, with transcriptional signatures related to structural support, stress response, and reproductive functions. Cross-species analysis indicates that capsule wall cells in liverworts are similar to tapetum cells. Notably, foot cells exhibit high expression of genes involved in sporopollenin biosynthesis and signaling pathways, serving as a central hub that mediates communication between the maternal gametophyte and the developing sporophyte. This study provides a comprehensive cellular and molecular map of M. polymorpha reproductive organs and sporophyte, establishing a framework for investigating the development and evolution of sexual reproduction in early land plants.

Cilia are vital subcellular organelles whose assembly is regulated by master transcription factors, such as Foxj1 and Rfx. However, the mechanisms of epigenetic regulation over cilia stability remain largely unclear. Here, we investigate the epigenetic control by manipulating chromatin remodeling genes in zebrafish. We demonstrate that the depletion of multiple components of the switch/sucrose non-fermentable (SWI/SNF) chromatin remodeling complex induces ciliopathy-like phenotypes in zebrafish embryos. Specifically, the loss of Actl6a, an essential component of the SWI/SNF complex, leads to cilia disassembly and cystic kidney defects, without affecting cilia motility. Our multi-omics analyses (RNA-seq, ATAC-seq, and FitCUT&RUN) consistently reveal that in Actl6a-depleted pronephros or embryos, a critical set of cilia genes, including the master regulators foxj1a and rfx2, exhibit concordant downregulation across the transcriptional level, chromatin accessibility, and SWI/SNF binding. Consistently, the depletion of foxj1a or rfx2 causes cilia assembly defects and cystic kidney formation in zebrafish. Furthermore, overexpression of either foxj1a or rfx2 mRNA substantially rescues the cystic kidney and cilia disassembly defects observed in actl6a^-/- mutant embryos. Collectively, these findings reveal that the SWI/SNF complex maintains cilia stability and kidney homeostasis by directly modulating the expression of the key ciliogenesis transcription factors foxj1a and rfx2.

Patients affected by monogenic diseases impose a substantial burden on both themselves and their families. The primary preventive measure, i.e., invasive prenatal diagnosis, carries a risk of miscarriage and cannot be performed early in pregnancy. Hence, there is a need for non-invasive prenatal testing (NIPT) for monogenic diseases. By utilizing enriched cell-free fetal DNA (cffDNA) from maternal plasma, we refine the NIPT method, which combines targeted region capture technology, haplotyping, and analysis of informative site frequency. We apply this method to 93 clinical families at genetic risk for thalassemia, encompassing various genetic variant types, to establish a workflow and evaluate its efficiency. Our approach requires only 3 ng of DNA input to generate 0.1 GB of informative target genomic data and leverages a minimum of 3% cffDNA. This method has a 98.16% success rate and 100% concordance with conventional invasive methods. Furthermore, we demonstrate the ability to analyze fetal genotypes as early as eight weeks of gestation. This study establishes an optimized NIPT method for the early detection of various thalassemia disorders during pregnancy. This technique demonstrates high accuracy and potential for clinical application in prenatal diagnosis.

The genetic basis of early-stage salt tolerance in alfalfa (Medicago sativa L.), a key factor limiting its productivity, remains poorly understand. To dissect this complex trait, we integrate genome-wide association study (GWAS) and transcriptomics (RNA-seq) from 176 accessions within a machine learning based genomic prediction framework. Analysis reveals weak genetic correlations among four salt-tolerance traits and a gradual decline in performance under increasing salt stress. GWAS identify 60 significant associated SNPs, with the highest number detected under 100 mM salt stress. Salt tolerance exhibits an additive effect from favorable haplotypes, which are most abundant in Chinese accessions. GWAS-associated genes are related to key regulators of hormone signaling and osmotic adjustment, while transcriptome analysis indicates a global repression of stress-responsive transcription factors. Integrating these multi-omics datasets allow us to identify 14 candidate genes, including MsHSD1 (seed dormancy) and MsMTATP6 (energy metabolism). Crucially, incorporating these markers into genomic prediction models improve cross-population predictive accuracy to an average of 54.4%. This study provides insights into the genetic architecture of salt tolerance in alfalfa and offers valuable markers to facilitate molecular breeding.

Hereditary cardiomyopathies and arrhythmias are major contributors to cardiovascular morbidity and mortality. The advent of next-generation sequencing (NGS) has made genetic testing more accessible, which is crucial for precise diagnosis and targeted therapeutic strategies. The aim of this study is to explore the landscape of genetic variants, the relationship between specific variants and clinical phenotypes, and the impact on clinical decision-making in China. A total of 1,536 probands (median age, 37 years; 1,025 males [66.7%]) with suspected hereditary cardiomyopathy or arrhythmia (covering 15 clinical phenotypes) are recruited from 146 hospitals across 30 provinces and cities in China. Positive results are confirmed in 390 of 1536 probands, leading to a diagnostic yield of 25.4%. Forty-two and three tenths percent (n = 169) of family members carry the same variants as positive probands. Hypertrophic cardiomyopathy (HCM) and dilated cardiomyopathy (DCM) are the predominant phenotypes, with MYBPC3 variants having the highest frequency in HCM and TTN variants in DCM. In 76.9% of the positive probands, the identified variants are helpful in clinical management, family screening and fertility. This large-scale study provides significant insights into the genetic landscape of hereditary cardiomyopathies and arrhythmias in China.

Phytochrome-interacting factors (PIFs) have been established as negative regulators of vascular patterning and xylem differentiation in the herbaceous plant Arabidopsis thaliana, however, PIFs’ regulatory role in secondary growth in woody species remains unclear. Here, we examined the expression patterns and involvement of PtoPIF3.1 and PtoPIF3.2 during stem growth and secondary xylem development in Populus tomentosa. Overexpression of either PtoPIF3.1 or PtoPIF3.2 significantly enhances both longitudinal stem growth and radial wood development. Conversely, CRISPR-generated Ptopif3.1 and Ptopif3.2 mutants exhibit reciprocal phenotypic defects. Exogenous auxin application partially restores the phenotypes of Ptopif3.1 and Ptopif3.2 mutants, and the auxin biosynthesis-deficient mutant Ptoyuc8 exhibits developmental abnormalities similar to those observed in Ptopif3 mutants. Further analysis revealed that PtoPIF3s directly bind to and activate expression of PtoYUC8 and cell expansion-related genes PtoEXPA1.1/1.2, while modulating cambial division and secondary xylem development marker genes (PtoWOX4, PtoANT, PtoCYCD3s, and PtoHB7/8) through auxin-mediated signaling. Together, our findings establish PtoPIF3.1/3.2 as key regulators that coordinate stem elongation and secondary growth in Populus, highlighting the functional divergence of PIF homologs between herbaceous and woody species.

Penetrance is a crucial indicator for accurately assessing disease risk and plays a vital role in disease research, gene therapy, and genetic counseling. However, with penetrance data dispersed across various sources, efficiently accessing and consolidating this information becomes a challenge. A comprehensive platform that integrates penetrance is urgently needed. Here, we present PenCards, a global, community-contributed public archive of variant penetrance, by first collecting penetrance data from all published literature and then using large international cohorts to specifically calculate the penetrance of autism-related variants. PenCards contains a total of 244,531 variants—including 239,244 single nucleotide variants, 4,994 insertions and deletions, and 293 copy number variants, covering approximately 300 phenotypes. We also provide a submission portal for the dynamic updating of penetrance. Additionally, to help users efficiently access genetic information, we comprehensively integrate over 150 variant- and gene-level resources. In summary, PenCards is a powerful platform designed to advance genetic research and diagnostics. PenCards is publicly available at https://genemed.tech/pencards/.

Cilia are microtubule-based organelles projecting from the cell surface with important sensory and motility functions. Ciliary defects are associated with diverse diseases collectively known as ciliopathies. However, the molecular mechanisms that govern ciliogenesis remain not fully understood. Herein, we demonstrate that ubiquitin-specific protease 21 (USP21) is indispensable for cilium formation through its deubiquitinating activity. Usp21 knockout mice exhibit ciliary defects in multiple organs, such as the kidney, liver, and trachea. Our data also reveal a constant localization of USP21 at the centrosome and basal body during ciliogenesis. Mechanistically, USP21 interacts with dihydropyrimidinase-like 2 (DYPSL2) at the centrosome and removes lysine 48-linked ubiquitination from DYPSL2. Loss of USP21 leads to the proteasomal degradation of DPYSL2 and causes a significant reduction in its centrosome abundance, ultimately resulting in ciliary defects. These findings thus identify a critical role for the USP21–DPYSL2 axis in ciliogenesis and have important implications for health and disease.

Leaf angle is a pivotal agronomic trait that significantly influences crop architecture and yield. Plant hormones, such as auxin, play a critical role in regulating leaf angle; however, the underlying molecular mechanisms remain to be fully elucidated. Here, we reveal that the auxin response factor gene, OsARF12, which is highly expressed in the leaf lamina joint, negatively regulates rice (Oryza sativa) leaf angle via affecting shoot gravitropism. Overexpression lines of OsARF12 exhibit more erect leaf angles, while the osarf12 mutants display enlarged leaf angles compared to the wild type. Further studies demonstrate that OsARF12 directly activates the expression of Loose Plant Architecture1 (LPA1) and LAZY1 by binding to their promoters. The osarf12 mutant presents impaired shoot gravitropism, a phenotype consistent with that of the lpa1 or lazy1 mutant. Collectively, we elucidate the biological functions of OsARF12, which modulates leaf angle through its impact on shoot gravitropism by regulating the expression levels of LPA1 and LAZY1. This study provides insight into the role of auxin in determining rice leaf angle, potentially holding profound effects for the optimization of crop architecture.

Modern agriculture relies heavily on resource-intensive and environmentally harmful inputs, while the increasing global population and decreasing arable land demand new strategies to improve sustainable productivity of cereal crops, particularly reducing inorganic nitrogen fertilizer use while simultaneously increasing photosynthesis and grain yield in rice. To improve rice productivity, it is essential to improve photosynthetic nitrogen assimilation and optimize the translocation of carbon and nitrogen products from source to sink tissues. In this review, we first summarize recent advances in the genetic basis for improving grain yield by enhancing photosynthetic carbon and nitrogen assimilation. We then discuss progress in modulating the source-sink relationships to achieve higher yield and improved harvest index. Finally, we explore the necessary optimizations for adapting rice to high-density planting. These advancements are driving the development of sustainable green revolution varieties through the rational design of multi-gene pyramids and artificial intelligence (AI)-driven protein engineering.

The frequency of aneuploid gamete formation increases with maternal age, yet the effects of genetic variants on meiotic chromosome segregation accuracy during aging remain poorly understood. Using the multicellular organism Caenorhabditis elegans, we investigate the impact of mutations in the conserved cohesin complex on age-associated meiotic errors. Point mutations in the head domain of the cohesin component SMC-1, which alter local hydrophobicity, cause meiotic defects that vary with age. A severe mutation causes incomplete synapsis and defective crossover formation, and a minor one causes age-related diakinesis bivalent abnormalities. Notably, while the mild mutation causes defects only in aged worms, worms with the severe mutation exhibit significantly alleviated phenotypes with age. Genetic and cytological analyses suggest that this alleviation results from a slowed meiotic progression during early prophase, which restores impaired cohesin loading. These findings reveal that cohesin variants, meiotic progression speed during early prophase, and the overall duration of meiosis collectively shape the accuracy of meiotic chromosome segregation.

Existing quantitative trait locus (QTL) mapping had low efficiency in identifying small-effect and closely linked QTL-by-environment interactions (QEIs) in RILs, especially in the era of global climate change. To address this challenge, here we integrate the compressed variance component mixed model with our GCIM to propose 3vGCIM for identifying QEIs in RILs, and extend 3vGCIM-random to 3vGCIM-fixed. 3vGCIM integrates genome-wide scanning with machine learning, significantly improving power. In the mixed full model, we consider all possible effects and control for all possible polygenic backgrounds. In simulation studies, 3vGCIM exhibits higher power (∼92.00%), higher accuracy of the estimates for QTL position (∼1.900 cM²) and effect (∼0.050), and lower false positive rate (∼0.48‰) and false negative rate (<8.10%) in three environments of 300 RILs each than ICIM (47.57%; 3.607 cM², 0.583; 2.81‰; 52.43%) and MCIM (60.30%; 5.279 cM², 0.274; 2.17‰; 39.70%). In the real data analysis of yield-related in 240 RILs, 3vGCIM mines more known genes (57∼60) and known gene-by-environment interactions (GEIs) (14∼19) and candidate GEIs (21∼23) than ICIM (27, 2, and 7), and MCIM (21, 1, and 3), especially in small-effect and linked QTLs and QEIs. This makes 3vGCIM a powerful and sensitive tool for QTL mapping and molecular QTL mapping.

Common wheat (Triticum aestivum L.) is a staple hexaploid crop with numerous wild relatives. However, domestication and modern breeding have significantly narrowed its genetic diversity, diminishing its capacity to adapt to climate change. Wild relatives of wheat serve as a vital reservoir of genetic diversity, offering traits that enhance its resistance to various biotic and abiotic stresses. Over recent decades, remarkable progress has been made in utilizing superior genes from wild relatives to bolster wheat's defenses against diseases and pests, though the exploration of genes conferring abiotic stress tolerance has lagged behind. In this review, we summarize key advancements in the utilization of wild relatives for wheat enhancement over the past century, emphasizing both theoretical and technological innovations. Furthermore, we evaluate the potential contributions of wild relatives to address production challenges posed by climate changes. We also explore strategies for isolating superior genes and developing pre-breeding germplasm to support the future development of climate-resilient wheat varieties.

Down syndrome (DS) is caused by an extra copy of chromosome 21 (Hsa21). Children with DS have an increased frequency of respiratory tract infections, impaired alveolar and vascular development, and pulmonary hypertension. How trisomy 21 causes lung diseases remains poorly understood. In this study, we use the Dp16 mouse model, which contains a segmental chromosomal duplication of the entire Hsa21 syntenic region on mouse chromosome 16, to explore the gene dosage effects on DS-related lung diseases. The Dp16 mice present impaired alveolar development and inflammatory-like pathological changes. Single-cell transcriptome (scRNA-seq) analysis highlights increased APP-related interactions among male Dp16 lung cells. Specifically, altered antigen processing and presentation with increased MHC-II signaling are found in Dp16 immune cells. Reduced angiogenesis and altered inflammatory responses of Dp16 endothelial cells are also suggested. Moreover, scRNA-seq indicates hyperplasia of Dp16 vascular smooth muscle cells, which is validated by tissue immunofluorescence assessment. Transthoracic echocardiography further shows the existence of pulmonary hypertension in young Dp16 mice. Independent scRNA-seq analysis of the female lung cells recapitulates the majority of key findings identified in male mice, confirming the reproducibility of the results. Collectively, our results provide important clues for the further development of therapeutic approaches for DS-related lung diseases.

Amyotrophic lateral sclerosis (ALS), one of the most prevalent neurodegenerative disorders, is pathologically characterized by the progressive degeneration of both upper and lower motor neurons, leading to muscle weakness, paralysis, and death within 2-4 years post-diagnosis. ALS is categorized into familial ALS (FALS) and sporadic ALS, with FALS accounting for approximately 10% of ALS cases. As a genetically heterogeneous disease, ALS exhibits diverse inheritance patterns, including autosomal dominant, autosomal recessive, and X-linked transmission, and genetic factors play pivotal roles in disease pathogenesis. To date, at least 34 disease-causing loci and 32 genes for ALS have been identified. The investigations of mutant protein products and the establishment of animal models have unraveled potential pathogenic pathways, offering insights into the mechanisms of neurodegeneration in ALS. This review focuses on ALS clinical characteristics, neuropathological features, causative loci/genes, genetic susceptibility factors, animal models, and pathogenic mechanisms, with particular attention to recent advances in genetic findings and pathogenic pathways of ALS. Elucidation of the genetic basis of ALS could provide the scientific foundation for personalized treatments to address this recalcitrant disease.

Body size control is fundamental to development and requires proper energy engagement. One of the key energy sensing factors is AMP-activated protein kinase (AMPK), which regulates glucose uptake to ensure ATP production and nutrition supply during development. Here, we identify that the mutation of xgr, a gene encoding an ATPase, results in a reduced body size in Drosophila. Xgr is primarily expressed in the epithelial cells of the Malpighian tubules and the midguts. Loss of xgr leads to the inactivation of the AMPK signaling due to an increased ATP level. Glucose reabsorption in the Malpighian tubules is significantly reduced, as the Glut1 translocation to the plasma membrane is significantly disrupted in the absence of Xgr function. Our results suggest that Xgr function in the Malpighian tubules is essential to systemic glucose supply and energy homeostasis at the organismal level, thereby impacting body size. Our findings provide a mechanistic connection between energy homeostasis and animal size control during development.

Anther is a key male reproductive organ that is essential for the plant life cycle, from the sporophyte to the gametophyte generation. To explore the isoform-level transcriptional landscape of developing anthers in maize (Zea mays L.), we analyzed Iso-Seq data from anthers collected at 10 developmental stages, together with strand-specific RNA-seq, CAGE-seq, and PAS-seq data. Of the 152,026 high-confidence full-length isoforms identified, 68.8% have not been described; these include 22,365 isoforms that originate from previously unannotated loci and 82,167 novel isoforms that originate from annotated protein-coding genes. Using our newly developed strategy to detect dynamic expression patterns of isoforms, we identified 13,899 differentially variable regions (DVRs); surprisingly, 1275 genes contain more than two DVRs, revealing highly efficient utilization of limited genic regions. We identified 7876 long non-coding RNAs (lncRNAs) from 4098 loci, most of which were preferentially expressed during cell differentiation and meiosis. We also detected 371 long-range interactions involving intergenic lncRNAs (lincRNAs); interestingly, 243 were lincRNA–gene ones, and the interacting genes were highly expressed in anthers, suggesting that many potential lncRNA regulators of key genes are required for anther development. This study provides valuable resources and fundamental information for studying the essential transcripts of key genes during anther development.

Potassium channels regulate diverse biological processes, ranging from cell proliferation to immune responses. However, the functions of potassium homeostasis and its regulatory mechanisms in adult stem cells and tumors remain poorly characterized. Here, we identify Sandman, a two-pore-domain potassium channel in Drosophila, as an essential regulator for the proliferation of intestinal stem cells and malignant tumors, while dispensable for the normal development processes. Mechanistically, loss of sandman elevates intracellular K⁺ concentration, leading to growth inhibition. This phenotype is rescued by pharmacological reduction of intracellular K⁺ levels using the K⁺ ionophore. Conversely, overexpression of sandman triggers stem cell death in most regions of the midgut, inhibits tumor growth, and induces a Notch loss-of-function phenotype in the posterior midgut. These effects are mediated predominantly via the induction of endoplasmic reticulum (ER) stress, as demonstrated by the complete rescue of phenotypes through the co-expression of Ire1 or Xbp1s. Additionally, human homologs of Sandman demonstrated similar ER stress-inducing capabilities, suggesting an evolutionarily conserved relationship between this channel and ER stress. Together, our findings identify Sandman as a shared regulatory node that governs Drosophila adult stem cell dynamics and tumorigenesis through bioelectric homeostasis, and reveal a link between the two-pore potassium channel and ER stress signaling.

In vitro maturation (IVM) of human oocytes offers cost efficiency and minimal invasiveness, serving as a valuable supplementary tool in assisted reproduction for fertility preservation, ovarian hyperstimulation syndrome prevention, and other reproductive strategies. Despite its availability for three decades, the clinical use of IVM remains limited due to efficacy and safety concerns. This study examines the DNA methylation profile of IVM oocytes collected during laparoscopic/hysteroscopic surgeries compared to in vivo matured oocytes via reduced representation bisulfite sequencing. Results indicate IVM oocytes exhibit a higher global methylation level. Differentially methylated regions (DMR) analysis reveals that the in vitro group displays more hypermethylated and fewer hypomethylated DMRs compared to the in vivo group. Additionally, the in vitro group exhibits a higher level of non-CpG methylation than the in vivo group. However, no significant correlation between methylation levels and transcriptional activity in these oocytes is found, especially for those specific imprinted genes or genes related to embryonic development. These findings shed light on the epigenetic landscape of IVM oocytes, contributing to the ongoing assessment of their clinical feasibility and safety in assisted reproduction.

While methodology for determining the mode of evolution in coding sequences has been well established, evaluation of adaptation events in emerging types of phenotype data needs further development. Here we propose an analysis framework (expression variance decomposition, EVaDe) for comparative single-cell expression data based on phenotypic evolution theory. After decomposing the gene expression variance into separate components, we use two strategies to identify genes exhibiting large between-taxon expression divergence and small within-cell-type expression noise in certain cell types, attributing this pattern to putative adaptive evolution. In a dataset of primate prefrontal cortex, we find that such human-specific key genes enrich with neurodevelopment-related functions, while most other genes exhibit neutral evolution patterns. Specific neuron types are found to harbor more of these key genes than other cell types, thus likely to have experienced more extensive adaptation. Reassuringly, at molecular sequence level, the key genes are significantly associated with the rapidly evolving conserved non-coding elements. An additional case analysis comparing the naked mole-rat (NMR) with the mouse suggests that innate-immunity-related genes and cell types have undergone putative expression adaptation in NMR. Overall, the EVaDe framework may effectively probe adaptive evolution mode in single-cell expression data.

Understanding how species diverge and adapt is fundamental to unraveling biodiversity. While environmental impacts on species evolution are well-documented, the roles of intrinsic life-history traits remain underexplored. The Qinghai-Tibet Plateau, with its harsh conditions and unique biodiversity, offers a natural laboratory for such investigations. Here, we examined two sympatric small mammals—the solitary, low-dispersal plateau zokor (Eospalax baileyi) and the social, high-dispersal plateau pika (Ochotona curzoniae)—to elucidate how life-history traits shape population structures and adaptive strategies. Through whole-genome sequencing and cardiac-blood phenotype analyses, we reveal striking differences in their evolutionary trajectories. Despite enduring similar environmental pressures, plateau zokor populations exhibit pronounced genetic subdivisions, high inbreeding, and distinct local adaptations. In contrast, plateau pika populations display genetic panmixia, widespread diversity, and adaptive uniformity. Demographic inference highlights plateau zokors experienced severe population bottlenecks and restricted gene flow during glacial periods, underscoring the impact of dispersal capacity on evolutionary outcomes. Our findings demonstrate that intrinsic biological traits, particularly dispersal ability, fundamentally influence genetic architecture, population connectivity, and local adaptation. This study not only provides empirical evidence of how life-history traits shape evolutionary dynamics but also offers a framework for integrating intrinsic and extrinsic factors in understanding biodiversity formation.