A pollen selection system links self and interspecific incompatibility in the Brassicaceae.

Self-incompatibility and recurrent transitions to self-compatibility have shaped the extant mating systems underlying the nonrandom mating critical for speciation in angiosperms. Linkage between self-incompatibility and speciation is illustrated by the shared pollen rejection pathway between self-incompatibility and interspecific unilateral incompatibility (UI) in the Brassicaceae. However, the pollen discrimination system that activates this shared pathway for heterospecific pollen rejection remains unknown. Here we show that Stigma UI3.1, the genetically identified stigma determinant of UI in Arabidopsis lyrata × Arabidopsis arenosa crosses, encodes the S-locus-related glycoprotein 1 (SLR1). Heterologous expression of A. lyrata or Capsella grandiflora SLR1 confers on some Arabidopsis thaliana accessions the ability to discriminate against heterospecific pollen. Acquisition of this ability also requires a functional S-locus receptor kinase (SRK), whose ligand-induced dimerization activates the self-pollen rejection pathway in the stigma. SLR1 interacts with SRK and interferes with SRK homomer formation. We propose a pollen discrimination system based on competition between basal or ligand-induced SLR1-SRK and SRK-SRK complex formation. The resulting SRK homomer levels would be sensed by the common pollen rejection pathway, allowing discrimination among conspecific self- and cross-pollen as well as heterospecific pollen. Our results establish a mechanistic link at the pollen recognition phase between self-incompatibility and interspecific incompatibility.

Transposable element-initiated enhancer-like elements generate the subgenome-biased spike specificity of polyploid wheat.

Transposable elements (TEs) comprise ~85% of the common wheat genome, which are highly diverse among subgenomes, possibly contribute to polyploid plasticity, but the causality is only assumed. Here, by integrating data from gene expression cap analysis and epigenome profiling via hidden Markov model in common wheat, we detect a large proportion of enhancer-like elements (ELEs) derived from TEs producing nascent noncoding transcripts, namely ELE-RNAs, which are well indicative of the regulatory activity of ELEs. Quantifying ELE-RNA transcriptome across typical developmental stages reveals that TE-initiated ELE-RNAs are mainly from RLG_famc7.3 specifically expanded in subgenome A. Acquisition of spike-specific transcription factor binding likely confers spike-specific expression of RLG_famc7.3-initiated ELE-RNAs. Knockdown of RLG_famc7.3-initiated ELE-RNAs resulted in global downregulation of spike-specific genes and abnormal spike development. These findings link TE expansion to regulatory specificity and polyploid developmental plasticity, highlighting the functional impact of TE-driven regulatory innovation on polyploid evolution.

Phase separation of S-RNase promotes self-incompatibility in Petunia hybrida.

Self-incompatibility (SI) is an intraspecific reproductive barrier widely present in angiosperms. The SI system with the broadest occurrence in angiosperms is based on an S-RNase linked to a cluster of multiple S-locus F-box (SLF) genes found in the Solanaceae, Plantaginaceae, Rosaceae, and Rutaceae. Recent studies reveal that non-self S-RNase is degraded by the Skip Cullin F-box (SCF)SLF -mediated ubiquitin-proteasome system in a collaborative manner in Petunia, but how self-RNase functions largely remains mysterious. Here, we show that S-RNases form S-RNase condensates (SRCs) in the self-pollen tube cytoplasm through phase separation and the disruption of SRC formation breaks SI in self-incompatible Petunia hybrida. We further find that the pistil SI factors of a small asparagine-rich protein HT-B and thioredoxin h together with a reduced state of the pollen tube all promote the expansion of SRCs, which then sequester several actin-binding proteins, including the actin polymerization factor PhABRACL, the actin polymerization activity of which is reduced by S-RNase in vitro. Meanwhile, we find that S-RNase variants lacking condensation ability fail to recruit PhABRACL and are unable to induce actin foci formation required for pollen tube growth inhibition. Taken together, our results demonstrate that phase separation of S-RNase promotes SI response in P. hybrida, revealing a new mode of S-RNase action.

