Transcriptional co-regulation of evolutionarily conserved microRNA/cone opsin gene pairs: Implications for photoreceptor subtype specification

Partitioning of gene expression among zebrafish photoreceptor subtypes

Vertebrate photoreceptors are categorized into two broad classes, rods and cones, responsible for dim- and bright-light vision, respectively. While many molecular features that distinguish rods and cones are known, gene expression differences among cone subtypes remain poorly understood. Teleost fishes are renowned for the diversity of their photoreceptor systems. Here, we used single-cell RNA-seq to profile adult photoreceptors in zebrafish, a teleost. We found that in addition to the four canonical zebrafish cone types, there exist subpopulations of green and red cones (previously shown to be located in the ventral retina) that express red-shifted opsin paralogs (opn1mw4 or opn1lw1) as well as a unique combination of cone phototransduction genes. Furthermore, the expression of many paralogous phototransduction genes is partitioned among cone subtypes, analogous to the partitioning of the phototransduction paralogs between rods and cones seen across vertebrates. The partitioned cone-gene pairs arose via the teleost-specific whole-genome duplication or later clade-specific gene duplications. We also discovered that cone subtypes express distinct transcriptional regulators, including many factors not previously implicated in photoreceptor development or differentiation. Overall, our work suggests that partitioning of paralogous gene expression via the action of differentially expressed transcriptional regulators enables diversification of cone subtypes in teleosts.

Independent recurrent evolution of MICRORNA genes converging onto similar non-canonical organisation across green plant lineages is driven by local and segmental duplication events in species, family and lineages

The relationship between evolutionary history, organisation and transcriptional regulation of genes are intrinsically linked. These have been well studied in canonically organised protein-coding genes but not of MIRNAs. In the present study, we investigated the non-canonical arrangement of MIRNAs across taxonomic boundaries from algae to angiosperms employing a combination of genome organization, phylogeny and synteny. We retrieved the complete dataset of MIRNA from twenty-five species to identify and classify based on organisational patterns. The median size of cluster was between 2-5 kb and between 1-20 % of all MIRNAs are organized in head-to-head (with bidirectional promoter), head-to-tail (tandem), and overlapping manner. Although majority of the clusters are composed of MIRNA homologs, 25% of all clusters comprises of non-homologous genes with a potential of generating functional and regulatory complexity. A comparison of phylogeny and organizational patterns revealed that multiple independent events, some of which are species-specific, and some ancient, in different lineages, are responsible for non-canonical organization. Detailed investigation of MIR395 family across the plants revealed a complex origin of non-canonical arrangement through ancient and recent, segmental and local duplications; analysis of MIR399 family revealed major expansion occurred prior to monocot-dicot split, with few lineage-specific events. Evolution of "convergent" organization pattern of non-canonical arrangement originating from independent loci through recurrent event highlights our poor understanding of evolutionary process of MIRNA genes. The present investigation thus paves way for comparative functional genomics to understand the role of non-canonical organization on transcriptional regulation and regulatory diversity in MIRNA gene families.

Comparative microRNAs expression profiles analysis during embryonic development of common carp, Cyprinus carpio

MicroRNAs (miRNAs) play important roles in biological processes by regulating specific gene expression. Limited miRNAs information is available on embryonic development in common carp (Cyprinus carpio) so far. In this study, six important embryonic development stages of C.carpio were collected to perform a times-series of small RNA-seq experiments from cleavage, blastocyst, gastrulation, organ formation, hatching stage to 1 day post-hatching larva. The expression profiles of miRNAs were identified and differentially expressed miRNAs (DEMs) were screened out based on pairwise comparison. A mean of 12,744,989 raw reads and 9,888,123 clean reads were obtained from each library. A total of 2565 miRNAs were identified. 68 of 204 DEMs were overlapped with stage-specific miRNAs, in which 15 were known miRNAs and seemed to play a key role in embryogenesis. Additionally, time-course expression reveals several intriguing fluctuations during embryogenesis. Numerous signaling pathways were identified in embryonic development, including the phototransduction, hippo signaling pathway, Wnt, melanogenesis, histidine metabolism and fatty acid biosynthesis. The results would provide new insight into the roles of miRNAs in embryonic development, and would help us to advance the understanding of miRNA-mediated mechanisms in embryonic development of fish.

