CTCF blocks anti-sense transcription initiation at divergent gene promoters

Transcription at most promoters is divergent, initiating at closely spaced oppositely oriented core promoters to produce sense transcripts along with often unstable upstream antisense (uasTrx). How antisense transcription is regulated and to what extent it is coordinated with sense transcription is largely unknown. Here by combining acute degradation of the multi-functional transcription factor CTCF and nascent transcription measurements, we find that CTCF specifically suppresses antisense but not sense transcription at hundreds of divergent promoters, the great majority of which bear proximal CTCF binding sites. Genome editing, chromatin conformation studies, and high-resolution transcript mapping revealed that precisely positioned CTCF directly suppresses the initiation of uasTrx, in a manner independent of its chromatin architectural function. Primary transcript RNA FISH revealed co-bursting of sense and anti-sense transcripts is disfavored, suggesting CTCF-regulated competition for transcription initiation. In sum, CTCF shapes the transcriptional landscape in part by suppressing upstream antisense transcription.

RNA-Seq reveals divergent gene expression between larvae with contrasting trophic modes in the poecilogonous polychaete Boccardia wellingtonensis

The polychaete Boccardia wellingtonensis is a poecilogonous species that produces different larval types. Females may lay Type I capsules, in which only planktotrophic larvae are present, or Type III capsules that contain planktotrophic and adelphophagic larvae as well as nurse eggs. While planktotrophic larvae do not feed during encapsulation, adelphophagic larvae develop by feeding on nurse eggs and on other larvae inside the capsules and hatch at the juvenile stage. Previous works have not found differences in the morphology between the two larval types; thus, the factors explaining contrasting feeding abilities in larvae of this species are still unknown. In this paper, we use a transcriptomic approach to study the cellular and genetic mechanisms underlying the different larval trophic modes of B. wellingtonensis. By using approximately 624 million high-quality reads, we assemble the de novo transcriptome with 133,314 contigs, coding 32,390 putative proteins. We identify 5221 genes that are up-regulated in larval stages compared to their expression in adult individuals. The genetic expression profile differed between larval trophic modes, with genes involved in lipid metabolism and chaetogenesis over expressed in planktotrophic larvae. In contrast, up-regulated genes in adelphophagic larvae were associated with DNA replication and mRNA synthesis.

The Bradyrhizobium Sp. LmicA16 Type VI Secretion System is Required for Efficient Nodulation of Lupinus Spp.

Many bacteria of the genus Bradyrhizobium are capable of inducing nodules in legumes. In this work, the importance of a type VI secretion system (T6SS) in a symbiotic strain of the genus Bradyrhizobium is described. T6SS of Bradyrhizobium sp. LmicA16 (A16) is necessary for efficient nodulation with Lupinus micranthus and L. angustifolius . A mutant in the gene vgrG, coding for a component of the T6SS nanostructure, induced less nodules and smaller plants than the wild type strain (wt) and was less competitive when co-inoculated with the wt strain. A16 T6SS genes are organized in a 26 kb DNA region in two divergent gene clusters of nine genes each. One of these genes codes for a protein (Tsb1) of unknown function but containing a methyltransferase domain. A tsb1 mutant showed an intermediate symbiotic phenotype regarding vgrG mutant and higher mucoidy and motility than the wt strain in free living conditions. T6SS promoter fusions to the lacZ reporter indicate expression in nodules but not in free living cells grown in different media and conditions. The analysis of nodule structure revealed that the level of nodule colonization was significantly reduced in the mutants with respect to the wt strain.

Tracing the evolution of human gene regulation and its association with shifts in environment

As humans spread throughout the world, they adapted to variation in many environmental factors, including climate, diet, and pathogens. Because many of these adaptations were likely mediated by multiple non-coding variants with small effects on gene regulation, it has been difficult to link genomic signals of selection to specific genes, and to describe the regulatory response to selection. To overcome this challenge, we adapted PrediXcan, a machine learning method for imputing gene regulation from genotype data, to analyze low-coverage ancient human DNA (aDNA). First, we used simulated genomes to benchmark strategies for adapting gene regulatory prediction to increase robustness to incomplete aDNA data. Applying the resulting models to 490 ancient Eurasians, we found that genes with the strongest divergent regulation among ancient populations with hunter-gatherer, pastoralist, and agricultural lifestyles are enriched for metabolic and immune functions. Next, we explored the contribution of divergent gene regulation to two traits with strong evidence of recent adaptation: dietary metabolism and skin pigmentation. We found enrichment for divergent regulation among genes previously proposed to be involved in diet-related local adaptation, and in many cases, the predicted effects on regulation provide explanations for previously observed signals of selection, e.g., at FADS1, GPX1, and LEPR. For skin pigmentation, we applied new models trained in melanocytes to a time series of 2999 ancient Europeans spanning ~38,000 years BP. In contrast to diet, skin pigmentation genes show little regulatory change over time, suggesting that adaptation mainly involved large-effect coding variants. This work demonstrates how aDNA can be combined with present-day genomes to shed light on the biological differences among ancient populations, the role of gene regulation in adaptation, and the relationship between ancient genetic diversity and the present-day distribution of complex traits.

