The Roles of CircRNAs in Mitochondria

Mitochondria participate in varieties of cellular events. It is widely accepted that human mitochondrial genome encodes 13 proteins, 2 rRNAs, and 22 tRNAs. Gene variation derived from human nuclear genome cannot completely explain mitochondrial diseases. The advent of high-throughput sequencing coupled with novel bioinformatic analyses decode the complexity of mitochondria-derived transcripts. Recently, circular RNAs (circRNAs) from both human mitochondrial genome and nuclear genome have been found to be located at mitochondria. Studies about the roles and molecular mechanisms underlying trafficking of the nucleus encoded circRNAs to mitochondria and mitochondria encoded circRNAs to the nucleus or cytoplasm in mammals are only beginning to emerge. These circRNAs have been associated with a variety of diseases, especially cancers. Here, we discuss the emerging field of mitochondria-located circRNAs by reviewing their identification, expression patterns, regulatory roles, and functional mechanisms. Mitochondria-located circRNAs have regulatory roles in cellular physiology and pathology. We also highlight future perspectives and challenges in studying mitochondria-located circRNAs, as well as their potential biomedical applications.

Large variance in inbreeding within the Iberian wolf population

The gray wolf (Canis lupus) population on the Iberian Peninsula was the largest in western and central Europe during most of the 20th century, with its size apparently never under a few hundred individuals. After partial legal protection in the 1970s in Spain, the northwest Iberian population increased to about 300-350 packs and then stabilized. In contrast to many current European wolf populations, which have been connected through gene flow, the Iberian wolf population has been isolated for decades. Here we measured changes on genomic diversity and inbreeding through the last decades in a geographic context. We find that the level of genomic diversity in Iberian wolves is low compared to other Eurasian wolf populations. Despite population expansion in the last 50 years, some modern wolves had very high inbreeding, especially in the recently recolonized and historical edge areas. These individuals contrast with others with low inbreeding within the same population. The high variance in inbreeding despite population expansion seems associated with small-scale fragmentation of the range that is revealed by the genetic similarity between modern and historical samples from close localities despite being separated by decades, remaining differentiated from other individuals that are just over 100 km away, a small distance for a species with great dispersal capacity inhabiting a continuous range. This illustrates that, despite its demographically stable condition, the population would probably benefit from favoring connectivity within the population as well as genetic exchange with other European wolf populations to avoid excessive fragmentation and local inbreeding depression.

Bridging the Gap: Combining Genomics and Transcriptomics Approaches to Understand Stylosanthes scabra, an Orphan Legume from the Brazilian Caatinga

Stylosanthes scabra is a scientifically orphaned legume found in the Brazilian Caatinga biome (a semi-arid environment). This work utilized omics approaches to investigate some ecophysiological aspects of stress tolerance/resistance in S. scabra, study its genomic landscape, and predict potential metabolic pathways. Considering its high-confidence conceptual proteome, 1694 (~2.6%) proteins were associated with resistance proteins, some of which were found in soybean QTL regions that confer resistance to Asian soybean rust. S. scabra was also found to be a potential source of terpenes, as biosynthetic gene clusters associated with terpene biosynthesis were identified in its genome. The analysis revealed that mobile elements comprised approximately 59% of the sequenced genome. In the remaining 41% of the sections, some of the 22,681 protein-coding gene families were categorized into two informational groups: those that were specific to S. scabra and those that expanded significantly compared to their immediate ancestor. Biological process enrichment analyses indicated that these gene families play fundamental roles in the adaptation of S. scabra to extreme environments. Additionally, phylogenomic analysis indicated a close evolutionary relationship between the genera Stylosanthes and Arachis. Finally, this study found a high number (57) of aquaporin-encoding loci in the S. scabra genome. RNA-Seq and qPCR data suggested that the PIP subfamily may play a key role in the species' adaptation to water deficit conditions. Overall, these results provide valuable insights into S. scabra biology and a wealth of gene/transcript information for future legume omics studies.

Characteristics of the spiny dogfish (Squalus acanthias) nuclear genome

Sequenced shark nuclear genomes are underrepresented, with reference genomes available for only four out of nine orders so far. Here, we present the nuclear genome, with annotations, of the spiny dogfish (Squalus acanthias), a shark of interest to biomedical and conservation efforts, and the first representative of the second largest order of sharks (Squaliformes) with nuclear genome annotations available. Using Pacific Biosciences Continuous Long Read data in combination with Illumina paired-end and Hi-C sequencing, we assembled the genome de novo, followed by RNA-Seq-supported annotation. The final chromosome-level assembly is 3.7 Gb in size, has a BUSCO completeness score of 91.6%, and an error rate of less than 0.02%. Annotation predicted 33,283 gene models in the spiny dogfish's genome, of which 31,979 are functionally annotated.

