Extensive allele mining discovers novel genetic diversity in the loci controlling frost tolerance in barley.

KEY MESSAGE: Exome sequencing-based allele mining for frost tolerance suggests HvCBF14 rather than CNV at Fr-H2 locus is the main responsible of frost tolerance in barley. Wild relatives, landraces and old cultivars of barley represent a reservoir of untapped and potentially important genes for crop improvement, and the recent sequencing technologies provide the opportunity to mine the existing genetic diversity and to identify new genes/alleles for the traits of interest. In the present study, we use frost tolerance and vernalization requirement as case studies to demonstrate the power of allele mining carried out on exome sequencing data generated from > 400 barley accessions. New deletions in the first intron of VRN-H1 were identified and linked to a reduced vernalization requirement, while the allelic diversity of HvCBF2a, HvCBF4b and HvCBF14 was investigated by combining the analysis of SNPs and read counts. This approach has proven very effective to identify gene paralogs and copy number variants of HvCBF2 and the HvCBF4b-HvCBF2a segment. A multiple linear regression model which considers allelic variation at these genes suggests a major involvement of HvCBF14, rather than copy number variation of HvCBF4b-HvCBF2a, in controlling frost tolerance in barley. Overall, the present study provides powerful resource and tools to discover novel alleles at relevant genes in barley.

Exome sequences and multi-environment field trials elucidate the genetic basis of adaptation in barley.

Broadening the genetic base of crops is crucial for developing varieties to respond to global agricultural challenges such as climate change. Here, we analysed a diverse panel of 371 domesticated lines of the model crop barley to explore the genetics of crop adaptation. We first collected exome sequence data and phenotypes of key life history traits from contrasting multi-environment common garden trials. Then we applied refined statistical methods, including some based on exomic haplotype states, for genotype-by-environment (G×E) modelling. Sub-populations defined from exomic profiles were coincident with barley's biology, geography and history, and explained a high proportion of trial phenotypic variance. Clear G×E interactions indicated adaptation profiles that varied for landraces and cultivars. Exploration of circadian clock-related genes, associated with the environmentally adaptive days to heading trait (crucial for the crop's spread from the Fertile Crescent), illustrated complexities in G×E effect directions, and the importance of latitudinally based genic context in the expression of large-effect alleles. Our analysis supports a gene-level scientific understanding of crop adaption and leads to practical opportunities for crop improvement, allowing the prioritisation of genomic regions and particular sets of lines for breeding efforts seeking to cope with climate change and other stresses.

Variability of mitochondrial ORFans hints at possible differences in the system of doubly uniparental inheritance of mitochondria among families of freshwater mussels (Bivalvia: Unionida).

BACKGROUND: Supernumerary ORFan genes (i.e., open reading frames without obvious homology to other genes) are present in the mitochondrial genomes of gonochoric freshwater mussels (Bivalvia: Unionida) showing doubly uniparental inheritance (DUI) of mitochondria. DUI is a system in which distinct female-transmitted and male-transmitted mitotypes coexist in a single species. In families Unionidae and Margaritiferidae, the transition from dioecy to hermaphroditism and the loss of DUI appear to be linked, and this event seems to affect the integrity of the ORFan genes. These observations led to the hypothesis that the ORFans have a role in DUI and/or sex determination. Complete mitochondrial genome sequences are however scarce for most families of freshwater mussels, therefore hindering a clear localization of DUI in the various lineages and a comprehensive understanding of the influence of the ORFans on DUI and sexual systems. Therefore, we sequenced and characterized eleven new mitogenomes from poorly sampled freshwater mussel families to gather information on the evolution and variability of the ORFan genes and their protein products. RESULTS: We obtained ten complete plus one almost complete mitogenome sequence from ten representative species (gonochoric and hermaphroditic) of families Margaritiferidae, Hyriidae, Mulleriidae, and Iridinidae. ORFan genes are present only in DUI species from Margaritiferidae and Hyriidae, while non-DUI species from Hyriidae, Iridinidae, and Mulleriidae lack them completely, independently of their sexual system. Comparisons among the proteins translated from the newly characterized ORFans and already known ones provide evidence of conserved structures, as well as family-specific features. CONCLUSIONS: The ORFan proteins show a comparable organization of secondary structures among different families of freshwater mussels, which supports a conserved physiological role, but also have distinctive family-specific features. Given this latter observation and the fact that the ORFans can be either highly mutated or completely absent in species that secondarily lost DUI depending on their respective family, we hypothesize that some aspects of the connection among ORFans, sexual systems, and DUI may differ in the various lineages of unionids.

