Search

Populations and species are threatened by human pressure, but their fate is variable. Some depleted populations, such as that of the northern elephant seal (Mirounga angustirostris), recover rapidly even when the surviving population was small. The northern elephant seal was hunted extensively and taken by collectors between the early 1800s and 1892, suffering an extreme population bottleneck as a consequence. Recovery was rapid and now there are over 200,000 individuals. We sequenced 260 modern and 8 historical northern elephant seal nuclear genomes to assess the impact of the population bottleneck on individual northern elephant seals and to better understand their recovery. Here we show that inbreeding, an increase in the frequency of alleles compromised by lost function, and allele frequency distortion, reduced the fitness of breeding males and females, as well as the performance of adult females on foraging migrations. We provide a detailed investigation of the impact of a severe bottleneck on fitness at the genomic level and report on the role of specific gene systems.

The genus Acacia is a large group of woody legumes containing an enormous amount of morphological diversity in leaf shape. This diversity is at least in part the result of an innovation in leaf development where many Acacia species are capable of developing leaves of both bifacial and unifacial morphologies. While not unique in the plant kingdom, unifaciality is most commonly associated with monocots, and its developmental genetic mechanisms have yet to be explored beyond this group. In this study, we identify an accession of Acacia crassicarpa with high regeneration rates and isolate a clone for genome sequencing. We generate a chromosome-level assembly of this readily transformable clone, and using comparative analyses, confirm a whole-genome duplication unique to Caesalpinoid legumes. This resource will be important for future work examining genome evolution in legumes and the unique developmental genetic mechanisms underlying unifacial morphogenesis in Acacia.

Pathogen lineage nomenclature systems are a key component of effective communication and collaboration for researchers and public health workers. Since February 2021, the Pango dynamic lineage nomenclature for SARS-CoV-2 has been sustained by crowdsourced lineage proposals as new isolates were sequenced. This approach is vulnerable to time-critical delays as well as regional and personal bias. Here we developed a simple heuristic approach for dividing phylogenetic trees into lineages, including the prioritization of key mutations or genes. Our implementation is efficient on extremely large phylogenetic trees consisting of millions of sequences and produces similar results to existing manually curated lineage designations when applied to SARS-CoV-2 and other viruses including chikungunya virus, Venezuelan equine encephalitis virus complex and Zika virus. This method offers a simple, automated and consistent approach to pathogen nomenclature that can assist researchers in developing and maintaining phylogeny-based classifications in the face of ever-increasing genomic datasets.

Macrophages clear infections by engulfing and digesting pathogens within phagolysosomes. Pathogens escape this fate by engaging in a molecular arms race; they use WxxxE motif-containing "effector" proteins to subvert the host cells they invade and seek refuge within protective vacuoles. Here, we define the host component of the molecular arms race as an evolutionarily conserved polar "hot spot" on the PH domain of ELMO1 (Engulfment and Cell Motility protein 1), which is targeted by diverse WxxxE effectors. Using homology modeling and site-directed mutagenesis, we show that a lysine triad within the "patch" directly binds all WxxxE effectors tested: SifA (Salmonella), IpgB1 and IpgB2 (Shigella), and Map (enteropathogenic Escherichia coli). Using an integrated SifA-host protein-protein interaction network, in silico network perturbation, and functional studies, we show that the major consequences of preventing SifA-ELMO1 interaction are reduced Rac1 activity and microbial invasion. That multiple effectors of diverse structure, function, and sequence bind the same hot spot on ELMO1 suggests that the WxxxE effector(s)-ELMO1 interface is a convergence point of intrusion detection and/or host vulnerability. We conclude that the interface may represent the fault line in coevolved molecular adaptations between pathogens and the host, and its disruption may serve as a therapeutic strategy.

The Steller's jay is a familiar bird of western forests from Alaska south to Nicaragua. Here, we report a draft reference assembly for the species generated from PacBio HiFi long-read and Omni-C chromatin-proximity sequencing data as part of the California Conservation Genomics Project (CCGP). Sequenced reads were assembled into 352 scaffolds totaling 1.16 Gb in length. Assembly metrics indicate a highly contiguous and complete assembly with a contig N50 of 7.8 Mb, scaffold N50 of 25.8 Mb, and BUSCO completeness score of 97.2%. Repetitive elements span 16.6% of the genome including nearly 90% of the W chromosome. Compared with high-quality assemblies from other members of the family Corvidae, the Steller's jay genome contains a larger proportion of repetitive elements than 4 crow species (Corvus), but a lower proportion of repetitive elements than the California scrub-jay (Aphelocoma californica). This reference genome will serve as an essential resource for future studies on speciation, local adaptation, phylogeography, and conservation genetics in this species of significant biological interest.

Violacein is a water-insoluble violet pigment produced by various Gram-negative bacteria. The compound and the bacteria that produce it have been gaining attention due to the antimicrobial and proposed antitumour properties of violacein and the possibility that strains producing it may have broad industrial uses. Bacteria that produce violacein have been isolated from diverse environments including fresh and ocean waters, glaciers, tropical soils, trees, fish and the skin of amphibians. We report here the isolation and characterization of six violacein-producing bacterial strains and three non-violacein-producing close relatives, each isolated from either an aquatic environment or moist food materials in northern California, USA. For each isolate, we characterized traditional phenotypes, generated and analysed draft genome sequences, and carried out multiple types of taxonomic, phylogenetic and phylogenomic analyses. Based on these analyses we assign putative identifications to the nine isolates, which include representatives of the genera Chromobacterium, Aquitalea, Iodobacter, Duganella, Massilia and Janthinobacterium. In addition, we discuss the utility of various metrics for taxonomic assignment in these groups including average nucleotide identity, whole genome phylogenetic analysis and extent of recent homologous recombination using the software program PopCOGenT.

While oral desensitization is capable of alleviating peanut allergen anaphylaxis, long-term immune tolerance is the sought-after goal. We developed a liver-targeting lipid nanoparticle (LNP) platform to deliver mRNA-encoded peanut allergen epitopes to liver sinusoidal endothelial cells (LSECs), which function as robust tolerogenic antigen-presenting cells that induce FoxP3+ regulatory T-cells (Tregs). The mRNA strand was constructed by including nucleotide sequences encoding for nonallergenic MHC-II binding T-cell epitopes, identified in the dominant peanut allergen, Ara h2. These epitopes were inserted in the mRNA strand downstream of an MHC-II targeting sequence, further endowed in vitro with 5 and 3 capping sequences, a PolyA tail, and uridine substitution. Codon-optimized mRNA was used for microfluidics synthesis of LNPs with an ionizable cationic lipid, also decorated with a lipid-anchored mannose ligand for LSEC targeting. Biodistribution to the liver was confirmed by in vivo imaging, while ELISpot assays demonstrated an increase in IL-10-producing Tregs in the spleen. Prophylactic administration of tandem-repeat or a combination of encapsulated Ara h2 epitopes induced robust tolerogenic effects in C3H/HeJ mice, sensitized to and subsequently challenged with crude peanut allergen extract. In addition to alleviating physical manifestations of anaphylaxis, there was suppression of Th2-mediated cytokine production, IgE synthesis, and mast cell release, accompanied by increased IL-10 and TGF-β production in the peritoneum. Similar efficacy was demonstrated during LNP administration postsensitization. While nondecorated particles had lesser but significant effects, PolyA/LNP-Man lacked protective effects. These results demonstrate an exciting application of mRNA/LNP for treatment of food allergen anaphylaxis, with the promise to be widely applicable to the allergy field.

