Intraspecific variation in the mitochondrial genome among local populations of Medaka Oryzias latipes

Exploring the mitochondrial genomes and phylogenetic relationships of trans-Andean Bryconidae species (Actinopterygii: Ostariophysi: Characiformes)

Comparative mitogenomics and its evolutionary relationships within Bryconidae remains largely unexplored. To bridge this gap, this study assembled 15 mitogenomes from 11 Bryconidae species, including five newly sequenced. Salminus mitogenomes, exceeding 17,700 bp, exhibited the largest size, contrasting with a median size of 16,848 bp in the remaining species (Brycon and Chilobrycon). These mitogenomes encode 37 typical mitochondrial genes, including 13 protein-coding, 2 ribosomal RNA, and 22 transfer RNA genes, and exhibit the conserved gene arrangement found in most fish species. Phylogenetic relationships, based on the maximum-likelihood method, revealed that the trans-Andean species (found in northwestern South America) clustered into two main sister clades. One clade comprised the trans-Andean species from the Pacific slope, Brycon chagrensis and Chilobrycon deuterodon. The other clade grouped the trans-Andean species from the Magdalena-Cauca Basin Brycon moorei and Salminus affinis, with their respective cis-Andean congeners (found in eastern South America), with Brycon rubricauda as its sister clade. Since the current members of Brycon are split in three separated lineages, the systematic classification of Bryconidae requires further examination. This study provides novel insights into mitogenome characteristics and evolutionary pathways within Bryconidae, standing as crucial information for prospective phylogenetic and taxonomic studies, molecular ecology, and provides a valuable resource for environmental DNA applications.

Assembly and comparative analysis of the complete multichromosomal mitochondrial genome of Cymbidium ensifolium, an orchid of high economic and ornamental value

BACKGROUND

Orchidaceae is one of the largest groups of angiosperms, and most species have high economic value and scientific research value due to their ornamental and medicinal properties. In China, Chinese Cymbidium is a popular ornamental orchid with high economic value and a long history. However, to date, no detailed information on the mitochondrial genome of any species of Chinese Cymbidium has been published.

RESULTS

Here, we present the complete assembly and annotation of the mitochondrial genome of Cymbidium ensifolium (L.) Sw. The mitogenome of C. ensifolium was 560,647 bp in length and consisted of 19 circular subgenomes ranging in size from 21,995 bp to 48,212 bp. The genome encoded 35 protein-coding genes, 36 tRNAs, 3 rRNAs, and 3405 ORFs. Repeat sequence analysis and prediction of RNA editing sites revealed a total of 915 dispersed repeats, 162 simple repeats, 45 tandem repeats, and 530 RNA editing sites. Analysis of codon usage showed a preference for codons ending in A/T. Interorganellar DNA transfer was identified in 13 of the 19 chromosomes, with plastid-derived DNA fragments representing 6.81% of the C. ensifolium mitochondrial genome. The homologous fragments of the mitochondrial genome and nuclear genome were also analysed. Comparative analysis showed that the GC content was conserved, but the size, structure, and gene content of the mitogenomes varied greatly among plants with multichromosomal mitogenome structure. Phylogenetic analysis based on the mitogenomes reflected the evolutionary and taxonomic statuses of C. ensifolium. Interestingly, compared with the mitogenomes of Cymbidium lancifolium Hook. and Cymbidium macrorhizon Lindl., the mitogenome of C. ensifolium lost 8 ribosomal protein-coding genes.

CONCLUSION

In this study, we assembled and annotated the mitogenome of C. ensifolium and compared it with the mitogenomes of other Liliidae and plants with multichromosomal mitogenome structures. Our findings enrich the mitochondrial genome database of orchid plants and reveal the rapid structural evolution of Cymbidium mitochondrial genomes, highlighting the potential for mitochondrial genes to help decipher plant evolutionary history.

Population genetic characteristics of Hainan medaka with whole-genome resequencing

The DMY gene is deleted in all males of the Sanya population (SY-medaka) of the Hainan medaka, Oryzias curvinotus, as recently reported by us. However, due to limited knowledge regarding their population genetic background, it is difficult to explore the possible evolutionary pathway. Herein, we resequenced the whole genome of four populations, including SY-medaka. A total of 56 mitogenomes and 32,826,105 SNPs were identified. We found that the genetic differentiation is highest between SY-medaka and the other populations. The results of the population history of the O. curvinotus suggest that the SY-medaka has been in a bottleneck period recently. Further analysis shows that SY-medaka are the most strongly affected by environmental selection. Moreover, we screened some potential genomic regions, and the genes contained in these regions may explain the potential mechanism of the selection process of the SY-medaka. In conclusion, our study can provide new clues for the adaptation process of medaka in the new environment of Sanya.