Grapevine comparative early transcriptomic profiling suggests that Flavescence dorée phytoplasma represses plant responses induced by vector feeding in susceptible varieties.

BACKGROUND: Flavescence dorée is the most serious grapevine yellows disease in Europe. It is caused by phytoplasmas which are transmitted from grapevine to grapevine by the leafhopper Scaphoideus titanus. Differences in susceptibility among grapevine varieties suggest the existence of specific genetic features associated with resistance to the phytoplasma and/or possibly with its vector. In this work, RNA-Seq was used to compare early transcriptional changes occurring during the three-trophic interaction between the phytoplasma, its vector and the grapevine, represented by two different cultivars, one very susceptible to the disease and the other scarcely susceptible. RESULTS: The comparative analysis of the constitutive transcriptomic profiles suggests the existence of passive defense strategies against the insect and/or the phytoplasma in the scarcely-susceptible cultivar. Moreover, the attack by the infective vector on the scarcely-susceptible variety prompted immediate and substantial transcriptomic changes that led to the rapid erection of further active defenses. On the other hand, in the most susceptible variety the response was delayed and mainly consisted of the induction of phytoalexin synthesis. Surprisingly, the jasmonic acid- and ethylene-mediated defense reactions, activated by the susceptible cultivar following FD-free insect feeding, were not detected in the presence of the phytoplasma-infected vector. CONCLUSIONS: The comparison of the transcriptomic response in two grapevine varieties with different levels of susceptibility to Flavescence dorèe highlighted both passive and active defense mechanisms against the vector and/or the pathogen in the scarcely-susceptible variety, as well as the capacity of the phytoplasmas to repress the defense reaction against the insect in the susceptible variety.

Mitochondrial gene order evolution in Mollusca: Inference of the ancestral state from the mtDNA of Chaetopleura apiculata (Polyplacophora, Chaetopleuridae).

The mitochondrial genome architecture of polyplacophorans has been usually regarded as being very ancient in comparison to all mollusks. However, even if some complete chiton mtDNAs have been recently sequenced, thorough studies of their evolution are lacking. To further expand the set of complete chiton mtDNAs and perform such analysis, we sequenced the mitochondrial genome of the Eastern beaded chiton Chaetopleura apiculata (Chaetopleuridae) using next-generation sequencing. With mitochondrial sequences from all available chiton mtDNAs, we also built a phylogeny on which we reconstructed the evolution of gene arrangement in this class. The arrangement of C. apiculata proved to be the most primitive known so far for polyplacophorans. Comparing this gene order to those of other molluscan classes, we found that it most probably is the original gene order of the last common ancestor to all extant Mollusca. The ancient mitochondrial genome organization of C. apiculata is an important information that may help reconstructing the phylogeny of Mollusca and their relationship with other lophotrochozoans.

A resourceful genome: updating the functional repertoire and evolutionary role of animal mitochondrial DNAs.

Recent data from mitochondrial genomics and proteomics research demonstrate the existence of several atypical mitochondrial protein-coding genes (other than the standard set of 13) and the involvement of mtDNA-encoded proteins in functions other than energy production in several animal species including humans. These results are of considerable importance for evolutionary and cellular biology because they indicate that animal mtDNAs have a larger functional repertoire than previously believed. This review summarizes recent studies on animal species with a non-standard mitochondrial functional repertoire and discusses how these genetic novelties represent promising candidates for studying the role of the mitochondrial genome in speciation.

