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.

Elucidating the ecophysiology of soybean pod-sucking stinkbug Riptortus pedestris (Hemiptera: Alydidae) based on de novo genome assembly and transcriptome analysis.

Food security is important for the ever-growing global population. Soybean, Glycine max (L.) Merr., is cultivated worldwide providing a key source of food, protein and oil. Hence, it is imperative to maintain or to increase its yield under different conditions including challenges caused by abiotic and biotic stresses. In recent years, the soybean pod-sucking stinkbug Riptortus pedestris has emerged as an important agricultural insect pest in East, South and Southeast Asia. Here, we present a genomics resource for R. pedestris including its genome assembly, messenger RNA (mRNA) and microRNA (miRNA) transcriptomes at different developmental stages and from different organs. As insect hormone biosynthesis genes (genes involved in metamorphosis) and their regulators such as miRNAs are potential targets for pest control, we analyzed the sesquiterpenoid (juvenile) and ecdysteroid (molting) hormone biosynthesis pathway genes including their miRNAs and relevant neuropeptides. Temporal gene expression changes of these insect hormone biosynthesis pathways were observed at different developmental stages. Similarly, a diet-specific response in gene expression was also observed in both head and salivary glands. Furthermore, we observed that microRNAs (bantam, miR-14, miR-316, and miR-263) of R. pedestris fed with different types of soybeans were differentially expressed in the salivary glands indicating a diet-specific response. Interestingly, the opposite arms of miR-281 (-5p and -3p), a miRNA involved in regulating development, were predicted to target Hmgs genes of R. pedestris and soybean, respectively. These observations among others highlight stinkbug's responses as a function of its interaction with soybean. In brief, the results of this study not only present salient findings that could be of potential use in pest management and mitigation but also provide an invaluable resource for R. pedestris as an insect model to facilitate studies on plant-pest interactions.

Self-supervised contrastive learning for integrative single cell RNA-seq data analysis.

We present a novel self-supervised Contrastive LEArning framework for single-cell ribonucleic acid (RNA)-sequencing (CLEAR) data representation and the downstream analysis. Compared with current methods, CLEAR overcomes the heterogeneity of the experimental data with a specifically designed representation learning task and thus can handle batch effects and dropout events simultaneously. It achieves superior performance on a broad range of fundamental tasks, including clustering, visualization, dropout correction, batch effect removal, and pseudo-time inference. The proposed method successfully identifies and illustrates inflammatory-related mechanisms in a COVID-19 disease study with 43 695 single cells from peripheral blood mononuclear cells.

stVAE deconvolves cell-type composition in large-scale cellular resolution spatial transcriptomics.

