Benefits and limitations of a new genome‐based PCR‐RFLP genotyping assay (GB‐RFLP): A SNP‐based detection method for identification of species in extremely young adaptive radiations

High‐throughput DNA sequencing technologies make it possible now to sequence entire genomes relatively easily. Complete genomic information obtained by whole‐genome resequencing (WGS) can aid in identifying and delineating species even if they are extremely young, cryptic, or morphologically difficult to discern and closely related. Yet, for taxonomic or conservation biology purposes, WGS can remain cost‐prohibitive, too time‐consuming, and often constitute a “data overkill.” Rapid and reliable identification of species (and populations) that is also cost‐effective is made possible by species‐specific markers that can be discovered by WGS. Based on WGS data, we designed a PCR restriction fragment length polymorphism (PCR‐RFLP) assay for 19 Neotropical Midas cichlid populations (Amphilophus cf. citrinellus), that includes all 13 described species of this species complex. Our work illustrates that identification of species and populations (i.e., fish from different lakes) can be greatly improved by designing genetic markers using available “high resolution” genomic information. Yet, our work also shows that even in the best‐case scenario, when whole‐genome resequencing information is available, unequivocal assignments remain challenging when species or populations diverged very recently, or gene flow persists. In summary, we provide a comprehensive workflow on how to design RFPL markers based on genome resequencing data, how to test and evaluate their reliability, and discuss the benefits and pitfalls of our approach.     

Drosophila Interspecific Hybrids Phenocopy piRNA-Pathway Mutants

The Piwi-interacting RNA (piRNA) pathway defends the germline of animals from the deleterious activity of selfish transposable elements (TEs) through small-RNA mediated silencing. Adaptation to novel invasive TEs is proposed to occur by incorporating their sequences into the piRNA pool that females produce and deposit into their eggs, which then propagates immunity against specific TEs to future generations. In support of this model, the F1 offspring of crosses between strains of the same Drosophila species sometimes suffer from germline derepression of paternally inherited TE families, caused by a failure of the maternal strain to produce the piRNAs necessary for their regulation. However, many protein components of the Drosophila piRNA pathway exhibit signatures of positive selection, suggesting that they also contribute to the evolution of host genome defense. Here we investigate piRNA pathway function and TE regulation in the F1 hybrids of interspecific crosses between D. melanogaster and D. simulans and compare them with intraspecific control crosses of D. melanogaster. We confirm previous reports showing that intraspecific crosses are characterized by derepression of paternally inherited TE families that are rare or absent from the maternal genome and piRNA pool, consistent with the role of maternally deposited piRNAs in shaping TE silencing. In contrast to the intraspecific cross, we discover that interspecific hybrids are characterized by widespread derepression of both maternally and paternally inherited TE families. Furthermore, the pattern of derepression of TE families in interspecific hybrids cannot be attributed to their paucity or absence from the piRNA pool of the maternal species. Rather, we demonstrate that interspecific hybrids closely resemble piRNA effector-protein mutants in both TE misregulation and aberrant piRNA production. We suggest that TE derepression in interspecific hybrids largely reflects adaptive divergence of piRNA pathway genes rather than species-specific differences in TE-derived piRNAs.     

Structural and mechanistic insight into stem‐loop RNA processing by yeast Pichia stipitis Dicer

Dicer is a member of the ribonuclease III enzyme family and processes double‐stranded RNA into small functional RNAs. The variation in the domain architecture of Dicer among different species whilst preserving its biological dicing function is intriguing. Here, we describe the structure and function of a novel catalytically active RNase III protein, a non‐canonical Dicer (PsDCR1), found in budding yeast Pichia stipitis. The structure of the catalytically active region (the catalytic RNase III domain and double‐stranded RNA‐binding domain 1 [dsRBD1]) of DCR1 showed that RNaseIII domain is structurally similar to yeast RNase III (Rnt1p) but uniquely presents dsRBD1 in a diagonal orientation, forming a catalytic core made of homodimer and large RNA‐binding surface. The second dsRNA binding domain at C‐terminus, which is absent in Rnt1, enhances the RNA cleavage activity. Although the cleavage pattern of PsDCR1 anchors an apical loop similar to Rnt1, the cleavage activity depended on the sequence motif at the lower stem, not the apical loop, of hairpin RNA. Through RNA sequencing and RNA mutations, we showed that RNA cleavage by PsDCR1 is determined by the stem‐loop structure of the RNA substrate, suggesting the possibility that stem‐loop RNA‐guided gene silencing pathway exists in budding yeast.     

DEAD‐box polypeptide 43 facilitates piRNA amplification by actively liberating RNA from Ago3‐piRISC