Improved sequencing technologies have profoundly altered global views of fungal diversity and evolution. High-throughput sequencing methods are critical for studying fungi due to the cryptic, symbiotic nature of many species, particularly those that are difficult to culture. However, the low coverage genome sequencing (LCGS) approach to phylogenomic inference has not been widely applied to fungi. Here we analyzed 171 Kickxellomycotina fungi using LCGS methods to obtain hundreds of marker genes for robust phylogenomic reconstruction. Additionally, we mined our LCGS data for a set of nine rDNA and protein coding genes to enable analyses across species for which no LCGS data were obtained. The main goals of this study were to: 1) evaluate the quality and utility of LCGS data for both phylogenetic reconstruction and functional annotation, 2) test relationships among clades of Kickxellomycotina, and 3) perform comparative functional analyses between clades to gain insight into putative trophic modes. In opposition to previous studies, our nine-gene analyses support two clades of arthropod gut dwelling species and suggest a possible single evolutionary event leading to this symbiotic lifestyle. Furthermore, we resolve the mycoparasitic Dimargaritales as the earliest diverging clade in the subphylum and find four major clades of Coemansia species. Finally, functional analyses illustrate clear variation in predicted carbohydrate active enzymes and secondary metabolites (SM) based on ecology, that is biotroph versus saprotroph. Saprotrophic Kickxellales broadly lack many known pectinase families compared with saprotrophic Mucoromycota and are depauperate for SM but have similar numbers of predicted chitinases as mycoparasitic.

Aberrant insulin-like growth factor 1 (IGF-1) signaling has been proposed as a contributing factor to the development of neurodegenerative disorders including diabetic neuropathy, and delivery of exogenous IGF-1 has been explored as a treatment for Alzheimer's disease and amyotrophic lateral sclerosis. However, the role of autocrine/paracrine IGF-1 in neuroprotection has not been well established. We therefore used in vitro cell culture systems and animal models of diabetic neuropathy to characterize endogenous IGF-1 in sensory neurons and determine the factors regulating IGF-1 expression and/or affecting neuronal health. Single-cell RNA sequencing (scRNA-Seq) and in situ hybridization analyses revealed high expression of endogenous IGF-1 in non-peptidergic neurons and satellite glial cells (SGCs) of dorsal root ganglia (DRG). Brain cortex and DRG had higher IGF-1 gene expression than sciatic nerve. Bidirectional transport of IGF-1 along sensory nerves was observed. Despite no difference in IGF-1 receptor levels, IGF-1 gene expression was significantly (P

A global international initiative, such as the Earth BioGenome Project (EBP), requires both agreement and coordination on standards to ensure that the collective effort generates rapid progress toward its goals. To this end, the EBP initiated five technical standards committees comprising volunteer members from the global genomics scientific community: Sample Collection and Processing, Sequencing and Assembly, Annotation, Analysis, and IT and Informatics. The current versions of the resulting standards documents are available on the EBP website, with the recognition that opportunities, technologies, and challenges may improve or change in the future, requiring flexibility for the EBP to meet its goals. Here, we describe some highlights from the proposed standards, and areas where additional challenges will need to be met.

Life on Earth has evolved from initial simplicity to the astounding complexity we experience today. Bacteria and archaea have largely excelled in metabolic diversification, but eukaryotes additionally display abundant morphological innovation. How have these innovations come about and what constraints are there on the origins of novelty and the continuing maintenance of biodiversity on Earth? The history of life and the code for the working parts of cells and systems are written in the genome. The Earth BioGenome Project has proposed that the genomes of all extant, named eukaryotes-about 2 million species-should be sequenced to high quality to produce a digital library of life on Earth, beginning with strategic phylogenetic, ecological, and high-impact priorities. Here we discuss why we should sequence all eukaryotic species, not just a representative few scattered across the many branches of the tree of life. We suggest that many questions of evolutionary and ecological significance will only be addressable when whole-genome data representing divergences at all of the branchings in the tree of life or all species in natural ecosystems are available. We envisage that a genomic tree of life will foster understanding of the ongoing processes of speciation, adaptation, and organismal dependencies within entire ecosystems. These explorations will resolve long-standing problems in phylogenetics, evolution, ecology, conservation, agriculture, bioindustry, and medicine.

Methods to deliver gene editing agents in vivo as ribonucleoproteins could offer safety advantages over nucleic acid delivery approaches. We report the development and application of engineered DNA-free virus-like particles (eVLPs) that efficiently package and deliver base editor or Cas9 ribonucleoproteins. By engineering VLPs to overcome cargo packaging, release, and localization bottlenecks, we developed fourth-generation eVLPs that mediate efficient base editing in several primary mouse and human cell types. Using different glycoproteins in eVLPs alters their cellular tropism. Single injections of eVLPs into mice support therapeutic levels of base editing in multiple tissues, reducing serum Pcsk9 levels 78% following 63% liver editing, and partially restoring visual function in a mouse model of genetic blindness. In vitro and in vivo off-target editing from eVLPs was virtually undetected, an improvement over AAV or plasmid delivery. These results establish eVLPs as promising vehicles for therapeutic macromolecule delivery that combine key advantages of both viral and nonviral delivery.

Evolutionary biologists have long sought to disentangle phylogenetic relationships among taxa spanning the tree of life, an increasingly important task as anthropogenic influences accelerate population declines and species extinctions, particularly in insects. Phylogenetic analyses are commonly used to identify unique evolutionary lineages, to clarify taxonomic designations of the focal taxa, and to inform conservation decisions. Advances in DNA sequencing techniques have increasingly facilitated the ability of researchers to apply genomic methods to phylogenetic analyses, even for non-model organisms. Stoneflies are non-model insects that are important bioindicators of the quality of freshwater habitats and landscape disturbance as they spend the immature stages of their life cycles in fresh water, and the adult stages in terrestrial environments. Phylogenetic relationships within the stonefly genus Suwallia (Insecta: Plecoptera: Chloroperlidae) are poorly understood, and have never been assessed using molecular data. We used DNA sequence data from genome-wide ultraconserved element loci to generate the first molecular phylogeny for the group and assess its monophyly. We found that Palearctic and Nearctic Suwallia do not form reciprocally monophyletic clades, and that a biogeographic history including dispersal, vicariance, and founder event speciation via jump dispersal best explains the geographic distribution of this group. Our results also strongly suggest that Neaviperla forcipata (Neave, 1929) is nested within Suwallia, and the concept of the genus Suwallia should be revised to include it. Thus, we formally propose a new taxonomic combination wherein Neaviperla forcipata (Neave, 1929) is reclassified as Suwallia forcipata (Neave, 1929). Moreover, some Suwallia species (e.g., S. amoenacolens, S. kerzhneri, S. marginata, S. pallidula, and S. starki) exhibit pronounced cryptic diversity that is worthy of further investigation. These findings provide a first glimpse into the evolutionary history of Suwallia, improve our understanding of stonefly diversity in the tribe Suwallini, and highlight areas where additional research is needed.

