The bowfin genome illuminates the developmental evolution of ray-finned fishes

The bowfin (Amia calva) is a ray-finned fish that possesses a unique suite of ancestral and derived phenotypes, which are key to understanding vertebrate evolution. The phylogenetic position of bowfin as a representative of neopterygian fishes, its archetypical body plan and its unduplicated and slowly evolving genome make bowfin a central species for the genomic exploration of ray-finned fishes. Here we present a chromosome-level genome assembly for bowfin that enables gene-order analyses, settling long-debated neopterygian phylogenetic relationships. We examine chromatin accessibility and gene expression through bowfin development to investigate the evolution of immune, scale, respiratory and fin skeletal systems and identify hundreds of gene-regulatory loci conserved across vertebrates. These resources connect developmental evolution among bony fishes, further highlighting the bowfin’s importance for illuminating vertebrate biology and diversity in the genomic era.

The genome of the bowfin (Amia calva) illuminates the developmental evolution of ray-finned fishes

The bowfin fish (Amia calva) diverged before the genome duplication in teleost fishes, and its archetypical body plan and slow rate of molecular evolution make it a key species for genomic exploration as a basal representative of the neopterygian fishes. To investigate the evolution and development of ray-finned fishes, we generated a chromosome-level genome assembly for bowfin that enables gene-order analyses which settle its long-debated, phylogenetic relationship with gars. We analyze the genomic underpinnings of the bowfin’s unique combination of derived and ancestral phenotypes involving the immune system as well as scale, respiratory organ, and skeletal development. By detailing chromatin accessibility and gene expression through bowfin development, we connect developmental gene regulatory loci across vertebrates. We illustrate the utility of these genomic resources to connect developmental evolution across bony fishes, showing the importance of bowfin in understanding vertebrate biology and diversity.

Genetic architecture of gene regulation in Indonesian populations identifies QTLs associated with local ancestry and archaic introgression

Lack of diversity in human genomics limits our understanding of the genetic underpinnings of complex traits, hinders precision medicine, and contributes to health disparities. To map genetic effects on gene regulation in the underrepresented Indonesian population, we have integrated genotype, gene expression, and CpG methylation data from 115 participants across three island populations that capture the major sources of genomic diversity on the region. In a comparison with a European dataset, we identify 166 uniquely Indonesia-specific eQTLs, highlighting the benefits of performing association studies on non-European populations. By combining local ancestry and archaic introgression inference eQTLs and methylQTLs, we identify regulatory loci driven by modern Papuan ancestry as well as introgressed Denisovan and Neanderthal variation. GWAS colocalization connects QTLs detected here to hematological traits. Our findings illustrate how local ancestry and archaic introgression drive variation in gene regulation across genetically distinct and in admixed populations.

Co-Expression Network Analyses of Anthocyanin Biosynthesis Genes in Ruellia (Wild Petunias; Acanthaceae)

Background. Anthocyanins are major pigments contributing to flower coloration and as such knowledge of molecular architecture underlying the anthocyanin biosynthetic pathway (ABP) is key to understanding flower color diversification. To identify ABP structural genes and associated regulatory networks, we sequenced 16 transcriptomes generated from 10 species of Ruellia and then conducted co-expression analyses among resulting data. Results. Complete coding sequences for 12 candidate structural loci representing eight genes plus nine candidate regulatory loci were assembled. Analysis of non-synonymous/synonymous (dn/ds) mutation rates indicated all identified loci are under purifying selection, suggesting overall selection to prevent the accumulation of deleterious mutations. Additionally, upstream enzymes have lower rates of molecular evolution compared to downstream enzymes. However, site-specific tests of selection yielded evidence for positive selection at several sites, including four in F3'H2 and five in DFR3, and these sites are located in protein binding regions. A species-level phylogenetic tree constructed using a newly implemented hybrid transcriptome–RADseq approach implicates numerous flower color transitions among the 10 species. We found evidence of both regulatory and structural mutations to F3'5'H in helping to explain the evolution of red flowers from purple-flowered ancestors.Conclusions. Sequence comparisons and co-expression analyses of ABP loci revealed that mutations in regulatory loci are likely to play a greater role in flower color transitions in Ruellia compared to mutations in underlying structural genes.

Identification and analysis of splicing quantitative trait loci across multiple tissues in the human genome