The piRNA amplification pathway in Bombyx is operated by Ago3 and Siwi in their piRISC form. The DEAD‐box protein, Vasa, facilitates Ago3‐piRISC production by liberating cleaved RNAs from Siwi‐piRISC in an ATP hydrolysis‐dependent manner. However, the Vasa‐like factor facilitating Siwi‐piRISC production along this pathway remains unknown. Here, we identify DEAD‐box polypeptide 43 (DDX43) as the Vasa‐like protein functioning in Siwi‐piRISC production. DDX43 belongs to the helicase superfamily II along with Vasa, and it contains a similar helicase core. DDX43 also contains a K‐homology (KH) domain, a prevalent RNA‐binding domain, within its N‐terminal region. Biochemical analyses show that the helicase core is responsible for Ago3‐piRISC interaction and ATP hydrolysis, while the KH domain enhances the ATPase activity of the helicase core. This enhancement is independent of the RNA‐binding activity of the KH domain. For maximal DDX43 RNA‐binding activity, both the KH domain and helicase core are required. This study not only provides new insight into the piRNA amplification mechanism but also reveals unique collaborations between the two domains supporting DDX43 function within the pathway.     

Innate,translation‐dependent silencing of an invasive transposon in Arabidopsis

Co‐evolution between hosts’ and parasites’ genomes shapes diverse pathways of acquired immunity based on silencing small (s)RNAs. In plants, sRNAs cause heterochromatinization, sequence degeneration, and, ultimately, loss of autonomy of most transposable elements (TEs). Recognition of newly invasive plant TEs, by contrast, involves an innate antiviral‐like silencing response. To investigate this response’s activation, we studied the single‐copy element EVADÉ (EVD), one of few representatives of the large Ty1/Copia family able to proliferate in Arabidopsis when epigenetically reactivated. In Ty1/Copia elements, a short subgenomic mRNA (shGAG) provides the necessary excess of structural GAG protein over the catalytic components encoded by the full‐length genomic flGAG‐POL. We show here that the predominant cytosolic distribution of shGAG strongly favors its translation over mostly nuclear flGAG‐POL. During this process, an unusually intense ribosomal stalling event coincides with mRNA breakage yielding unconventional 5’OH RNA fragments that evade RNA quality control. The starting point of sRNA production by RNA‐DEPENDENT‐RNA‐POLYMERASE‐6 (RDR6), exclusively on shGAG, occurs precisely at this breakage point. This hitherto‐unrecognized “translation‐dependent silencing” (TdS) is independent of codon usage or GC content and is not observed on TE remnants populating the Arabidopsis genome, consistent with their poor association, if any, with polysomes. We propose that TdS forms a primal defense against EVD de novo invasions that underlies its associated sRNA pattern.     

Behavioral “bycatch” from camera trap surveys yields insights on prey responses to human‐mediated predation risk

Human disturbance directly affects animal populations and communities, but indirect effects of disturbance on species behaviors are less well understood. For instance, disturbance may alter predator activity and cause knock‐on effects to predator‐sensitive foraging in prey. Camera traps provide an emerging opportunity to investigate such disturbance‐mediated impacts to animal behaviors across multiple scales. We used camera trap data to test predictions about predator‐sensitive behavior in three ungulate species (caribou Rangifer tarandus; white‐tailed deer, Odocoileus virginianus; moose, Alces alces) across two western boreal forest landscapes varying in disturbance. We quantified behavior as the number of camera trap photos per detection event and tested its relationship to inferred human‐mediated predation risk between a landscape with greater industrial disturbance and predator activity and a “control” landscape with lower human and predator activity. We also assessed the finer‐scale influence on behavior of variation in predation risk (relative to habitat variation) across camera sites within the more disturbed landscape. We predicted that animals in areas with greater predation risk (e.g., more wolf activity, less cover) would travel faster past cameras and generate fewer photos per detection event, while animals in areas with less predation risk would linger (rest, forage, investigate), generating more photos per event. Our predictions were supported at the landscape‐level, as caribou and moose had more photos per event in the control landscape where disturbance‐mediated predation risk was lower. At a finer‐scale within the disturbed landscape, no prey species showed a significant behavioral response to wolf activity, but the number of photos per event decreased for white‐tailed deer with increasing line of sight (m) along seismic lines (i.e., decreasing visual cover), consistent with a predator‐sensitive response. The presence of juveniles was associated with shorter behavioral events for caribou and moose, suggesting greater predator sensitivity for females with calves. Only moose demonstrated a positive behavioral association (i.e., longer events) with vegetation productivity (16‐day NDVI), suggesting that for other species bottom‐up influences of forage availability were generally weaker than top‐down influences from predation risk. Behavioral insights can be gleaned from camera trap surveys and provide complementary information about animal responses to predation risk, and thus about the indirect impacts of human disturbances on predator–prey interactions.     