Surrounded by Kindred: Spermophilus major Hybridization with Other Spermophilus Species in Space and Time

Among the numerous described cases of hybridization in mammals, the most intriguing are (a) cases of introgressive hybridization deeply affecting the evolutionary history of species, and (b) models involving not a pair of species but a multi-species complex. Therefore, the hybridization history of the russet ground squirrel Spermophilus major, whose range has repeatedly changed due to climatic fluctuations and now borders the ranges of four related species, is of great interest. The main aims of this study were to determine the direction and intensity of gene introgression, the spatial depth of the infiltration of extraneous genes into the S. major range, and to refine the hypothesis of the hybridogenic replacement of mitochondrial genomes in the studied group. Using phylogenetic analysis of the variability of mitochondrial (CR, cytb) and nuclear (SmcY, BGN, PRKCI, c-myc, i6p53) markers, we determined the contribution of neighboring species to the S. major genome. We showed that 36% of S. major individuals had extraneous alleles. All peripheral species that were in contact with S. major contributed towards its genetic variability. We also proposed a hypothesis for the sequence and localization of serial hybridization events. Our assessment of the S. major genome implications of introgression highlights the importance of implementing conservation measures to protect this species.

Combining methods for non-invasive fecal DNA enables whole genome and metagenomic analyses in wildlife biology

Non-invasive biological samples benefit studies that investigate rare, elusive, endangered, or dangerous species. Integrating genomic techniques that use non-invasive biological sampling with advances in computational approaches can benefit and inform wildlife conservation and management. Here, we used non-invasive fecal DNA samples to generate low- to medium-coverage genomes (e.g., >90% of the complete nuclear genome at six X-fold coverage) and metagenomic sequences, combining widely available and accessible DNA collection cards with commonly used DNA extraction and library building approaches. DNA preservation cards are easy to transport and can be stored non-refrigerated, avoiding cumbersome or costly sample methods. The genomic library construction and shotgun sequencing approach did not require enrichment or targeted DNA amplification. The utility and potential of the data generated was demonstrated through genome scale and metagenomic analyses of zoo and free-ranging African savanna elephants (Loxodonta africana). Fecal samples collected from free-ranging individuals contained an average of 12.41% (5.54-21.65%) endogenous elephant DNA. Clustering of these elephants with others from the same geographic region was demonstrated by a principal component analysis of genetic variation using nuclear genome-wide SNPs. Metagenomic analyses identified taxa that included Loxodonta, green plants, fungi, arthropods, bacteria, viruses and archaea, showcasing the utility of this approach for addressing complementary questions based on host-associated DNA, e.g., pathogen and parasite identification. The molecular and bioinformatic analyses presented here contributes towards the expansion and application of genomic techniques to conservation science and practice.

De Novo Long-Read Whole-Genome Assemblies and the Comparative Pan-Genome Analysis of Ascochyta Blight Pathogens Affecting Field Pea

Ascochyta Blight (AB) is a major disease of many cool-season legumes globally. In field pea, three fungal pathogens have been identified to be responsible for this disease in Australia, namely Peyronellaea pinodes, Peyronellaea pinodella and Phoma koolunga. Limited genomic resources for these pathogens have been generated, which has hampered the implementation of effective management strategies and breeding for resistant cultivars. Using Oxford Nanopore long-read sequencing, we report the first high-quality, fully annotated, near-chromosome-level nuclear and mitochondrial genome assemblies for 18 isolates from the Australian AB complex. Comparative genome analysis was performed to elucidate the differences and similarities between species and isolates using phylogenetic relationships and functional diversity. Our data indicated that P. pinodella and P. koolunga are heterothallic, while P. pinodes is homothallic. More homology and orthologous gene clusters are shared between P. pinodes and P. pinodella compared to P. koolunga. The analysis of the repetitive DNA content showed differences in the transposable repeat composition in the genomes and their expression in the transcriptomes. Significant repeat expansion in P. koolunga’s genome was seen, with strong repeat-induced point mutation (RIP) activity being evident. Phylogenetic analysis revealed that genetic diversity can be exploited for species marker development. This study provided the much-needed genetic resources and characterization of the AB species to further drive research in key areas such as disease epidemiology and host–pathogen interactions.

Nuclear genome of a pedinophyte pinpoints genomic innovation and streamlining in the green algae