BackgroundSeahorses, seadragons, pygmy pipehorses, and pipefishes (Syngnathidae, Syngnathiformes) are among the most recognizable groups of fishes because of their derived morphology, unusual life history, and worldwide distribution. Despite previous phylogenetic studies and recent new species descriptions of syngnathids, the evolutionary relationships among several major groups within this family remain unresolved.ResultsHere, we provide a reconstruction of syngnathid phylogeny based on genome-wide sampling of 1314 ultraconserved elements (UCEs) and expanded taxon sampling to assess the current taxonomy and as a basis for macroevolutionary insights. We sequenced a total of 244 new specimens across 117 species and combined with published UCE data for a total of 183 species of Syngnathidae, about 62% of the described species diversity, to compile the most data-rich phylogeny to date. We estimated divergence times using 14 syngnathiform fossils, including nine fossils with newly proposed phylogenetic affinities, to better characterize current and historical biogeographical patterns, and to reconstruct diversification through time. We present a phylogenetic hypothesis that is well-supported and provides several notable insights into syngnathid evolution. We found nine non-monophyletic genera, evidence for seven cryptic species, five potentially invalid synonyms, and identified a novel sister group to the seahorses, the Indo-Pacific pipefishes Halicampus macrorhynchus and H. punctatus. In addition, the morphologically distinct southwest Pacific seahorse Hippocampus jugumus was recovered as the sister to all other non-pygmy seahorses. As found in many other groups, a high proportion of syngnathid lineages appear to have originated in the Central Indo-Pacific and subsequently dispersed to adjoining regions. Conversely, we also found an unusually high subsequent return of lineages from southern Australasia to the Central Indo-Pacific. Diversification rates rose abruptly during the Middle Miocene Climate Transition and peaked after the closure of the Tethys Sea.ConclusionsOur results reveal a previously underappreciated diversity of syngnathid lineages. The observed biogeographic patterns suggest a significant role of the southern Australasian region as a source and sink of lineages. Shifts in diversification rates imply possible links to declining global temperatures, the separation of the Atlantic and Pacific faunas, and the environmental changes associated with these events.

BackgroundTranscriptomic data has demonstrated utility to advance the study of physiological diversity and organisms' responses to environmental stressors. However, a lack of genomic resources and challenges associated with collecting high-quality RNA can limit its application for many wild populations. Minimally invasive blood sampling combined with de novo transcriptomic approaches has great potential to alleviate these barriers. Here, we advance these goals for marine turtles by generating high quality de novo blood transcriptome assemblies to characterize functional diversity and compare global transcriptional profiles between tissues, species, and foraging aggregations.ResultsWe generated high quality blood transcriptome assemblies for hawksbill (Eretmochelys imbricata), loggerhead (Caretta caretta), green (Chelonia mydas), and leatherback (Dermochelys coriacea) turtles. The functional diversity in assembled blood transcriptomes was comparable to those from more traditionally sampled tissues. A total of 31.3% of orthogroups identified were present in all four species, representing a core set of conserved genes expressed in blood and shared across marine turtle species. We observed strong species-specific expression of these genes, as well as distinct transcriptomic profiles between green turtle foraging aggregations that inhabit areas of greater or lesser anthropogenic disturbance.ConclusionsObtaining global gene expression data through non-lethal, minimally invasive sampling can greatly expand the applications of RNA-sequencing in protected long-lived species such as marine turtles. The distinct differences in gene expression signatures between species and foraging aggregations provide insight into the functional genomics underlying the diversity in this ancient vertebrate lineage. The transcriptomic resources generated here can be used in further studies examining the evolutionary ecology and anthropogenic impacts on marine turtles.

BackgroundGenes with multiple co-active promoters appear common in brain, yet little is known about functional requirements for these potentially redundant genomic regulatory elements. SCN1A, which encodes the NaV1.1 sodium channel alpha subunit, is one such gene with two co-active promoters. Mutations in SCN1A are associated with epilepsy, including Dravet syndrome (DS). The majority of DS patients harbor coding mutations causing SCN1A haploinsufficiency; however, putative causal non-coding promoter mutations have been identified.MethodsTo determine the functional role of one of these potentially redundant Scn1a promoters, we focused on the non-coding Scn1a 1b regulatory region, previously described as a non-canonical alternative transcriptional start site. We generated a transgenic mouse line with deletion of the extended evolutionarily conserved 1b non-coding interval and characterized changes in gene and protein expression, and assessed seizure activity and alterations in behavior.ResultsMice harboring a deletion of the 1b non-coding interval exhibited surprisingly severe reductions of Scn1a and NaV1.1 expression throughout the brain. This was accompanied by electroencephalographic and thermal-evoked seizures, and behavioral deficits.ConclusionsThis work contributes to functional dissection of the regulatory wiring of a major epilepsy risk gene, SCN1A. We identified the 1b region as a critical disease-relevant regulatory element and provide evidence that non-canonical and seemingly redundant promoters can have essential function.

Nullomers and nullpeptides are short DNA or amino acid sequences that are absent from a genome or proteome, respectively. One potential cause for their absence could be their having a detrimental impact on an organism. RESULTS: Here, we identify all possible nullomers and nullpeptides in the genomes and proteomes of thirty eukaryotes and demonstrate that a significant proportion of these sequences are under negative selection. We also identify nullomers that are unique to specific functional categories: coding sequences, exons, introns, 5'UTR, 3'UTR, promoters, and show that coding sequence and promoter nullomers are most likely to be selected against. By analyzing all protein sequences across the tree of life, we further identify 36,081 peptides up to six amino acids in length that do not exist in any known organism, termed primes. We next characterize all possible single base pair mutations that can lead to the appearance of a nullomer in the human genome, observing a significantly higher number of mutations than expected by chance for specific nullomer sequences in transposable elements, likely due to their suppression. We also annotate nullomers that appear due to naturally occurring variants and show that a subset of them can be used to distinguish between different human populations. Analysis of nullomers and nullpeptides across vertebrate evolution shows they can also be used as phylogenetic classifiers. CONCLUSIONS: We provide a catalog of nullomers and nullpeptides in distinct functional categories, develop methods to systematically study them, and highlight the use of variability in these sequences in other analyses.

Retrotransposons mediate gene regulation in important developmental and pathological processes. Here, we characterized the transient retrotransposon induction during preimplantation development of eight mammals. Induced retrotransposons exhibit similar preimplantation profiles across species, conferring gene regulatory activities, particularly through long terminal repeat (LTR) retrotransposon promoters. A mouse-specific MT2B2 retrotransposon promoter generates an N-terminally truncated Cdk2ap1ΔN that peaks in preimplantation embryos and promotes proliferation. In contrast, the canonical Cdk2ap1 peaks in mid-gestation and represses cell proliferation. This MT2B2 promoter, whose deletion abolishes Cdk2ap1ΔN production, reduces cell proliferation and impairs embryo implantation, is developmentally essential. Intriguingly, Cdk2ap1ΔN is evolutionarily conserved in sequence and function yet is driven by different promoters across mammals. The distinct preimplantation Cdk2ap1ΔN expression in each mammalian species correlates with the duration of its preimplantation development. Hence, species-specific transposon promoters can yield evolutionarily conserved, alternative protein isoforms, bestowing them with new functions and species-specific expression to govern essential biological divergence.