Molecular and structural mechanisms of ZZ domain‐mediated cargo selection by Nbr1

In selective autophagy, cargo selectivity is determined by autophagy receptors. However, it remains scarcely understood how autophagy receptors recognize specific protein cargos. In the fission yeast Schizosaccharomyces pombe, a selective autophagy pathway termed Nbr1‐mediated vacuolar targeting (NVT) employs Nbr1, an autophagy receptor conserved across eukaryotes including humans, to target cytosolic hydrolases into the vacuole. Here, we identify two new NVT cargos, the mannosidase Ams1 and the aminopeptidase Ape4, that bind competitively to the first ZZ domain of Nbr1 (Nbr1‐ZZ1). High‐resolution cryo‐EM analyses reveal how a single ZZ domain recognizes two distinct protein cargos. Nbr1‐ZZ1 not only recognizes the N‐termini of cargos via a conserved acidic pocket, similar to other characterized ZZ domains, but also engages additional parts of cargos in a cargo‐specific manner. Our findings unveil a single‐domain bispecific mechanism of autophagy cargo recognition, elucidate its underlying structural basis, and expand the understanding of ZZ domain‐mediated protein–protein interactions.     

Non‐zero‐sum neutrality test for the tropical rain forest community using long‐term between‐census data

For community ecologists, “neutral or not?” is a fundamental question, and thus, rejecting neutrality is an important first step before investigating the deterministic processes underlying community dynamics. Hubbell''s neutral model is an important contribution to the exploration of community dynamics, both technically and philosophically. However, the neutrality tests for this model are limited by a lack of statistical power, partly because the zero‐sum assumption of the model is unrealistic. In this study, we developed a neutrality test for local communities that implements non‐zero‐sum community dynamics and determines the number of new species (N _sp) between observations. For the non‐zero‐sum neutrality test, the model distributed the expected N _sp, as calculated by extensive simulations, which allowed us to investigate the neutrality of the observed community by comparing the observed N _sp with distributions of the expected N _sp derived from the simulations. For this comparison, we developed a new “non‐zero‐sum N _sp test,” which we validated by running multiple neutral simulations using different parameter settings. We found that the non‐zero‐sum N _sp test rejected neutrality at a near‐significance level, which justified the validity of our approach. For an empirical test, the non‐zero‐sum N _sp test was applied to real tropical tree communities in Panama and Malaysia. The non‐zero‐sum N _sp test rejected neutrality in both communities when the observation interval was long and N _sp was large. Hence, the non‐zero‐sum N _sp test is an effective way to examine neutrality and has reasonable statistical power to reject the neutral model, especially when the observed N _sp is large. This unique and simple approach is statistically powerful, even though it only employs two temporal sequences of community data. Thus, this test can be easily applied to existing datasets. In addition, application of the test will provide significant benefits for detecting changing biodiversity under climate change and anthropogenic disturbance.     

Endogenous piRNA-guided slicing triggers responder and trailer piRNA production from viral RNA in Aedes aegypti mosquitoes

In the germline of animals, PIWI interacting (pi)RNAs protect the genome against the detrimental effects of transposon mobilization. In Drosophila, piRNA-mediated cleavage of transposon RNA triggers the production of responder piRNAs via ping-pong amplification. Responder piRNA 3′ end formation by the nuclease Zucchini is coupled to the production of downstream trailer piRNAs, expanding the repertoire of transposon piRNA sequences. In Aedes aegypti mosquitoes, piRNAs are generated from viral RNA, yet, it is unknown how viral piRNA 3′ ends are formed and whether viral RNA cleavage gives rise to trailer piRNA production. Here we report that in Ae. aegypti, virus- and transposon-derived piRNAs have sharp 3′ ends, and are biased for downstream uridine residues, features reminiscent of Zucchini cleavage of precursor piRNAs in Drosophila. We designed a reporter system to study viral piRNA 3′ end formation and found that targeting viral RNA by abundant endogenous piRNAs triggers the production of responder and trailer piRNAs. Using this reporter, we identified the Ae. aegypti orthologs of Zucchini and Nibbler, two nucleases involved in piRNA 3′ end formation. Our results furthermore suggest that autonomous piRNA production from viral RNA can be triggered and expanded by an initial cleavage event guided by genome-encoded piRNAs.     

Piwi/piRNAs control food intake by promoting neuropeptide F expression in locusts

Animal feeding, which directly affects growth and metabolism, is an important physiological process. However, the contribution of PIWI proteins and PIWI‐interacting RNAs (piRNAs) to the regulatory mechanism of animal feeding is unknown. Here, we report a novel function of Piwi and piRNAs in regulating food intake in locusts. Our study shows that the locust can serve as a representative species for determining PIWI function in insects. Knockdown of Piwi1 expression suppresses anabolic processes and reduces food consumption and body weight. The reduction in food intake by knockdown of Piwi1 expression results from decreased expression of neuropeptide NPF1 in a piRNA‐dependent manner. Mechanistically, intronic piRNAs might enhance RNA splicing of NPF1 by preventing hairpin formation at the branch point sites. These results suggest a novel nuclear PIWI/piRNA‐mediated mechanism that controls food intake in the locust nervous system.     