The genomic diversity underpinning high ecological and species diversity in the green algae (Chlorophyta) remains little known. Here, we aimed to track genome evolution in the Chlorophyta, focusing on loss and gain of homologous genes, and lineage-specific innovations of the core Chlorophyta. We generated a high-quality nuclear genome for pedinophyte YPF701, a sister lineage to others in the core Chlorophyta and incorporated this genome in a comparative analysis with 25 other genomes from diverse Viridiplantae taxa. The nuclear genome of pedinophyte YPF701 has an intermediate size and gene number between those of most prasinophytes and the remainder of the core Chlorophyta. Our results suggest positive selection for genome streamlining in the Pedinophyceae, independent from genome minimisation observed among prasinophyte lineages. Genome expansion was predicted along the branch leading to the UTC clade (classes Ulvophyceae, Trebouxiophyceae and Chlorophyceae) after divergence from their last common ancestor with pedinophytes, with genomic novelty implicated in a range of basic biological functions. Results emphasise multiple independent signals of genome minimisation within the Chlorophyta, as well as the genomic novelty arising before diversification in the UTC clade, which may underpin the success of this species-rich clade in a diversity of habitats.

De novo chromosome level assembly of a plant genome from long read sequence data

Recent advances in the sequencing and assembly of plant genomes have allowed the generation of genomes with increasing contiguity and sequence accuracy. Chromosome level genome assemblies using sequence contigs generated from long read sequencing have involved the use of proximity analysis (Hi-C) or traditional genetic maps to guide the placement of sequence contigs within chromosomes. The development of highly accurate long reads by repeated sequencing of circularized DNA (HiFi; PacBio) has greatly increased the size of contigs. We now report the use of HiFiasm to assemble the genome of Macadamia jansenii, a genome that has been used as a model to test sequencing and assembly. This achieved almost complete chromosome level assembly from the sequence data alone without the need for higher level chromosome map information. Eight of the 14 chromosomes were represented by a single large contig (six with telomere repeats at both ends) and the other six assembled from two to four main contigs. The small number of chromosome breaks appears to be the result of highly repetitive regions including ribosomal genes that cannot be assembled by these approaches. De novo assembly of near complete chromosome level plant genomes now appears possible using these sequencing and assembly tools. Further targeted strategies might allow these remaining gaps to be closed.

The SNP-Based Profiling of Montecristo Feral Goat Populations Reveals a History of Isolation, Bottlenecks, and the Effects of Management

The Montecristo wild goat is an endangered feral population that has been on the homonymous island in the Tuscan Archipelago since ancient times. The origins of Montecristo goats are still debated, with authors dating their introduction either back to Neolithic times or between the 6th and 13th century of the Common Era. To investigate the evolutionary history and relationships of this population we assembled a 50K SNP dataset including 55 Mediterranean breeds and two nuclei of Montecristo goats sampled on the island and from an ex situ conservation project. Diversity levels, gene flow, population structure, and genetic relationships were assessed through multiple approaches. The insular population scored the lowest values of both observed and expected heterozygosity, highlighting reduced genetic variation, while the ex situ nucleus highlighted a less severe reduction. Multivariate statistics, network, and population structure analyses clearly separated the insular nucleus from all other breeds, including the population of Montecristo goats from the mainland. Moreover, admixture and gene flow analyses pinpointed possible genetic inputs received by the two Montecristo goat nuclei from different sources, while Runs of Homozygosity (ROHs) indicated an ancient bottleneck/founder effect in the insular population and recent extensive inbreeding in the ex situ one. Overall, our results suggest that Montecristo goats experienced several demographic fluctuations combined with admixture events over time and highlighted a noticeable differentiation between the two nuclei.

[Risk of Multiple Sclerosis: Analysis of Interactions Between Variants of Nuclear and Mitochondrial Genomes]

There is increasing evidence that the interaction of the mitochondrial and nuclear genomes substantially affects the risk of neurodegenerative diseases. The role of mitonuclear interactions in the development of multiple sclerosis, a severe chronic neurodegenerative disease of a polygenic nature, is poorly understood. In this work, we analyzed the association of multiple sclerosis with two-component mitonuclear combinations that include each of seven polymorphic variants of the nuclear genome localized in the region of the UCP2, and KIF1B genes and in the PVT1 locus (MYC, PVT1, and MIR1208 genes) and each often polymorphisms of the mitochondrial genome, as well as individual genetic variants that make up these combinations. Association of the individual components of these combinations with multiple sclerosis was also evaluated. 507 patients with multiple sclerosis and 321 healthy individuals were enrolled in the study, all participants were ethnic Russians. Two mitonuclear combinations associated with multiple sclerosis were identified: the UCP2 (rs660339)*A + MT-ATP6 (rs193303045)*G combination was characterized by p-value = 0.015 and OR= 1.39 [95% CI 1.05-1.87], and the PVT1 (rs2114358)*G + MT-ND1 (rs1599988)*С combination - by p-value = 0.012 and OR = 1.77 [95% CI 1.10-2.84]. Only one of the individual components of these combinations, allele rs660339*A of the nuclear gene UCP2 encoding uncoupling protein 2 of the mitochondrial anion carrier family, was independently associated with multiple sclerosis (p = 0.028; OR = 1.36 [95% CI 1.01-1.84]). This study expands the current understanding of the role of mitonuclear interactions and variance of nuclear genes, whose products function in mitochondria, and in risk of MS.