In silico analyses of mitochondrial ORFans in freshwater mussels (Bivalvia: Unionoida) provide a framework for future studies of their origin and function.

BACKGROUND: Many species of bivalves exhibit a unique system of mtDNA transmission named Doubly Uniparental Inheritance (DUI). Under this system, species have two distinct, sex-linked mitochondrial genomes: the M-type mtDNA, which is transmitted by males to male offspring and found in spermatozoa, and the F-type mtDNA, which is transmitted by females to all offspring, and found in all tissues of females and in somatic tissues of males. Bivalves with DUI also have sex-specific mitochondrial ORFan genes, (M-orf in the M mtDNA, F-orf in the F mtDNA), which are open reading frames having no detectable homology and no known function. DUI ORFan proteins have previously been characterized in silico in a taxonomically broad array of bivalves including four mytiloid, one veneroid and one unionoid species. However, the large evolutionary distance among these taxa prevented a meaningful comparison of ORFan properties among these divergent lineages. The present in silico study focuses on a suite of more closely-related Unionoid freshwater mussel species to provide more reliably interpretable information on patterns of conservation and properties of DUI ORFans. Unionoid species typically have separate sexes, but hermaphroditism also occurs, and hermaphroditic species lack the M-type mtDNA and possess a highly mutated version of the F-orf in their maternally transmitted mtDNA (named H-orf in these taxa). In this study, H-orfs and their respective proteins are analysed for the first time. RESULTS: Despite a rapid rate of evolution, strong structural and functional similarities were found for M-ORF proteins compared among species, and among the F-ORF and H-ORF proteins across the studied species. In silico analyses suggest that M-ORFs have a role in transport and cellular processes such as signalling, cell cycle and division, and cytoskeleton organisation, and that F-ORFs may be involved in cellular traffic and transport, and in immune response. H-ORFs appear to be structural glycoproteins, which may be involved in signalling, transport and transcription. Our results also support either a viral or a mitochondrial origin for the ORFans. CONCLUSIONS: Our findings reveal striking structural and functional similarities among proteins encoded by mitochondrial ORFans in freshwater mussels, and strongly support a role for these genes in the DUI mechanism. Our analyses also support the possibility of DUI systems with elements of different sources/origins and different mechanisms of action in the distantly-related DUI taxa. Parallel situations to the novel mitochondrially-encoded functions of freshwater mussel ORFans present in some other eukaryotes are also discussed.

Identification of a protein network interacting with TdRF1, a wheat RING ubiquitin ligase with a protective role against cellular dehydration.

Plants exploit ubiquitination to modulate the proteome with the final aim to ensure environmental adaptation and developmental plasticity. Ubiquitination targets are specifically driven to degradation through the action of E3 ubiquitin ligases. Genetic analyses have indicated wide functions of ubiquitination in plant life; nevertheless, despite the large number of predicted E3s, only a few of them have been characterized so far, and only a few ubiquitination targets are known. In this work, we characterized durum wheat (Triticum durum) RING Finger1 (TdRF1) as a durum wheat nuclear ubiquitin ligase. Moreover, its barley (Hordeum vulgare) homolog was shown to protect cells from dehydration stress. A protein network interacting with TdRF1 has been defined. The transcription factor WHEAT BEL1-TYPE HOMEODOMAIN1 (WBLH1) was degraded in a TdRF1-dependent manner through the 26S proteasome in vivo, the mitogen-activated protein kinase TdWNK5 [for Triticum durum WITH NO LYSINE (K)5] was able to phosphorylate TdRF1 in vitro, and the RING-finger protein WHEAT VIVIPAROUS-INTERACTING PROTEIN2 (WVIP2) was shown to have a strong E3 ligase activity. The genes coding for the TdRF1 interactors were all responsive to cold and/or dehydration stress, and a negative regulative function in dehydration tolerance was observed for the barley homolog of WVIP2. A role in the control of plant development was previously known, or predictable based on homology, for wheat BEL1-type homeodomain1(WBLH1). Thus, TdRF1 E3 ligase might act regulating the response to abiotic stress and remodeling plant development in response to environmental constraints.