Molecular basis of cross‐species ACE2 interactions with SARS‐CoV‐2‐like viruses of pangolin origin

Correction to: The EMBO Journal (2021) 40: e107786. DOI 10.15252/embj.2021107786 | Published online 8 June 2021The authors would like to add three references to the paper: Starr et al and Zahradník et al also reported that the Q498H or Q498R mutation has enhanced binding affinity to ACE2; and Liu et al reported on the binding of bat coronavirus to ACE2.Starr et al and Zahradník et al have now been cited in the Discussion section, and the following sentence has been corrected from:“According to our data, the SARS‐CoV‐2 RBD with Q498H increases the binding strength to hACE2 by 5‐fold, suggesting the Q498H mutant is more ready to interact with human receptor than the wildtype and highlighting the necessity for more strict control of virus and virus‐infected animals”.to“Here, according to our data and two recently published papers, the SARS‐CoV‐2 RBD with Q498H or Q498R increases the binding strength to hACE2 (Starr et al, 2020; Zahradník et al, 2021), suggesting the mutant with Q498H or Q498R is more ready to interact with human receptor than the wild type and highlighting the necessity for more strict control of virus and virus‐infected animals”.The Liu et al citation has been added to the following sentence:“In another paper published by our group recently, RaTG13 RBD was found to bind to hACE2 with much lower binding affinity than SARS‐CoV‐2 though RaTG13 displays the highest whole‐genome sequence identity (96.2%) with the SARS‐CoV‐2 (Liu et al, 2021)”.Additionally, the authors have added the GISAID accession IDs to the sequence names of the SARS‐CoV‐2 in two human samples (Discussion section). To make identification unambiguous, the sequence names have been updated from “SA‐lsf‐27 and SA‐lsf‐37” to “GISAID accession ID: EPI_ISL_672581 and EPI_ISL_672589”.Lastly, the authors declare in the Materials and Methods section that all experiments employed SARS‐CoV‐2 pseudovirus in cultured cells. These experiments were performed in a BSL‐2‐level laboratory and approved by Science and Technology Conditions Platform Office, Institute of Microbiology, Chinese Academy of Sciences.These changes are herewith incorporated into the paper.     

Neutral theory reveals the challenge of bending the curve for the post‐2020 Global Biodiversity Framework

In October, nations of the world will begin negotiations for the post‐2020 Global Biodiversity Framework under the Convention on Biological Diversity. An influential ambition is “bending the curve of biodiversity loss,” which aims to reverse the decline of global biodiversity indicators. A second relevant, yet less prominent, milestone is the 20th anniversary of the publication of The Unified Neutral Theory of Biodiversity and Biogeography. Here, I apply neutral theory to show how global biodiversity indicators for population size (Living Planet Index) and extinction threat (Red List Index) decline under neutral ecological drift. This demonstrates that declining indicators are not solely caused by deterministic species‐specific or geographical patterns of biodiversity loss. Instead, indicators are sensitive to nondirectional stochasticity. Thus, “bending the curve” could be assessed relative to a counterfactual based on neutral theory, rather than static baselines. If used correctly, the 20‐year legacy of neutral theory can be extended to make a valuable contribution to the post‐2020 Global Biodiversity Framework.     

SARS‐CoV‐2 nucleocapsid protein phase‐separates with RNA and with human hnRNPs

Tightly packed complexes of nucleocapsid protein and genomic RNA form the core of viruses and assemble within viral factories, dynamic compartments formed within the host cells associated with human stress granules. Here, we test the possibility that the multivalent RNA‐binding nucleocapsid protein (N) from severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) condenses with RNA via liquid–liquid phase separation (LLPS) and that N protein can be recruited in phase‐separated forms of human RNA‐binding proteins associated with SG formation. Robust LLPS with RNA requires two intrinsically disordered regions (IDRs), the N‐terminal IDR and central‐linker IDR, as well as the folded C‐terminal oligomerization domain, while the folded N‐terminal domain and the C‐terminal IDR are not required. N protein phase separation is induced by addition of non‐specific RNA. In addition, N partitions in vitro into phase‐separated forms of full‐length human hnRNPs (TDP‐43, FUS, hnRNPA2) and their low‐complexity domains (LCs). These results provide a potential mechanism for the role of N in SARS‐CoV‐2 viral genome packing and in host‐protein co‐opting necessary for viral replication and infectivity.