Medaka Population Genome Structure and Demographic History Described via Genotyping-by-Sequencing

Medaka is a model organism in medicine, genetics, developmental biology and population genetics. Lab stocks composed of more than 100 local wild populations are available for research in these fields. Thus, medaka represents a potentially excellent bioresource for screening disease-risk- and adaptation-related genes in genome-wide association studies. Although the genetic population structure should be known before performing such an analysis, a comprehensive study on the genome-wide diversity of wild medaka populations has not been performed. Here, we performed genotyping-by-sequencing (GBS) for 81 and 12 medakas captured from a bioresource and the wild, respectively. Based on the GBS data, we evaluated the genetic population structure and estimated the demographic parameters using an approximate Bayesian computation (ABC) framework. The genome-wide data confirmed that there were substantial differences between local populations and supported our previously proposed hypothesis on medaka dispersal based on mitochondrial genome (mtDNA) data. A new finding was that a local group that was thought to be a hybrid between the northern and the southern Japanese groups was actually an origin of the northern Japanese group. Thus, this paper presents the first population-genomic study of medaka and reveals its population structure and history based on chromosomal genetic diversity.

Comparative mitogenomic analysis of three species of periwinkles: Littorina fabalis, L. obtusata and L. saxatilis

The flat periwinkles, Littorina fabalis and L. obtusata, offer an interesting system for local adaptation and ecological speciation studies. In order to provide genomic resources for these species, we sequenced their mitogenomes together with that of the rough periwinkle L. saxatilis by means of next-generation sequencing technologies. The three mitogenomes present the typical repertoire of 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes and a putative control region. Although the latter could not be fully recovered in flat periwinkles using short-reads due to a highly repetitive fragment, in L. saxatilis this problem was overcome with additional long-reads and we were able to assemble the complete mitogenome. Both gene order and nucleotide composition are similar between the three species as well as compared to other Littorinimorpha. A large variance in divergence was observed across mitochondrial regions, with six- to ten-fold difference between the highest and the lowest divergence rates. Based on nucleotide changes on the whole molecule and assuming a molecular clock, L. fabalis and L. obtusata started to diverge around 0.8 Mya (0.4-1.1 Mya). The evolution of the mitochondrial protein-coding genes in the three Littorina species appears mainly influenced by purifying selection as revealed by phylogenetic tests based on d_N/d_S ratios that did not detect any evidence for positive selection, although some caution is required given the limited power of the dataset and the implemented approaches.

Genetic control of startle behavior in medaka fish

Genetic polymorphisms are thought to generate intraspecific behavioral diversities, both within and among populations. The mechanisms underlying genetic control of behavioral properties, however, remain unclear in wild-type vertebrates, including humans. To explore this issue, we used diverse inbred strains of medaka fish (Oryzias latipes) established from the same and different local populations. Medaka exhibit a startle response to a visual stimulus (extinction of illumination) by rapidly bending their bodies (C-start) 20-ms after the stimulus presentation. We measured the rates of the response to repeated stimuli (1-s interval, 40 times) among four inbred strains, HNI-I, HNI-II, HO5, and Hd-rR-II1, and quantified two properties of the startle response: sensitivity (response rate to the first stimulus) and attenuation of the response probability with repeated stimulus presentation. Among the four strains, the greatest differences in these properties were detected between HNI-II and Hd-rR-II1. HNI-II exhibited high sensitivity (approximately 80%) and no attenuation, while Hd-rR-II1 exhibited low sensitivity (approximately 50%) and almost complete attenuation after only five stimulus presentations. Our findings suggested behavioral diversity of the startle response within a local population as well as among different populations. Linkage analysis with F2 progeny between HNI-II and Hd-rR-II1 detected quantitative trait loci (QTL) highly related to attenuation, but not to sensitivity, with a maximum logarithm of odds score of 11.82 on linkage group 16. The three genotypes (homozygous for HNI-II and Hd-rR-II1 alleles, and heterozygous) at the marker nearest the QTL correlated with attenuation. Our findings are the first to suggest that a single genomic region might be sufficient to generate individual differences in startle behavior between wild-type strains. Further identification of genetic polymorphisms that define the behavioral trait will contribute to our understanding of the neural mechanisms underlying behavioral diversity, allowing us to investigate the adaptive significance of intraspecific behavioral polymorphisms of the startle response.

Regular heartbeat rhythm at the heartbeat initiation stage is essential for normal cardiogenesis at low temperature