Polygenic routes lead to parallel altitudinal adaptation in Heliosperma pusillum (Caryophyllaceae)

Understanding how organisms adapt to the environment is a major goal of modern biology. Parallel evolution - the independent evolution of similar phenotypes in different populations - provides a powerful framework to explore this question. Here, we quantified the degree of gene expression and functional parallelism across replicated ecotype formation in Heliosperma pusillum (Caryophyllaceae) and gained insights into the architecture of adaptive traits. Population structure analyses and demographic modelling confirm the previously formulated hypothesis of parallel polytopic divergence of montane and alpine ecotypes. We detect a large proportion of differentially expressed genes (DEGs) underlying adaptation of each replicate ecotype pair, with a strikingly low amount of shared DEGs across pairs. Functional enrichment of DEGs reveals that the traits affected by divergent gene expression are the same across ecotype pairs, in strong contrast to the non-shared genetic basis. The remarkable redundancy of differential gene expression indicates that diverged adaptive traits are highly polygenic. We conclude that polygenic traits appear key to opening multiple routes for adaptation, widening the adaptive potential of organisms.

Independent pseudogenizations and losses of sox15 during amniote diversification following asymmetric ohnolog evolution

Background Four ohnologous genes (sox1, sox2, sox3, and sox15) were generated by two rounds of whole-genome duplication in a vertebrate ancestor. In eutherian mammals, Sox1, Sox2, and Sox3 participate in central nervous system (CNS) development. Sox15 has a function in skeletal muscle regeneration and has little functional overlap with the other three ohnologs. In contrast, the frog Xenopus laevis and zebrafish orthologs of sox15 as well as sox1-3 function in CNS development. We previously reported that Sox15 is involved in mouse placental development as neofunctionalization, but is pseudogenized in the marsupial opossum. These findings suggest that sox15 might have evolved with divergent gene fates during vertebrate evolution. However, knowledge concerning sox15 in other vertebrate lineages than therian mammals, anuran amphibians, and teleost fish is scarce. Our purpose in this study was to clarify the fate and molecular evolution of sox15 during vertebrate evolution. Results We searched for sox15 orthologs in all vertebrate classes from agnathans to mammals by significant sequence similarity and synteny analyses using vertebrate genome databases. Interestingly, sox15 was independently pseudogenized at least twice during diversification of the marsupial mammals. Moreover, we observed independent gene loss of sox15 at least twice during reptile evolution in squamates and crocodile-bird diversification. Codon-based phylogenetic tree and selective analyses revealed an increased dN/dS ratio for sox15 compared to the other three ohnologs during jawed vertebrate evolution. Conclusions The findings revealed an asymmetric evolution of sox15 among the four ohnologs during vertebrate evolution, which was supported by the increased dN/dS values in cartilaginous fishes, anuran amphibians, and amniotes. The increased dN/dS value of sox15 may have been caused mainly by relaxed selection. Notably, independent pseudogenizations and losses of sox15 were observed during marsupial and reptile evolution, respectively. Both might have been caused by strong relaxed selection. The drastic gene fates of sox15, including neofunctionalization and pseudogenizations/losses during amniote diversification, might be caused by a release from evolutionary constraints.

Diversity of Conopeptides and Conoenzymes from the Venom Duct of the Marine Cone Snail Conus bayani as Determined from Transcriptomic and Proteomic Analyses

