Evolutionary changes of the target sites of two microRNAs encoded in the Hox gene cluster of Drosophila and other insect species

High-resolution chromosome-level genome of Scylla paramamosain provides molecular insights into adaptive evolution in crabs

BACKGROUND

Evolutionary adaptation drives organismal adjustments to environmental pressures, exemplified in the diverse morphological and ecological adaptations seen in Decapoda crustaceans, particularly brachyuran crabs. Crabs thrive in diverse ecosystems, from coral reefs to hydrothermal vents and terrestrial habitats. Despite their ecological importance, the genetic mechanisms underpinning their developmental processes, reproductive strategies, and nutrient acquisition remain poorly understood.

RESULTS

Here, we report a comprehensive genomic analysis of the green mud crab Scylla paramamosain using ultralong sequencing technologies, achieving a high-quality chromosome-level assembly. The refined 1.21 Gb genome, with an impressive contig N50 of 11.45 Mb, offers a valuable genomic resource. The genome exhibits 33,662 protein-coding genes, enriched in various pathways related to development and environmental adaptation. Gene family analysis shows expansion in development-related pathways and contraction in metabolic pathways, indicating niche adaptations. Notably, investigation into Hox gene regulation sheds light on their role in pleopod development, with the Abd-A gene identified as a linchpin. Post-transcriptional regulation involving novel-miR1317 negatively regulates Abd-A levels. Furthermore, the potential role of fru gene in ovarian development and the identification of novel-miR35 as a regulator of Spfru2 add complexity to gene regulatory networks. Comparative functional analysis across Decapoda species reveals neo-functionalization of the elovl6 gene in the synthesis of long-chain polyunsaturated fatty acids (LC-PUFA), suggesting its importance in environmental adaptation.

CONCLUSIONS

Our findings shed light on various aspects of crab biology, including genome sequencing, assembly, and annotation, as well as gene family expansion, contraction, and regulatory mechanisms governing crucial developmental processes such as metamorphosis, reproductive strategies, and fatty acid metabolism.

Annotation of Hox cluster and Hox cofactor genes in the Asian citrus psyllid, Diaphorina citri, reveals novel features

Hox genes and their cofactors are essential developmental genes specifying regional identity in animals. Hox genes have a conserved arrangement in clusters in the same order in which they specify identity along the anterior-posterior axis. A few insect species have breaks in the cluster, but these are exceptions. We annotated the 10 Hox genes of the Asian citrus psyllid Diaphorina citri, and found a split in its Hox cluster between the Deformed and Sex combs reduced genes - the first time a break at this position has been observed in an insect Hox cluster. We also annotated D. citri orthologs of the Hox cofactor genes homothorax, PKNOX and extradenticle and found an additional copy of extradenticle in D. citri that appears to be a retrogene. Expression data and sequence conservation suggest that the extradenticle retrogene may have retained the original extradenticle function and allowed divergence of the parental extradenticle gene.

The Chinese mitten crab genome provides insights into adaptive plasticity and developmental regulation

The infraorder Brachyura (true or short-tailed crabs) represents a successful group of marine invertebrates yet with limited genomic resources. Here we report a chromosome-anchored reference genome and transcriptomes of the Chinese mitten crab Eriocheir sinensis, a catadromous crab and invasive species with wide environmental tolerance, strong osmoregulatory capacity and high fertility. We show the expansion of specific gene families in the crab, including F-ATPase, which enhances our knowledge on the adaptive plasticity of this successful invasive species. Our analysis of spatio-temporal transcriptomes and the genome of E. sinensis and other decapods shows that brachyurization development is associated with down-regulation of Hox genes at the megalopa stage when tail shortening occurs. A better understanding of the molecular mechanism regulating sexual development is achieved by integrated analysis of multiple omics. These genomic resources significantly expand the gene repertoire of Brachyura, and provide insights into the biology of this group, and Crustacea in general.

The Evolution of Imprinted microRNAs and Their RNA Targets

Mammalian genomes contain many imprinted microRNAs. When an imprinted miRNA targets an unimprinted mRNA their interaction may have different fitness consequences for the loci encoding the miRNA and mRNA. In one possible outcome, the mRNA sequence evolves to evade regulation by the miRNA by a simple change of target sequence. Such a response is unavailable if the targeted sequence is strongly constrained by other functions. In these cases, the mRNA evolves to accommodate regulation by the imprinted miRNA. These evolutionary dynamics are illustrated using the examples of the imprinted C19MC cluster of miRNAs in primates and C2MC cluster in mice that are paternally expressed in placentas. The 3′ UTR of PTEN, a gene with growth-related and metabolic functions, appears to be an important target of miRNAs from both clusters.