BACKGROUND

The development of blood flow in the heart is crucial for heart function and embryonic survival. Recent studies have revealed the importance of the extracellular matrix and the mechanical stress applied to the valve cushion that controls blood flow to the formation of the cardiac valve during embryogenesis. However, the events that trigger such valve formation and mechanical stress, and their temperature dependence have not been explained completely. Medaka (Oryzias latipes) inhabits a wide range of East Asia and adapts to a wide range of climates. We used medaka embryos from different genomic backgrounds and analyzed heartbeat characteristics including back-and-forth blood flow and bradyarrhythmia in embryos incubated at low temperature. We also used high-speed imaging analysis to examine the heartbeat of these animals after transient exposure to low temperature.

RESULTS

Embryos of the Hd-rR medaka strain exhibited back-and-forth blood flow in the heart (blood regurgitation) after incubation at 15 °C. This regurgitation was induced by exposure to low temperature around the heartbeat initiation period and was related to abnormalities in the maintenance or pattern of contraction of the atrium or the atrioventricular canal. The Odate strain from the northern Japanese group exhibited normal blood flow after incubation at 15 °C. High-speed time-lapse analysis of the heartbeat revealed that bradyarrhythmia occurred only in Hd-rR embryos incubated at 15 °C. The coefficient of contraction, defined as the quotient of the length of the atrium at systole divided by its length at diastole, was not affected in either strain. The average heart rate after removing the effect of arrhythmia did not differ significantly between the two strains, suggesting that the mechanical stress of individual myocardial contractions and the total mechanical stress could be equivalent, regardless of the presence of arrhythmia or the heart rate. Test-cross experiments suggested that this circulation phenotype was caused by a single major genomic locus.

CONCLUSIONS

These results suggest that cardiogenesis at low temperature requires a constant heartbeat. Abnormal contraction rhythms at the stage of heartbeat initiation may cause regurgitation at later stages. From the evolutionary viewpoint, strains that exhibit normal cardiogenesis during development at low temperature inhabit northern environments.

A comprehensive description and evolutionary analysis of 22 grouper (perciformes, epinephelidae) mitochondrial genomes with emphasis on two novel genome organizations

Groupers of the family Epinephelidae are a diverse and economically valuable group of reef fishes. To investigate the evolution of their mitochondrial genomes we characterized and compared these genomes among 22 species, 17 newly sequenced. Among these fishes we identified three distinct genome organizations, two of them never previously reported in vertebrates. In 19 of these species, mitochondrial genomes followed the typical vertebrate canonical organization with 13 protein-coding genes, 22 tRNAs, two rRNAs, and a non-coding control region. Differing from this, members of genus Variola have an extra tRNA-Ile between tRNA-Val and 16S rRNA. Evidence suggests that this evolved from tRNA-Val via a duplication event due to slipped strand mispairing during replication. Additionally, Cephalopholisargus has an extra tRNA-Asp in the midst of the control region, likely resulting from long-range duplication of the canonical tRNA-Asp through illicit priming of mitochondrial replication by tRNAs. Along with their gene contents, we characterized the regulatory elements of these mitochondrial genomes' control regions, including putative termination-associated sequences and conserved sequence blocks. Looking at the mitochondrial genomic constituents, rRNA and tRNA are the most conserved, followed by protein-coding genes, and non-coding regions are the most divergent. Divergence rates vary among the protein-coding genes, and the three cytochrome oxidase subunits (COI, II, III) are the most conserved, while NADH dehydrogenase subunit 6 (ND6) and the ATP synthase subunit 8 (ATP8) are the most divergent. We then tested the phylogenetic utility of this new mt genome data using 12 protein-coding genes of 48 species from the suborder Percoidei. From this, we provide further support for the elevation of the subfamily Epinephelinae to family Epinephelidae, the resurrection of the genus Hyporthodus, and the combination of the monotypic genera Anyperodon and Cromileptes to genus Epinephelus, and Aethaloperca to genus Cephalopholis.

Structure and organization of the mitochondrial DNA control region with tandemly repeated sequence in the Amazon ornamental fish

Tandemly repeated sequences are a common feature of vertebrate mitochondrial DNA control regions. However, questions still remain about their mode of evolution and function. To better understand patterns of variation in length and to explore the existence of previously described domain, we have characterized the control region structure of the Amazonian ornamental fish Nannostomus eques and Nannostomus unifasciatus. The control region ranged from 1121 to 1142 bp in length and could be separated into three domains: the domain associated with the extended terminal associated sequences, the central conserved domain, and the conserved sequence blocks domain. In the first domain, we encountered a sequence repeated 10 times in tandem (variable number tandem repeat (VNTR)) that could adopt an "inverted repetitions" type structural conformation. The results suggest that the VNTR pattern encountered in both N. eques and N. unifasciatus is consistent with the prerequisites of the illegitimate elongation model in which the unequal pairing of the chains near the 5'-end of the control region favors the formation of repetitions.