MOTIVATION: Recent rapid developments in spatial transcriptomic techniques at cellular resolution have gained increasing attention. However, the unique characteristics of large-scale cellular resolution spatial transcriptomic datasets, such as the limited number of transcripts captured per spot and the vast number of spots, pose significant challenges to current cell-type deconvolution methods. RESULTS: In this study, we introduce stVAE, a method based on the variational autoencoder framework to deconvolve the cell-type composition of cellular resolution spatial transcriptomic datasets. To assess the performance of stVAE, we apply it to five datasets across three different biological tissues. In the Stereo-seq and Slide-seqV2 datasets of the mouse brain, stVAE accurately reconstructs the laminar structure of the pyramidal cell layers in the cortex, which are mainly organized by the subtypes of telencephalon projecting excitatory neurons. In the Stereo-seq dataset of the E12.5 mouse embryo, stVAE resolves the complex spatial patterns of osteoblast subtypes, which are supported by their marker genes. In Stereo-seq and Pixel-seq datasets of the mouse olfactory bulb, stVAE accurately delineates the spatial distributions of known cell types. In summary, stVAE can accurately identify spatial patterns of cell types and their relative proportions across spots for cellular resolution spatial transcriptomic data. It is instrumental in understanding the heterogeneity of cell populations and their interactions within tissues. AVAILABILITY AND IMPLEMENTATION: stVAE is available in GitHub (https://github.com/lichen2018/stVAE) and Figshare (https://figshare.com/articles/software/stVAE/23254538).

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.

The genome and sex-dependent responses to temperature in the common yellow butterfly, Eurema hecabe.

BACKGROUND: Lepidoptera (butterflies and moths) is one of the most geographically widespread insect orders in the world, and its species play important and diverse ecological and applied roles. Climate change is one of the biggest challenges to biodiversity this century, and lepidopterans are vulnerable to climate change. Temperature-dependent gene expression differences are of relevance under the ongoing climate crisis. However, little is known about how climate affects gene expression in lepidopterans and the ecological consequences of this, particularly with respect to genes with biased expression in one of the sexes. The common yellow butterfly, Eurema hecabe (Family Pieridae), is one of the most geographically widespread lepidopterans that can be found in Asia, Africa, and Australia. Nevertheless, what temperature-dependent effects there may be and whether the effects differ between the sexes remain largely unexplored. RESULTS: Here, we generated high-quality genomic resources for E. hecabe along with transcriptomes from eight developmental stages. Male and female butterflies were subjected to varying temperatures to assess sex-specific gene expression responses through mRNA and microRNA transcriptomics. We find that there are more temperature-dependent sex-biased genes in females than males, including genes that are involved in a range of biologically important functions, highlighting potential ecological impacts of increased temperatures. Further, by considering available butterfly data on sex-biased gene expression in a comparative genomic framework, we find that the pattern of sex-biased gene expression identified in E. hecabe is highly species-specific, rather than conserved across butterfly species, suggesting that sex-biased gene expression responses to climate change are complex in butterflies. CONCLUSIONS: Our study lays the foundation for further understanding of differential responses to environmental stress in a widespread lepidopteran model and demonstrates the potential complexity of sex-specific responses of lepidopterans to climate change.

Millipede genomes reveal unique adaptations during myriapod evolution.

The Myriapoda, composed of millipedes and centipedes, is a fascinating but poorly understood branch of life, including species with a highly unusual body plan and a range of unique adaptations to their environment. Here, we sequenced and assembled 2 chromosomal-level genomes of the millipedes Helicorthomorpha holstii (assembly size = 182 Mb; shortest scaffold/contig length needed to cover 50% of the genome [N50] = 18.11 Mb mainly on 8 pseudomolecules) and Trigoniulus corallinus (assembly size = 449 Mb, N50 = 26.78 Mb mainly on 17 pseudomolecules). Unique genomic features, patterns of gene regulation, and defence systems in millipedes, not observed in other arthropods, are revealed. Both repeat content and intron size are major contributors to the observed differences in millipede genome size. Tight Hox and the first loose ecdysozoan ParaHox homeobox clusters are identified, and a myriapod-specific genomic rearrangement including Hox3 is also observed. The Argonaute (AGO) proteins for loading small RNAs are duplicated in both millipedes, but unlike in insects, an AGO duplicate has become a pseudogene. Evidence of post-transcriptional modification in small RNAs-including species-specific microRNA arm switching-providing differential gene regulation is also obtained. Millipedes possesses a unique ozadene defensive gland unlike the venomous forcipules found in centipedes. We identify sets of genes associated with the ozadene that play roles in chemical defence as well as antimicrobial activity. Macro-synteny analyses revealed highly conserved genomic blocks between the 2 millipedes and deuterostomes. Collectively, our analyses of millipede genomes reveal that a series of unique adaptations have occurred in this major lineage of arthropod diversity. The 2 high-quality millipede genomes provided here shed new light on the conserved and lineage-specific features of millipedes and centipedes. These findings demonstrate the importance of the consideration of both centipede and millipede genomes-and in particular the reconstruction of the myriapod ancestral situation-for future research to improve understanding of arthropod evolution, and animal evolutionary genomics more widely.

Effects of In Utero PFOS Exposure on Epigenetics and Metabolism in Mouse Fetal Livers.

Prenatal exposure to perfluorooctanesulfonate (PFOS) increases fetus' metabolic risk; however, the investigation of the underlying mechanism is limited. In this study, pregnant mice in the gestational days (GD, 4.5-17.5) were exposed to PFOS (0.3 and 3 µg/g of body weight). At GD 17.5, PFOS perturbed maternal lipid metabolism and upregulated metabolism-regulating hepatokines (Angptl4, Angptl8, and Selenop). Mass-spectrometry imaging and whole-genome bisulfite sequencing revealed, respectively, selective PFOS localization and deregulation of gene methylation in fetal livers, involved in inflammation, glucose, and fatty acid metabolism. PCR and Western blot analysis of lipid-laden fetal livers showed activation of AMPK signaling, accompanied by significant increases in the expression of glucose transporters (Glut2/4), hexose-phosphate sensors (Retsat and ChREBP), and the key glycolytic enzyme, pyruvate kinase (Pk) for glucose catabolism. Additionally, PFOS modulated the expression levels of PPARα and PPARγ downstream target genes, which simultaneously stimulated fatty acid oxidation (Cyp4a14, Acot, and Acox) and lipogenesis (Srebp1c, Acaca, and Fasn). Using human normal hepatocyte (MIHA) cells, the underlying mechanism of PFOS-elicited nuclear translocation of ChREBP, associated with a fatty acid synthesizing pathway, was revealed. Our finding implies that in utero PFOS exposure altered the epigenetic landscape associated with dysregulation of fetal liver metabolism, predisposing postnatal susceptibility to metabolic challenges.

RNA G-quadruplexes (rG4s): genomics and biological functions.

G-quadruplexes (G4s) are non-classical DNA or RNA secondary structures that have been first observed decades ago. Over the years, these four-stranded structural motifs have been demonstrated to have significant regulatory roles in diverse biological processes, but challenges remain in detecting them globally and reliably. Compared to DNA G4s (dG4s), the study of RNA G4s (rG4s) has received less attention until recently. In this review, we will summarize the innovative high-throughput methods recently developed to detect rG4s on a transcriptome-wide scale, highlight the many novel and important functions of rG4 being discovered in vivo across the tree of life, and discuss the key biological questions to be addressed in the near future.

G-quadruplex RNA motifs influence gene expression in the malaria parasite Plasmodium falciparum.

G-quadruplexes are non-helical secondary structures that can fold in vivo in both DNA and RNA. In human cells, they can influence replication, transcription and telomere maintenance in DNA, or translation, transcript processing and stability of RNA. We have previously showed that G-quadruplexes are detectable in the DNA of the malaria parasite Plasmodium falciparum, despite a very highly A/T-biased genome with unusually few guanine-rich sequences. Here, we show that RNA G-quadruplexes can also form in P. falciparum RNA, using rG4-seq for transcriptome-wide structure-specific RNA probing. Many of the motifs, detected here via the rG4seeker pipeline, have non-canonical forms and would not be predicted by standard in silico algorithms. However, in vitro biophysical assays verified formation of non-canonical motifs. The G-quadruplexes in the P. falciparum transcriptome are frequently clustered in certain genes and associated with regions encoding low-complexity peptide repeats. They are overrepresented in particular classes of genes, notably those that encode PfEMP1 virulence factors, stress response genes and DNA binding proteins. In vitro translation experiments and in vivo measures of translation efficiency showed that G-quadruplexes can influence the translation of P. falciparum mRNAs. Thus, the G-quadruplex is a novel player in post-transcriptional regulation of gene expression in this major human pathogen.

Enhanced transcriptome-wide RNA G-quadruplex sequencing for low RNA input samples with rG4-seq 2.0.

BACKGROUND: RNA G-quadruplexes (rG4s) are non-canonical structural motifs that have diverse functional and regulatory roles, for instance in transcription termination, alternative splicing, mRNA localization and stabilization, and translational process. We recently developed the RNA G-quadruplex structure sequencing (rG4-seq) technique and described rG4s in both eukaryotic and prokaryotic transcriptomes. However, rG4-seq suffers from a complicated gel purification step and limited PCR product yield, thus requiring a high amount of RNA input, which limits its applicability in more physiologically or clinically relevant studies often characterized by the limited availability of biological material and low RNA abundance. Here, we redesign and enhance the workflow of rG4-seq to address this issue. RESULTS: We developed rG4-seq 2.0 by introducing a new ssDNA adapter containing deoxyuridine during library preparation to enhance library quality with no gel purification step, less PCR amplification cycles and higher yield of PCR products. We demonstrate that rG4-seq 2.0 produces high-quality cDNA libraries that support reliable and reproducible rG4 identification at varying RNA inputs, including RNA mounts as low as 10 ng. rG4-seq 2.0 also improved the rG4-seq calling outcome and nucleotide bias in rG4 detection persistent in rG4-seq 1.0. We further provide in vitro mapping of rG4 in the HEK293T cell line, and recommendations for assessing RNA input and sequencing depth for individual rG4 studies based on transcript abundance. CONCLUSIONS: rG4-seq 2.0 can improve the identification and study of rG4s in low abundance transcripts, and our findings can provide insights to optimize cDNA library preparation in other related methods.

Increased copy number of gibberellin 2-oxidase 8 genes reduced trailing growth and shoot length during soybean domestication.

Copy number variations (CNVs) play important roles in crop domestication. However, there is only very limited information on the involvement of CNVs in soybean domestication. Trailing growth and long shoots are soybean adaptations for natural habitats but cause lodging that hampers yield in cultivation. Previous studies have focused on Dt1/2 affecting the indeterminate/determinate growth habit, whereas the possible role of the gibberellin pathway remained unclear. In the present study, quantitative trait locus (QTL) mapping of a recombinant inbred population of 460 lines revealed a trailing-growth-and-shoot-length QTL. A CNV region within this QTL was identified, featuring the apical bud-expressed gibberellin 2-oxidase 8A/B, the copy numbers of which were positively correlated with expression levels and negatively with trailing growth and shoot length, and their effects were demonstrated by transgenic soybean and Arabidopsis thaliana. Based on the fixation index, this CNV region underwent intense selection during the initial domestication process.

Differential microRNA expression, microRNA arm switching, and microRNA:long noncoding RNA interaction in response to salinity stress in soybean.

BACKGROUND: Soybean is a major legume crop with high nutritional and environmental values suitable for sustainable agriculture. Noncoding RNAs (ncRNAs), including microRNAs (miRNAs) and long noncoding RNAs (lncRNAs), are important regulators of gene functions in eukaryotes. However, the interactions between these two types of ncRNAs in the context of plant physiology, especially in response to salinity stress, are poorly understood. RESULTS: Here, we challenged a cultivated soybean accession (C08) and a wild one (W05) with salt treatment and obtained their small RNA transcriptomes at six time points from both root and leaf tissues. In addition to thoroughly analyzing the differentially expressed miRNAs, we also documented the first case of miRNA arm-switching (miR166m), the swapping of dominant miRNA arm expression, in soybean in different tissues. Two arms of miR166m target different genes related to salinity stress (chloroplastic beta-amylase 1 targeted by miR166m-5p and calcium-dependent protein kinase 1 targeted by miR166m-3p), suggesting arm-switching of miR166m play roles in soybean in response to salinity stress. Furthermore, two pairs of miRNA:lncRNA interacting partners (miR166i-5p and lncRNA Gmax_MSTRG.35921.1; and miR394a-3p and lncRNA Gmax_MSTRG.18616.1) were also discovered in reaction to salinity stress. CONCLUSIONS: This study demonstrates how ncRNA involves in salinity stress responses in soybean by miRNA arm switching and miRNA:lncRNA interactions. The behaviors of ncRNAs revealed in this study will shed new light on molecular regulatory mechanisms of stress responses in plants, and hence provide potential new strategies for crop improvement.

Uncovering the functions of plasma proteins in ulcerative colitis and identifying biomarkers for BPA-induced severe ulcerative colitis: A plasma proteome analysis.

Ulcerative colitis (UC), a long-term inflammation of the colon, is a worldwide disease. Accumulating reports have suggested the contribution of environmental pollutants to UC development. As such, the identification of biomarkers to evaluate pollutant-induced UC could provide a better assessment on the world's pollution problem. In the present study, we applied the plasma proteome to profile the plasma protein changes in three models: dextran sulfate sodium (DSS)-induced colitis, bisphenol A (BPA), and BPA-severe colitis. We aimed to investigate the functional roles of plasma proteins related to colitis development and further understand the synergistic effect of BPA on colitis. In addition, we aimed to identify novel biomarkers for UC non-invasive diagnosis and assessment of BPA-induced colitis. Our results showed a significant dysregulation of plasma proteins in these three models. Bioinformatics analysis, including gene ontology, Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis, and Ingenuity Pathway Analysis, highlighted the important effects of these dysregulated plasma proteins in immune and inflammatory responses through the regulation of CCR3 signaling in eosinophils, PI3K signaling in B lymphocytes, CD28 signaling in T helper cells, and leukocyte extravasation signaling in DSS-induced colitis model. Furthermore, our data suggested that BPA exposure altered the plasma proteins involved in lipid-related metabolic processes, leukocyte cell-cell adhesion and cytokine response. More importantly, we identified plasma proteins, ALB, APOA4, C3, CFB, DPEP1, HP, LTF, and Retnlg as biomarkers for assessing BPA-induced colitis.

Single-Cell Atlas of the Drosophila Leg Disc Identifies a Long Non-Coding RNA in Late Development.

The Drosophila imaginal disc has been an excellent model for the study of developmental gene regulation. In particular, long non-coding RNAs (lncRNAs) have gained widespread attention in recent years due to their important role in gene regulation. Their specific spatiotemporal expressions further support their role in developmental processes and diseases. In this study, we explored the role of a novel lncRNA in Drosophila leg development by dissecting and dissociating w1118 third-instar larval third leg (L3) discs into single cells and single nuclei, and performing single-cell RNA-sequencing (scRNA-seq) and single-cell assays for transposase-accessible chromatin (scATAC-seq). Single-cell transcriptomics analysis of the L3 discs across three developmental timepoints revealed different cell types and identified lncRNA:CR33938 as a distal specific gene with high expression in late development. This was further validated by fluorescence in-situ hybridization (FISH). The scATAC-seq results reproduced the single-cell transcriptomics landscape and elucidated the distal cell functions at different timepoints. Furthermore, overexpression of lncRNA:CR33938 in the S2 cell line increased the expression of leg development genes, further elucidating its potential role in development.

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.

Identification of Smad3-related transcriptomes in type-2 diabetic nephropathy by whole transcriptome RNA sequencing.

Smad3 deficiency prevents the development of type 2 diabetic nephropathy; however, the underlying molecular mechanisms remain unknown. In this study, we aimed to identify Smad3-related genes involved in the pathogenesis of diabetic kidney disease. High-throughput RNA sequencing was performed to profile the whole transcriptome in the diabetic kidney of Smad3 WT-db/db, Smad3 KO-db/db, Smad3+/- db/db and their littermate control db/m mice at 20 weeks. The gene ontology, pathways and alternative splicing of differentially expressed protein-coding genes and long non-coding RNAs related to Smad3 in diabetic kidney were analysed. Compared to Smad3 WT-db/db mice, Smad3 KO-db/db mice exhibited an alteration of genes associated with RNA splicing and metabolism, whereas heterozygosity deletion of Smad3 (Smad3+/- db/db mice) significantly altered genes related to cell division and cell cycle. Notably, three protein-coding genes (Upk1b, Psca and Gdf15) and two lncRNAs (NONMMUG023520.2 and NONMMUG032975.2) were identified to be Smad3-dependent and to be associated with the development of diabetic nephropathy. By using whole transcriptome RNA sequencing, we identified novel Smad3 transcripts related to the development of diabetic nephropathy. Thus, targeting these transcripts may represent a novel and effective therapy for diabetic nephropathy.

Comparative genomics of the coconut crab and other decapod crustaceans: exploring the molecular basis of terrestrial adaptation.

BACKGROUND: The complex life cycle of the coconut crab, Birgus latro, begins when an obligate terrestrial adult female visits the intertidal to hatch zoea larvae into the surf. After drifting for several weeks in the ocean, the post-larval glaucothoes settle in the shallow subtidal zone, undergo metamorphosis, and the early juveniles then subsequently make their way to land where they undergo further physiological changes that prevent them from ever entering the sea again. Here, we sequenced, assembled and analyzed the coconut crab genome to shed light on its adaptation to terrestrial life. For comparison, we also assembled the genomes of the long-tailed marine-living ornate spiny lobster, Panulirus ornatus, and the short-tailed marine-living red king crab, Paralithodes camtschaticus. Our selection of the latter two organisms furthermore allowed us to explore parallel evolution of the crab-like form in anomurans. RESULTS: All three assembled genomes are large, repeat-rich and AT-rich. Functional analysis reveals that the coconut crab has undergone proliferation of genes involved in the visual, respiratory, olfactory and cytoskeletal systems. Given that the coconut crab has atypical mitochondrial DNA compared to other anomurans, we argue that an abundance of kif22 and other significantly proliferated genes annotated with mitochondrial and microtubule functions, point to unique mechanisms involved in providing cellular energy via nuclear protein-coding genes supplementing mitochondrial and microtubule function. We furthermore detected in the coconut crab a significantly proliferated HOX gene, caudal, that has been associated with posterior development in Drosophila, but we could not definitively associate this gene with carcinization in the Anomura since it is also significantly proliferated in the ornate spiny lobster. However, a cuticle-associated coatomer gene, gammacop, that is significantly proliferated in the coconut crab, may play a role in hardening of the adult coconut crab abdomen in order to mitigate desiccation in terrestrial environments. CONCLUSION: The abundance of genomic features in the three assembled genomes serve as a source of hypotheses for future studies of anomuran environmental adaptations such as shell-utilization, perception of visual and olfactory cues in terrestrial environments, and cuticle sclerotization. We hypothesize that the coconut crab exhibits gene proliferation in lieu of alternative splicing as a terrestrial adaptation mechanism and propose life-stage transcriptomic assays to test this hypothesis.

Micro-RNA Clusters Integrate Evolutionary Constraints on Expression and Target Affinities: The miR-6/5/4/286/3/309 Cluster in Drosophila.

A striking feature of micro-RNAs is that they are often clustered in the genomes of animals. The functional and evolutionary consequences of this clustering remain obscure. Here, we investigated a micro-RNA cluster miR-6/5/4/286/3/309 that is conserved across drosophilid lineages. Small RNA sequencing revealed expression of this micro-RNA cluster in Drosophila melanogaster leg discs, and conditional overexpression of the whole cluster resulted in leg appendage shortening. Transgenic overexpression lines expressing different combinations of micro-RNA cluster members were also constructed. Expression of individual micro-RNAs from the cluster resulted in a normal wild-type phenotype, but either the expression of several ancient micro-RNAs together (miR-5/4/286/3/309) or more recently evolved clustered micro-RNAs (miR-6-1/2/3) can recapitulate the phenotypes generated by the whole-cluster overexpression. Screening of transgenic fly lines revealed downregulation of leg-patterning gene cassettes in generation of the leg-shortening phenotype. Furthermore, cell transfection with different combinations of micro-RNA cluster members revealed a suite of downstream genes targeted by all cluster members, as well as complements of targets that are unique for distinct micro-RNAs. Considered together, the micro-RNA targets and the evolutionary ages of each micro-RNA in the cluster demonstrate the importance of micro-RNA clustering, where new members can reinforce and modify the selection forces on both the cluster regulation and the gene regulatory network of existing micro-RNAs. Key words: micro-RNA, cluster, evolution.