Long-Read Genome Sequencing and Assembly of Leptopilina boulardi: A Specialist Drosophila Parasitoid

Leptopilinaboulardi (Hymenoptera: Figitidae) is a specialist parasitoid of Drosophila The Drosophila-Leptopilina system has emerged as a suitable model for understanding several aspects of host-parasitoid biology. However, a good quality genome of the wasp counterpart was lacking. Here, we report a whole-genome assembly of L. boulardi to bring it in the scope of the applied and fundamental research on Drosophila parasitoids with access to epigenomics and genome editing tools. The 375Mb draft genome has an N50 of 275Kb with 6315 scaffolds >500bp and encompasses >95% complete BUSCOs. Using a combination of ab-initio and RNA-Seq based methods, 25259 protein-coding genes were predicted and 90% (22729) of them could be annotated with at least one function. We demonstrate the quality of the assembled genome by recapitulating the phylogenetic relationship of L. boulardi with other Hymenopterans. The key developmental regulators like Hox genes and sex determination genes are well conserved in L. boulardi, and so is the basic toolkit for epigenetic regulation. The search for epigenetic regulators has also revealed that L. boulardi genome possesses DNMT1 (maintenance DNA methyltransferase), DNMT2 (tRNA methyltransferase) but lacks the de novo DNA methyltransferase (DNMT3). Also, the heterochromatin protein 1 family appears to have expanded as compared to other hymenopterans. The draft genome of L. boulardi (Lb17) will expedite the research on Drosophila parasitoids. This genome resource and early indication of epigenetic aspects in its specialization make it an interesting system to address a variety of questions on host-parasitoid biology.

A homeotic shift late in development drives mimetic color variation in a bumble bee

Natural phenotypic radiations, with their high diversity and convergence, are well-suited for informing how genomic changes translate to natural phenotypic variation. New genomic tools enable discovery in such traditionally nonmodel systems. Here, we characterize the genomic basis of color pattern variation in bumble bees (Hymenoptera, Apidae, Bombus), a group that has undergone extensive convergence of setal color patterns as a result of Müllerian mimicry. In western North America, multiple species converge on local mimicry patterns through parallel shifts of midabdominal segments from red to black. Using genome-wide association, we establish that a cis-regulatory locus between the abdominal fate-determining Hox genes, abd-A and Abd-B, controls the red-black color switch in a western species, Bombus melanopygus Gene expression analysis reveals distinct shifts in Abd-B aligned with the duration of setal pigmentation at the pupal-adult transition. This results in atypical anterior Abd-B expression, a late developmental homeotic shift. Changing expression of Hox genes can have widespread effects, given their important role across segmental phenotypes; however, the late timing reduces this pleiotropy, making Hox genes suitable targets. Analysis of this locus across mimics and relatives reveals that other species follow independent genetic routes to obtain the same phenotypes.

Coping with living in the soil: the genome of the parthenogenetic springtail Folsomia candida

Background

Folsomia candida is a model in soil biology, belonging to the family of Isotomidae, subclass Collembola. It reproduces parthenogenetically in the presence of Wolbachia, and exhibits remarkable physiological adaptations to stress. To better understand these features and adaptations to life in the soil, we studied its genome in the context of its parthenogenetic lifestyle.

Results

We applied Pacific Bioscience sequencing and assembly to generate a reference genome for F. candida of 221.7 Mbp, comprising only 162 scaffolds. The complete genome of its endosymbiont Wolbachia, was also assembled and turned out to be the largest strain identified so far. Substantial gene family expansions and lineage-specific gene clusters were linked to stress response. A large number of genes (809) were acquired by horizontal gene transfer. A substantial fraction of these genes are involved in lignocellulose degradation. Also, the presence of genes involved in antibiotic biosynthesis was confirmed. Intra-genomic rearrangements of collinear gene clusters were observed, of which 11 were organized as palindromes. The Hox gene cluster of F. candida showed major rearrangements compared to arthropod consensus cluster, resulting in a disorganized cluster.

Conclusions

The expansion of stress response gene families suggests that stress defense was important to facilitate colonization of soils. The large number of HGT genes related to lignocellulose degradation could be beneficial to unlock carbohydrate sources in soil, especially those contained in decaying plant and fungal organic matter. Intra- as well as inter-scaffold duplications of gene clusters may be a consequence of its parthenogenetic lifestyle. This high quality genome will be instrumental for evolutionary biologists investigating deep phylogenetic lineages among arthropods and will provide the basis for a more mechanistic understanding in soil ecology and ecotoxicology.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3852-x) contains supplementary material, which is available to authorized users.

Composition and genomic organization of arthropod Hox clusters

Background