Estrogen receptor α (ERα) is a hormone receptor and key driver for over 70% of breast cancers that has been studied for decades as a transcription factor. Unexpectedly, we discover that ERα is a potent non-canonical RNA-binding protein. We show that ERα RNA binding function is uncoupled from its activity to bind DNA and critical for breast cancer progression. Employing genome-wide cross-linking immunoprecipitation (CLIP) sequencing and a functional CRISPRi screen, we find that ERα-associated mRNAs sustain cancer cell fitness and elicit cellular responses to stress. Mechanistically, ERα controls different steps of RNA metabolism. In particular, we demonstrate that ERα RNA binding mediates alternative splicing of XBP1 and translation of the eIF4G2 and MCL1 mRNAs, which facilitates survival upon stress conditions and sustains tamoxifen resistance of cancer cells. ERα is therefore a multifaceted RNA-binding protein, and this activity transforms our knowledge of post-transcriptional regulation underlying cancer development and drug response.

The PD-L1 overexpression is an important event of immune escape and metastasis in triple-negative breast cancer (TNBC), but the molecular mechanism remains to be determined. Interferon gamma (IFNγ) represents a major driving force behind PD-L1 expression in tumor microenvironment, and histone deacetylase 2 (HDAC2) is required for IFN signaling. Here, we investigated the regulation of HDAC2 on the IFNγ-induced PD-L1 expression in TNBC cells. We found the HDAC2 and PD-L1 expression in TNBC was significantly higher than that in non-TNBC, and HDAC2 was positively correlated with PD-L1 expression. HDAC2 promoted PD-L1 induction by upregulating the phosphorylation of JAK1, JAK2, and STAT1, as well as the translocation of STAT1 to the nucleus and the recruitment of STAT1 to the PD-L1 promoter. Meanwhile, HDAC2 was recruited to the PD-L1 promoter by STAT1, and HDAC2 knockout compromised IFNγ-induced upregulation of H3K27, H3K9 acetylation, and the BRD4 recruitment in PD-L1 promoter. In addition, significant inhibition of proliferation, colony formation, migration, and cell cycle of TNBC cells were observed following knockout of HDAC2 in vitro. Furthermore, HDAC2 knockout reduced IFNγ-induced PD-L1 expression, lymphocyte infiltration, and retarded tumor growth and metastasis in the breast cancer mouse models. This study may provide evidence that HDAC2 promotes IFNγ-induced PD-L1 expression, suggesting a way for enhanced antitumor immunity when targeting the HDAC2 in TNBC.

Spermatogonial stem cells (SSCs) are essential for male fertility. Here, we report that mouse SSC generation is driven by a transcription factor (TF) cascade controlled by the homeobox protein, RHOX10, which acts by driving the differentiation of SSC precursors called pro-spermatogonia (ProSG). We identify genes regulated by RHOX10 in ProSG in vivo and define direct RHOX10-target genes using several approaches, including a rapid temporal induction assay: iSLAMseq. Together, these approaches identify temporal waves of RHOX10 direct targets, as well as RHOX10 secondary-target genes. Many of the RHOX10-regulated genes encode proteins with known roles in SSCs. Using an in vitro ProSG differentiation assay, we find that RHOX10 promotes mouse ProSG differentiation through a conserved transcriptional cascade involving the key germ-cell TFs DMRT1 and ZBTB16. Our study gives important insights into germ cell development and provides a blueprint for how to define TF cascades.

Given the global impact and severity of COVID-19, there is a pressing need for a better understanding of the SARS-CoV-2 genome and mutations. Multi-strain sequence alignments of coronaviruses (CoV) provide important information for interpreting the genome and its variation. We apply a comparative genomics method, ConsHMM, to the multi-strain alignments of CoV to annotate every base of the SARS-CoV-2 genome with conservation states based on sequence alignment patterns among CoV. The learned conservation states show distinct enrichment patterns for genes, protein domains, and other regions of interest. Certain states are strongly enriched or depleted of SARS-CoV-2 mutations, which can be used to predict potentially consequential mutations. We expect the conservation states to be a resource for interpreting the SARS-CoV-2 genome and mutations.

Using the Cap Analysis of Gene Expression (CAGE) technology, the FANTOM5 consortium provided one of the most comprehensive maps of transcription start sites (TSSs) in several species. Strikingly, ~72% of them could not be assigned to a specific gene and initiate at unconventional regions, outside promoters or enhancers. Here, we probe these unassigned TSSs and show that, in all species studied, a significant fraction of CAGE peaks initiate at microsatellites, also called short tandem repeats (STRs). To confirm this transcription, we develop Cap Trap RNA-seq, a technology which combines cap trapping and long read MinION sequencing. We train sequence-based deep learning models able to predict CAGE signal at STRs with high accuracy. These models unveil the importance of STR surrounding sequences not only to distinguish STR classes, but also to predict the level of transcription initiation. Importantly, genetic variants linked to human diseases are preferentially found at STRs with high transcription initiation level, supporting the biological and clinical relevance of transcription initiation at STRs. Together, our results extend the repertoire of non-coding transcription associated with DNA tandem repeats and complexify STR polymorphism.

Dysregulated mTORC1 signaling alters a wide range of cellular processes, contributing to metabolic disorders and cancer. Defining the molecular details of downstream effectors is thus critical for uncovering selective therapeutic targets. We report that mTORC1 and its downstream kinase S6K enhance eIF4A/4B-mediated translation of Wilms' tumor 1-associated protein (WTAP), an adaptor for the N6-methyladenosine (m6A) RNA methyltransferase complex. This regulation is mediated by 5' UTR of WTAP mRNA that is targeted by eIF4A/4B. Single-nucleotide-resolution m6A mapping revealed that MAX dimerization protein 2 (MXD2) mRNA contains m6A, and increased m6A modification enhances its degradation. WTAP induces cMyc-MAX association by suppressing MXD2 expression, which promotes cMyc transcriptional activity and proliferation of mTORC1-activated cancer cells. These results elucidate a mechanism whereby mTORC1 stimulates oncogenic signaling via m6A RNA modification and illuminates the WTAP-MXD2-cMyc axis as a potential therapeutic target for mTORC1-driven cancers.

Most eukaryotes harbor two distinct pre-mRNA splicing machineries: the major spliceosome, which removes >99% of introns, and the minor spliceosome, which removes rare, evolutionarily conserved introns. Although hypothesized to serve important regulatory functions, physiologic roles of the minor spliceosome are not well understood. For example, the minor spliceosome component ZRSR2 is subject to recurrent, leukemia-associated mutations, yet functional connections among minor introns, hematopoiesis and cancers are unclear. Here, we identify that impaired minor intron excision via ZRSR2 loss enhances hematopoietic stem cell self-renewal. CRISPR screens mimicking nonsense-mediated decay of minor intron-containing mRNA species converged on LZTR1, a regulator of RAS-related GTPases. LZTR1 minor intron retention was also discovered in the RASopathy Noonan syndrome, due to intronic mutations disrupting splicing and diverse solid tumors. These data uncover minor intron recognition as a regulator of hematopoiesis, noncoding mutations within minor introns as potential cancer drivers and links among ZRSR2 mutations, LZTR1 regulation and leukemias.