Identification of SUMO targets by a novel proteomic approach in plants(F).

Post-translational modifications (PTMs) chemically and physically alter the properties of proteins, including their folding, subcellular localization, stability, activity, and consequently their function. In spite of their relevance, studies on PTMs in plants are still limited. Small Ubiquitin-like Modifier (SUMO) modification regulates several biological processes by affecting protein-protein interactions, or changing the subcellular localizations of the target proteins. Here, we describe a novel proteomic approach to identify SUMO targets that combines 2-D liquid chromatography, immunodetection, and mass spectrometry (MS) analyses. We have applied this approach to identify nuclear SUMO targets in response to heat shock. Using a bacterial SUMOylation system, we validated that some of the targets identified here are, in fact, labeled with SUMO1. Interestingly, we found that GIGANTEA (GI), a photoperiodic-pathway protein, is modified with SUMO in response to heat shock both in vitro and in vivo.

QTL mapping for kernel-related traits in a durum wheat x T. dicoccum segregating population.

Durum wheat breeding relies on grain yield improvement to meet its upcoming demand while coping with climate change. Kernel size and shape are the determinants of thousand kernel weight (TKW), which is a key component of grain yield, and the understanding of the genetic control behind these traits supports the progress in yield potential. The present study aimed to dissect the genetic network responsible for kernel size components (length, width, perimeter, and area) and kernel shape traits (width-to-length ratio and formcoefficient) as well as their relationships with kernel weight, plant height, and heading date in durum wheat. Quantitative Trait Locus (QTL) mapping was performed on a segregating population of 110 recombinant inbred lines, derived from a cross between the domesticated emmer wheat accession MG5323 and the durum wheat cv. Latino, evaluated in four different environments. A total of 24 QTLs stable across environments were found and further grouped in nine clusters on chromosomes 2A, 2B, 3A, 3B, 4B, 6B, and 7A. Among them, a QTL cluster on chromosome 4B was associated with kernel size traits and TKW, where the parental MG5323 contributed the favorable alleles, highlighting its potential to improve durum wheat germplasm. The physical positions of the clusters, defined by the projection on the T. durum reference genome, overlapped with already known genes (i.e., BIG GRAIN PROTEIN 1 on chromosome 4B). These results might provide genome-based guidance for the efficient exploitation of emmer wheat diversity in wheat breeding, possibly through yield-related molecular markers.

The complete mitogenome of the inarticulate brachiopod Glottidia pyramidata reveals insights into gene order variation, deviant ATP8 and mtORFans in the Brachiopoda.

Brachiopods are a clade of marine organisms with a tremendously diverse and abundant fossil record but with fewer than 500 species extant today. Even if a better understanding of their biology and genetics could help to test hypotheses about their impressive decline, knowledge of genetics and evolutionary genomics in extant brachiopods is very poor. Here, we present the complete mitochondrial genome sequence of the inarticulate Glottidia pyramidata, an eastern North American extant representative of the phylum Brachiopoda. Besides the general characteristics of the sequenced mitogenome, we present its unusual features such as deviant ATP8 protein sequence and supernumerary ORFs, and also unique gene order, considering the available genome sequences of other brachiopod species.

Unorthodox features in two venerid bivalves with doubly uniparental inheritance of mitochondria.