Determining a multimodal aging clock in a cohort of Chinese women.

BACKGROUND: Translating aging rejuvenation strategies into clinical practice has the potential to address the unmet needs of the global aging population. However, to successfully do so requires precise quantification of aging and its reversal in a way that encompasses the complexity and variation of aging. METHODS: Here, in a cohort of 113 healthy women, tiled in age from young to old, we identified a repertoire of known and previously unknown markers associated with age based on multimodal measurements, including transcripts, proteins, metabolites, microbes, and clinical laboratory values, based on which an integrative aging clock and a suite of customized aging clocks were developed. FINDINGS: A unified analysis of aging-associated traits defined four aging modalities with distinct biological functions (chronic inflammation, lipid metabolism, hormone regulation, and tissue fitness), and depicted waves of changes in distinct biological pathways peak around the third and fifth decades of life. We also demonstrated that the developed aging clocks could measure biological age and assess partial aging deceleration by hormone replacement therapy, a prevalent treatment designed to correct hormonal imbalances. CONCLUSIONS: We established aging metrics that capture systemic physiological dysregulation, a valuable framework for monitoring the aging process and informing clinical development of aging rejuvenation strategies. FUNDING: This work was supported by the National Natural Science Foundation of China (32121001), the National Key Research and Development Program of China (2022YFA1103700 and 2020YFA0804000), the National Natural Science Foundation of China (81502304), and the Quzhou Technology Projects (2022K46).

fastHaN: a fast and scalable program for constructing haplotype network for large-sample sequences.

Haplotype networks can be used to demonstrate the genealogical relationships of DNA sequences within species, and thus are widely applied in population genetics, molecular ecology, epidemiology and evolutionary studies. However, existing programs become computationally infeasible as the sample size increases. Here, we present fastHaN, an efficient and scalable program suitable for constructing haplotype networks for large samples. On a data set with the haplotype length of 30 kb, the Median Joining Network (MJN) algorithm implemented by fastHaN is 3000 times faster than PopART and 70 times faster than NETWORK in single-threaded mode. The implementation of the Templeton-Crandall-Sing (TCS) algorithm is 100 times faster than PopART and 5800 times faster than the TCS software. Moreover, fastHaN also enables multi-threaded mode with scalability. The source code is freely available on https://github.com/ChenHuaLab/fastHaN/. A web-based version is also available on https://ngdc.cncb.ac.cn/haplotype/.

Brassinosteroid coordinates cell layer interactions in plants via cell wall and tissue mechanics.

Growth coordination between cell layers is essential for development of most multicellular organisms. Coordination may be mediated by molecular signaling and/or mechanical connectivity between cells, but how genes modify mechanical interactions between layers is unknown. Here we show that genes driving brassinosteroid synthesis promote growth of internal tissue, at least in part, by reducing mechanical epidermal constraint. We identified a brassinosteroid-deficient dwarf mutant in the aquatic plant Utricularia gibba with twisted internal tissue, likely caused by mechanical constraint from a slow-growing epidermis. We tested this hypothesis by showing that a brassinosteroid mutant in Arabidopsis enhances epidermal crack formation, indicative of increased tissue stress. We propose that by remodeling cell walls, brassinosteroids reduce epidermal constraint, showing how genes can control growth coordination between layers by means of mechanics.

McAN: a novel computational algorithm and platform for constructing and visualizing haplotype networks.

