A High-Throughput Screening System for Populus Wood-Associated Transcription Factors and Its Application to Lignin Regulation

Wood formation of trees is a complex and costly developmental process, whose regulatory network is involved in the protein-protein and protein-DNA interactions. To detect such interactions in wood development, we developed a high-throughput screening system with 517 Gal4-AD-wood-associated transcription factors (TFs) library from Populus alba × P. glandulosa cv “84K.” This system can be used for screening the upstream regulators and interacting proteins of targets by mating-based yeast-one hybrid (Y1H) and yeast-two-hybrid (Y2H) method, respectively. Multiple regulatory modules of lignin biosynthesis were identified based on this Populus system. Five TFs interacted with the 500-bp promoter fragment of PHENYLALANINE AMMONIA-LYASE 2 (PAL2), the first rate-limiting enzyme gene in the lignin biosynthesis pathway, and 10 TFs interacted with PaMYB4/LTF1, a key regulator of lignin biosynthesis. Some of these interactions were further validated by EMSA and BiFC assays. The TF-PaPAL2 promoter interaction and TF-PaMYB4 interaction revealed a complex mechanism governing the regulation of lignin synthesis in wood cells. Our high-throughput Y1H/Y2H screening system may be an efficient tool for studying regulatory network of wood formation in tree species.

Identification of miRNA–mRNA Regulatory Modules Involved in Lipid Metabolism and Seed Development in a Woody Oil Tree (Camellia oleifera)

Tea oil camellia (Camellia oleifera), an important woody oil tree, is a source of seed oil of high nutritional and medicinal value that is widely planted in southern China. However, there is no report on the identification of the miRNAs involved in lipid metabolism and seed development in the high- and low-oil cultivars of tea oil camellia. Thus, we explored the roles of miRNAs in the key periods of oil formation and accumulation in the seeds of tea oil camellia and identified miRNA–mRNA regulatory modules involved in lipid metabolism and seed development. Sixteen small RNA libraries for four development stages of seed oil biosynthesis in high- and low-oil cultivars were constructed. A total of 196 miRNAs, including 156 known miRNAs from 35 families, and 40 novel miRNAs were identified, and 55 significantly differentially expressed miRNAs were found, which included 34 upregulated miRNAs, and 21 downregulated miRNAs. An integrated analysis of the miRNA and mRNA transcriptome sequence data revealed that 10 miRNA–mRNA regulatory modules were related to lipid metabolism; for example, the regulatory modules of ath-miR858b–MYB82/MYB3/MYB44 repressed seed oil biosynthesis, and a regulation module of csi-miR166e-5p–S-ACP-DES6 was involved in the formation and accumulation of oleic acid. A total of 23 miRNA–mRNA regulatory modules were involved in the regulation of the seed size, such as the regulatory module of hpe-miR162a_L-2–ARF19, involved in early seed development. A total of 12 miRNA–mRNA regulatory modules regulating growth and development were identified, such as the regulatory modules of han-miR156a_L+1–SPL4/SBP2, promoting early seed development. The expression changes of six miRNAs and their target genes were validated using quantitative real-time PCR, and the targeting relationship of the cpa-miR393_R-1–AFB2 regulatory module was verified by luciferase assays. These data provide important theoretical values and a scientific basis for the genetic improvement of new cultivars of tea oil camellia in the future.

At the Root of Nodule Organogenesis: Conserved Regulatory Pathways Recruited by Rhizobia

The interaction between legume plants and soil bacteria rhizobia results in the formation of new organs on the plant roots, symbiotic nodules, where rhizobia fix atmospheric nitrogen. Symbiotic nodules represent a perfect model to trace how the pre-existing regulatory pathways have been recruited and modified to control the development of evolutionary “new” organs. In particular, genes involved in the early stages of lateral root development have been co-opted to regulate nodule development. Other regulatory pathways, including the players of the KNOX-cytokinin module, the homologues of the miR172-AP2 module, and the players of the systemic response to nutrient availability, have also been recruited to a unique regulatory program effectively governing symbiotic nodule development. The role of the NIN transcription factor in the recruitment of such regulatory modules to nodulation is discussed in more details.