Egg-laying mammals (monotremes) are the only extant mammalian outgroup to therians (marsupial and eutherian animals) and provide key insights into mammalian evolution1,2. Here we generate and analyse reference genomes of the platypus (Ornithorhynchus anatinus) and echidna (Tachyglossus aculeatus), which represent the only two extant monotreme lineages. The nearly complete platypus genome assembly has anchored almost the entire genome onto chromosomes, markedly improving the genome continuity and gene annotation. Together with our echidna sequence, the genomes of the two species allow us to detect the ancestral and lineage-specific genomic changes that shape both monotreme and mammalian evolution. We provide evidence that the monotreme sex chromosome complex originated from an ancestral chromosome ring configuration. The formation of such a unique chromosome complex may have been facilitated by the unusually extensive interactions between the multi-X and multi-Y chromosomes that are shared by the autosomal homologues in humans. Further comparative genomic analyses unravel marked differences between monotremes and therians in haptoglobin genes, lactation genes and chemosensory receptor genes for smell and taste that underlie the ecological adaptation of monotremes.

Ultraconserved enhancer sequences show perfect conservation between human and rodent genomes, suggesting that their functions are highly sensitive to mutation. However, current models of enhancer function do not sufficiently explain this extreme evolutionary constraint. We subjected 23 ultraconserved enhancers to different levels of mutagenesis, collectively introducing 1,547 mutations, and examined their activities in transgenic mouse reporter assays. Overall, we find that the regulatory properties of ultraconserved enhancers are robust to mutation. Upon mutagenesis, nearly all (19/23, 83%) still functioned as enhancers at one developmental stage, as did most of those tested again later in development (5/9, 56%). Replacement of endogenous enhancers with mutated alleles in mice corroborated results of transgenic assays, including the functional resilience of ultraconserved enhancers to mutation. Our findings show that the currently known activities of ultraconserved enhancers do not necessarily require the perfect conservation observed in evolution and suggest that additional regulatory or other functions contribute to their sequence constraint.

RNAcentral is a comprehensive database of non-coding RNA (ncRNA) sequences that provides a single access point to 44 RNA resources and >18 million ncRNA sequences from a wide range of organisms and RNA types. RNAcentral now also includes secondary (2D) structure information for >13 million sequences, making RNAcentral the world's largest RNA 2D structure database. The 2D diagrams are displayed using R2DT, a new 2D structure visualization method that uses consistent, reproducible and recognizable layouts for related RNAs. The sequence similarity search has been updated with a faster interface featuring facets for filtering search results by RNA type, organism, source database or any keyword. This sequence search tool is available as a reusable web component, and has been integrated into several RNAcentral member databases, including Rfam, miRBase and snoDB. To allow for a more fine-grained assignment of RNA types and subtypes, all RNAcentral sequences have been annotated with Sequence Ontology terms. The RNAcentral database continues to grow and provide a central data resource for the RNA community. RNAcentral is freely available at https://rnacentral.org.

Advances in sequencing technology have allowed researchers to sequence DNA with greater ease and at decreasing costs. Main developments have focused on either sequencing many short sequences or fewer large sequences. Methods for sequencing mid-sized sequences of 600-5,000 bp are currently less efficient. For example, the PacBio Sequel I system yields ~ 100,000-300,000 reads with an accuracy per base pair of 90-99%. We sought to sequence several DNA populations of ~ 870 bp in length with a sequencing accuracy of 99% and to the greatest depth possible. We optimised a simple, robust method to concatenate genes of ~ 870 bp five times and then sequenced the resulting DNA of ~ 5,000 bp by PacBioSMRT long-read sequencing. Our method improved upon previously published concatenation attempts, leading to a greater sequencing depth, high-quality reads and limited sample preparation at little expense. We applied this efficient concatenation protocol to sequence nine DNA populations from a protein engineering study. The improved method is accompanied by a simple and user-friendly analysis pipeline, DeCatCounter, to sequence medium-length sequences efficiently at one-fifth of the cost.

BACKGROUND: Duck viral enteritis, commonly known as duck plague (DP), is an acute and contagious fatal disease in ducks, geese, and swans caused by the DP virus (DPV). It poses a serious threat to the growth of duck farming in the Haor (wetland) areas of Bangladesh. AIM: This study aimed to detect the circulating DPV by molecular characterization, followed by phylogenetic analysis, targeting the UL30 gene in infected ducks from five Haor districts in Bangladesh and to observe the variation in the genome sequence between the field virus and vaccine strain of DPV. METHODS: A total of 150 samples (liver, 50; intestine, 50; and oropharyngeal tissue, 50) were collected from DP-suspected sick/dead ducks from 50 affected farms in Kishoreganj, Netrokona, B. Baria, Habiganj, and Sunamganj districts in Bangladesh. For the identification of DPV in collected samples, polymerase chain reaction (PCR) was utilized. Nucleotide sequences of the amplified UL30 gene were compared with those of other DPV strains available in GenBank. RESULTS: Of the 150 samples, 90 (60%) were found to be positive for DPV, as confirmed by PCR. Organ-wise prevalence was higher in the liver (72%), followed by the intestine (64%) and oropharyngeal tissue (44%). Regarding areas, the highest and lowest prevalence in the liver and intestine was observed in Habiganj and B. Baria, respectively, whereas the highest and lowest prevalence in the oropharyngeal tissue was observed in B. Baria and Habiganj, respectively. Two isolates, BAU/KA/DPV(B1)/2014 from Kishoreganj and BAU/KA/DPV(B4)/2014 from Sunamganj were sequenced, and phylogenetic analysis revealed that these isolates are evolutionarily closely related to Chinese isolates of DPV. Additionally, the isolates of DPV BAU/KA/DPV(B1)/2014 and BAU/KA/DPV(B4)/2014 showed the highest (98%) similarity to each other. The nucleotide sequence of the isolate BAU/KA/DPV(B1)/2014 exhibited higher nucleotide variability (246 nucleotides) than that of the vaccine strain (accession no. EU082088), which may affect protein function and additional drug sensitivity. CONCLUSION: Based on the findings of the molecular study, it can be assumed that the Bangladeshi isolates and all Chinese isolates of DPV may have a common ancestry.

