Genomic analysis and characterization of new loci associated with seed protein and oil content in soybeans.

Breeding for increased protein without a reduction in oil content in soybeans [Glycine max (L.) Merr.] is a challenge for soybean breeders but an expected goal. Many efforts have been made to develop new soybean varieties with high yield in combination with desirable protein and/or oil traits. An elite line, R05-1415, was reported to be high yielding, high protein, and low oil. Several significant quantitative trait loci (QTL) for protein and oil were reported in this line, but many of them were unstable across environments or genetic backgrounds. Thus, a new study under multiple field environments using the Infinium BARCSoySNP6K BeadChips was conducted to detect and confirm stable genomic loci for these traits. Genetic analyses consistently detected a single major genomic locus conveying these two traits with remarkably high phenotypic variation explained (R2 ), varying between 24.2% and 43.5%. This new genomic locus is located between 25.0 and 26.7 Mb, distant from the previously reported QTL and did not overlap with other commonly reported QTL and the recently cloned gene Glyma.20G085100. Homolog analysis indicated that this QTL did not result from the paracentric chromosome inversion with an adjacent genomic fragment that harbors the reported QTL. The pleiotropic effect of this QTL could be a challenge for improving protein and oil simultaneously; however, a further study of four candidate genes with significant expressions in the seed developmental stages coupled with haplotype analysis may be able to pinpoint causative genes. The functionality and roles of these genes can be determined and characterized, which lay a solid foundation for the improvement of protein and oil content in soybeans.

Integrated network findings reveal ubiquitin-specific protease 44 overexpression suppresses tumorigenicity of liver cancer.

Hepatocellular carcinoma (HCC) is the sixth most common cancer and third leading cause of cancer-related deaths worldwide. HCC is a multistep disease marked by various signaling alterations. A better understanding of the new molecular drivers of HCC could therefore provide an opportunity to develop effective diagnostic and therapeutic targets. Ubiquitin-specific protease 44 (USP44), a member of the cysteine protease family, has been reported to play a role in many cancer types. However, its contribution to HCC development remains unknown. In the present study, we observed suppression of USP44 expression in HCC tissue. Clinicopathologic analysis further showed that low USP44 expression correlated with poorer survival and a later tumor stage in HCC, suggesting that USP44 could be a predictor of poor prognosis in patients with HCC. Gain-of-function analysis in vitro demonstrated the importance of USP44 in HCC cell growth and G0/G1 cell cycle arrest. To investigate the downstream targets of USP44 and the molecular mechanisms underlying its regulation of cell proliferation in HCC, we conducted a comparative transcriptomic analysis and identified a cluster of proliferation-related genes, including CCND2, CCNG2, and SMC3. Ingenuity Pathway Analysis further delineated the gene networks controlled by USP44 through the regulation of membrane proteins and receptors, enzymes, transcriptional factors, and cyclins involved in the control of cell proliferation, metastasis, and apoptosis in HCC. To summarize, our results highlight, for the first time, the tumor-suppression role of USP44 in HCC and suggest a new prognostic biomarker in this disease.

Cell-free circulating tumor RNAs in plasma as the potential prognostic biomarkers in colorectal cancer.

Background: Cell free RNA (cfRNA) contains transcript fragments from multiple cell types, making it useful for cancer detection in clinical settings. However, the pathophysiological origins of cfRNAs in plasma from colorectal cancer (CRC) patients remain unclear. Methods: To identify the tissue-specific contributions of cfRNAs transcriptomic profile, we used a published single-cell transcriptomics profile to deconvolute cell type abundance among paired plasma samples from CRC patients who underwent tumor-ablative surgery. We further validated the differentially expressed cfRNAs in 5 pairs of CRC tumor samples and adjacent tissue samples as well as 3 additional CRC tumor samples using RNA-sequencing. Results: The transcriptomic component from intestinal secretory cells was significantly decreased in the in-house post-surgical cfRNA. The HPGD, PACS1, and TDP2 expression was consistent across cfRNA and tissue samples. Using the Cancer Genome Atlas (TCGA) CRC datasets, we were able to classify the patients into two groups with significantly different survival outcomes. Conclusions: The three-gene signature holds promise in applying minimal residual disease (MRD) testing, which involves profiling remnants of cancer cells after or during treatment. Biomarkers identified in the present study need to be validated in a larger cohort of samples in order to ascertain their possible use in early diagnosis of CRC.