In animals, strictly maternal inheritance (SMI) of mitochondria is the rule, but one exception (doubly uniparental inheritance or DUI), marked by the transmission of sex-specific mitogenomes, has been reported in bivalves. Associated with DUI is a frequent modification of the mitochondrial cox2 gene, as well as additional sex-specific mitochondrial genes not involved in oxidative phosphorylation. With the exception of freshwater mussels (for 3 families of the order Unionida), these DUI-associated features have only been shown in few species [within Mytilidae (order Mytilida) and Veneridae (order Venerida)] because of the few complete sex-specific mitogenomes published for these orders. Here, we present the complete sex-specific mtDNAs of two recently-discovered DUI species in two families of the order Venerida, Scrobicularia plana (Semelidae) and Limecola balthica (Tellinidae). These species display the largest differences in genome size between sex-specific mitotypes in DUI species (>10 kb), as well as the highest mtDNA divergences (sometimes reaching >50%). An important in-frame insertion (>3.5 kb) in the male cox2 gene is partly responsible for the differences in genome size. The S. plana cox2 gene is the largest reported so far in the Kingdom Animalia. The mitogenomes may be carrying sex-specific genes, indicating that general mitochondrial features are shared among DUI species.

Genomic Regions From an Iranian Landrace Increase Kernel Size in Durum Wheat.

Kernel size and shape are important parameters determining the wheat profitability, being main determinants of yield and its technological quality. In this study, a segregating population of 118 recombinant inbred lines, derived from a cross between the Iranian durum landrace accession "Iran_249" and the Iranian durum cultivar "Zardak", was used to investigate durum wheat kernel morphology factors and their relationships with kernel weight, and to map the corresponding QTLs. A high density genetic map, based on wheat 90k iSelect Infinium SNP assay, comprising 6,195 markers, was developed and used to perform the QTL analysis for kernel length and width, traits related to kernel shape and weight, and heading date, using phenotypic data from three environments. Overall, a total of 31 different QTLs and 9 QTL interactions for kernel size, and 21 different QTLs and 5 QTL interactions for kernel shape were identified. The landrace Iran_249 contributed the allele with positive effect for most of the QTLs related to kernel length and kernel weight suggesting that the landrace might have considerable potential toward enhancing the existing gene pool for grain shape and size traits and for further yield improvement in wheat. The correlation among traits and co-localization of corresponding QTLs permitted to define 11 clusters suggesting causal relationships between simplest kernel size trait, like kernel length and width, and more complex secondary trait, like kernel shape and weight related traits. Lastly, the recent release of the T. durum reference genome sequence allowed to define the physical interval of our QTL/clusters and to hypothesize novel candidate genes inspecting the gene content of the genomic regions associated to target traits.

Genome Survey of the Freshwater Mussel Venustaconcha ellipsiformis (Bivalvia: Unionida) Using a Hybrid De Novo Assembly Approach.

Freshwater mussels (Bivalvia: Unionida) serve an important role as aquatic ecosystem engineers but are one of the most critically imperilled groups of animals. Here, we used a combination of sequencing strategies to assemble and annotate a draft genome of Venustaconcha ellipsiformis, which will serve as a valuable genomic resource given the ecological value and unique "doubly uniparental inheritance" mode of mitochondrial DNA transmission of freshwater mussels. The genome described here was obtained by combining high-coverage short reads (65× genome coverage of Illumina paired-end and 11× genome coverage of mate-pairs sequences) with low-coverage Pacific Biosciences long reads (0.3× genome coverage). Briefly, the final scaffold assembly accounted for a total size of 1.54 Gb (366,926 scaffolds, N50 = 6.5 kb, with 2.3% of "N" nucleotides), representing 86% of the predicted genome size of 1.80 Gb, while over one third of the genome (37.5%) consisted of repeated elements and >85% of the core eukaryotic genes were recovered. Given the repeated genetic bottlenecks of V. ellipsiformis populations as a result of glaciations events, heterozygosity was also found to be remarkably low (0.6%), in contrast to most other sequenced bivalve species. Finally, we reassembled the full mitochondrial genome and found six polymorphic sites with respect to the previously published reference. This resource opens the way to comparative genomics studies to identify genes related to the unique adaptations of freshwater mussels and their distinctive mitochondrial inheritance mechanism.