The functional annotation of livestock genomes is crucial for understanding the molecular mechanisms that underpin complex traits of economic importance, adaptive evolution and comparative genomics. Here, we provide the most comprehensive catalogue to date of regulatory elements in the pig (Sus scrofa) by integrating 223 epigenomic and transcriptomic data sets, representing 14 biologically important tissues. We systematically describe the dynamic epigenetic landscape across tissues by functionally annotating 15 different chromatin states and defining their tissue-specific regulatory activities. We demonstrate that genomic variants associated with complex traits and adaptive evolution in pig are significantly enriched in active promoters and enhancers. Furthermore, we reveal distinct tissue-specific regulatory selection between Asian and European pig domestication processes. Compared with human and mouse epigenomes, we show that porcine regulatory elements are more conserved in DNA sequence, under both rapid and slow evolution, than those under neutral evolution across pig, mouse, and human. Finally, we provide biological insights on tissue-specific regulatory conservation, and by integrating 47 human genome-wide association studies, we demonstrate that, depending on the traits, mouse or pig might be more appropriate biomedical models for different complex traits and diseases.

RNAcentral is a comprehensive database of non-coding RNA (ncRNA) sequences that provides a single access point to 44 RNA resources and >18 million ncRNA sequences from a wide range of organisms and RNA types. RNAcentral now also includes secondary (2D) structure information for >13 million sequences, making RNAcentral the world's largest RNA 2D structure database. The 2D diagrams are displayed using R2DT, a new 2D structure visualization method that uses consistent, reproducible and recognizable layouts for related RNAs. The sequence similarity search has been updated with a faster interface featuring facets for filtering search results by RNA type, organism, source database or any keyword. This sequence search tool is available as a reusable web component, and has been integrated into several RNAcentral member databases, including Rfam, miRBase and snoDB. To allow for a more fine-grained assignment of RNA types and subtypes, all RNAcentral sequences have been annotated with Sequence Ontology terms. The RNAcentral database continues to grow and provide a central data resource for the RNA community. RNAcentral is freely available at https://rnacentral.org.

Axon injury is a hallmark of many neurodegenerative diseases, often resulting in neuronal cell death and functional impairment. Dual leucine zipper kinase (DLK) has emerged as a key mediator of this process. However, while DLK inhibition is robustly protective in a wide range of neurodegenerative disease models, it also inhibits axonal regeneration. Indeed, there are no genetic perturbations that are known to both improve long-term survival and promote regeneration. To identify such a neuroprotective target, we conducted a set of complementary high-throughput screens using a protein kinase inhibitor library in human stem cell-derived retinal ganglion cells (hRGCs). Overlapping compounds that promoted both neuroprotection and neurite outgrowth were bioinformatically deconvoluted to identify specific kinases that regulated neuronal death and axon regeneration. This work identified the role of germinal cell kinase four (GCK-IV) kinases in cell death and additionally revealed their unexpected activity in suppressing axon regeneration. Using an adeno-associated virus (AAV) approach, coupled with genome editing, we validated that GCK-IV kinase knockout improves neuronal survival, comparable to that of DLK knockout, while simultaneously promoting axon regeneration. Finally, we also found that GCK-IV kinase inhibition also prevented the attrition of RGCs in developing retinal organoid cultures without compromising axon outgrowth, addressing a major issue in the field of stem cell-derived retinas. Together, these results demonstrate a role for the GCK-IV kinases in dissociating the cell death and axonal outgrowth in neurons and their druggability provides for therapeutic options for neurodegenerative diseases.

Fluid intake is an essential innate behaviour that is mainly caused by two distinct types of thirst1-3. Increased blood osmolality induces osmotic thirst that drives animals to consume pure water. Conversely, the loss of body fluid induces hypovolaemic thirst, in which animals seek both water and minerals (salts) to recover blood volume. Circumventricular organs in the lamina terminalis are critical sites for sensing both types of thirst-inducing stimulus4-6. However, how different thirst modalities are encoded in the brain remains unknown. Here we employed stimulus-to-cell-type mapping using single-cell RNA sequencing to identify the cellular substrates that underlie distinct types of thirst. These studies revealed diverse types of excitatory and inhibitory neuron in each circumventricular organ structure. We show that unique combinations of these neuron types are activated under osmotic and hypovolaemic stresses. These results elucidate the cellular logic that underlies distinct thirst modalities. Furthermore, optogenetic gain of function in thirst-modality-specific cell types recapitulated water-specific and non-specific fluid appetite caused by the two distinct dipsogenic stimuli. Together, these results show that thirst is a multimodal physiological state, and that different thirst states are mediated by specific neuron types in the mammalian brain.

Annotations of evolutionary sequence constraint based on multi-species genome alignments and genome-wide maps of epigenomic marks and transcription factor binding provide important complementary information for understanding the human genome and genetic variation. Here we developed the Constrained Non-Exonic Predictor (CNEP) to quantify the evidence of each base in the genome being in an evolutionarily constrained non-exonic element from an input of over 60,000 epigenomic and transcription factor binding features. We find that the CNEP score outperforms baseline and related existing scores at predicting evolutionarily constrained non-exonic bases from such data. However, a subset of them are still not well predicted by CNEP. We developed a complementary Conservation Signature Score by CNEP (CSS-CNEP) that is predictive of those bases. We further characterize the nature of constrained non-exonic bases with low CNEP scores using additional types of information. CNEP and CSS-CNEP are resources for analyzing constrained non-exonic bases in the genome.

Genetic sequencing, or DNA sequencing, using the Sanger technique has become widely used in the veterinary diagnostic community. This technology plays a role in verification of PCR results and is used to provide the genetic sequence data needed for phylogenetic analysis, epidemiologic studies, and forensic investigations. The Laboratory Technology Committee of the American Association of Veterinary Laboratory Diagnosticians has prepared guidelines for sample preparation, submission to sequencing facilities or instrumentation, quality assessment of nucleic acid sequence data performed, and for generating basic sequencing data and phylogenetic analysis for diagnostic applications. This guidance is aimed at assisting laboratories in providing consistent, high-quality, and reliable sequence data when using Sanger-based genetic sequencing as a component of their laboratory services.

The mammary epithelial cell (MEC) system is a bilayered ductal epithelium of luminal and basal cells, maintained by a lineage of stem and progenitor populations. Here, we used integrated single-cell transcriptomics and chromatin accessibility analysis to reconstruct the cell types of the mouse MEC system and their underlying gene regulatory features in an unbiased manner. We define differentiation states within the secretory type of luminal cells, which forms a continuous spectrum of general luminal progenitor and lactation-committed progenitor cells. By integrating single-cell transcriptomics and chromatin accessibility landscapes, we identify cis- and trans-regulatory elements that are differentially activated in the specific epithelial cell types and our newly defined luminal differentiation states. Our work provides a resource to reveal cis/trans-regulatory elements associated with MEC identity and differentiation that will serve as a reference to determine how the chromatin accessibility landscape changes during breast cancer.

Blood vascular endothelial cells (BECs) control the immune response by regulating blood flow and immune cell recruitment in lymphoid tissues. However, the diversity of BEC and their origins during immune angiogenesis remain unclear. Here we profile transcriptomes of BEC from peripheral lymph nodes and map phenotypes to the vasculature. We identify multiple subsets, including a medullary venous population whose gene signature predicts a selective role in myeloid cell (vs lymphocyte) recruitment to the medulla, confirmed by videomicroscopy. We define five capillary subsets, including a capillary resident precursor (CRP) that displays stem cell and migratory gene signatures, and contributes to homeostatic BEC turnover and to neogenesis of high endothelium after immunization. Cell alignments show retention of developmental programs along trajectories from CRP to mature venous and arterial populations. Our single cell atlas provides a molecular roadmap of the lymph node blood vasculature and defines subset specialization for leukocyte recruitment and vascular homeostasis.

