SMARTer single cell total RNA sequencing.

Single cell RNA sequencing methods have been increasingly used to understand cellular heterogeneity. Nevertheless, most of these methods suffer from one or more limitations, such as focusing only on polyadenylated RNA, sequencing of only the 3' end of the transcript, an exuberant fraction of reads mapping to ribosomal RNA, and the unstranded nature of the sequencing data. Here, we developed a novel single cell strand-specific total RNA library preparation method addressing all the aforementioned shortcomings. Our method was validated on a microfluidics system using three different cancer cell lines undergoing a chemical or genetic perturbation and on two other cancer cell lines sorted in microplates. We demonstrate that our total RNA-seq method detects an equal or higher number of genes compared to classic polyA[+] RNA-seq, including novel and non-polyadenylated genes. The obtained RNA expression patterns also recapitulate the expected biological signal. Inherent to total RNA-seq, our method is also able to detect circular RNAs. Taken together, SMARTer single cell total RNA sequencing is very well suited for any single cell sequencing experiment in which transcript level information is needed beyond polyadenylated genes.

Unraveling the KNOTTED1 regulatory network in maize meristems.

KNOTTED1 (KN1)-like homeobox (KNOX) transcription factors function in plant meristems, self-renewing structures consisting of stem cells and their immediate daughters. We defined the KN1 cistrome in maize inflorescences and found that KN1 binds to several thousand loci, including 643 genes that are modulated in one or multiple tissues. These KN1 direct targets are strongly enriched for transcription factors (including other homeobox genes) and genes participating in hormonal pathways, most significantly auxin, demonstrating that KN1 plays a key role in orchestrating the upper levels of a hierarchical gene regulatory network that impacts plant meristem identity and function.

Gene regulatory interactions at lateral organ boundaries in maize.

Maize leaves have distinct tissues that serve specific purposes. The blade tilts back to photosynthesize and the sheath wraps around the stem to provide structural support and protect young leaves. At the junction between blade and sheath are the ligule and auricles, both of which are absent in the recessive liguleless1 (lg1) mutant. Using an antibody against LG1, we reveal LG1 accumulation at the site of ligule formation and in the axil of developing tassel branches. The dominant mutant Wavy auricle in blade1 (Wab1-R) produces ectopic auricle tissue in the blade and increases the domain of LG1 accumulation. We determined that wab1 encodes a TCP transcription factor by positional cloning and revertant analysis. Tassel branches are few and upright in the wab1 revertant tassel and have an increased branch angle in the dominant mutant. wab1 mRNA is expressed at the base of branches in the inflorescence and is necessary for LG1 expression. wab1 is not expressed in leaves, except in the dominant mutant. The domain of wab1 expression in the Wab1-R leaf closely mirrors the accumulation of LG1. Although wab1 is not needed to induce lg1 expression in the leaf, LG1 is needed to counteract the severe phenotype of the dominant Wab1-R mutant. The regulatory interaction of LG1 and WAB1 reveals a link between leaf shape and tassel architecture, and suggests the ligule is a boundary similar to that at the base of lateral organs.

A comparative study of ChIP-seq sequencing library preparation methods.

BACKGROUND: ChIP-seq is the primary technique used to investigate genome-wide protein-DNA interactions. As part of this procedure, immunoprecipitated DNA must undergo "library preparation" to enable subsequent high-throughput sequencing. To facilitate the analysis of biopsy samples and rare cell populations, there has been a recent proliferation of methods allowing sequencing library preparation from low-input DNA amounts. However, little information exists on the relative merits, performance, comparability and biases inherent to these procedures. Notably, recently developed single-cell ChIP procedures employing microfluidics must also employ library preparation reagents to allow downstream sequencing. RESULTS: In this study, seven methods designed for low-input DNA/ChIP-seq sample preparation (Accel-NGS® 2S, Bowman-method, HTML-PCR, SeqPlex™, DNA SMART™, TELP and ThruPLEX®) were performed on five replicates of 1 ng and 0.1 ng input H3K4me3 ChIP material, and compared to a "gold standard" reference PCR-free dataset. The performance of each method was examined for the prevalence of unmappable reads, amplification-derived duplicate reads, reproducibility, and for the sensitivity and specificity of peak calling. CONCLUSIONS: We identified consistent high performance in a subset of the tested reagents, which should aid researchers in choosing the most appropriate reagents for their studies. Furthermore, we expect this work to drive future advances by identifying and encouraging use of the most promising methods and reagents. The results may also aid judgements on how comparable are existing datasets that have been prepared with different sample library preparation reagents.