Movement of transposable elements contributes to cichlid diversity

African cichlid fishes are a prime model for studying speciation mechanisms. Despite the development of extensive genomic resources, it has been difficult to determine which sources of genetic variation are responsible for cichlid phenotypic variation. One of their most variable phenotypes is visual sensitivity, with some of the largest spectral shifts among vertebrates. These shifts arise primarily from differential expression of seven cone opsin genes. By mapping expression quantitative trait loci (eQTL) in intergeneric crosses of Lake Malawi cichlids, we previously identified four causative genetic variants that correspond to indels in the promoters of either key transcription factors or an opsin gene. In this comprehensive study, we show that these indels are the result of the movement of transposable elements (TEs) that correlate with opsin expression variation across the Malawi flock. In tracking the evolutionary history of these particular indels, we found they are endemic to Lake Malawi, suggesting that these TEs are recently active and are segregating within the Malawi cichlid lineage. However, an independent indel has arisen at a similar genomic location in one locus outside of the Malawi flock. The convergence in TE movement suggests these loci are primed for TE insertion and subsequent deletions. Increased TE mobility may be associated with interspecific hybridization, which disrupts mechanisms of TE suppression. This might provide a link between cichlid hybridization and accelerated regulatory variation. Overall, our study suggests that TEs may be an important driver of key regulatory changes, facilitating rapid phenotypic change and possibly speciation in African cichlids.

Axes of visual adaptation in the ecologically diverse family Cichlidae

The family Cichlidae contains approximately 2000 species that live in diverse freshwater habitats including murky lakes, turbid rivers, and clear lakes from both the Old and New Worlds. Their visual systems are similarly diverse and have evolved specific sensitivities that differ along several axes of variation. Variation in cornea and lens transmission affect which wavelengths reach the retina. Variation in photoreceptor number and distribution affect brightness sensitivity, spectral sensitivity and resolution. Probably their most dynamic characteristic is the variation in visual pigment peak sensitivities. Visual pigments can be altered through changes in chromophore, opsin sequence and opsin expression. Opsin expression varies by altering which of the seven available cone opsins in their genomes are turned on. These opsins can even be coexpressed to produce seemingly infinitely tunable cone sensitivities. Both chromophore and opsin expression can vary on either rapid (hours or days), slower (seasonal or ontogenetic) or evolutionary timescales. Such visual system shifts have enabled cichlids to adapt to different habitats and foraging styles. Through both short term plasticity and longer evolutionary adaptations, cichlids have proven to be ecologically successful and an excellent model for studying organismal adaptation.

Thyroid hormone receptors mediate two distinct mechanisms of long-wavelength vision

Thyroid hormone (TH) signaling plays an important role in the regulation of long-wavelength vision in vertebrates. In the retina, thyroid hormone receptor β (thrb) is required for expression of long-wavelength-sensitive opsin (lws) in red cone photoreceptors, while in retinal pigment epithelium (RPE), TH regulates expression of a cytochrome P450 enzyme, cyp27c1, that converts vitamin A₁ into vitamin A₂ to produce a red-shifted chromophore. To better understand how TH controls these processes, we analyzed the phenotype of zebrafish with mutations in the three known TH nuclear receptor transcription factors (thraa, thrab, and thrb). We found that no single TH nuclear receptor is required for TH-mediated induction of cyp27c1 but that deletion of all three (thraa ^-/- ;thrab ^-/- ;thrb ^-/- ) completely abrogates its induction and the resulting conversion of A₁- to A₂-based retinoids. In the retina, loss of thrb resulted in an absence of red cones at both larval and adult stages without disruption of the underlying cone mosaic. RNA-sequencing analysis revealed significant down-regulation of only five genes in adult thrb ^-/- retina, of which three (lws1, lws2, and miR-726) occur in a single syntenic cluster. In the thrb ^-/- retina, retinal progenitors destined to become red cones were transfated into ultraviolet (UV) cones and horizontal cells. Taken together, our findings demonstrate cooperative regulation of cyp27c1 by TH receptors and a requirement for thrb in red cone fate determination. Thus, TH signaling coordinately regulates both spectral sensitivity and sensory plasticity.