[The possibility of using multiple sclerosis-associated variants of the mitochondrial genome to predict the development of multiple sclerosis]

Multiple sclerosis (MS) is a chronic disease of the central nervous system characterized by autoimmune inflammation, demyelination, and neurodegeneration. MS is a complex disease that develops under the influence of environmental factors in genetically predisposed individuals. Currently, more than 200 genetic loci associated with MS have been identified by various methods. Some of them are located in the mitochondrial DNA. This paper collects data on mtDNA variants associated with MS in the Russian ethnic group, and shows the possibility of using this information to construct and refine models for predicting the development of MS.

Morphological stasis masks ecologically divergent coral species on tropical reefs

Coral reefs are the epitome of species diversity, yet the number of described scleractinian coral species, the framework-builders of coral reefs, remains moderate by comparison. DNA sequencing studies are rapidly challenging this notion by exposing a wealth of undescribed diversity, but the evolutionary and ecological significance of this diversity remains largely unclear. Here, we present an annotated genome for one of the most ubiquitous corals in the Indo-Pacific (Pachyseris speciosa) and uncover, through a comprehensive genomic and phenotypic assessment, that it comprises morphologically indistinguishable but ecologically divergent lineages. Demographic modeling based on whole-genome resequencing indicated that morphological crypsis (across micro- and macromorphological traits) was due to ancient morphological stasis rather than recent divergence. Although the lineages occur sympatrically across shallow and mesophotic habitats, extensive genotyping using a rapid molecular assay revealed differentiation of their ecological distributions. Leveraging "common garden" conditions facilitated by the overlapping distributions, we assessed physiological and quantitative skeletal traits and demonstrated concurrent phenotypic differentiation. Lastly, spawning observations of genotyped colonies highlighted the potential role of temporal reproductive isolation in the limited admixture, with consistent genomic signatures in genes related to morphogenesis and reproduction. Overall, our findings demonstrate the presence of ecologically and phenotypically divergent coral species without substantial morphological differentiation and provide new leads into the potential mechanisms facilitating such divergence. More broadly, they indicate that our current taxonomic framework for reef-building corals may be scratching the surface of the ecologically relevant diversity on coral reefs, consequently limiting our ability to protect or restore this diversity effectively.

Genomic adaptations to an endolithic lifestyle in the coral-associated alga Ostreobium

The green alga Ostreobium is an important coral holobiont member, playing key roles in skeletal decalcification and providing photosynthate to bleached corals that have lost their dinoflagellate endosymbionts. Ostreobium lives in the coral's skeleton, a low-light environment with variable pH and O₂ availability. We present the Ostreobium nuclear genome and a metatranscriptomic analysis of healthy and bleached corals to improve our understanding of Ostreobium's adaptations to its extreme environment and its roles as a coral holobiont member. The Ostreobium genome has 10,663 predicted protein-coding genes and shows adaptations for life in low and variable light conditions and other stressors in the endolithic environment. This alga presents a rich repertoire of light-harvesting complex proteins but lacks many genes for photoprotection and photoreceptors. It also has a large arsenal of genes for oxidative stress response. An expansion of extracellular peptidases suggests that Ostreobium may supplement its energy needs by feeding on the organic skeletal matrix, and a diverse set of fermentation pathways allows it to live in the anoxic skeleton at night. Ostreobium depends on other holobiont members for vitamin B12, and our metatranscriptomes identify potential bacterial sources. Metatranscriptomes showed Ostreobium becoming a dominant agent of photosynthesis in bleached corals and provided evidence for variable responses among coral samples and different Ostreobium genotypes. Our work provides a comprehensive understanding of the adaptations of Ostreobium to its extreme environment and an important genomic resource to improve our comprehension of coral holobiont resilience, bleaching, and recovery.

Some Molecular and Cellular Stress Mechanisms Associated with Neurodegenerative Diseases and Atherosclerosis