Unequal redundancy in maize knotted1 homeobox genes.

The knotted1 (kn1) homeobox (knox) gene family was first identified through gain-of-function dominant mutants in maize (Zea mays). Class I knox members are expressed in meristems but excluded from leaves. In maize, a loss-of-function phenotype has only been characterized for kn1. To assess the function of another knox member, we characterized a loss-of-function mutation of rough sheath1 (rs1). rs1-mum1 has no phenotype alone but exacerbates several aspects of the kn1 phenotype. In permissive backgrounds in which kn1 mutants grow to maturity, loss of a single copy of rs1 enhances the tassel branch reduction phenotype, while loss of both copies results in limited shoots. In less introgressed lines, double mutants can grow to maturity but are shorter. Using a KNOX antibody, we demonstrate that RS1 binds in vivo to some of the KN1 target genes, which could partially explain why KN1 binds many genes but modulates few. Our results demonstrate an unequal redundancy between knox genes, with a role for rs1 only revealed in the complete absence of kn1.

SARS-CoV-2 can infect human embryos.

The spread of SARS-CoV-2 has led to a devastating pandemic, with infections resulting in a range of symptoms collectively known as COVID-19. The full repertoire of human tissues and organs susceptible to infection is an area of active investigation, and some studies have implicated the reproductive system. The effects of COVID-19 on human reproduction remain poorly understood, and particularly the impact on early embryogenesis and establishment of a pregnancy are not known. In this work, we explore the susceptibility of early human embryos to SARS-CoV-2 infection. By using RNA-seq and immunofluorescence, we note that ACE2 and TMPRSS2, two canonical cell entry factors for SARS-CoV-2, are co-expressed in cells of the trophectoderm in blastocyst-stage preimplantation embryos. For the purpose of viral entry studies, we used fluorescent reporter virions pseudotyped with Spike (S) glycoprotein from SARS-CoV-2, and we observe robust infection of trophectoderm cells. This permissiveness could be attenuated with blocking antibodies targeting S or ACE2. When exposing human blastocysts to the live, fully infectious SARS-CoV-2, we detected cases of infection that compromised embryo health. Therefore, we identify a new human target tissue for SARS-CoV-2 with potential medical implications for reproductive health during the COVID-19 pandemic and its aftermath.

Compound leaf development in the palm Chamaedorea elegans is KNOX-independent.

PREMISE OF THE STUDY: How a leaf acquires its shape is a major and largely unresolved question in plant biology. This problem is particularly complex in the case of compound leaves, where the leaf blade is subdivided into leaflets. In many eudicots with compound leaves, class I KNOTTED1-LIKE HOMEOBOX (KNOX) genes are upregulated in the leaf primordium and promote leaflet initiation, while KNOX genes are restricted to the shoot apical meristem in simple-leaved plants. In monocots, however, little is known about the extent of KNOX contribution to compound leaf development, and we aimed to address this issue in the palm Chamaedorea elegans. METHODS: We investigated the accumulation pattern of KNOX proteins in shoot apical meristems and leaf primordia of the palm C. elegans using immunolocalization experiments. KEY RESULTS: KNOX proteins accumulated in vegetative and inflorescence apical meristems and in the subtending stem tissue, but not in the plicated regions of the leaf primordia. These plicated areas form during primary morphogenesis and are the only meristematic tissue in the developing primordium. In addition, KNOX proteins did not accumulate in any region of the developing leaf during secondary morphogenesis, when leaflets separate to create the final pinnately compound leaf. CONCLUSIONS: The compound leaf character in palms, C. elegans in particular and likely other pinnately compound palms, does not depend on the activities of KNOX proteins.