Reassessment of the function of somatolactin alpha in lipid metabolism using medaka mutant and transgenic strains

Background

Somatolactin alpha (SLa) is a fish-specific peptide hormone secreted from the pituitary. In medaka, SLa functions to darken the skin color and lack of SLa makes it pale. Transcription of SLa is enhanced or suppressed when fish are kept in dark or bright conditions, respectively, indicating SLa’s important role in background acclimation of the skin color. Bizarrely, however, the lack of SLa seems to cause the additional defect of increased triglycerides in organs, which could not be rescued (decreased) by its overexpression.

Results

To assess this enigmatic result, we investigated genetic (the SLa, Slc45a2, r, and Y genes) and nongenetic (age, fasting, water temperature, and background color) effects on hepatic triglycerides. These experiments found that percent hepatic triglycerides quickly change in response to external/internal environments. Effects of SLa seemed to be much less obvious, although it may increase the proportion of hepatic triglycerides at least during certain breeding conditions or under certain genetic backgrounds.

Conclusions

The present results do not exclude the possibility that SLa takes part in lipid metabolism or other physiological processes. However, we suggest that skin-color regulation is the only definite role of SLa so far demonstrated in this species.

Complete mitochondrial genome of the marine medaka Oryzias melastigma (Beloniformes, Adrianichthyidae)

The complete mitochondrial genome was obtained from the assembled genome data sequenced by next-generation sequencer from the marine medaka Oryzias melastigma. The mitochondrial genome sequence was 16,864 bp in size, and the gene order and contents were identical with those of previously reported fish mitochondrial genomes. Of 13 protein-coding genes (PCGs), 4 genes (CO3, ND3, ND4, and Cytb) had incomplete stop codons. The base composition of O. melastigma mitogenome showed high A+T (59.65%) and anti-G bias (8.73%) on the 3rd position of PCGs.

Genetic divergence among invasive and native populations of Plagioscion squamosissimus (Perciformes, Sciaenidae) in Neotropical regions

The genetic divergence among invasive and native populations of Plagioscion squamosissimus from four Neotropical hydrographic basins was assessed using the hypervariable domain of the mitochondrial DNA (mtDNA) control region. Plagioscion squamosissimus is native to the neighbouring hydrographic basins of the Parnaíba and Amazon Rivers, and the latter includes the Araguaia-Tocantins drainage, but it is invasive in other basins due to introductions. The mtDNA nucleotide polymorphism supported the hypothesis that the Amazon and Parnaíba populations constitute the same species and are separated into two independent evolutionary lineages. Absence of nucleotide polymorphism was observed within and among P. squamosissimus populations invasive to the uppper and middle Paraná River basins. Nucleotide divergence was null or low comparing the Paraná invasive populations with the populations native to the Parnaíba River basin, whereas it was significantly high compared to Tocantins populations. These results ascertain that P. squamosissimus populations invasive to the upper Paraná River basin and to the middle Paraná River basin downstream of the Itaipu dam are derived from the Parnaíba River basin. The genetic data presented are potentially useful to assist further studies on P. squamosissimus taxonomic and geographic distribution, development of ecological guidelines for managing populations invasive to the upper Paraná River basin and for preservation of native fish diversity.

Variation in repeat length and heteroplasmy of the mitochondrial DNA control region along a core-edge gradient in the eastern spadefoot toad (Pelobates syriacus)

Peripheral populations are those situated at the distribution margins of a species and are often subjected to more extreme abiotic and biotic conditions than those at the core. Here, we hypothesized that shorter repeat length and fewer heteroplasmic mitochondrial DNA (mtDNA) copies, which are associated with more efficient mitochondrial function, may be related to improved survival under extreme environmental conditions. We sampled eastern spadefoot toads (mostly as tadpoles) from 43 rain pools distributed along a 300-km gradient from core to edge of the species' distribution. We show that mean pool tandem repeat length and heteroplasmy increase from edge to core, even after controlling for body size. We evaluate several alternative hypotheses and propose the Fisher hypothesis as the most likely explanation. However, additional sequential sampling and experimental studies are required to determine whether selection under extreme conditions, or alternative mechanisms, could account for the gradient in heteroplasmy and repeat length in the mtDNA control region.