We have developed an efficient and reproducible pipeline for the discovery of genetic variants affecting splicing (sQTLs), based on an approach that captures the intrinsically multivariate nature of this phenomenon. We employed it to analyze the multi-tissue transcriptome GTEx dataset, generating a comprehensive catalogue of sQTLs in the human genome. A core set of these sQTLs is shared across multiple tissues. Downstream analyses of this catalogue contribute to the understanding of the mechanisms underlying splicing regulation. We found that sQTLs often target the global splicing pattern of genes, rather than individual splicing events. Many of them also affect gene expression, but not always of the same gene, potentially uncovering regulatory loci that act on different genes through different mechanisms. sQTLs tend to be preferentially located in introns that are post-transcriptionally spliced, which would act as hotspots for splicing regulation. While many variants affect splicing patterns by directly altering the sequence of splice sites, many more modify the binding of RNA-binding proteins (RBPs) to target sequences within the transcripts. Genetic variants affecting splicing can have a phenotypic impact comparable or even stronger than variants affecting expression, with those that alter RBP binding playing a prominent role in disease.

Regulators of an ancient polyphenism evolved through episodic protein divergence and parallel gene radiations

Polyphenism is a form of developmental plasticity that transduces environmental cues into discontinuous, often disparate phenotypes. In some cases, polyphenism has been attributed to facilitating morphological diversification and even the evolution of novel traits. However, this process is predicated on the origins and evolutionary maintenance of genetic mechanisms that specify alternate developmental networks. When and how regulatory loci arise and change, specifically before and throughout the history of a polyphenism, is little understood. Here, we establish a phylogenetic and comparative molecular context for two dynamically evolving genes, eud-1 and seud-1 , which regulate polyphenism in the nematode Pristionchus pacificus . This species is dimorphic in its adult feeding-structures, allowing individuals to become microbivores or facultative predators depending on the environment. Although polyphenism regulation is increasingly well understood in P. pacificus , the polyphenism is far older than this species and has diversified morphologically to enable an array of ecological functions across polyphenic lineages. To bring this taxonomic diversity into a comparative context, we reconstructed the histories of eud-1 and seud-1 relative to the origin and diversification of polyphenism, finding that homologues of both genes have undergone lineage-specific radiations across polyphenic taxa. Further, we detected signatures of episodic diversifying selection on eud-1 , particularly in early diplogastrid lineages. Lastly, transgenic rescue experiments suggest that the gene's product has functionally diverged from its orthologue's in a non-polyphenic outgroup. In summary, we provide a comparative framework for the molecular components of a plasticity switch, enabling studies of how polyphenism, its regulation, and ultimately its targets evolve.

Dissecting the regulatory activity and sequence content of loci with exceptional numbers of transcription factor associations

DNA associated proteins (DAPs) classically regulate gene expression by binding to regulatory loci such as enhancers or promoters. As expanding catalogs of genome-wide DAP binding maps reveal thousands of loci that, unlike the majority of conventional enhancers and promoters, associate with dozens of different DAPs with apparently little regard for motif preference, an understanding of DAP association and coordination at such regulatory loci is essential to deciphering how these regions contribute to normal development and disease. In this study, we aggregated publicly available ChIP-seq data from 469 human DAPs assayed in three cell lines and integrated these data with an orthogonal dataset of 352 non-redundant, in vitro-derived motifs mapped to the genome within DNase hypersensitivity footprints in an effort to characterize regions of the genome that have exceptionally high numbers of DAP associations. We subsequently performed a massively parallel mutagenesis assay to search for sequence elements driving transcriptional activity at such loci and explored plausible biological mechanisms underlying their formation. We establish a generalizable definition for High Occupancy Target (HOT) loci and identify putative driver DAP motifs in HEPG2 cells, including HNF4A, SP1, SP5, and ETV4, that are highly prevalent and exhibit sequence conservation at HOT loci. The number of different DAPs associated with an element is positively associated with evidence of regulatory activity and, by systematically mutating 245 HOT loci, we localized regulatory activity to a central core region that depends on the motif sequences of our previously nominated driver DAPs. In sum, this work leverages the increasingly large number of DAP motif and ChIP-seq data publicly available to explore how DAP associations contribute to genome-wide transcriptional regulation.