Chronic stress is a combination of nonspecific adaptive reactions of the body to the influence of various adverse stress factors which disrupt its homeostasis, and it is also a corresponding state of the organism's nervous system (or the body in general). We hypothesized that chronic stress may be one of the causes occurence of several molecular and cellular types of stress. We analyzed literary sources and considered most of these types of stress in our review article. We examined genes and mutations of nuclear and mitochondrial genomes and also molecular variants which lead to various types of stress. The end result of chronic stress can be metabolic disturbance in humans and animals, leading to accumulation of reactive oxygen species (ROS), oxidative stress, energy deficiency in cells (due to a decrease in ATP synthesis) and mitochondrial dysfunction. These changes can last for the lifetime and lead to severe pathologies, including neurodegenerative diseases and atherosclerosis. The analysis of literature allowed us to conclude that under the influence of chronic stress, metabolism in the human body can be disrupted, mutations of the mitochondrial and nuclear genome and dysfunction of cells and their compartments can occur. As a result of these processes, oxidative, genotoxic, and cellular stress can occur. Therefore, chronic stress can be one of the causes forthe occurrence and development of neurodegenerative diseases and atherosclerosis. In particular, chronic stress can play a large role in the occurrence and development of oxidative, genotoxic, and cellular types of stress.

Aquatic Insects Are Dramatically Underrepresented in Genomic Research

Aquatic insects comprise 10% of all insect diversity, can be found on every continent except Antarctica, and are key components of freshwater ecosystems. However, aquatic insect genome biology lags dramatically behind that of terrestrial insects. If genomic effort was spread evenly, one aquatic insect genome would be sequenced for every ~9 terrestrial insect genomes. Instead, ~24 terrestrial insect genomes have been sequenced for every aquatic insect genome. This discrepancy is even more dramatic if the quality of genomic resources is considered; for instance, while no aquatic insect genome has been assembled to the chromosome level, 29 terrestrial insect genomes spanning four orders have. We argue that a lack of aquatic insect genomes is not due to any underlying difficulty (e.g., small body sizes or unusually large genomes), yet it is severely hampering aquatic insect research at both fundamental and applied scales. By expanding the availability of aquatic insect genomes, we will gain key insight into insect diversification and empower future research for a globally important taxonomic group.

Cellular Mechanisms of Human Atherogenesis: Focus on Chronification of Inflammation and Mitochondrial Mutations

Atherosclerosis is one of the most common diseases of the cardiovascular system that leads to the development of life-threatening conditions, such as heart attack and stroke. Arthrosclerosis affects various arteries in the human body, but is especially dangerous in the arteries alimenting heart and brain, aorta, and arteries of the lower limbs. By its pathophysiology, atherosclerosis is an inflammatory disease. During the pathological process, lesions of arterial intima in the form of focal thickening are observed, which form atherosclerotic plaques as the disease progresses further. Given the significance of atherosclerosis for the global health, the search for novel effective therapies is highly prioritized. However, despite the constant progress, our understanding of the mechanisms of atherogenesis is still incomplete. One of the remaining puzzles in atherosclerosis development is the focal distribution of atherosclerotic lesions in the arterial wall. It implies the existence of certain mosaicism within the tissue, with some areas more susceptible to disease development than others, which may prove to be important for novel therapy development. There are many hypotheses explaining this phenomenon, for example, the influence of viruses, and the spread in the endothelium of the vessel multinucleated giant endothelial cells. We suggest the local variations of the mitochondrial genome as a possible explanation of this mosaicism. In this review, we discuss the role of genetic variations in the nuclear and mitochondrial genomes that influence the development of atherosclerosis. Changes in the mitochondrial and nuclear genome have been identified as independent factors for the development of the disease, as well as potential diagnostic markers.

Subchromosome-Scale Nuclear and Complete Mitochondrial Genome Characteristics of Morchella crassipes

Morchella crassipes (Vent.) Pers., a typical yellow morel species with high economic value, is mainly distributed in the low altitude plains of Eurasia. However, rare research has been performed on its genomics and polarity, thus limiting its research and development. Here, we reported a fine physical map of the nuclear genome at the subchromosomal-scale and the complete mitochondrial genome of M. crassipes. The complete size of the nuclear genome was 56.7 Mb, and 23 scaffolds were assembled, with eight of them being complete chromosomes. A total of 11,565 encoding proteins were predicted. The divergence time analysis showed that M. crassipes representing yellow morels differentiated with black morels at ~33.98 Mya (million years), with 150 gene families contracted and expanded in M. crassipes versus the two black morels (M. snyderi and M. importuna). Furthermore, 409 CAZYme genes were annotated in M. crassipes, containing almost all plant cell wall degrading enzymes compared with the mycorrhizal fungi (truffles). Genomic annotation of mating type loci and amplification of the mating genes in the monospore population was conducted, the results indicated that M. crassipes is a heterothallic fungus. Additionally, a complete circular mitochondrial genome of M. crassipes was assembled, the size reached as large as 531,195 bp. It can be observed that the strikingly large size was the biggest up till now, coupled with 14 core conserved mitochondrial protein-coding genes, two rRNAs, 31 tRNAs, 51 introns, and 412 ncORFs. The total length of intron sequences accounted for 53.67% of the mitochondrial genome, with 19 introns having a length over 5 kb. Particularly, 221 of 412 ncORFs were distributed within 51 introns, and the total length of the ncORFs sequence accounted for 40.83% of the mitochondrial genome, and 297 ncORFs had expression activity in the mycelium stage, suggesting their potential functions in M. crassipes. Meanwhile, there was a high degree of repetition (51.31%) in the mitochondria of M. crassipes. Thus, the large number of introns, ncORFs and internal repeat sequences may contribute jointly to the largest fungal mitochondrial genome to date. The fine physical maps of nuclear genome and mitochondrial genome obtained in this study will open a new door for better understanding of the mysterious species of M. crassipes.