Unexpected allelic heterogeneity and spectrum of mutations in Fowler syndrome revealed by next-generation exome sequencing.

Protein coding genes constitute approximately 1% of the human genome but harbor 85% of the mutations with large effects on disease-related traits. Therefore, efficient strategies for selectively sequencing complete coding regions (i.e., "whole exome") have the potential to contribute our understanding of human diseases. We used a method for whole-exome sequencing coupling Agilent whole-exome capture to the Illumina DNA-sequencing platform, and investigated two unrelated fetuses from nonconsanguineous families with Fowler Syndrome (FS), a stereotyped phenotype lethal disease. We report novel germline mutations in feline leukemia virus subgroup C cellular-receptor-family member 2, FLVCR2, which has recently been shown to cause FS. Using this technology, we identified three types of genetic abnormalities: point-mutations, insertions-deletions, and intronic splice-site changes (first pathogenic report using this technology), in the fetuses who both were compound heterozygotes for the disease. Although revealing a high level of allelic heterogeneity and mutational spectrum in FS, this study further illustrates the successful application of whole-exome sequencing to uncover genetic defects in rare Mendelian disorders. Of importance, we show that we can identify genes underlying rare, monogenic and recessive diseases using a limited number of patients (n=2), in the absence of shared genetic heritage and in the presence of allelic heterogeneity.

Distal expression of knotted1 in maize leaves leads to reestablishment of proximal/distal patterning and leaf dissection.

Maize (Zea mays) leaves provide a useful system to study how proximal/distal patterning is established because of the distinct tissues found in the distal blade and the proximal sheath. Several mutants disrupt this pattern, including the dominant knotted1-like homeobox (knox) mutants. knox genes encode homeodomain proteins of the TALE superclass of transcription factors. Class I knox genes are expressed in the meristem and down-regulated as leaves initiate. Gain-of-function phenotypes result from misexpression in leaves. We identified a new dominant allele of maize knotted1, Kn1-DL, which contains a transposon insertion in the promoter in addition to a tandem duplication of the kn1 locus. In situ hybridization shows that kn1 is misexpressed in two different parts of the blade that correlate with the different phenotypes observed. When kn1 is misexpressed along the margins, flaps of sheath-like tissue form along the margins. Expression in the distal tip leads to premature termination of the midrib into a knot and leaf bifurcation. The gain-of-function phenotypes suggest that kn1 establishes proximal/distal patterning when expressed in distal locations and lead to the hypothesis that kn1 normally participates in the establishment of proximal/distal polarity in the incipient leaf.

Modulation of Bax Inhibitor-1 and cytosolic Ca2+ by cytokinins in Nicotiana tabacum cells.

The protein Bax Inhibitor-1 (BI-1) has recently emerged as a negative regulator of plant programmed cell death (PCD), but how it functions at the biochemical level remains unknown. To elucidate its regulation and mode of action, we used suspension cells of Nicotiana tabacum to study the effects of cytokinins (CKs) on the expression level of NtBI-1 via western analysis. We found that the NtBI-1 protein is up-regulated following treatments with CKs at concentrations inducing a stress response (determined by growth reduction and PR1a accumulation), but not at PCD-inducing concentrations. These data point toward a role for NtBI-1 in the stress response to CKs. Application of CKs was also accompanied by a rapid cytosolic Ca(2+) pulse, and inhibition of this pulse with La(3+) or EGTA partially restored viability, indicating a signaling role for Ca(2+) in CK-induced cell death. However, CK-induced NtBI-1 accumulation was not altered by pretreatment with La(3+), nor by treatment with several modulators of intracellular Ca(2+) homeostasis and signaling, suggesting that CK-dependent regulation of NtBI-1 accumulation is not directly mediated by Ca(2+).

Preparation of Low-Input and Ligation-Free ChIP-seq Libraries Using Template-Switching Technology.