Marine cone snails are predatory gastropods characterized by a well-developed venom apparatus and highly evolved hunting strategies that utilize toxins to paralyze prey and defend against predators. The venom of each species of cone snail has a large number of pharmacologically active peptides known as conopeptides or conotoxins that are usually unique in each species. Nevertheless, venoms of only very few species have been characterized so far by transcriptomic approaches. In this study, we used transcriptome sequencing technologies and mass spectrometric methods to describe the diversity of venom components expressed by a worm-hunting species, Conus bayani. A total of 82 conotoxin sequences were retrieved from transcriptomic data that contain 54 validated conotoxin sequences clustered into 21 gene superfamilies including divergent gene family, 17 sequences clustered to 6 different conotoxin classes, and 11 conotoxins classified as unassigned gene family. Seven new conotoxin sequences showed unusual cysteine patterns. We were also able to identify 19 peptide sequences using mass spectrometry that completely overlapped with the conotoxin sequences obtained from transcriptome analysis. Importantly, herein we document the presence of 16 proteins that include five post-translational modifying enzymes obtained from transcriptomic data. Our results revealed diverse and novel conopeptides of an unexplored species that could be used extensively in biomedical research due to their therapeutic potentials.

Divergent and Overlapping Roles for Selected Phytochemicals in the Regulation of Pathological Cardiac Hypertrophy

Plant-based foods, like fruits, vegetables, whole grains, legumes, nuts, seeds and other foodstuffs, have been deemed as heart healthy. The chemicals within these plant-based foods, i.e., phytochemicals, are credited with protecting the heart. However, the mechanistic actions of phytochemicals, which prevent clinical endpoints, such as pathological cardiac hypertrophy, are still being elucidated. We sought to characterize the overlapping and divergent mechanisms by which 18 selected phytochemicals prevent phenylephrine- and phorbol 12-myristate 13-acetate-mediated cardiomyocyte enlargement. Of the tested 18 compounds, six attenuated PE- and PMA-mediated enlargement of neonatal rat ventricular myocytes. Cell viability assays showed that apigenin, baicalein, berberine hydrochloride, emodin, luteolin and quercetin dihydrate did not reduce cell size through cytotoxicity. Four of the six phytochemicals, apigenin, baicalein, berberine hydrochloride and emodin, robustly inhibited stress-induced hypertrophy and were analyzed further against intracellular signaling and genome-wide changes in mRNA expression. The four phytochemicals differentially regulated mitogen-activated protein kinases and protein kinase D. RNA-sequencing further showed divergence in gene regulation, while pathway analysis demonstrated overlap in the regulation of inflammatory pathways. Combined, this study provided a comprehensive analysis of cardioprotective phytochemicals. These data highlight two defining observations: (1) that these compounds predominantly target divergent gene pathways within cardiac myocytes and (2) that regulation of overlapping signaling and gene pathways may be of particular importance for the anti-hypertrophic actions of these phytochemicals. Despite these new findings, future works investigating rodent models of heart failure are still needed to understand the roles for these compounds in the heart.

Comparative Transcriptomics Analyses in Livers of Mice, Humans, and Humanized Mice Define Human-Specific Gene Networks

Mouse is the most widely used animal model in biomedical research, but it remains unknown what causes the large number of differentially regulated genes between human and mouse livers identified in recent years. In this report, we aim to determine whether these divergent gene regulations are primarily caused by environmental factors or some of them are the result of cell-autonomous differences in gene regulation in human and mouse liver cells. The latter scenario would suggest that many human genes are subject to human-specific regulation and can only be adequately studied in a human or humanized system. To understand the similarity and divergence of gene regulation between human and mouse livers, we performed stepwise comparative analyses in human, mouse, and humanized livers with increased stringency to gradually remove the impact of factors external to liver cells, and used bioinformatics approaches to retrieve gene networks to ascertain the regulated biological processes. We first compared liver gene regulation by fatty liver disease in human and mouse under the condition where the impact of genetic and gender biases was minimized, and identified over 50% of all commonly regulated genes, that exhibit opposite regulation by fatty liver disease in human and mouse. We subsequently performed more stringent comparisons when a single specific transcriptional or post-transcriptional event was modulated in vitro or vivo or in liver-specific humanized mice in which human and mouse hepatocytes colocalize and share a common circulation. Intriguingly and strikingly, the pattern of a high percentage of oppositely regulated genes persists under well-matched conditions, even in the liver of the humanized mouse model, which represents the most closely matched in vivo condition for human and mouse liver cells that is experimentally achievable. Gene network analyses further corroborated the results of oppositely regulated genes and revealed substantial differences in regulated biological processes in human and mouse cells. We also identified a list of regulated lncRNAs that exhibit very limited conservation and could contribute to these differential gene regulations. Our data support that cell-autonomous differences in gene regulation might contribute substantially to the divergent gene regulation between human and mouse livers and there are a significant number of biological processes that are subject to human-specific regulation and need to be carefully considered in the process of mouse to human translation.