Haplotype networks are graphs used to represent evolutionary relationships between a set of taxa and are characterized by intuitiveness in analyzing genealogical relationships of closely related genomes. We here propose a novel algorithm termed McAN that considers mutation spectrum history (mutations in ancestry haplotype should be contained in descendant haplotype), node size (corresponding to sample count for a given node) and sampling time when constructing haplotype network. We show that McAN is two orders of magnitude faster than state-of-the-art algorithms without losing accuracy, making it suitable for analysis of a large number of sequences. Based on our algorithm, we developed an online web server and offline tool for haplotype network construction, community lineage determination, and interactive network visualization. We demonstrate that McAN is highly suitable for analyzing and visualizing massive genomic data and is helpful to enhance the understanding of genome evolution. Availability: Source code is written in C/C++ and available at https://github.com/Theory-Lun/McAN and https://ngdc.cncb.ac.cn/biocode/tools/BT007301 under the MIT license. Web server is available at https://ngdc.cncb.ac.cn/bit/hapnet/. SARS-CoV-2 dataset are available at https://ngdc.cncb.ac.cn/ncov/. Contact: songshh@big.ac.cn (Song S), zhaowm@big.ac.cn (Zhao W), baoym@big.ac.cn (Bao Y), zhangzhang@big.ac.cn (Zhang Z), ybxue@big.ac.cn (Xue Y).

The Snapdragon Genomes Reveal the Evolutionary Dynamics of the S-Locus Supergene.

The genus Antirrhinum has been used as a model to study self-incompatibility extensively. The multi-allelic S-locus, carrying a pistil S-RNase and dozens of S-locus F-box (SLF) genes, underlies the genetic control of self-incompatibility (SI) in Antirrhinum hispanicum. However, there have been limited studies on the genomic organization of the S-locus supergene due to a lack of high-quality genomic data. Here, we present the chromosome-level reference and haplotype-resolved genome assemblies of a self-incompatible A. hispanicum line, AhS7S8. For the first time, 2 complete A. hispanicum S-haplotypes spanning ∼1.2 Mb and containing a total of 32 SLFs were reconstructed, whereas most of the SLFs derived from retroelement-mediated proximal or tandem duplication ∼122 Mya. Back then, the S-RNase gene and incipient SLFs came into linkage to form the pro-type of type-1 S-locus in the common ancestor of eudicots. Furthermore, we detected a pleiotropic cis-transcription factor (TF) associated with regulating the expression of SLFs, and two miRNAs may control the expression of this TF. Interspecific S-locus and intraspecific S-haplotype comparisons revealed the dynamic nature and polymorphism of the S-locus supergene mediated by continuous gene duplication, segmental translocation or loss, and TE-mediated transposition events. Our data provide an excellent resource for future research on the evolutionary studies of the S-RNase-based self-incompatibility system.

Genomic decoding of breeding history to guide breeding-by-design in rice.

Deciphering the intrinsic molecular logic of empirical crop breeding from a genomic perspective is a decisive prerequisite for breeding-by-design (BbD), but remains not well established. Here, we decoded the historical features of past rice breeding by phenotyping and haplotyping 546 accessions covering the majority of cultivars bred in the history of Northeast China (NEC). We revealed that three groups founded the genetic diversities in NEC rice with distinct evolution patterns and traced and verified the breeding footprints to known or genome-wide association study (GWAS)-detected quantitative trait loci (QTLs), or introgressions from indica sub-species with chronological changes in allele frequencies. Then we summarized a rice breeding trend/principle in NEC, and combined with the successful example in breeding and application of Zhongkefa5 to demonstrate the guiding value of our conclusion for BbD in practice. Our study provides a paradigm for decoding the breeding history of a specific crop to guide BbD, which may have implications in different crop breeding.

Brassinosteroid signals cooperate with katanin-mediated microtubule severing to control stamen filament elongation.