Chromatin immunoprecipitation (ChIP) followed by high-throughput sequencing (ChIP-seq) has become the gold standard for mapping of transcription factors and histone modifications throughout the genome. However, for ChIP experiments involving few cells or targeting low-abundance transcription factors, the small amount of DNA recovered makes ligation of adapters very challenging. In this unit, we describe a ChIP-seq workflow that can be applied to small cell numbers, including a robust single-tube and ligation-free method for preparation of sequencing libraries from sub-nanogram amounts of ChIP DNA. An example ChIP protocol is first presented, resulting in selective enrichment of DNA-binding proteins and cross-linked DNA fragments immobilized on beads via an antibody bridge. This is followed by a protocol for fast and easy cross-linking reversal and DNA recovery. Finally, we describe a fast, ligation-free library preparation protocol, featuring DNA SMART technology, resulting in samples ready for Illumina sequencing. © 2016 by John Wiley & Sons, Inc.

Strand-Specific Transcriptome Sequencing Using SMART Technology.

Next-generation sequencing is empowering a deeper understanding of biology by enabling RNA expression analysis over the entire transcriptome with high sensitivity and a wide dynamic range. One powerful application within this field is stranded RNA sequencing (RNA-seq), which is necessary to distinguish overlapping genes and to conduct comprehensive annotation and quantification of long non-coding RNAs. Commonly used methods for generating strand-specific RNA-seq libraries are often complicated by protocols that require several rounds of enzymatic treatments and clean-up steps, making them time-intensive, insensitive, and unsuitable for processing several samples simultaneously. An additional challenge in the generation of RNA-seq libraries from total RNA involves the high amount of ribosomal RNA (rRNA) in the starting material. This unit presents streamlined workflows for generating strand-specific RNA-seq libraries from 10 ng to 1 µg total RNA, representing a minimum of 1000 cells, in less than 7 hr with minimal carryover rRNA. These methods allow scientists to evaluate the expression of all transcripts, including non-polyadenylated long non-coding RNAs, even in limited biological samples. Combination of the RNase H-based RiboGone rRNA removal system and SMARTer Stranded RNA-seq technology enables depletion of over 95% of rRNA from mammalian samples, and direct production of Illumina-ready libraries that maintain strand-of-origin information. An alternate method for low input of highly degraded samples is also presented. © 2016 by John Wiley & Sons, Inc.

Antisense down regulation of NtBI-1 in tobacco BY-2 cells induces accelerated cell death upon carbon starvation.

Bax inhibitor-1 (BI-1) protein is proposed to be a conserved programmed cell death suppressor. In this report, we investigate the anti-apoptotic function of plant BI-1 by antisense (AS) down regulation of NtBI-1 in Nicotiana tabacum cv. BY-2 cells. We observed that AS cell lines were more susceptible to autophagy, internucleosomal DNA fragmentation and death than control cells when subjected to sucrose starvation and hypo-osmotic shock, in agreement with a role of BI-1 as a death inhibitor.

The maize transcription factor KNOTTED1 directly regulates the gibberellin catabolism gene ga2ox1.

KNOTTED1 (KN1)-like homeobox (KNOX) transcription factors are involved in the establishment and maintenance of plant meristems; however, few direct targets of KNOX proteins have been recognized. Using a combination of double mutant analysis and biochemistry, we found that in maize (Zea mays), KN1 negatively modulates the accumulation of gibberellin (GA) through the control of ga2ox1, which codes for an enzyme that inactivates GA. The ga2ox1 mRNA level is elevated in immature leaves of dominant KNOX mutants and downregulated in reproductive meristems of the null allele kn1-e1. KN1 binds in vivo to an intron of ga2ox1 through a cis-regulatory element containing two TGAC motifs. VP16-KN1 activates transcription in planta from a chimeric promoter containing this binding site. The domains of expression of kn1 and ga2ox1 mRNAs overlap at the base of the shoot apical meristem and the base of newly initiated leaves, suggesting that KN1-mediated activation of ga2ox1 maintains a boundary between meristem cell identity and rapidly elongating cells of the shoot. The KN1 binding site is conserved in ga2ox1 genes of different grasses, suggesting that the local regulation of bioactive GA levels through KNOX proteins is a common theme in grasses.