SOX4 is a novel phenotypic regulator of endothelial cells in atherosclerosis revealed by single-cell analysis.

INTRODUCTION: Atherosclerotic complications represent the leading cause of cardiovascular mortality globally. Dysfunction of endothelial cells (ECs) often initiates the pathological events in atherosclerosis. OBJECTIVES: In this study, we sought to investigate the transcriptional profile of atherosclerotic aortae, identify novel regulator in dysfunctional ECs and hence provide mechanistic insights into atherosclerotic progression. METHODS: We applied single-cell RNA sequencing (scRNA-seq) on aortic cells from Western diet-fed apolipoprotein E-deficient (ApoE-/-) mice to explore the transcriptional landscape and heterogeneity of dysfunctional ECs. In vivo validation of SOX4 upregulation in ECs were performed in atherosclerotic tissues, including mouse aortic tissues, human coronary arteries, and human renal arteries. Single-cell analysis on human aortic aneurysmal tissue was also performed. Downstream vascular abnormalities induced by EC-specific SOX4 overexpression, and upstream modulators of SOX4 were revealed by biochemical assays, immunostaining, and wire myography. Effects of shear stress on endothelial SOX4 expression was investigated by in vitro hemodynamic study. RESULTS: Among the compendium of aortic cells, mesenchymal markers in ECs were significantly enriched. Two EC subsets were subsequently distinguished, as the 'endothelial-like' and 'mesenchymal-like' subsets. Conventional assays consistently identified SOX4 as a novel atherosclerotic marker in mouse and different human arteries, additional to a cancer marker. EC-specific SOX4 overexpression promoted atherogenesis and endothelial-to-mesenchymal transition (EndoMT). Importantly, hyperlipidemia-associated cytokines and oscillatory blood flow upregulated, whereas the anti-diabetic drug metformin pharmacologically suppressed SOX4 level in ECs. CONCLUSION: Our study unravels SOX4 as a novel phenotypic regulator during endothelial dysfunction, which exacerbates atherogenesis. Our study also pinpoints hyperlipidemia-associated cytokines and oscillatory blood flow as endogenous SOX4 inducers, providing more therapeutic insights against atherosclerotic diseases.

Chromatin modifiers: A new class of pollutants with potential epigenetic effects revealed by in vitro assays and transcriptomic analyses.

A great variety of endocrine-disrupting chemicals (EDCs) have been used extensively and become widespread in the environment nowadays. Limited mammalian studies have shown that certain EDCs may target chromosome and epigenome of the germline, leading to adverse effects in subsequent generations, despite these progenies having never been exposed to the EDC before. However, the underlying mechanisms of chromosomal changes induced by these pollutants remain poorly known. Using the human ovarian granulosa tumor cell line COV434 as a model, we investigated and compared the transcriptomic changes induced by nine EDCs with diverse chemical structures (i.e. BDE-47, BPA, BP-3, DEHP, DHP, EE2, TCS, TDCPP and NP), to inquire if there is any common epigenetic modification associated with reproductive functions induced by these EDCs. Our results showed that COV434 cells were more responsive to BP-3, NP, DEHP and EE2, and more importantly, these four EDCs altered the expression of gene clusters related to DNA damage response, cell cycle, proliferation, and chromatin remodeling, which can potentially lead to epigenetic modifications and transgenerational inheritance. Furthermore, dysregulation of similar gene clusters was common in DEHP and NP treatments. Bioinformatics analysis further revealed that BP-3 disturbed signaling pathways associated with reproductive functions, whereas alterations in telomere-related pathways were highlighted upon EE2 exposure. Overall, this study highlighted chromatin modifications caused by a class of chemicals which that may potentially lead to epigenetic changes and transgenerational reproductive impairments.

LAFITE Reveals the Complexity of Transcript Isoforms in Subcellular Fractions.