Proper stamen filament elongation is essential for pollination and plant reproduction. Plant hormones are extensively involved in every stage of stamen development; however, the cellular mechanisms by which phytohormone signals couple with microtubule dynamics to control filament elongation remain unclear. Here, we screened a series of Arabidopsis thaliana mutants showing different microtubule defects and revealed that only those unable to sever microtubules, lue1 and ktn80.1234, displayed differential floral organ elongation with less elongated stamen filaments. Prompted by short stamen filaments and severe decrease in KTN1 and KTN80s expression in qui-2 lacking five BZR1-family transcription factors (BFTFs), we investigated the crosstalk between microtubule severing and brassinosteroid (BR) signaling. The BFTFs transcriptionally activate katanin-encoding genes, and the microtubule-severing frequency was severely reduced in qui-2. Taken together, our findings reveal how BRs can regulate cytoskeletal dynamics to coordinate the proper development of reproductive organs.

Transposable elements orchestrate subgenome-convergent and -divergent transcription in common wheat.

The success of common wheat as a global staple crop was largely attributed to its genomic diversity and redundancy due to the merge of different genomes, giving rise to the major question how subgenome-divergent and -convergent transcription is mediated and harmonized in a single cell. Here, we create a catalog of genome-wide transcription factor-binding sites (TFBSs) to assemble a common wheat regulatory network on an unprecedented scale. A significant proportion of subgenome-divergent TFBSs are derived from differential expansions of particular transposable elements (TEs) in diploid progenitors, which contribute to subgenome-divergent transcription. Whereas subgenome-convergent transcription is associated with balanced TF binding at loci derived from TE expansions before diploid divergence. These TFBSs have retained in parallel during evolution of each diploid, despite extensive unbalanced turnover of the flanking TEs. Thus, the differential evolutionary selection of paleo- and neo-TEs contribute to subgenome-convergent and -divergent regulation in common wheat, highlighting the influence of TE repertory plasticity on transcriptional plasticity in polyploid.

A legume kinesin controls vacuole morphogenesis for rhizobia endosymbiosis.

Symbioses between legumes and rhizobia require establishment of the plant-derived symbiosome membrane, which surrounds the rhizobia and accommodates the symbionts by providing an interface for nutrient and signal exchange. The host cytoskeleton and endomembrane trafficking systems play central roles in the formation of a functional symbiotic interface for rhizobia endosymbiosis; however, the underlying mechanisms remain largely unknown. Here we demonstrate that the nodulation-specific kinesin-like calmodulin-binding protein (nKCBP), a plant-specific microtubule-based kinesin motor, controls central vacuole morphogenesis in symbiotic cells in Medicago truncatula. Phylogenetic analysis further indicated that nKCBP duplication occurs solely in legumes of the clade that form symbiosomes. Knockout of nKCBP results in central vacuole deficiency, defective symbiosomes and abolished nitrogen fixation. nKCBP decorates linear particles along microtubules, and crosslinks microtubules with the actin cytoskeleton, to control central vacuole formation by modulating vacuolar vesicle fusion in symbiotic cells. Together, our findings reveal that rhizobia co-opted nKCBP to achieve symbiotic interface formation by regulating cytoskeletal assembly and central vacuole morphogenesis during nodule development.

Ongoing Positive Selection Drives the Evolution of SARS-CoV-2 Genomes.

SARS-CoV-2 is a new RNA virus affecting humans and spreads extensively throughout the world since its first outbreak in December, 2019. Whether the transmissibility and pathogenicity of SARS-CoV-2 in humans after zoonotic transfer are actively evolving, and driven by adaptation to the new host and environments is still under debate. Understanding the evolutionary mechanism underlying epidemiological and pathological characteristics of COVID-19 is essential for predicting the epidemic trend, and providing guidance for disease control and treatments. Interrogating novel strategies for identifying natural selection using within-species polymorphisms and 3,674,076 SARS-CoV-2 genome sequences of 169 countries as of December 30, 2021, we demonstrate with population genetic evidence that during the course of SARS-CoV-2 pandemic in humans, 1) SARS-CoV-2 genomes are overall conserved under purifying selection, especially for the 14 genes related to viral RNA replication, transcription, and assembly; 2) ongoing positive selection is actively driving the evolution of 6 genes (e.g., S, ORF3a, and N) that play critical roles in molecular processes involving pathogen-host interactions, including viral invasion into and egress from host cells, and viral inhibition and evasion of host immune response, possibly leading to high transmissibility and mild symptom in SARS-CoV-2 evolution. According to an established haplotype phylogenetic relationship of 138 viral clusters, a spatial and temporal landscape of 556 critical mutations is constructed based on their divergence among viral haplotype clusters or repeatedly increase in frequency within at least 2 clusters, of which multiple mutations potentially conferring alterations in viral transmissibility, pathogenicity, and virulence of SARS-CoV-2 are highlighted, warranting attention.