Dual functions of the KNOTTED1 homeodomain: sequence-specific DNA binding and regulation of cell-to-cell transport.

Homeodomain proteins are well-characterized developmental regulators that control expression of target genes through sequence-specific DNA binding. The homeodomain forms a trihelical structure, with the third helix conferring specific interactions with the DNA major groove. A specific class of plant homeodomain proteins, called KNOX [KNOTTED1 (KN1)-like homeobox], also has the ability to signal between cells by directly trafficking through intercellular channels called plasmodesmata. Trafficking is mediated by a signal that is also contained within the homeodomain. Movement protein binding protein 2C was identified as a protein that interacts with the KN1 homeodomain and regulates the cell-to-cell trafficking of KN1 by sequestering the protein on microtubules. Therefore, KN1 has multiple potential cellular addresses, each of which is conferred by its homeodomain.

Molecular characterization of two plant BI-1 homologues which suppress Bax-induced apoptosis in human 293 cells.

To date, few homologues of animal programmed cell death (PCD) regulators have been identified in plants. Among these is the plant Bax Inhibitor-1 (BI-1) protein, which possesses, like its human counterpart, the ability to suppress Bax-induced lethality in yeast cells. As the role of BI-1 in the regulation of plant PCD remains to be elucidated, we cloned BnBI-1 and NtBI-1 from cDNA libraries of oilseed rape ( Brassica napus L.) and tobacco ( Nicotiana tabacum L.). The analysis of the deduced amino acid sequences of BnBI-1 and NtBI-1 indicated that these proteins share a relatively high level of identity with other plant BI-1 proteins (73-95%) as well as with animal BI-1 proteins (26-42%). Comparative analysis with other available plant BI-1 proteins allowed the establishment of a structural model presenting seven transmembrane domains. Moreover, transient co-transfection of Bax with BnBI-1 or NtBI-1 in human embryonic kidney 293 cells revealed that both proteins can substantially inhibit apoptosis induced by Bax overexpression. Localization studies were also conducted using stable transformation of tobacco BY-2 cells and Saccharomyces cerevisiae, or transient expression in tobacco leaves, with the fusion protein BnBI-1GFP under control of the cauliflower mosaic virus 35S promoter. All transformants showed a fluorescence pattern of distribution typical of an endoplasmic reticulum (ER) protein. Results from differential permeabilization experiments in BY-2 cells expressing BnBI-1GFP also showed that the C-terminus is located on the cytosolic side of the ER. Taken altogether, our results suggest that BI-1 is evolutionarily conserved and could act as a key regulator of a death pathway common to plants and animals.

The two isoforms of the cAMP-dependent protein kinase catalytic subunit are involved in the control of dimorphism in the human fungal pathogen Candida albicans.

We have cloned the Candida albicans TPK2 gene encoding a cAMP-dependent protein kinase (PKA) catalytic subunit and generated a tpk2 homozygous null mutant to assess its ability to germinate in liquid media. N-acetylglucosamine (GlcNAc)-induced germ-tube formation was attenuated in the tpk2 strain and enhanced by compounds that are known to increase the PKA activity in situ. Germination was completely blocked in the presence of the myristoylated derivative of the heat-stable PKA inhibitor (MyrPKI). These results indicate that TPK1 acts positively in regulating the morphogenetic transition in C. albicans in the absence of the TPK2 gene. We were able to identify an mRNA from this second form of PKA in both wild-type and tpk2 null mutant cells. We found that PKA activity measured in the mutant lacking the TPK2 gene was about 10% of that displayed by the wild-type. The finding that the germinative response of tpk2 null mutant to serum was severely diminished at low serum concentrations indicates that the level of PKA is an important determinant of filamentous growth at low serum concentrations. The extent of germination attained at higher serum concentrations (5%) was similar in the wild-type and in the tpk2 null mutant strains suggesting that under these conditions germination was triggered through a PKA-independent pathway.