Regulatory modules of human thermogenic adipocytes: functional genomics of large cohort and Meta-analysis derived marker-genes

Background Recently, ProFAT and BATLAS studies identified brown and white adipocytes marker genes based on analysis of large databases. They offered scores to determine the thermogenic status of adipocytes using the gene-expression data of these markers. In this work, we investigated the functional context of these genes. Results Gene Set Enrichment Analyses (KEGG, Reactome) of the BATLAS and ProFAT marker-genes identified pathways deterministic in the formation of brown and white adipocytes. The collection of the annotated proteins of the defined pathways resulted in expanded white and brown characteristic protein-sets, which theoretically contain all functional proteins that could be involved in the formation of adipocytes. Based on our previously obtained RNA-seq data, we visualized the expression profile of these proteins coding genes and found patterns consistent with the two adipocyte phenotypes. The trajectory of the regulatory processes could be outlined by the transcriptional profile of progenitor and differentiated adipocytes, highlighting the importance of suppression processes in browning. Protein interaction network-based functional genomics by STRING, Cytoscape and R-Igraph platforms revealed that different biological processes shape the brown and white adipocytes and highlighted key regulatory elements and modules including GAPDH-CS, DECR1, SOD2, IL6, HRAS, MTOR, INS-AKT, ERBB2 and 4-NFKB, and SLIT-ROBO-MAPK. To assess the potential role of a particular protein in shaping adipocytes, we assigned interaction network location-based scores (betweenness centrality, number of bridges) to them and created a freely accessible platform, the AdipoNET (https//adiponet.com), to conveniently use these data. The Eukaryote Promoter Database predicted the response elements in the UCP1 promoter for the identified, potentially important transcription factors (HIF1A, MYC, REL, PPARG, TP53, AR, RUNX, and FoxO1). Conclusion Our integrative approach-based results allowed us to investigate potential regulatory elements of thermogenesis in adipose tissue. The analyses revealed that some unique biological processes form the brown and white adipocyte phenotypes, which presumes the existence of the transitional states. The data also suggests that the two phenotypes are not mutually exclusive, and differentiation of thermogenic adipocyte requires induction of browning as well as repressions of whitening. The recognition of these simultaneous actions and the identified regulatory modules can open new direction in obesity research.

Discrete regulatory modules instruct hematopoietic lineage commitment and differentiation

Lineage commitment and differentiation is driven by the concerted action of master transcriptional regulators at their target chromatin sites. Multiple efforts have characterized the key transcription factors (TFs) that determine the various hematopoietic lineages. However, the temporal interactions between individual TFs and their chromatin targets during differentiation and how these interactions dictate lineage commitment remains poorly understood. Here we perform dense, daily, temporal profiling of chromatin accessibility (DNase I-seq) and gene expression changes (total RNA-seq) along ex vivo human erythropoiesis to comprehensively define developmentally regulated DNase I hypersensitive sites (DHSs) and transcripts. We link both distal DHSs to their target gene promoters and individual TFs to their target DHSs, revealing that the regulatory landscape is organized in distinct sequential regulatory modules that regulate lineage restriction and maturation. Finally, direct comparison of transcriptional dynamics (bulk and single-cell) and lineage potential between erythropoiesis and megakaryopoiesis uncovers differential fate commitment dynamics between the two lineages as they exit the stem and progenitor stage. Collectively, these data provide insights into the temporally regulated synergy of the cis- and the trans-regulatory components underlying hematopoietic lineage commitment and differentiation.

Single-sample proteome enrichment enables missing protein recovery and phenotype association

Proteomic studies characterize the protein composition of complex biological samples. Despite recent developments in mass spectrometry instrumentation and computational tools, low proteome coverage remains a challenge. To address this, we present Proteome Support Vector Enrichment (PROSE), a fast, scalable, and effective pipeline for scoring protein identifications based on gene co-expression matrices. Using a simple set of observed proteins as input, PROSE gauges the relative importance of proteins in the phenotype. The resultant enrichment scores are interpretable and stable, corresponding well to the source phenotype, thus enabling reproducible recovery of missing proteins. We further demonstrate its utility via reanalysis of the Cancer Cell Line Encyclopedia (CCLE) proteomic data, with prediction of oncogenic dependencies and identification of well-defined regulatory modules. PROSE is available as a user-friendly Python module from https://github.com/bwbio/PROSE.

Identification of limb-specific Lmx1b auto-regulatory modules with Nail-patella syndrome pathogenicity

LMX1B haploinsufficiency causes Nail-patella syndrome (NPS; MIM 161200), characterized by nail dysplasia, absent/hypoplastic patellae, chronic kidney disease, and glaucoma. Accordingly in mice, Lmx1b has been shown to play crucial roles in the development of the limb, kidney and eye. Although one functional allele of Lmx1b appears adequate for development, Lmx1b null mice display ventral-ventral distal limbs with abnormal kidney, eye and cerebellar development, more disruptive, but fully concordant with NPS. In Lmx1b functional knockouts (KOs), Lmx1b transcription in the limb is decreased nearly 6-fold, indicating autoregulation. Herein, we report on two conserved Lmx1b-associated cis-regulatory modules (LARM1 and LARM2) that are bound by Lmx1b, amplify Lmx1b expression with unique spatial modularity in the limb, and are necessary for Lmx1b-mediated limb dorsalization. These enhancers, being conserved across vertebrates (including coelacanth, but not other fish species), and required for normal locomotion, provide a unique opportunity to study the role of dorsalization in the fin to limb transition. We also report on two NPS patient families with normal LMX1B coding sequence, but with loss-of-function variations in the LARM1/2 region, stressing the role of regulatory modules in disease pathogenesis.

Large-scale miRNA-Target Data Analysis to Discover miRNA Co-regulation Network of Abiotic Stress Tolerance in Soybeans

Although growing evidence shows that microRNA (miRNA) regulates plant growth and development, miRNA regulatory networks in plants are not well understood. Current experimental studies cannot characterize miRNA regulatory networks on a large scale. This information gap provides a good opportunity to employ computational methods for global analysis and to generate useful models and hypotheses. To address this opportunity, we collected miRNA-target interactions (MTIs) and used MTIs from Arabidopsis thaliana and Medicago truncatula to predict homologous MTIs in soybeans, resulting in 80,235 soybean MTIs in total. A multi-level iterative bi-clustering method was developed to identify 483 soybean miRNA-target regulatory modules (MTRMs). Furthermore, we collected soybean miRNA expression data and corresponding gene expression data in response to abiotic stresses. By clustering these data, 37 MTRMs related to abiotic stresses were identified including stress-specific MTRMs and shared MTRMs. These MTRMs have gene ontology (GO) enrichment in resistance response, iron transport, positive growth regulation, etc. Our study predicts soybean miRNA-target regulatory modules with high confidence under different stresses, constructs miRNA-GO regulatory networks for MTRMs under different stresses and provides miRNA targeting hypotheses for experimental study. The method can be applied to other biological processes and other plants to elucidate miRNA co-regulation mechanisms.

Dog colour patterns explained by modular promoters of ancient canid origin

Distinctive colour patterns in dogs are an integral component of canine diversity. Colour pattern differences are thought to have arisen from mutation and artificial selection during and after domestication from wolves but important gaps remain in understanding how these patterns evolved and are genetically controlled. In other mammals, variation at the ASIP gene controls both the temporal and spatial distribution of yellow and black pigments. Here, we identify independent regulatory modules for ventral and hair cycle ASIP expression, and we characterize their action and evolutionary origin. Structural variants define multiple alleles for each regulatory module and are combined in different ways to explain five distinctive dog colour patterns. Phylogenetic analysis reveals that the haplotype combination for one of these patterns is shared with Arctic white wolves and that its hair cycle-specific module probably originated from an extinct canid that diverged from grey wolves more than 2 million years ago. Natural selection for a lighter coat during the Pleistocene provided the genetic framework for widespread colour variation in dogs and wolves.