Recognition of highly degenerate mammalian splice sites by the core spliceosomal machinery is regulated by several protein factors that predominantly bind exonic splicing motifs. These are postulated to be single-stranded in order to be functional, yet knowledge of secondary structural features that regulate the exposure of exonic splicing motifs across the transcriptome is not currently available. Using transcriptome-wide RNA structural information we show that retained introns in mouse are commonly flanked by a short (≲70 nucleotide), highly base-paired segment upstream and a predominantly single-stranded exonic segment downstream. Splicing assays with select pre-mRNA substrates demonstrate that loops immediately upstream of the introns contain pre-mRNA-specific splicing enhancers, the substitution or hybridization of which impedes splicing. Additionally, the exonic segments flanking the retained introns appeared to be more enriched in a previously identified set of hexameric exonic splicing enhancer (ESE) sequences compared to their spliced counterparts, suggesting that base-pairing in the exonic segments upstream of retained introns could be a means for occlusion of ESEs. The upstream exonic loops of the test substrate promoted recruitment of splicing factors and consequent pre-mRNA structural remodeling, leading up to assembly of the early spliceosome. These results suggest that disruption of exonic stem-loop structures immediately upstream (but not downstream) of the introns regulate alternative splicing events, likely through modulating accessibility of splicing factors.

Pediatric cataract is clinically and genetically heterogeneous and is the most common cause of childhood blindness worldwide. In this study, we aimed to identify disease-causing variants in three large British families and one isolated case with autosomal dominant congenital cataract, using whole exome sequencing. We identified four different heterozygous variants, three in the large families and one in the isolated case. Family A, with a novel missense variant (c.178G>C, p.Gly60Arg) in GJA8 with lamellar cataract; family B, with a recurrent variant in GJA8 (c.262C>T, p.Pro88Ser) associated with nuclear cataract; and family C, with a novel variant in GJA3 (c.771dupC, p.Ser258GlnfsTer68) causing a lamellar phenotype. Individual D had a novel variant in GJA3 (c.82G>T, p.Val28Leu) associated with congenital cataract. Each sequence variant was found to co-segregate with disease. Here, we report three novel and one recurrent disease-causing sequence variant in the gap junctional protein encoding genes causing autosomal dominant congenital cataract. Our study further extends the mutation spectrum of these genes and further facilitates clinical diagnosis. A recurrent p.P88S variant in GJA8 causing isolated nuclear cataract provides evidence of further phenotypic heterogeneity associated with this variant.

Developmental progression depends on temporally defined changes in gene expression mediated by transient exposure of lineage intermediates to signals in the progenitor niche. To determine whether cell-intrinsic epigenetic mechanisms contribute to signal-induced transcriptional responses, here we manipulate the signalling environment and activity of the histone demethylase LSD1 during differentiation of hESC-gut tube intermediates into pancreatic endocrine cells. We identify a transient requirement for LSD1 in endocrine cell differentiation spanning a short time-window early in pancreas development, a phenotype we reproduced in mice. Examination of enhancer and transcriptome landscapes revealed that LSD1 silences transiently active retinoic acid (RA)-induced enhancers and their target genes. Furthermore, prolonged RA exposure phenocopies LSD1 inhibition, suggesting that LSD1 regulates endocrine cell differentiation by limiting the duration of RA signalling. Our findings identify LSD1-mediated enhancer silencing as a cell-intrinsic epigenetic feedback mechanism by which the duration of the transcriptional response to a developmental signal is limited.

Mucopolysaccharidosis (MPS) is a metabolic storage disorder caused by the deficiency of any lysosomal enzyme required for the breakdown of glycosaminoglycans. A 15-month-old Boston Terrier presented with clinical signs consistent with lysosomal storage disease including corneal opacities, multifocal central nervous system disease and progressively worsening clinical course. Diagnosis was confirmed at necropsy based on histopathologic evaluation of multiple organs demonstrating accumulation of mucopolysaccharides. Whole genome sequencing was used to uncover a frame-shift insertion affecting the alpha-L-iduronidase (IDUA) gene (c.19_20insCGGCCCCC), a mutation confirmed in another Boston Terrier presented 2 years later with a similar clinical picture. Both dogs were homozygous for the IDUA mutation and shared coat colors not recognized as normal for the breed by the American Kennel Club. In contrast, the mutation was not detected in 120 unrelated Boston Terriers as well as 202 dogs from other breeds. Recent inbreeding to select for recessive and unusual coat colors may have concentrated this relatively rare allele in the breed. The identification of the variant enables ante-mortem diagnosis of similar cases and selective breeding to avoid the spread of this disease in the breed. Boston Terriers carrying this variant represent a promising model for MPS I with neurological abnormalities in humans.

Higher expression of the human DNA repair enzyme MUTYH has previously been shown to be strongly associated with reduced survival in a panel of 24 human lymphoblastoid cell lines exposed to the alkylating agent N-methyl-N'-nitro-N-nitrosoguanidine (MNNG). The molecular mechanism of MUTYH-enhanced MNNG cytotoxicity is unclear, because MUTYH has a well-established role in the repair of oxidative DNA lesions. Here, we show in mouse embryonic fibroblasts (MEFs) that this MNNG-dependent phenotype does not involve oxidative DNA damage and occurs independently of both O6-methyl guanine adduct cytotoxicity and MUTYH-dependent glycosylase activity. We found that blocking of abasic (AP) sites abolishes higher survival of Mutyh-deficient (Mutyh-/-) MEFs, but this blockade had no additive cytotoxicity in WT MEFs, suggesting the cytotoxicity is due to MUTYH interactions with MNNG-induced AP sites. We found that recombinant mouse MUTYH tightly binds AP sites opposite all four canonical undamaged bases and stimulated apurinic/apyrimidinic endonuclease 1 (APE1)-mediated DNA incision. Consistent with these observations, we found that stable expression of WT, but not catalytically-inactive MUTYH, enhances MNNG cytotoxicity in Mutyh-/- MEFs and that MUTYH expression enhances MNNG-induced genomic strand breaks. Taken together, these results suggest that MUTYH enhances the rapid accumulation of AP-site intermediates by interacting with APE1, implicating MUTYH as a factor that modulates the delicate process of base-excision repair independently of its glycosylase activity.

A single enzyme active site that catalyzes multiple reactions is a well-established biochemical theme, but how one nuclease site cleaves both DNA strands of a double helix has not been well understood. In analyzing site-specific DNA cleavage by the mammalian RAG1-RAG2 recombinase, which initiates V(D)J recombination, we find that the active site is reconfigured for the two consecutive reactions and the DNA double helix adopts drastically different structures. For initial nicking of the DNA, a locally unwound and unpaired DNA duplex forms a zipper via alternating interstrand base stacking, rather than melting as generally thought. The second strand cleavage and formation of a hairpin-DNA product requires a global scissor-like movement of protein and DNA, delivering the scissile phosphate into the rearranged active site.