The ancestral arthropod is believed to have had a clustered arrangement of ten Hox genes. Within arthropods, Hox gene mutations result in transformation of segment identities. Despite the fact that variation in segment number/character was common in the diversification of arthropods, few examples of Hox gene gains/losses have been correlated with morphological evolution. Furthermore, a full appreciation of the variation in the genomic arrangement of Hox genes in extant arthropods has not been recognized, as genome sequences from each major arthropod clade have not been reported until recently. Initial genomic analysis of the chelicerate Tetranychusurticae suggested that loss of Hox genes and Hox gene clustering might be more common than previously assumed. To further characterize the genomic evolution of arthropod Hox genes, we compared the genomic arrangement and general characteristics of Hox genes from representative taxa from each arthropod subphylum.

Results

In agreement with others, we find arthropods generally contain ten Hox genes arranged in a common orientation in the genome, with an increasing number of sampled species missing either Hox3 or abdominal-A orthologs. The genomic clustering of Hox genes in species we surveyed varies significantly, ranging from 0.3 to 13.6 Mb. In all species sampled, arthropod Hox genes are dispersed in the genome relative to the vertebrate Mus musculus. Differences in Hox cluster size arise from variation in the number of intervening genes, intergenic spacing, and the size of introns and UTRs. In the arthropods surveyed, Hox gene duplications are rare and four microRNAs are, in general, conserved in similar genomic positions relative to the Hox genes.

Conclusions

The tightly clustered Hox complexes found in the vertebrates are not evident within arthropods, and differential patterns of Hox gene dispersion are found throughout the arthropods. The comparative genomic data continue to support an ancestral arthropod Hox cluster of ten genes with a shared orientation, with four Hox gene-associated miRNAs, although the degree of dispersion between genes in an ancestral cluster remains uncertain. Hox3 and abdominal-A orthologs have been lost in multiple, independent lineages, and current data support a model in which inversions of the Abdominal-B locus that result in the loss of abdominal-A correlate with reduced trunk segmentation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13227-016-0048-4) contains supplementary material, which is available to authorized users.

miR-1279, miR-548j, miR-548m, and miR-548d-5p binding sites in CDSs of paralogous and orthologous PTPN12, MSH6, and ZEB1 Genes

Only PTPN12, MSH6, and ZEB1 have significant miR-1279 binding sites among paralogous genes of human tyrosine phosphatase family, DNA mismatch repair family, and zinc finger family, respectively. All miRNA binding sites are located within CDSs of studied mRNAs. Nucleotide sequences of hsa-miR-1279 binding sites with mRNAs of human PTPN12, MSH6, and ZEB1 genes encode TKEQYE, EGSSDE, and GEKPYE oligopeptides, respectively. The conservation of miRNA binding sites encoding oligopeptides has been revealed. MRNAs of many paralogs of zinc finger gene family have from 1 to 12 binding sites coding the same GEKPYE hexapeptide. MRNAs of PTPN12, MSH6, and ZEB1 orthologous genes from different animal species have binding sites for hsa-miR-1279 which consist of homologous oligonucleotides encoding similar human oligopeptides TKEQYE, EGSSDE, and GEKPYE. MiR-548j, miR-548m, and miR-548d-5p have homologous binding sites in the mRNA of PTPN12 orthologous genes which encode PRTRSC, TEATDI, and STASAT oligopeptides, respectively. All regions of miRNA are important for binding with the mRNA.

miR156- and miR171-binding sites in the protein-coding sequences of several plant genes

We identified the interaction sites of several miRNAs with the mRNAs from paralogs and orthologs of the SPL and HAM genes in A. thaliana. miRNAs from the miR156 and miR157 families in A. thaliana are shown to have binding sites within the mRNAs of SPL genes. The ath-miR156a–j binding sites located in the mRNAs of the SPL paralogs contain the sequence GUGCUCUCUCUCUUCUGUCA. This sequence encodes the ALSLLS motif. miR157a–d bind to mRNAs of the SPL family at the same site. We suggest merging the miR156 and miR157 families into one family. Several SPL genes in eight plants contain conserved miR156 binding sites. GUGCUCUCUCUCUUCUGUCA polynucleotide is homologous in its binding sites. The ALSLLS hexapeptide is also conserved in the SPL proteins from these plants. Binding sites for ath-miR171a–c and ath-miR170 in HAM1, HAM2, and HAM3 paralog mRNAs are located in the CDSs. The conserved miRNA binding sequence GAUAUUGGCGCGGCUCAAUCA encodes the ILARLN hexapeptide. Nucleotides within the HAM1, HAM2, and HAM3 miRNA binding sites are conserved in the mRNAs of 37 orthologs from 13 plants. The miR171- and miR170-binding sites within the ortholog mRNAs were conserved and encode the ILARLN motif. We suggest that the ath-miR170 and ath-miR171a–c families should be in one family.