The twin-beginnings of COVID-19 in Asia and Europe-one prevails quickly.

In the spread of SARS-CoV-2, there have been multiple waves of replacement between strains, each of which having a distinct set of mutations. The first wave is a group of four mutations (C241T, C3037T, C14408T and A23403G [this being the amino acid change D614G]; all designated 0 to 1 below). This DG (D614G) group, fixed at the start of the pandemic, is the foundation of all subsequent waves of strains. Curiously, the DG group is absent in early Asian samples but present (and likely common) in Europe from the beginning. European data show that the high fitness of DG1111 requires the synergistic effect of all four mutations. However, the European strains would have had no time to evolve the four DG mutations (0 to 1), had they come directly from the early Asian DG0000 strain. Very likely, the European DG1111 strain had acquired the highly adaptive DG mutations in pre-pandemic Europe and had been spreading in parallel with the Asian strains. Two recent reports further support this twin-beginning interpretation. There was a period of two-way spread between Asia and Europe but, by May 2020, the European strains had supplanted the Asian strains globally. This large-scale replacement of one set of mutations for another has since been replayed many times as COVID-19 progresses.

Origin, loss, and regain of self-incompatibility in angiosperms.

The self-incompatibility (SI) system with the broadest taxonomic distribution in angiosperms is based on multiple S-locus F-box genes (SLFs) tightly linked to an S-RNase termed type-1. Multiple SLFs collaborate to detoxify nonself S-RNases while being unable to detoxify self S-RNases. However, it is unclear how such a system evolved, because in an ancestral system with a single SLF, many nonself S-RNases would not be detoxified, giving low cross-fertilization rates. In addition, how the system has been maintained in the face of whole-genome duplications (WGDs) or lost in other lineages remains unclear. Here we show that SLFs from a broad range of species can detoxify S-RNases from Petunia with a high detoxification probability, suggestive of an ancestral feature enabling cross-fertilization and subsequently modified as additional SLFs evolved. We further show, based on its genomic signatures, that type-1 was likely maintained in many lineages, despite WGD, through deletion of duplicate S-loci. In other lineages, SI was lost either through S-locus deletions or by retaining duplications. Two deletion lineages regained SI through type-2 (Brassicaceae) or type-4 (Primulaceae), and one duplication lineage through type-3 (Papaveraceae) mechanisms. Thus, our results reveal a highly dynamic process behind the origin, maintenance, loss, and regain of SI.

Metagenomic evidence for the co-existence of SARS and H1N1 in patients from 2007-2012 flu seasons in France.

By re-analzying public metagenomic data from 101 patients infected with influenza A virus during the 2007-2012 H1N1 flu seasons in France, we identified 22 samples with SARS-CoV sequences. In 3 of them, the SARS genome sequences could be fully assembled out of each. These sequences are highly similar (99.99% and 99.7%) to the artificially constructed recombinant SARS-CoV (SARSr-CoV) strains generated by the J. Craig Venter Institute in the USA. Moreover, samples from different flu seasons have different SARS-CoV strains, and the divergence between these strains cannot be explained by natural evolution. Our study also shows that retrospective studies using public metagenomic data from past major epidemic outbreaks serve as a genomic strategy for researching the origins or spread of infectious diseases.