Genome editing followed by reproductive cloning was previously used to produce two hornless dairy bulls. We crossed one genome-edited dairy bull, homozygous for the dominant PC Celtic POLLED allele, with horned cows (pp) and obtained six heterozygous (PCp) polled calves. The calves had no horns and were otherwise healthy and phenotypically unremarkable. We conducted whole-genome sequencing of all animals using an Illumina HiSeq4000 to achieve ~20× coverage. Bioinformatics analyses revealed the bull was a compound heterozygote, carrying one naturally occurring PC Celtic POLLED allele and an allele containing an additional introgression of the homology-directed repair donor plasmid along with the PC Celtic allele. These alleles segregated in the offspring of this bull, and inheritance of either allele produced polled calves. No other unintended genomic alterations were observed. These data can be used to inform conversations in the scientific community, with regulatory authorities and with the public around 'intentional genomic alterations' and future regulatory actions regarding genome-edited animals.

Nuclear factor κB (NF-κB) RelA is the potent transcriptional activator of inflammatory response genes. We stringently defined a list of direct RelA target genes by integrating physical (chromatin immunoprecipitation sequencing [ChIP-seq]) and functional (RNA sequencing [RNA-seq] in knockouts) datasets. We then dissected each gene's regulatory strategy by testing RelA variants in a primary-cell genetic-complementation assay. All endogenous target genes require RelA to make DNA-base-specific contacts, and none are activatable by the DNA binding domain alone. However, endogenous target genes differ widely in how they employ the two transactivation domains. Through model-aided analysis of the dynamic time-course data, we reveal the gene-specific synergy and redundancy of TA1 and TA2. Given that post-translational modifications control TA1 activity and intrinsic affinity for coactivators determines TA2 activity, the differential TA logics suggests context-dependent versus context-independent control of endogenous RelA-target genes. Although some inflammatory initiators appear to require co-stimulatory TA1 activation, inflammatory resolvers are a part of the NF-κB RelA core response.

DNA from formalin-preserved tissue could unlock a vast repository of genetic information stored in museums worldwide. However, formaldehyde crosslinks proteins and DNA, and prevents ready amplification and DNA sequencing. Formaldehyde acylation also fragments the DNA. Treatment with proteinase K proteolyzes crosslinked proteins to rescue the DNA, though the process is quite slow. To reduce processing time and improve rescue efficiency, we applied the mechanical energy of a vortex fluidic device (VFD) to drive the catalytic activity of proteinase K and recover DNA from American lobster tissue (Homarus americanus) fixed in 3.7% formalin for >1-year. A scan of VFD rotational speeds identified the optimal rotational speed for recovery of PCR-amplifiable DNA and while 500+ base pairs were sequenced, shorter read lengths were more consistently obtained. This VFD-based method also effectively recovered DNA from formalin-preserved samples. The results provide a roadmap for exploring DNA from millions of historical and even extinct species.

CRISPR-based gene drives can spread through wild populations by biasing their own transmission above the 50% value predicted by Mendelian inheritance. These technologies offer population-engineering solutions for combating vector-borne diseases, managing crop pests, and supporting ecosystem conservation efforts. Current technologies raise safety concerns for unintended gene propagation. Herein, we address such concerns by splitting the drive components, Cas9 and gRNAs, into separate alleles to form a trans-complementing split-gene-drive (tGD) and demonstrate its ability to promote super-Mendelian inheritance of the separate transgenes. This dual-component configuration allows for combinatorial transgene optimization and increases safety by restricting escape concerns to experimentation windows. We employ the tGD and a small-molecule-controlled version to investigate the biology of component inheritance and resistant allele formation, and to study the effects of maternal inheritance and impaired homology on efficiency. Lastly, mathematical modeling of tGD spread within populations reveals potential advantages for improving current gene-drive technologies for field population modification.

Sequences encoding Olduvai protein domains (formerly DUF1220) show the greatest human lineage-specific increase in copy number of any coding region in the genome and have been associated, in a dosage-dependent manner, with brain size, cognitive aptitude, autism, and schizophrenia. Tandem intragenic duplications of a three-domain block, termed the Olduvai triplet, in four NBPF genes in the chromosomal 1q21.1-0.2 region, are primarily responsible for the striking human-specific copy number increase. Interestingly, most of the Olduvai triplets are adjacent to, and transcriptionally coregulated with, three human-specific NOTCH2NL genes that have been shown to promote cortical neurogenesis. Until now, the underlying genomic events that drove the Olduvai hyperamplification in humans have remained unexplained. Here, we show that the presence or absence of an alternative first exon of the Olduvai triplet perfectly discriminates between amplified (58/58) and unamplified (0/12) triplets. We provide sequence and breakpoint analyses that suggest the alternative exon was produced by an nonallelic homologous recombination-based mechanism involving the duplicative transposition of an existing Olduvai exon found in the CON3 domain, which typically occurs at the C-terminal end of NBPF genes. We also provide suggestive in vitro evidence that the alternative exon may promote instability through a putative G-quadraplex (pG4)-based mechanism. Lastly, we use single-molecule optical mapping to characterize the intragenic structural variation observed in NBPF genes in 154 unrelated individuals and 52 related individuals from 16 families and show that the presence of pG4-containing Olduvai triplets is strongly correlated with high levels of Olduvai copy number variation. These results suggest that the same driver of genomic instability that allowed the evolutionarily recent, rapid, and extreme human-specific Olduvai expansion remains highly active in the human genome.

Accurate methods for tracking individuals are crucial to the success of fisheries and aquaculture management. Management of migratory salmonid populations, which are important for the health of many economies, ecosystems, and indigenous cultures, is particularly dependent on data gathered from tagged fish. However, the physical tagging methods currently used have many challenges including cost, variable marker retention, and information limited to tagged individuals. Genetic tracking methods combat many of the problems associated with physical tags, but have their own challenges including high cost, potentially difficult marker design, and incompatibility of markers across species. Here we show the feasibility of a new genotyping method for parent-based tagging (PBT), where individuals are tracked through the inherent genetic relationships with their parents. We found that Rapture sequencing, a combination of restriction-site associated DNA and capture sequencing, provides sufficient data for parentage assignment. Additionally, the same capture bait set, which targets specific restriction-site associated DNA loci, can be used for both Rainbow Trout Oncorhynchus mykiss and Chinook Salmon Oncorhynchus tshawytscha. We input 248 single nucleotide polymorphisms from 1,121 samples to parentage assignment software and compared parent-offspring relationships of the spawning pairs recorded in a hatchery. Interestingly, our results suggest sperm contamination during hatchery spawning occurred in the production of 14% of offspring, further confirming the need for genetic tagging in accurately tracking individuals. PBT with Rapture successfully assigned progeny to parents with a 98.86% accuracy with sufficient genetic data. Cost for this pilot study was approximately $3 USD per sample. As costs vary based on the number of markers used and individuals sequenced, we expect that when implemented at a large-scale, per sample costs could be further decreased. We conclude that Rapture PBT provides a cost-effective and accurate alternative to the physical coded wire tags, and other genetic-based methods.