Transposable element (TE) invasions have shaped vertebrate genomes over the course of evolution. They have contributed an extra layer of species-specific gene regulation by providing novel transcription factor binding sites. In humans, SINE-VNTR-Alu (SVA) elements are one of three still active TE families; approximately 2800 SVA insertions exist in the human genome, half of which are human-specific. TEs are often silenced by KRAB zinc finger (KZNF) proteins recruiting corepressor proteins that establish a repressive chromatin state. A number of KZNFs have been reported to bind SVAs, but their individual contribution to repressing SVAs and their roles in suppressing SVA-mediated gene-regulatory effects remains elusive. We analyzed the genome-wide binding profile for ZNF91 in human cells and found that ZNF91 interacts with the VNTR region of SVAs. Through CRISPR-Cas9-mediated deletion of ZNF91 in human embryonic stem cells, we established that loss of ZNF91 results in increased transcriptional activity of SVAs. In contrast, SVA activation was not observed upon genetic deletion of the ZNF611 gene encoding another strong SVA interactor. Epigenetic profiling confirmed the loss of SVA repression in the absence of ZNF91 and revealed that mainly evolutionary young SVAs gain gene activation-associated epigenetic modifications. Genes close to activated SVAs showed a mild up-regulation, indicating SVAs adopt properties of cis-regulatory elements in the absence of repression. Notably, genome-wide derepression of SVAs elicited the communal upregulation of KZNFs that reside in KZNF clusters. This phenomenon may provide new insights into the potential mechanisms used by the host genome to sense and counteract TE invasions.

SUMMARY Phosphorus (P) is a major element required for plant growth and development. To cope with P shortage, plants activate local and long‐distance signaling pathways, such as an increase in the production and exudation of strigolactones (SLs). The role of the latter in mitigating P deficiency is, however, still largely unknown. To shed light on this, we studied the transcriptional response to P starvation and replenishment in wild‐type rice and a SL mutant, dwarf10 (d10), and upon exogenous application of the synthetic SL GR24. P starvation resulted in major transcriptional alterations, such as the upregulation of P TRANSPORTER, SYG1/PHO81/XPR1 (SPX) and VACUOLAR PHOSPHATE EFFLUX TRANSPORTER. Gene Ontology (GO) analysis of the genes induced by P starvation showed enrichment in phospholipid catabolic process and phosphatase activity. In d10, P deficiency induced upregulation of genes enriched for sesquiterpenoid production, secondary shoot formation and metabolic processes, including lactone biosynthesis. Furthermore, several genes induced by GR24 treatment shared the same GO terms with P starvation‐induced genes, such as oxidation reduction, heme binding and oxidoreductase activity, hinting at the role that SLs play in the transcriptional reprogramming upon P starvation. Gene co‐expression network analysis uncovered a METHYL TRANSFERASE that displayed co‐regulation with known rice SL biosynthetic genes. Functional characterization showed that this gene encodes an enzyme catalyzing the conversion of carlactonoic acid to methyl carlactonoate. Our work provides a valuable resource to further studies on the response of crops to P deficiency and reveals a tool for the discovery of SL biosynthetic genes.

Genome-wide association studies (GWAS) have been highly informative in discovering disease-associated loci but are not designed to capture all structural variations in the human genome. Using long-read sequencing data, we discovered widespread structural variation within SINE-VNTR-Alu (SVA) elements, a class of great ape-specific transposable elements with gene-regulatory roles, which represents a major source of structural variability in the human population. We highlight the presence of structurally variable SVAs (SV-SVAs) in neurological disease–associated loci, and we further associate SV-SVAs to disease-associated SNPs and differential gene expression using luciferase assays and expression quantitative trait loci data. Finally, we genetically deleted SV-SVAs in the BIN1 and CD2AP Alzheimer's disease–associated risk loci and in the BCKDK Parkinson's disease–associated risk locus and assessed multiple aspects of their gene-regulatory influence in a human neuronal context. Together, this study reveals a novel layer of genetic variation in transposable elements that may contribute to identification of the structural variants that are the actual drivers of disease associations of GWAS loci.

Comparing the abundance of microbial communities between different groups or obtained under different experimental conditions using count sequence data is a challenging task due to various issues such as inflated zero counts, overdispersion, and non-normality. Several methods and procedures based on counts, their transformation and compositionality have been proposed in the literature to detect differentially abundant species in datasets containing hundreds to thousands of microbial species. Despite efforts to address the the large numbers of zeros present in microbiome datasets, even after careful data preprocessing, the performance of existing methods is impaired by the presence of inflated zero counts and group-wise structured zeros (i.e., all zero counts in a group). We propose and validate using extensive simulations an approach combining two differential abundance testing methods, namely DESeq2-ZINBWaVE and DESeq2, to address the issues of zero-inflation and group-wise structured zeros, respectively. This combined approach was subsequently successfully applied to two plant microbiome datasets that revealed a number of taxa as interesting candidates for further experimental validation.