Relationships between Iraqi Rice Varieties at the Nuclear and Plastid Genome Levels

Due to the importance of the rice crop in Iraq, this study was conducted to determine the origin of the major varieties and understand the evolutionary relationships between Iraqi rice varieties and other Asian rice accessions that could be significant in the improvement of this crop. Five varieties of Oryza sativa were obtained from Baghdad/Iraq, and the whole genomic DNA was sequenced, among these varieties, Amber33, Furat, Yasmin, Buhooth1 and Amber al-Baraka. Raw sequence reads of 33 domesticated Asian rice accessions were obtained from the Sequence Read Archive (SRA-NCBI). The sequence of the whole chloroplast-genome was assembled while only the sequence of 916 concatenated nuclear-genes was assembled. The phylogenetic analysis of both chloroplast and nuclear genomes showed that two main clusters, Indica and Japonica, and further five sub-clusters based upon their ecotype, indica, aus, tropical-japonica, temperate-japonica and basmati were created; moreover, Amber33, Furat, Yasmin and Buhooth1 belonged to the basmati, indica and japonica ecotypes, respectively, where Amber33 was placed in the basmati group as a sister of cultivars from Pakistan and India. This confirms the traditional story that Amber was transferred by a group of people who had migrated from India and settled in southern Iraq a long time ago.

Mitochondrial disease genetics update: recent insights into the molecular diagnosis and expanding phenotype of primary mitochondrial disease

PURPOSE OF REVIEW

Primary mitochondrial disease (PMD) is a genetically and phenotypically diverse group of inherited energy deficiency disorders caused by impaired mitochondrial oxidative phosphorylation (OXPHOS) capacity. Mutations in more than 350 genes in both mitochondrial and nuclear genomes are now recognized to cause primary mitochondrial disease following every inheritance pattern. Next-generation sequencing technologies have dramatically accelerated mitochondrial disease gene discovery and diagnostic yield. Here, we provide an up-to-date review of recently identified, novel mitochondrial disease genes and/or pathogenic variants that directly impair mitochondrial structure, dynamics, and/or function.

RECENT FINDINGS

A review of PubMed publications was performed from the past 12 months that identified 16 new PMD genes and/or pathogenic variants, and recognition of expanded phenotypes for a wide variety of mitochondrial disease genes.

SUMMARY

Broad-based exome sequencing has become the standard first-line diagnostic approach for PMD. This has facilitated more rapid and accurate disease identification, and greatly expanded understanding of the wide spectrum of potential clinical phenotypes. A comprehensive dual-genome sequencing approach to PMD diagnosis continues to improve diagnostic yield, advance understanding of mitochondrial physiology, and provide strong potential to develop precision therapeutics targeted to diverse aspects of mitochondrial disease pathophysiology.

Variants of Mitochondrial Genome and Risk of Multiple Sclerosis Development in Russians

For the first time in the history of ethnic Russians, an association analysis the development of multiple sclerosis (MS) was performed for the mitochondrial haplogroups H, J, K, and U, as well as for the individual mitochondrial DNA (mtDNA) polymorphisms discriminating these haplogroups (m.1719G > A, m. 7028C > T, m.9055G > A, m.10398A > G, m.12308A > G). A total of 283 unrelated patients with the relapsing-remitting form of MS and 290 healthy controls were enrolled in the study. Association of haplogroup J with MS was observed (P = 0.0055, OR = 2.00 [95% CI 1.21-3.41]). After gender stratification, the association remained significant in women (P = 0.0083, OR = 2.20 [95% CI 1.19-4.03]). A multilocus analysis of the association between combinations of mtDNA haplogroups with variants of 38 nuclear immune-related genes and MS risk was carried out. MS-associated biallelic combinations of haplogroup J with the alleles CCL5 rs2107538*A, PVT1 rs2114358*G, TNFSF14 rs1077667*C, and IL4 rs2243250*C, which were not associated with MS individually, were identified. For the combination of haplogroup J and the CCL5*A allele (P = 0.00043, OR = 5.47 [95% CI 1.85-16.15]), a epistatic (synergistic) interaction between the components was established using two statistical criteria: the PFLINT value in the Fisher-like interaction numeric test and the synergy factor, SF (PFLINT = 0.025, SF = 4.32 [95% CI 1.20-15.60]). The combination of haplogroup J and the PVT1*G allele is characterized by PFLINT = 0.084; SF = 3.05 [95% CI 1.00-9.31] and can also be epistatic. Thus, interaction between nuclear and mitochondrial genome components in the risk of developing MS was demonstrated for the first time.