Characterization of the subcellular distribution of RNA is essential for understanding the molecular basis of biological processes. Here, the subcellular nanopore direct RNA-sequencing (DRS) of four lung cancer cell lines (A549, H1975, H358, and HCC4006) is performed, coupled with a computational pipeline, Low-abundance Aware Full-length Isoform clusTEr (LAFITE), to comprehensively analyze the full-length cytoplasmic and nuclear transcriptome. Using additional DRS and orthogonal data sets, it is shown that LAFITE outperforms current methods for detecting full-length transcripts, particularly for low-abundance isoforms that are usually overlooked due to poor read coverage. Experimental validation of six novel isoforms exclusively identified by LAFITE further confirms the reliability of this pipeline. By applying LAFITE to subcellular DRS data, the complexity of the nuclear transcriptome is revealed in terms of isoform diversity, 3'-UTR usage, m6A modification patterns, and intron retention. Overall, LAFITE provides enhanced full-length isoform identification and enables a high-resolution view of the RNA landscape at the isoform level.

Accurate identification of structural variations from cancer samples.

Structural variations (SVs) are commonly found in cancer genomes. They can cause gene amplification, deletion and fusion, among other functional consequences. With an average read length of hundreds of kilobases, nano-channel-based optical DNA mapping is powerful in detecting large SVs. However, existing SV calling methods are not tailored for cancer samples, which have special properties such as mixed cell types and sub-clones. Here we propose the Cancer Optical Mapping for detecting Structural Variations (COMSV) method that is specifically designed for cancer samples. It shows high sensitivity and specificity in benchmark comparisons. Applying to cancer cell lines and patient samples, COMSV identifies hundreds of novel SVs per sample.

Transcriptomic Analysis in Marine Medaka Gill Reveals That the Hypo-Osmotic Stress Could Alter the Immune Response via the IL17 Signaling Pathway.

Fish gills are the major osmoregulatory tissue that contact the external water environment and have developed an effective osmoregulatory mechanism to maintain cellular function. Marine medaka (Oryzias melastigma) has the ability to live in both seawater and fresh water environments. The present study performed a seawater (SW) to 50% seawater (SFW) transfer, and the gill samples were used for comparative transcriptomic analysis to study the alteration of hypo-osmotic stress on immune responsive genes in this model organism. The result identified 518 differentiated expressed genes (DEGs) after the SW to SFW transfer. Various pathways such as p53 signaling, forkhead box O signaling, and the cell cycle were enriched. Moreover, the immune system was highlighted as one of the top altered biological processes in the enrichment analysis. Various cytokines, chemokines, and inflammatory genes that participate in the IL-17 signaling pathway were suppressed after the SW to SFW transfer. On the other hand, some immunoglobulin-related genes were up-regulated. The results were further validated by real-time qPCR. Taken together, our study provides additional gill transcriptome information in marine medaka; it also supports the notion that osmotic stress could influence the immune responses in fish gills.

An RNA G-Quadruplex Structure within the ADAR 5'UTR Interacts with DHX36 Helicase to Regulate Translation.

RNA G-quadruplex (rG4) structures in the 5' untranslated region (5'UTR) play crucial roles in fundamental cellular processes. ADAR is an important enzyme that binds to double-strand RNA and accounts for the conversion of Adenosine to Inosine in RNA editing. However, so far there is no report on the formation and regulatory role of rG4 on ADAR expression. Here, we identify and characterize a thermostable rG4 structure within the 5'UTR of the ADAR1 mRNA and demonstrate its formation and inhibitory role on translation in reporter gene and native gene constructs. We reveal rG4-specific helicase DHX36 interacts with this rG4 in vitro and in cells under knockdown and knockout conditions by GTFH (G-quadruplex-triggered fluorogenic hybridization) probes and modulates translation in an rG4-dependent manner. Our results further substantiate the rG4 structure-DHX36 protein interaction in cells and highlight rG4 to be a key player in controlling ADAR1 translation.

Corrigendum: Low-Dose Radiation Can Cause Epigenetic Alterations Associated With Impairments in Both Male and Female Reproductive Cell.

[This corrects the article DOI: 10.3389/fgene.2021.710143.].

Endothelium-targeted delivery of PPARδ by adeno-associated virus serotype 1 ameliorates vascular injury induced by hindlimb ischemia in obese mice.