Evolutionary history of the medaka long-wavelength sensitive genes and effects of artificial regression by gene loss on behavioural photosensitivity

Tandem gene duplication has led to an expansion of cone-opsin repertoires in many fish, but the resulting functional advantages have only been conjectured without empirical demonstration. Medaka (Oryzias latipes and O. sakaizumii) have eight (two red, three green, two blue, and one violet) cone opsin genes. Absorbance maxima (λ_max) of the proteins vary from 356 nm to 562 nm, but those of the red opsins (long-wavelength sensitive; LWS) are nearly identical, obscuring the necessity of their coexistence. Here, we compared the LWSa and LWSb loci of these sister species and found that the gene duplication occurred long before the latipes–sakaizumii speciation (4–18 million years ago), and the high sequence similarity between the paralogues is the result of at least two events of gene conversion. These repetitive gene conversions would indicate the importance for medaka of retaining two identical LWSs in the genome. However, a newly established medaka mutant with a single LWS showed no defect in LWS expression or behavioural red-light sensitivity, demonstrating functional redundancy of the paralogs. Thus, as with many other genes after whole-genome duplication, the redundant LWS might be on the way to being lost from the current cone opsin repertoire. Thus, non-allelic gene conversion may temporarily provide an easier and more frequent solution than gene loss for reducing genetic diversity, which should be considered when assessing history of gene evolution by phylogenetic analyses.

Multiple trans QTL and one cis-regulatory deletion are associated with the differential expression of cone opsins in African cichlids

Background

Dissecting the genetic basis of phenotypic diversity is one of the fundamental goals in evolutionary biology. Despite growing evidence for gene expression divergence being responsible for the evolution of complex traits, knowledge about the proximate genetic causes underlying these traits is still limited. African cichlids have diverse visual systems, with different species expressing different combinations of seven cone opsin genes. Using opsin expression variation in African cichlids as a model for gene expression evolution, this study aims to investigate the genetic architecture of opsin expression divergence in this group.

Results

Results from a genome-wide linkage mapping on the F₂ progeny of an intergeneric cross, between two species with differential opsin expression show that opsins in Lake Malawi cichlids are controlled by multiple quantitative trait loci (QTLs). Most of these QTLs are located in trans to the opsins except for one cis-QTL for SWS1 on LG17. A closer look at this major QTL revealed the presence of a 691 bp deletion in the promoter of the SWS1 opsin (located 751 bp upstream of the start site) that is associated with a decrease in its expression. Phylogenetic footprinting indicates that the region spanning the deletion harbors a microRNA miR-729 and a conserved non-coding element (CNE) that also occurs in zebrafish and other teleosts. This suggests that the deletion might contain ancestrally preserved regulators that have been tuned for SWS1 gene expression in Lake Malawi. While this deletion is not common, it does occur in several other species within the lake.

Conclusions

Differential expression of cichlid opsins is associated with multiple overlapping QTL, with all but one in trans to the opsins they regulate. The one cis-acting factor is a deletion in the promoter of the SWS1 opsin, suggesting that ancestral polymorphic deletions may contribute to cichlid’s visual diversity. In addition to expanding our understanding of the molecular landscape of opsin expression in African cichlids, this study sheds light on the molecular mechanisms underlying phenotypic variation in natural populations.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5328-z) contains supplementary material, which is available to authorized users.

Dynamic mRNA and miRNA expression analysis in response to hypoxia and reoxygenation in the blunt snout bream (Megalobrama amblycephala)