Northern Spotted Owl (Strix occidentalis caurina) Genome: Divergence with the Barred Owl (Strix varia) and Characterization of Light-Associated Genes

We report here the assembly of a northern spotted owl (Strix occidentalis caurina) genome. We generated Illumina paired-end sequence data at 90× coverage using nine libraries with insert lengths ranging from ∼250 to 9,600 nt and read lengths from 100 to 375 nt. The genome assembly is comprised of 8,108 scaffolds totaling 1.26 × 109 nt in length with an N50 length of 3.98 × 106 nt. We calculated the genome-wide fixation index (FST) of S. o. caurina with the closely related barred owl (Strix varia) as 0.819. We examined 19 genes that encode proteins with light-dependent functions in our genome assembly as well as in that of the barn owl (Tyto alba). We present genomic evidence for loss of three of these in S. o. caurina and four in T. alba. We suggest that most light-associated gene functions have been maintained in owls and their loss has not proceeded to the same extent as in other dim-light-adapted vertebrates.

Comparative analysis of single nucleotide polymorphisms in the nuclear, chloroplast, and mitochondrial genomes in identification of phylogenetic association among seven melon (Cucumis melo L.) cultivars

A variety of melons are cultivated worldwide, and their specific biological properties make them an attractive model for molecular studies. This study aimed to investigate the single nucleotide polymorphisms (SNPs) from the mitochondrial, chloroplast, and nuclear genomes of seven melon accessions (Cucumis melo L.) to identify the phylogenetic relationships among melon cultivars with the Illumina HiSeq 2000 platform and bioinformatical analyses. The data showed that there were a total of 658 mitochondrial SNPs (207-295 in each), while there were 0-60 chloroplast SNPs among these seven melon cultivars, compared to the reference genome. Bioinformatical analysis showed that the mitochondrial tree topology was unable to separate the melon features, whereas the maximum parsimony/neighbor joining (MP/NJ) tree of the chloroplast SNPs could define melon features such as seed length, width, thickness, 100-seed weight, and type. SNPs of the nuclear genome were better than the mitochondrial and chloroplast SNPs in the identification of melon features. The data demonstrated the usefulness of mitochondrial, chloroplast, and nuclear SNPs in identification of phylogenetic associations among these seven melon cultivars.

Molecular techniques to interrogate and edit the Chlamydomonas nuclear genome

The success of the green alga Chlamydomonas reinhardtii as a model organism is to a large extent due to the wide range of molecular techniques that are available for its characterization. Here, we review some of the techniques currently used to modify and interrogate the C. reinhardtii nuclear genome and explore several technologies under development. Nuclear mutants can be generated with ultraviolet (UV) light and chemical mutagens, or by insertional mutagenesis. Nuclear transformation methods include biolistic delivery, agitation with glass beads, and electroporation. Transforming DNA integrates into the genome at random sites, and multiple strategies exist for mapping insertion sites. A limited number of studies have demonstrated targeted modification of the nuclear genome by approaches such as zinc-finger nucleases and homologous recombination. RNA interference is widely used to knock down expression levels of nuclear genes. A wide assortment of transgenes has been successfully expressed in the Chlamydomonas nuclear genome, including transformation markers, fluorescent proteins, reporter genes, epitope tagged proteins, and even therapeutic proteins. Optimized expression constructs and strains help transgene expression. Emerging technologies such as the CRISPR/Cas9 system, high-throughput mutant identification, and a whole-genome knockout library are being developed for this organism. We discuss how these advances will propel future investigations.

Nuclear pore proteins and the control of genome functions