Diabetic vasculopathy is a major health problem worldwide. Peripheral arterial disease (PAD), and in its severe form, critical limb ischemia is a major form of diabetic vasculopathy with limited treatment options. Existing literature suggested an important role of PPARδ in vascular homeostasis. It remains elusive for using PPARδ as a potential therapeutic target due to mostly the side effects of PPARδ agonists. To explore the roles of PPARδ in endothelial homeostasis, endothelial cell (EC) selective Ppard knockout and controlled mice were subjected to hindlimb ischemia (HLI) injury. The muscle ECs were sorted for single-cell RNA sequencing (scRNA-seq) analysis. HLI was also performed in high fat diet (HFD)-induced obese mice to examine the function of PPARδ in obese mice with delayed vascular repair. Adeno-associated virus type 1 (AAV1) carrying ICAM2 promoter to target endothelium for overexpressing PPARδ was injected into the injured muscles of normal chow- and HFD-fed obese mice before HLI surgery was performed. scRNA-seq analysis of ECs in ischemic muscles revealed a pivotal role of PPARδ in endothelial homeostasis. PPARδ expression was diminished both after HLI injury, and also in obese mice, which showed further delayed vascular repair. Endothelium-targeted delivery of PPARδ by AAV1 improved perfusion recovery, increased capillary density, restored endothelial integrity, suppressed vascular inflammation, and promoted muscle regeneration in ischemic hindlimbs of both lean and obese mice. Our study indicated the effectiveness of endothelium-targeted PPARδ overexpression for restoring functional vasculature after ischemic injury, which might be a promising option of gene therapy to treat PAD and CLI.

Network Pharmacology and Comparative Transcriptome Reveals Biotargets and Mechanisms of Curcumol Treating Lung Adenocarcinoma Patients With COVID-19.

The coronavirus disease 2019 (COVID-19) pandemic has led to 4,255,892 deaths worldwide. Although COVID-19 vaccines are available, mutant forms of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have reduced the effectiveness of vaccines. Patients with cancer are more vulnerable to COVID-19 than patients without cancer. Identification of new drugs to treat COVID-19 could reduce mortality rate, and traditional Chinese Medicine(TCM) has shown potential in COVID-19 treatment. In this study, we focused on lung adenocarcinoma (LUAD) patients with COVID-19. We aimed to investigate the use of curcumol, a TCM, to treat LUAD patients with COVID-19, using network pharmacology and systematic bioinformatics analysis. The results showed that LUAD and patients with COVID-19 share a cluster of common deregulated targets. The network pharmacology analysis identified seven core targets (namely, AURKA, CDK1, CCNB1, CCNB2, CCNE1, CCNE2, and TTK) of curcumol in patients with COVID-19 and LUAD. Clinicopathological analysis of these targets demonstrated that the expression of these targets is associated with poor patient survival rates. The bioinformatics analysis further highlighted the involvement of this target cluster in DNA damage response, chromosome stability, and pathogenesis of LUAD. More importantly, these targets influence cell-signaling associated with the Warburg effect, which supports SARS-CoV-2 replication and inflammatory response. Comparative transcriptomic analysis on in vitro LUAD cell further validated the effect of curcumol for treating LUAD through the control of cell cycle and DNA damage response. This study supports the earlier findings that curcumol is a potential treatment for patients with LUAD and COVID-19.

A tRNA-derived fragment from Chinese yew suppresses ovarian cancer growth via targeting TRPA1.