Evolution of sex-dependent mtDNA transmission in freshwater mussels (Bivalvia: Unionida).

Doubly uniparental inheritance (DUI) describes a mode of mtDNA transmission widespread in gonochoric freshwater mussels (Bivalvia: Palaeoheterodonta: Unionida). In this system, both female- and male-transmitted mtDNAs, named F and M respectively, coexist in the same species. In unionids, DUI is strictly correlated to gonochorism and to the presence of the atypical open reading frames (ORFans) F-orf and M-orf, respectively inside F and M mtDNAs, which are hypothesized to participate in sex determination. However, DUI is not found in all three Unionida superfamilies (confirmed in Hyrioidea and Unionoidea but not in Etherioidea), raising the question of its origin in these bivalves. To reconstruct the co-evolution of DUI and of ORFans, we sequenced the mtDNAs of four unionids (two gonochoric with DUI, one gonochoric and one hermaphroditic without DUI) and of the related gonochoric species Neotrigonia margaritacea (Palaeoheterodonta: Trigoniida). Our analyses suggest that rearranged mtDNAs appeared early during unionid radiation, and that a duplicated and diverged atp8 gene evolved into the M-orf associated with the paternal transmission route in Hyrioidea and Unionoidea, but not in Etherioidea. We propose that novel mtDNA-encoded genes can deeply influence bivalve sex determining systems and the evolution of the mitogenomes in which they occur.

The extremely divergent maternally- and paternally-transmitted mitochondrial genomes are co-expressed in somatic tissues of two freshwater mussel species with doubly uniparental inheritance of mtDNA.

Freshwater mussel species with doubly uniparental inheritance (DUI) of mtDNA are unique because they are naturally heteroplasmic for two extremely divergent mtDNAs with ~50% amino acid differences for protein-coding genes. The paternally-transmitted mtDNA (or M mtDNA) clearly functions in sperm in these species, but it is still unknown whether it is transcribed when present in male or female soma. In the present study, we used PCR and RT-PCR to detect the presence and expression of the M mtDNA in male and female somatic and gonadal tissues of the freshwater mussel species Venustaconcha ellipsiformis and Utterbackia peninsularis (Unionidae). This is the first study demonstrating that the M mtDNA is transcribed not only in male gonads, but also in male and female soma in freshwater mussels with DUI. Because of the potentially deleterious nature of heteroplasmy, we suggest the existence of different mechanisms in DUI species to deal with this possibly harmful situation, such as silencing mechanisms for the M mtDNA at the transcriptional, post-transcriptional and/or post-translational levels. These hypotheses will necessitate additional studies in distantly-related DUI species that could possess different mechanisms of action to deal with heteroplasmy.

Post-transcriptional and post-translational regulations of drought and heat response in plants: a spider's web of mechanisms.

Drought and heat tolerance are complex quantitative traits. Moreover, the adaptive significance of some stress-related traits is more related to plant survival than to agronomic performance. A web of regulatory mechanisms fine-tunes the expression of stress-related traits and integrates both environmental and developmental signals. Both post-transcriptional and post-translational modifications contribute substantially to this network with a pivotal regulatory function of the transcriptional changes related to cellular and plant stress response. Alternative splicing and RNA-mediated silencing control the amount of specific transcripts, while ubiquitin and SUMO modify activity, sub-cellular localization and half-life of proteins. Interactions across these modification mechanisms ensure temporally and spatially appropriate patterns of downstream-gene expression. For key molecular components of these regulatory mechanisms, natural genetic diversity exists among genotypes with different behavior in terms of stress tolerance, with effects upon the expression of adaptive morphological and/or physiological target traits.

The largest unassigned regions of the male- and female-transmitted mitochondrial DNAs in Musculista senhousia (Bivalvia Mytilidae).