Nuclear pore complexes (NPCs) are composed of several copies of ∼30 different proteins called nucleoporins (Nups). NPCs penetrate the nuclear envelope (NE) and regulate the nucleocytoplasmic trafficking of macromolecules. Beyond this vital role, NPC components influence genome functions in a transport-independent manner. Nups play an evolutionarily conserved role in gene expression regulation that, in metazoans, extends into the nuclear interior. Additionally, in proliferative cells, Nups play a crucial role in genome integrity maintenance and mitotic progression. Here we discuss genome-related functions of Nups and their impact on essential DNA metabolism processes such as transcription, chromosome duplication, and segregation.

Detection and decay rates of prey and prey symbionts in the gut of a predator through metagenomics

DNA methods are useful to identify ingested prey items from the gut of predators, but reliable detection is hampered by low amounts of degraded DNA. PCR-based methods can retrieve minute amounts of starting material but suffer from amplification biases and cross-reactions with the predator and related species genomes. Here, we use PCR-free direct shotgun sequencing of total DNA isolated from the gut of the harlequin ladybird Harmonia axyridis at five time points after feeding on a single pea aphid Acyrthosiphon pisum. Sequence reads were matched to three reference databases: Insecta mitogenomes of 587 species, including H. axyridis sequenced here; A. pisum nuclear genome scaffolds; and scaffolds and complete genomes of 13 potential bacterial symbionts. Immediately after feeding, multicopy mtDNA of A. pisum was detected in tens of reads, while hundreds of matches to nuclear scaffolds were detected. Aphid nuclear DNA and mtDNA decayed at similar rates (0.281 and 0.11 h(-1) respectively), and the detectability periods were 32.7 and 23.1 h. Metagenomic sequencing also revealed thousands of reads of the obligate Buchnera aphidicola and facultative Regiella insecticola aphid symbionts, which showed exponential decay rates significantly faster than aphid DNA (0.694 and 0.80 h(-1) , respectively). However, the facultative aphid symbionts Hamiltonella defensa, Arsenophonus spp. and Serratia symbiotica showed an unexpected temporary increase in population size by 1-2 orders of magnitude in the predator guts before declining. Metagenomics is a powerful tool that can reveal complex relationships and the dynamics of interactions among predators, prey and their symbionts.

Polar bears exhibit genome-wide signatures of bioenergetic adaptation to life in the arctic environment

Polar bears (Ursus maritimus) face extremely cold temperatures and periods of fasting, which might result in more severe energetic challenges than those experienced by their sister species, the brown bear (U. arctos). We have examined the mitochondrial and nuclear genomes of polar and brown bears to investigate whether polar bears demonstrate lineage-specific signals of molecular adaptation in genes associated with cellular respiration/energy production. We observed increased evolutionary rates in the mitochondrial cytochrome c oxidase I gene in polar but not brown bears. An amino acid substitution occurred near the interaction site with a nuclear-encoded subunit of the cytochrome c oxidase complex and was predicted to lead to a functional change, although the significance of this remains unclear. The nuclear genomes of brown and polar bears demonstrate different adaptations related to cellular respiration. Analyses of the genomes of brown bears exhibited substitutions that may alter the function of proteins that regulate glucose uptake, which could be beneficial when feeding on carbohydrate-dominated diets during hyperphagia, followed by fasting during hibernation. In polar bears, genes demonstrating signatures of functional divergence and those potentially under positive selection were enriched in functions related to production of nitric oxide (NO), which can regulate energy production in several different ways. This suggests that polar bears may be able to fine-tune intracellular levels of NO as an adaptive response to control trade-offs between energy production in the form of adenosine triphosphate versus generation of heat (thermogenesis).

Conserved genetic regions across angiosperms as tools to develop single-copy nuclear markers in gymnosperms: an example using cycads

Several individuals of the Caribbean Zamia clade and other cycad genera were used to identify single-copy nuclear genes for phylogeographic and phylogenetic studies in Cycadales. Two strategies were employed to select target loci: (i) a tblastX search of Arabidopsis conserved ortholog sequence (COS) set and (ii) a tblastX search of Arabidopsis-Populus-Vitis-Oryza Shared Single-Copy genes (APVO SSC) against the EST Zamia databases in GenBank. From the first strategy, 30 loci were selected, and from the second, 16 loci. In both cases, the matching GenBank accessions of Zamia were used as a query for retrieving highly similar sequences from Cycas, Picea, Pinus species or Ginkgo biloba. After retrieving and aligning all the sequences in each locus, intron predictions were completed to assist in primer design. PCR was carried out in three rounds to detect paralogous loci. A total of 29 loci were successfully amplified as a single band of which 20 were likely single-copy loci. These loci showed different diversity and divergence levels. A preliminary screening allowed us to select 8 promising loci (40S, ATG2, BG, GroES, GTP, LiSH, PEX4 and TR) for the Zamia pumila complex and 4 loci (COS26, GroES, GTP and HTS) for all other cycad genera.