Drug discovery from plants usually focuses on small molecules rather than such biological macromolecules as RNAs. Although plant transfer RNA (tRNA)-derived fragment (tRF) has been associated with the developmental and defense mechanisms in plants, its regulatory role in mammals remains unclear. By employing a novel reverse small interfering RNA (siRNA) screening strategy, we show that a tRF mimic (antisense derived from the 5' end of tRNAHis(GUG) of Chinese yew) exhibits comparable anti-cancer activity with that of taxol on ovarian cancer A2780 cells, with a 16-fold lower dosage than that of taxol. A dual-luciferase reporter assay revealed that tRF-T11 directly targets the 3' UTR of oncogene TRPA1 mRNA. Furthermore, an Argonaute-RNA immunoprecipitation (AGO-RIP) assay demonstrated that tRF-T11 can interact with AGO2 to suppress TRPA1 via an RNAi pathway. This study uncovers a new role of plant-derived tRFs in regulating endogenous genes. This holds great promise for exploiting novel RNA drugs derived from nature and sheds light on the discovery of unknown molecular targets of therapeutics.

Type 2 innate immunity drives distinct neonatal immune profile conducive for heart regeneration.

Aims: Neonatal immunity is functionally immature and skewed towards a TH2-driven, anti-inflammatory profile. This neonatal immunotolerance is partly driven by the type 2 cytokines: interleukin-4 (IL-4) and interleukin-13 (IL-13). Studies on neonatal cardiac regeneration reveal the beneficial role of an anti-inflammatory response in restoring cardiac function after injury. However, the role of an imbalanced immune repertoire observed in neonates on tissue regeneration is poorly understood; specifically, whether IL-4 and IL-13 actively modulate neonatal immunity during cardiac injury. Methods and results: Neonatal mice lacking IL-4 and IL-13 (DKOs) examined at 2 days after birth exhibited reduced anti-inflammatory immune populations with basal cardiac immune populations like adult mice. Examination of neonates lacking IL-4 and IL-13 at 2 days post cardiac ischemic injury, induced on the second day after birth, showed impaired cardiac function compared to their control counterparts. Treatment with either IL-4 or IL-13 cytokine during injury restored both cardiac function and immune population profiles in knockout mice. Examination of IL-4/IL-13 downstream pathways revealed the role of STAT6 in mediating the regenerative response in neonatal hearts. As IL-4/IL-13 drives polarization of alternatively activated macrophages, we also examined the role of IL-4/IL-13 signaling within the myeloid compartment during neonatal cardiac regeneration. Injury of IL-4Rα myeloid specific knockout neonates 2 days after birth revealed that loss of IL-4/IL-13 signaling in macrophages alone was sufficient to impair cardiac regeneration. Conclusions: Our results confirm that the TH2 cytokines: IL-4 and IL-13, which skews neonatal immunity to a TH2 profile, are necessary for maintaining and mediating an anti-inflammatory response in the neonatal heart, in part through the activation of alternatively activated macrophages, thereby permitting a niche conducive for regeneration.

Detrimental effects of microplastic exposure on normal and asthmatic pulmonary physiology.

Concerns that airborne microplastics (MP) may be detrimental to human health are rising. However, research on the effects of MP on the respiratory system are limited. We tested the effect of MP exposure on both normal and asthmatic pulmonary physiology in mice. We show that MP exposure caused pulmonary inflammatory cell infiltration, bronchoalveolar macrophage aggregation, increased TNF-α level in bronchoalveolar lavage fluid (BALF), and increased plasma IgG1 production in normal mice. MP exposure also affected asthma symptoms by increasing mucus production and inflammatory cell infiltration with notable macrophage aggregation. Further, we found co-labeling of macrophage markers with MP incorporating fluorescence, which indicates phagocytosis of the MP by macrophages. A comparative transcriptomic analysis showed that MP exposure altered clusters of genes related to immune response, cellular stress response, and programmed cell death. A bioinformatics analysis further uncovered the molecular mechanism whereby MP stimulated production of tumor necrosis factor and immunoglobulins to activate a group of transmembrane B-cell antigens, leading to the modulation of cellular stress and programmed cell death in the asthma model. In summary, we show that MP exposure had detrimental effects on the respiratory system in both healthy and asthmatic mice, which calls for urgent discourse and action to mitigate environmental microplastic pollutants.

Low-Dose Radiation Can Cause Epigenetic Alterations Associated With Impairments in Both Male and Female Reproductive Cells.

Humans are regularly and continuously exposed to ionizing radiation from both natural and artificial sources. Cumulating evidence shows adverse effects of ionizing radiation on both male and female reproductive systems, including reduction of testis weight and sperm count and reduction of female germ cells and premature ovarian failure. While most of the observed effects were caused by DNA damage and disturbance of DNA repairment, ionizing radiation may also alter DNA methylation, histone, and chromatin modification, leading to epigenetic changes and transgenerational effects. However, the molecular mechanisms underlying the epigenetic changes and transgenerational reproductive impairment induced by low-dose radiation remain largely unknown. In this study, two different types of human ovarian cells and two different types of testicular cells were exposed to low dose of ionizing radiation, followed by bioinformatics analysis (including gene ontology functional analysis and Ingenuity Pathway Analysis), to unravel and compare epigenetic effects and pathway changes in male and female reproductive cells induced by ionizing radiation. Our findings showed that the radiation could alter the expression of gene cluster related to DNA damage responses through the control of MYC. Furthermore, ionizing radiation could lead to gender-specific reproductive impairment through deregulation of different gene networks. More importantly, the observed epigenetic modifications induced by ionizing radiation are mediated through the alteration of chromatin remodeling and telomere function. This study, for the first time, demonstrated that ionizing radiation may alter the epigenome of germ cells, leading to transgenerational reproductive impairments, and correspondingly call for research in this new emerging area which remains almost unknown.

CAG RNAs induce DNA damage and apoptosis by silencing NUDT16 expression in polyglutamine degeneration.

DNA damage plays a central role in the cellular pathogenesis of polyglutamine (polyQ) diseases, including Huntington's disease (HD). In this study, we showed that the expression of untranslatable expanded CAG RNA per se induced the cellular DNA damage response pathway. By means of RNA sequencing (RNA-seq), we found that expression of the Nudix hydrolase 16 (NUDT16) gene was down-regulated in mutant CAG RNA-expressing cells. The loss of NUDT16 function results in a misincorporation of damaging nucleotides into DNAs and leads to DNA damage. We showed that small CAG (sCAG) RNAs, species generated from expanded CAG transcripts, hybridize with CUG-containing NUDT16 mRNA and form a CAG-CUG RNA heteroduplex, resulting in gene silencing of NUDT16 and leading to the DNA damage and cellular apoptosis. These results were further validated using expanded CAG RNA-expressing mouse primary neurons and in vivo R6/2 HD transgenic mice. Moreover, we identified a bisamidinium compound, DB213, that interacts specifically with the major groove of the CAG RNA homoduplex and disfavors the CAG-CUG heteroduplex formation. This action subsequently mitigated RNA-induced silencing complex (RISC)-dependent NUDT16 silencing in both in vitro cell and in vivo mouse disease models. After DB213 treatment, DNA damage, apoptosis, and locomotor defects were rescued in HD mice. This work establishes NUDT16 deficiency by CAG repeat RNAs as a pathogenic mechanism of polyQ diseases and as a potential therapeutic direction for HD and other polyQ diseases.

Pan-cancer investigation reveals mechanistic insights of planar cell polarity gene Fuz in carcinogenesis.

The fuzzy planar cell polarity protein (Fuz) is an effector component of the planar cell polarity (PCP) signaling. Together with other core and effector proteins, the PCP pathway controls polarized cell movements. Fuz was also reported as a negative regulator of cell survival. In this study, we performed a pan-cancer survey to demonstrate the role of Fuz in multiple types of cancer. In head-neck squamous cell carcinoma and lung adenocarcinoma tumor samples, a reduction of Fuz transcript expression was detected. This coincides with the poor overall survival probabilities of these patients. We further showed that Fuz promoter hypermethylation contributes to its transcriptional downregulation. Meanwhile, we also identified a relatively higher mutation frequency at the 404th arginine amino acid residue in the coding sequence of Fuz locus, and further demonstrated that mutant Fuz proteins perturb the pro-apoptotic function of Fuz. In summary, our study unveiled an intriguing relationship between Fuz dysregulation and cancer prognosis, and further provides mechanistic insights of Fuz's involvement in carcinogenesis.

Horseshoe crab genomes reveal the evolution of genes and microRNAs after three rounds of whole genome duplication.

Whole genome duplication (WGD) has occurred in relatively few sexually reproducing invertebrates. Consequently, the WGD that occurred in the common ancestor of horseshoe crabs ~135 million years ago provides a rare opportunity to decipher the evolutionary consequences of a duplicated invertebrate genome. Here, we present a high-quality genome assembly for the mangrove horseshoe crab Carcinoscorpius rotundicauda (1.7 Gb, N50 = 90.2 Mb, with 89.8% sequences anchored to 16 pseudomolecules, 2n = 32), and a resequenced genome of the tri-spine horseshoe crab Tachypleus tridentatus (1.7 Gb, N50 = 109.7 Mb). Analyses of gene families, microRNAs, and synteny show that horseshoe crabs have undergone three rounds (3R) of WGD. Comparison of C. rotundicauda and T. tridentatus genomes from populations from several geographic locations further elucidates the diverse fates of both coding and noncoding genes. Together, the present study represents a cornerstone for improving our understanding of invertebrate WGD events on the evolutionary fates of genes and microRNAs, at both the individual and population level. We also provide improved genomic resources for horseshoe crabs, of applied value for breeding programs and conservation of this fascinating and unusual invertebrate lineage.