Musculista senhousia is a marine mussel with doubly uniparental inheritance (DUI) of mitochondria. In this study we analyzed the largest unassigned region (LUR) of its female- and male-transmitted mitochondrial genomes, described their fine characteristics and searched for shared features. Our results suggest that both LURs contain the control region of their respective mitochondrial genomes. The female-transmitted control region is duplicated in tandem, with the two copies evolving in concert. This makes the F-mtDNA of M. senhousia the first Bivalve mitochondrial genome with this feature. We also compared M. senhousia control regions to that of other Mytilidae, and demonstrated that signals for basic mtDNA functions are retained over evolutionary times even among the fast-evolving mitochondrial genomes of DUI species. Finally, we discussed how similarities between female and male LURs may be explained in the context of DUI evolution and if the duplicated female control region might have influenced the DUI system in this species.

A comparative analysis of mitochondrial ORFans: new clues on their origin and role in species with doubly uniparental inheritance of mitochondria.

Despite numerous comparative mitochondrial genomics studies revealing that animal mitochondrial genomes are highly conserved in terms of gene content, supplementary genes are sometimes found, often arising from gene duplication. Mitochondrial ORFans (ORFs having no detectable homology and unknown function) were found in bivalve molluscs with Doubly Uniparental Inheritance (DUI) of mitochondria. In DUI animals, two mitochondrial lineages are present: one transmitted through females (F-type) and the other through males (M-type), each showing a specific and conserved ORF. The analysis of 34 mitochondrial major Unassigned Regions of Musculista senhousia F- and M-mtDNA allowed us to verify the presence of novel mitochondrial ORFs in this species and to compare them with ORFs from other species with ascertained DUI, with other bivalves and with animals showing new mitochondrial elements. Overall, 17 ORFans from nine species were analyzed for structure and function. Many clues suggest that the analyzed ORFans arose from endogenization of viral genes. The co-option of such novel genes by viral hosts may have determined some evolutionary aspects of host life cycle, possibly involving mitochondria. The structure similarity of DUI ORFans within evolutionary lineages may also indicate that they originated from independent events. If these novel ORFs are in some way linked to DUI establishment, a multiple origin of DUI has to be considered. These putative proteins may have a role in the maintenance of sperm mitochondria during embryo development, possibly masking them from the degradation processes that normally affect sperm mitochondria in species with strictly maternal inheritance.

Structure, transcription, and variability of metazoan mitochondrial genome: perspectives from an unusual mitochondrial inheritance system.

Despite its functional conservation, the mitochondrial genome (mtDNA) presents strikingly different features among eukaryotes, such as size, rearrangements, and amount of intergenic regions. Nonadaptive processes such as random genetic drift and mutation rate play a fundamental role in shaping mtDNA: the mitochondrial bottleneck and the number of germ line replications are critical factors, and different patterns of germ line differentiation could be responsible for the mtDNA diversity observed in eukaryotes. Among metazoan, bivalve mollusc mtDNAs show unusual features, like hypervariable gene arrangements, high mutation rates, large amount of intergenic regions, and, in some species, an unique inheritance system, the doubly uniparental inheritance (DUI). The DUI system offers the possibility to study the evolutionary dynamics of mtDNAs that, despite being in the same organism, experience different genetic drift and selective pressures. We used the DUI species Ruditapes philippinarum to study intergenic mtDNA functions, mitochondrial transcription, and polymorphism in gonads. We observed: 1) the presence of conserved functional elements and novel open reading frames (ORFs) that could explain the evolutionary persistence of intergenic regions and may be involved in DUI-specific features; 2) that mtDNA transcription is lineage-specific and independent from the nuclear background; and 3) that male-transmitted and female-transmitted mtDNAs have a similar amount of polymorphism but of different kinds, due to different population size and selection efficiency. Our results are consistent with the hypotheses that mtDNA evolution is strongly dependent on the dynamics of germ line formation, and that the establishment of a male-transmitted mtDNA lineage can increase male fitness through selection on sperm function.