Structure, evolution and function of the bi-directionally transcribed iab-4/iab-8 microRNA locus in arthropods

In Drosophila melanogaster, the iab-4/iab-8 locus encodes bi-directionally transcribed microRNAs that regulate the function of flanking Hox transcription factors. We show that bi-directional transcription, temporal and spatial expression patterns and Hox regulatory function of the iab-4/iab-8 locus are conserved between fly and the beetle Tribolium castaneum. Computational predictions suggest iab-4 and iab-8 microRNAs can target common sites, and cell-culture assays confirm that iab-4 and iab-8 function overlaps on Hox target sites in both fly and beetle. However, we observe key differences in the way Hox genes are targeted. For instance, abd-A transcripts are targeted only by iab-8 in Drosophila, whereas both iab-4 and iab-8 bind to Tribolium abd-A. Our evolutionary and functional characterization of a bi-directionally transcribed microRNA establishes the iab-4/iab-8 system as a model for understanding how multiple products from sense and antisense microRNAs target common sites.

Building a robust a-p axis

Since the last common ancestor of Metazoa, animals have evolved complex body plans with specialized cells and spatial organization of tissues and organs. Arguably, one of the most significant innovations during animal evolutionary history was the establishment of a bilateral plane of symmetry on which morphological features (e.g. tissues, organs, appendages, skeleton) could be given specific coordinates within the animal along the anterior-posterior (A-P) and dorsal-ventral (D-V) axes. Hox genes are a known group of eumetazoan transcription factors central to regulating A-P patterning, but less well known and under current investigation is the broader regulatory landscape incorporating these genes, including microRNA (miRNA) regulation. The degree to which evolutionarily conserved targeting of Hox genes by Hox-embedded miRNAs contributes directly to A-P patterning is under investigation, yielding contrasting information dependent on the organism and miRNA of interest. The widespread A-P patterning defects observed in recent miR-196 loss-of-function studies solidifies the importance of miRNA regulation in Hox genetic hierarchies, and elucidating the developmental and evolutionary importance of all Hox-embedded miRNAs remains a challenge for the future.

Origins and evolution of microRNA genes in plant species

MicroRNAs (miRNAs) are among the most important regulatory elements of gene expression in animals and plants. However, their origin and evolutionary dynamics have not been studied systematically. In this paper, we identified putative miRNA genes in 11 plant species using the bioinformatic technique and examined their evolutionary changes. Our homology search indicated that no miRNA gene is currently shared between green algae and land plants. The number of miRNA genes has increased substantially in the land plant lineage, but after the divergence of eudicots and monocots, the number has changed in a lineage-specific manner. We found that miRNA genes have originated mainly by duplication of preexisting miRNA genes or protein-coding genes. Transposable elements also seem to have contributed to the generation of species-specific miRNA genes. The relative importance of these mechanisms in plants is quite different from that in Drosophila species, where the formation of hairpin structures in the genomes seems to be a major source of miRNA genes. This difference in the origin of miRNA genes between plants and Drosophila may be explained by the difference in the binding to target mRNAs between plants and animals. We also found that young miRNA genes are less conserved than old genes in plants as well as in Drosophila species. Yet, nearly half of the gene families in the ancestor of flowering plants have been lost in at least one species examined. This indicates that the repertoires of miRNA genes have changed more dynamically than previously thought during plant evolution.

Replication competent HIV-1 viruses that express intragenomic microRNA reveal discrete RNA-interference mechanisms that affect viral replication

BACKGROUND

It remains unclear whether retroviruses can encode and express an intragenomic microRNA (miRNA). Some have suggested that processing by the Drosha and Dicer enzymes might preclude the viability of a replicating retroviral RNA genome that contains a cis-embedded miRNA. To date, while many studies have shown that lentiviral vectors containing miRNAs can transduce mammalian cells and express the inserted miRNA efficiently, no study has examined the impact on the replication of a lentivirus such as HIV-1 after the deliberate intragenomic insertion of a bona fide miRNA.

RESULTS

We have constructed several HIV-1 molecular clones, each containing a discrete cellular miRNA positioned in Nef. These retroviral genomes express the inserted miRNA and are generally replication competent in T-cells. The inserted intragenomic miRNA was observed to elicit two different consequences for HIV-1 replication. First, the expression of miRNAs with predicted target sequences in the HIV-1 genome was found to reduce viral replication. Second, in one case, where an inserted miRNA was unusually well-processed by Drosha, this processing event inhibited viral replication.

CONCLUSION

This is the first study to examine in detail the replication competence of HIV-1 genomes that express cis-embedded miRNAs. The results indicate that a replication competent retroviral genome is not precluded from encoding and expressing a viral miRNA.