Adaptation to hypoxia is a complex process involving various pathways and regulation mechanisms. A better understanding of the genetic influence on these mechanisms could permit selection for hypoxia-sensitive fish. To aid this understanding, an integrated analysis of miRNA and mRNA expression was performed in Megalobrama amblycephala under four acute hypoxia and reoxygenation stages. A number of significantly differentially-expressed miRNAs and genes associated with oxidative stress were identified, and their functional characteristics were revealed by GO function and KEGG pathway analysis. They were found to be involved in HIF-1 pathways known to affect energy metabolism and apoptosis. MiRNA-mRNA interaction pairs were detected from comparison of expression between the four different stages. The function annotation results also showed that many miRNA-mRNA interaction pairs were likely to be involved in regulating hypoxia stress. As a unique resource for gene expression and regulation during hypoxia and reoxygenation, this study could provide a starting point for further studies to better understand the genetic background of hypoxia stress.

Genome-wide identification of novel ovarian-predominant miRNAs: new insights from the medaka (Oryzias latipes)

MicroRNAs (miRNAs) are small, highly conserved non-coding RNAs that play important roles in the regulation of many physiological processes. However, the role of miRNAs in vertebrate oocyte formation (i.e., oogenesis) remains poorly investigated. To gain new insights into the roles of miRNAs in oogenesis, we searched for ovarian-predominant miRNAs. Using a microarray displaying 3,800 distinct miRNAs originating from different vertebrate species, we identified 66 miRNAs that are expressed predominantly in the ovary. Of the miRNAs exhibiting the highest overabundance in the ovary, 20 were selected for further analysis. Using a combination of QPCR and in silico analyses, we identified 8 novel miRNAs that are predominantly expressed in the ovary, including 2 miRNAs (miR-4785 and miR-6352) that exhibit strict ovarian expression. Of these 8 miRNAs, 7 were previously uncharacterized in fish. The strict ovarian expression of miR-4785 and miR-6352 suggests an important role in oogenesis and/or early development, possibly involving a maternal effect. Together, these results indicate that, similar to protein-coding genes, a significant number of ovarian-predominant miRNA genes are found in fish.

Discovery and functional characterization of novel miRNAs in the marine medaka Oryzias melastigma

The marine medaka Oryzias melastigma has often been used as a marine fish model to investigate the biological responses to environmental stresses and pollutants in marine environments. miRNAs are post-transcriptional regulators of many biological processes in a variety of organisms, and have been shown to be affected by environmental stresses, but the novel miRNA profile of marine medaka has not been reported. Using both genome and small RNA sequencings coupled with different bioinformatics analyses, we have discovered 58, 82, 234, and 201 unannotated miRNAs in the brain, liver, ovary and testis tissues of marine medaka, respectively. Furthermore, these novel miRNAs were found to target genes with tissue-specific roles such as neuron development and synaptic transmission in the brain, glucose and fat metabolism in the liver and steroidogenesis in the gonads. We here report, for the first time, novel miRNA profile of marine medaka, which will provide a foundation for future biomarkers and transgenerational studies for the assessment of environmental stresses and pollutions in the marine environments. In a boarder context, our data will provide novel insight into our knowledge of miRNome and miR research.

Prediction and characterization of microRNAs from eleven fish species by computational methods

MicroRNAs (miRNAs) are a family of single-stranded RNA molecules about 22 nt in length, which can regulate protein-coding gene expression in various organisms by post-transcriptional repression of messenger. In this research, the potential miRNAs and their target genes were analyzed and predicted by computational methods from the EST and GSS databases of eleven fish species, 43 potential miRNAs were identified, they belong to 38 miRNA families, some miRNAs are highly conserved in animal kingdom, the predicted target genes are involved in development, signal transduction, response to environmental stress and pathogen invasion. Taken together, our data suggest that there are a plentiful of miRNAs in these eleven fish species, these miRNAs may play some important roles by regulating their target genes, and the data provide important information for further functional studies.