Ancestral Genome Research Articles

Restriction of AID activity and somatic hypermutation by PARP-1.

Affinity maturation of the humoral immune response depends on somatic hypermutation (SHM) of immunoglobulin (Ig) genes, which is initiated by targeted lesion introduction by activation-induced deaminase (AID), followed by error-prone DNA repair. Stringent regulation of this process is essential to prevent genetic instability, but no negative feedback control has been identified to date. Here we show that poly(ADP-ribose) polymerase-1 (PARP-1) is a key factor restricting AID activity during somatic hypermutation. Poly(ADP-ribose) (PAR) chains formed at DNA breaks trigger AID-PAR association, thus preventing excessive DNA damage induction at sites of AID action. Accordingly, AID activity and somatic hypermutation at the Ig variable region is decreased by PARP-1 activity. In addition, PARP-1 regulates DNA lesion processing by affecting strand biased A:T mutagenesis. Our study establishes a novel function of the ancestral genome maintenance factor PARP-1 as a critical local feedback regulator of both AID activity and DNA repair during Ig gene diversification.

Characterization of the mitochondrial genome of Analcellicampa xanthosoma gen. et sp. nov. (Hymenoptera: Tenthredinidae).

A new genus with a new species of the tribe Hoplocampini of Hoplocampinae was described from China: Analcellicampa xanthosoma Wei & Niu, gen. et sp. nov. Hoplocampa danfengensis G. Xiao 1994 was designated as the type species of the new genus. The characters of Analcellicampa danfengensis (G. Xiao) comb. nov. were briefly discussed. A key to the tribes and known genera of Hoplocampinae was provided. The nearly complete mitochondrial genome of A. xanthosoma was characterized as having a length of 15,512 bp and containing 37 genes (22 tRNAs, 13 protein-coding genes (PCGs), and 2 rRNAs). The gene order of this new specimen was the same as that in the inferred insect ancestral mitochondrial genome. All PCGs were initiated by ATN codons and ended with TAA or T stop codons. All tRNAs had a typical cloverleaf secondary structure, except for trnS1. Remarkably, the helices H991 of rrnS and H47 of rrnL were redundant, while helix H563 of rrnL was highly conserved. A phylogeny based on previously reported symphytan mitochondrial genomes showed that A. xanthosoma is a sister group to Monocellicampa pruni, with high support values. We suggest that A. xanthosoma and M. pruni belong to the tribe Hoplocampini of Hoplocampinae.

An empirical approach to demographic inference with genomic data

Inference with population genetic data usually treats the population pedigree as a nuisance parameter, the unobserved product of a past history of random mating. However, the history of genetic relationships in a given population is a fixed, unobserved object, and so an alternative approach is to treat this network of relationships as a complex object we wish to learn about, by observing how genomes have been noisily passed down through it. This paper explores this point of view, showing how to translate questions about population genetic data into calculations with a Poisson process of mutations on all ancestral genomes. This method is applied to give a robust interpretation to the f4 statistic used to identify admixture, and to design a new statistic that measures covariances in mean times to most recent common ancestor between two pairs of sequences. The method more generally interprets population genetic statistics in terms of sums of specific functions over ancestral genomes, thereby providing concrete, broadly interpretable interpretations for these statistics. This provides a method for describing demographic history without simplified demographic models. More generally, it brings into focus the population pedigree, which is averaged over in model-based demographic inference.

Rainbow trout (Oncorhynchus mykiss) urea cycle and polyamine synthesis gene families show dynamic expression responses to inflammation

The urea cycle is an endogenous source of arginine that also supports removal of nitrogenous waste following protein metabolism. This cycle is considered inefficient in salmonids, where only 10–15% of nitrogenous waste is excreted as urea. In rainbow trout, arginine is an essential amino acid that has attracted attention due to its many functional roles. These roles include the regulation of protein deposition, immune responses and polyamine synthesis; the latter is directly linked to the urea cycle and involved in tissue repair. The key enzymes used in the urea cycle, namely arginase, ornithine transcarbamylase, argininosuccinate synthase and argininosuccinate lyase, in addition to two rate limiting enzymes required for polyamine synthesis (ornithine decarboxylase and s-adenosylmethionine decarboxylase) are poorly studied in fishes, and their responses to inflammation remain unknown. To address this knowledge gap, we characterised these gene families using phylogenetics and comparative genomics, investigated their mRNA distribution among a panel of tissues and established their transcriptional responses to an acute inflammatory response caused by bacterial infection in liver and muscle. Gene duplicates (paralogues) were identified for arginase (ARG1a, 1b, 2a and 2b), ornithine decarboxylase (ODC1 and 2) and s-adenosylmethionine decarboxylase (SAMdc1 and 2), including paralogues retained from an ancestral salmonid-specific whole genome duplication. ARG2a and 2b were highly upregulated following bacterial infection in liver, whereas ARG1b was downregulated, while both paralogues of SAMdc and ODC were upregulated in liver and unchanged in muscle. Overall, these findings improve our understanding of the molecules supporting the urea cycle and polyamine synthesis in fish, highlighting major changes in the regulation of these systems during inflammation.

Reconstruction of ancestral genome reveals chromosome evolution history for selected legume species.

Reconstruction of an ancestral genome for a set of plant species has been a challenging task because of complex histories that may include whole-genome duplications, segmental duplications, independent gene duplications or losses, diploidization and rearrangement events. Here, we describe the reconstruction a hypothetical ancestral genome for the papilionoid legumes (the largest subfamily within the third largest family in flowering plants), and evaluate the results relative to phylogenetic and chromosomal count data for this group of legumes, spanning 294 diverse papilionoid genera. To reconstruct the ancestral genomes for nine legume species with sequenced genomes, we used a maximum likelihood approach combined with a novel method for identifying informative markers for this purpose. Analyzing genomes from four species within the Phaseoleae, two in Dalbergieae, two in the 'inverted repeat loss' clade, and one in the Robinieae, we infer a common ancestral genome with nine chromosomes. The reconstructed genome structural histories are consistent with chromosomal and phylogenetic histories, but we also infer that a common ancestor with nine chromosomes was probably intermediate to an earlier state of 14 chromosomes following a whole-genome duplication that pre-dated the radiation of the papilionoid legumes, evidence for which is found in early-diverging papilionoid lineages.

Perspectives on the clonal persistence of presumed \u2018ghost\u2019 genomes in unisexual or allopolyploid taxa arising via hybridization

Although hybridization between non-sibling species rarely results in viable or fertile offspring, it occasionally produces self-perpetuating or sexually-parasitic lineages in which ancestral genomes are inherited clonally and thus may persist as ‘ghost species’ after ancestor extinction. Ghost species have been detected in animals and plants, for polyploid and diploid organisms, and across clonal, semi-clonal, and even sexual reproductive modes. Here we use a detailed investigation of the evolutionary and taxonomic status of a newly-discovered, putative ghost lineage (HX) in the fish genus Hypseleotris to provide perspectives on several important issues not previously explored by other studies on ghost species, but relevant to ongoing discussions about their detection, conservation, and artificial re-creation. Our comprehensive genetic (allozymes, mtDNA) and genomic (SNPs) datasets successfully identified a threatened sexual population of HX in one tiny portion of the extensive distribution displayed by two hemi-clonal HX-containing lineages. We also discuss what confidence should be placed on any assertion that an ancestral species is actually extinct, and how to assess whether any putative sexual ancestor represents a pure remnant, as shown here, or a naturally-occurring resurrection via the crossing of compatible clones or hemi-clones.

Gram-Positive Bacteria-Like DNA Binding Machineries Involved in Replication Initiation and Termination Mechanisms of Mimivirus.

The detailed mechanisms of replication initiation, termination and segregation events were not yet known in Acanthamoeba polyphaga mimivirus (APMV). Here, we show detailed bioinformatics-based analyses of chromosomal replication in APMV from initiation to termination mediated by proteins bound to specific DNA sequences. Using GC/AT skew and coding sequence skew analysis, we estimated that the replication origin is located at 382 kb in the APMV genome. We performed homology-modeling analysis of the gamma domain of APMV-FtsK (DNA translocase coordinating chromosome segregation) related to FtsK-orienting polar sequences (KOPS) binding, suggesting that there was an insertion in the gamma domain which maintains the structure of the DNA binding motif. Furthermore, UvrD/Rep-like helicase in APMV was homologous to Bacillus subtilis AddA, while the chi-like quartet sequence 5′-CCGC-3′ was frequently found in the estimated ori region, suggesting that chromosomal replication of APMV is initiated via chi-like sequence recognition by UvrD/Rep-like helicase. Therefore, the replication initiation, termination and segregation of APMV are presumably mediated by DNA repair machineries derived from gram-positive bacteria. Moreover, the other frequently observed quartet sequence 5′-CGGC-3′ in the ori region was homologous to the mitochondrial signal sequence of replication initiation, while the comparison of quartet sequence composition in APMV/Rickettsia-genome showed significantly similar values, suggesting that APMV also conserves the mitochondrial replication system acquired from an ancestral genome of mitochondria during eukaryogenesis.

The story of promiscuous crucifers: origin and genome evolution of an invasive species, Cardamine occulta (Brassicaceae), and its relatives.

Cardamine occulta (Brassicaceae) is an octoploid weedy species (2n = 8x = 64) originated in Eastern Asia. It has been introduced to other continents including Europe and considered to be an invasive species. Despite its wide distribution, the polyploid origin of C. occulta remained unexplored. The feasibility of comparative chromosome painting (CCP) in crucifers allowed us to elucidate the origin and genome evolution in Cardamine species. We aimed to investigate the genome structure of C. occulta in comparison with its tetraploid (2n = 4x = 32, C. kokaiensis and C. scutata) and octoploid (2n = 8x = 64, C. dentipetala) relatives. Genomic in situ hybridization (GISH) and large-scale CCP were applied to uncover the parental genomes and chromosome composition of the investigated Cardamine species. All investigated species descended from a common ancestral Cardamine genome (n = 8), structurally resembling the Ancestral Crucifer Karyotype (n = 8), but differentiated by a translocation between chromosomes AK6 and AK8. Allotetraploid C. scutata originated by hybridization between two diploid species, C. parviflora and C. amara (2n = 2x = 16). By contrast, C. kokaiensis has an autotetraploid origin from a parental genome related to C. parviflora. Interestingly, octoploid C. occulta probably originated through hybridization between the tetraploids C. scutata and C. kokaiensis. The octoploid genome of C. dentipetala probably originated from C. scutata via autopolyploidization. Except for five species-specific centromere repositionings and one pericentric inversion post-dating the polyploidization events, the parental subgenomes remained stable in the tetra- and octoploids. Comparative genome structure, origin and evolutionary history was reconstructed in C. occulta and related species. For the first time, whole-genome cytogenomic maps were established for octoploid plants. Post-polyploid evolution in Asian Cardamine polyploids has not been associated with descending dysploidy and intergenomic rearrangements. The combination of different parental (sub)genomes adapted to distinct habitats provides an evolutionary advantage to newly formed polyploids by occupying new ecological niches.

Natural History of a Satellite DNA Family: From the Ancestral Genome Component to Species-Specific Sequences, Concerted and Non-Concerted Evolution.

Satellite DNA (satDNA) is the most variable fraction of the eukaryotic genome. Related species share a common ancestral satDNA library and changing of any library component in a particular lineage results in interspecific differences. Although the general developmental trend is clear, our knowledge of the origin and dynamics of satDNAs is still fragmentary. Here, we explore whole genome shotgun Illumina reads using the RepeatExplorer (RE) pipeline to infer satDNA family life stories in the genomes of Chenopodium species. The seven diploids studied represent separate lineages and provide an example of a species complex typical for angiosperms. Application of the RE pipeline allowed by similarity searches a determination of the satDNA family with a basic monomer of ~40 bp and to trace its transformation from the reconstructed ancestral to the species-specific sequences. As a result, three types of satDNA family evolutionary development were distinguished: (i) concerted evolution with mutation and recombination events; (ii) concerted evolution with a trend toward increased complexity and length of the satellite monomer; and (iii) non-concerted evolution, with low levels of homogenization and multidirectional trends. The third type is an example of entire repeatome transformation, thus producing a novel set of satDNA families, and genomes showing non-concerted evolution are proposed as a significant source for genomic diversity.

Organization and evolution of the chalcone synthase gene family in bread wheat and relative species

BackgroundFlavonoid compounds are secondary plant metabolites, having a functional importance in plant development, protection from pathogens and unfavorable environmental factors. Chalcone synthase (CHS) is a key enzyme in the biosynthesis of flavonoids; it is involved in biosynthesis of all classes of flavonoid compounds. Nevertheless, the Chs gene family in bread wheat (Triticum aestivum L.) has been not characterized yet. The aim of the current study was to investigate structural and functional organization of the Chs genes and evolution of this gene family in bread wheat and relative species.ResultsThe nucleotide sequences of the eight Chs copies in T. aestivum were identified. Among them, two homoeologous sets of the Chs genes were located on the short (Chs-A1, −B1, −D1) and the long (Chs-A4, −B4, −D4) arms of homoeologous group 2 chromosomes. Two paralogous gene copies in the B-genome (Chs-B2, −B3) were located in the distal regions of 2BS chromosome. To clarify the origin of Chs duplications in the B-genome the phylogenetic analysis with the Chs sequences of Triticum and Aegilops species carrying ancestral genomes was conducted. It was estimated that the first duplication event occurred in the genome of the common ancestor of Triticum and Aegilops genera about 10–12 million years ago (MYA), then another copy was formed in the ancestor of the B-genome about 6–7 MYA. A homology modeling revealed high sequence similarity of bread wheat CHS enzymes. A number of short deletions in coding regions of some Chs sequences are not expected to have any significant functional effects. Estimation of transcriptional activity of the Chs copies along with a comparative analysis of their promoters structure suggested their functional specialization, which likely contributed to the maintaining of the duplicated Chs genes in wheat genome.ConclusionsFrom possible ten Chs copies in bread wheat genome, eight members of this family retained their intact structure and activity, while two copies appear to be lost at the level of diploid and tetraploid ancestors. Transcriptional assay along with a comparative analysis of the cis-regulatory elements revealed their functional diversification. The multiple functions supported by the Chs family are assumed to be a driving force for duplications of the Chs gene and their retention in plant genome.

Comparative analysis of gene family size provides insight into the adaptive evolution of vertebrates

Copy numbers of homologous gene families vary greatly among different species, which is caused by the differences in the rates of gene gain and loss. It is well known that gene copy number variation can be responsible for the phenotypic novelties of particular species. In this study, 64 species that represent the main vertebrate groups spanning evolutionary period of about 600 million years were selected and the homology of gene families across these species were established, thereby revealing the evolutionary patterns of gene family size in vertebrates. The results show that among the 6857 gene families inferred to be present in the most recent common ancestor of the vertebrates, 6712 had changed their sizes in at least one lineage, and these gene families had contracted in most cases. Gene families in Choloepus hoffmanni and Danio rerio had undergone the greatest contraction and expansion, respectively. Based on the highly dynamic evolution of vertebrate gene family size, we sought to identify any genomic signals that might be related to the evolution of specific vertebrate populations from the perspective of the distinct gene family size changes. We observed a high proportion of gene family amplification occurred, probably due to genome-wide duplication in the recent common ancestral genome of teleosts, which was followed by contraction in the decedents due to the extensive gene fractionation. Furthermore, we found evidence that orphan genes in the bony fish might contribute to the adaptive evolution of fish in aquatic environment. For example, some orphan genes were involved in fin development, tail development and kidney physiology. Overall, our work provides novel insights into the evolution of vertebrate gene family size and provides several lines of evidence for understanding the relationship between the genome evolution and phenotypic diversity in vertebrates.

Constructing a High-Density Genetic Linkage Map for Large Yellow Croaker (Larimichthys crocea) and Mapping Resistance Trait Against Ciliate Parasite Cryptocaryon irritans.

The large yellow croaker (Larimichthys crocea) is the most economically important marine cage-farming fish in China in the past decade. However, the sustainable development of large yellow croaker aquaculture has been severely hampered by several diseases, of which, the white spot disease caused by ciliate protozoan parasite Cryptocaryon irritans ranks the most damaging disease in large yellow croaker cage farms. To better understand the genetic basis of parasite infection and disease resistance to C. irritans, it is vital to map the traits and localize the underlying candidate genes in L. crocea genome. Here, we constructed a high-density genetic linkage map using double-digest restriction-site associated DNA(ddRAD)-based high-throughput SNP genotyping data of a F1 mapping family, which had been challenged with C. irritans for resistant trait measure. A total of 5261 SNPs was grouped and oriented into 24 linkage groups (LGs), representing 24 chromosomes of L. crocea. The total genetic map length was 1885.67cM with an average inter-locus distance of 0.36cM. Quantitative trait loci (QTL) mapping identified seven significant QTLs in four LGs linked to C. irritans disease resistance. Candidate genes underlying disease resistance were identified from the reference genome, including ifnar1, ifngr2, ikbke, and CD112. Comparative genomic analysis between large yellow croaker and the four closely related species revealed high evolutionary conservation of chromosomes, though inter-chromosomal rearrangements do exist. Especially, the croaker genome structure was closer to the medaka genome than stickleback, indicating that the croaker genome might retain the teleost ancestral genome structure. The high-density genetic linkage map provides an important tool and resource for fine mapping, comparative genome analysis, and molecular selective breeding of large yellow croaker.

Phylogenetic analysis of S-RNase-like genes from Fragaria provides evidence for ancestral S-RNase duplication in the Rosaceae lineage

Self-incompatibility (SI) is associated with cultural and breeding difficulties in rosaceous fruit trees such as cherry, plum, apple, and pear. Many species in Rosaceae, Solanaceae, and Plantaginaceae have S-RNase-based gametophytic SI, with S-RNase and F-box proteins as pistil S and pollen S determinants, respectively. Although the S-RNase-based gametophytic SI system is generally believed to have evolved once in a eudicot ancestor before the divergence of asterids and rosids, we previously found that S-RNase and its homologs were duplicated in the ancestral genome of Rosaceae. Here, we reinvestigated the evolutionary paths of S-RNase-like genes using homologs from the diploid Fragaria species F. nipponica (SI), F. nubicola (SI), and F. iinumae (SI phenotype undetermined). Our phylogenetic analysis and estimation of proxy ages for the establishment of S-RNase and its homologs in rosaceous species supports the previously inferred ancestral S-RNase duplication in the Rosaceae lineage. The divergence of S-RNases in subtribe Malinae and Prunus predated the most recent common ancestor of Rosaceae, which includes genes from three rosaceous species (Fragaria, Prunus, and Malus). The results of the present study strengthen the idea that the genus Prunus and subtribe Malinae have evolved different S-RNases for SI recognition. These results may explain the genus- or subtribe-specific SI recognition mechanism operating in rosaceous fruit trees and will be valuable for artificial SI control and self-compatible breeding.

The complete mitochondrial genome of Monolepta occifuvis Gressitt & Kimoto (Coleoptera: Chrysomelidae: Galerucinae)

The complete mitochondrial (mt) genome of Monolepta occifuvis Gressitt & Kimoto (Coleoptera: Chrysomelidae: Galerucinae) was sequenced with 15,998 bp long. The sequenced genome has the standard metazoan complement of 37 genes and an A + T-rich region which is consistent with insect ancestral genome arrangement. All protein-coding genes (PCGs) initiate with ATN except that nad1 initiates with TTG. The control region occupied 1,407 bp with two conserved blocks. The phylogenetic results support Alticinae as a separate subfamily and being the sister group to Galerucinae while the monophyly of Galerucinae was not recovered.

Characterization of five complete Cyrtodactylus mitogenome structures reveals low structural diversity and conservation of repeated sequences in the lineage.

Mitochondrial genomes (mitogenomes) of five Cyrtodactylus were determined. Their compositions and structures were similar to most of the available gecko lizard mitogenomes as 13 protein-coding, two rRNA and 22 tRNA genes. The non-coding control region (CR) of almost all Cyrtodactylus mitogenome structures contained a repeated sequence named the 75-bp box family, except for C. auribalteatus which contained the 225-bp box. Sequence similarities indicated that the 225-bp box resulted from the duplication event of 75-bp boxes, followed by homogenization and fixation in C. auribalteatus. The 75-bp box family was found in most gecko lizards with high conservation (55–75% similarities) and could form secondary structures, suggesting that this repeated sequence family played an important role under selective pressure and might involve mitogenome replication and the likelihood of rearrangements in CR. The 75-bp box family was acquired in the common ancestral genome of the gecko lizard, evolving gradually through each lineage by independent nucleotide mutation. Comparison of gecko lizard mitogenomes revealed low structural diversity with at least six types of mitochondrial gene rearrangements. Cyrtodactylus mitogenome structure showed the same gene rearrangement as found in most gecko lizards. Advanced mitogenome information will enable a better understanding of structure evolution mechanisms.

Genetic ancestry, admixture and health determinants in Latin America

BackgroundModern Latin American populations were formed via genetic admixture among ancestral source populations from Africa, the Americas and Europe. We are interested in studying how combinations of genetic ancestry in admixed Latin American populations may impact genomic determinants of health and disease. For this study, we characterized the impact of ancestry and admixture on genetic variants that underlie health- and disease-related phenotypes in population genomic samples from Colombia, Mexico, Peru, and Puerto Rico.ResultsWe analyzed a total of 347 admixed Latin American genomes along with 1102 putative ancestral source genomes from Africans, Europeans, and Native Americans. We characterized the genetic ancestry, relatedness, and admixture patterns for each of the admixed Latin American genomes, finding a spectrum of ancestry proportions within and between populations. We then identified single nucleotide polymorphisms (SNPs) with anomalous ancestry-enrichment patterns, i.e. SNPs that exist in any given Latin American population at a higher frequency than expected based on the population’s genetic ancestry profile. For this set of ancestry-enriched SNPs, we inspected their phenotypic impact on disease, metabolism, and the immune system. All four of the Latin American populations show ancestry-enrichment for a number of shared pathways, yielding evidence of similar selection pressures on these populations during their evolution. For example, all four populations show ancestry-enriched SNPs in multiple genes from immune system pathways, such as the cytokine receptor interaction, T cell receptor signaling, and antigen presentation pathways. We also found SNPs with excess African or European ancestry that are associated with ancestry-specific gene expression patterns and play crucial roles in the immune system and infectious disease responses. Genes from both the innate and adaptive immune system were found to be regulated by ancestry-enriched SNPs with population-specific regulatory effects.ConclusionsAncestry-enriched SNPs in Latin American populations have a substantial effect on health- and disease-related phenotypes. The concordant impact observed for same phenotypes across populations points to a process of adaptive introgression, whereby ancestry-enriched SNPs with specific functional utility appear to have been retained in modern populations by virtue of their effects on health and fitness.

Discovery of highly divergent lineages of plant-associated astro-like viruses sheds light on the emergence of potyviruses

RNA viruses are believed to have originated from a common ancestor, but how this ancestral genome evolved into the large variety of genomic architectures and viral proteomes we see today remains largely unknown. Tackling this question is hindered by the lack of universally conserved proteins other than the RNA-dependent RNA polymerase (RdRp) as well as a limited RNA virus sampling. The latter is still heavily biased towards relatively few viral lineages from a non-representative collection of hosts, which complicates studies aiming to reveal possible trajectories during the evolution of RNA virus genomes that are favored over others.We report the discovery of 11 highly divergent lineages of viruses with genomic architectures that resemble those of the astroviruses. These genomes were initially identified through a sequence homology search in more than 6600 plant transcriptome projects from the Sequence Read Archive (SRA) using astrovirus representatives as query. Seed-based viral genome assembly of unprocessed SRA data for several dozens of the most promising hits resulted in two viral genome sequences with full-length coding regions, nine partial genomes and a much larger number of short sequence fragments. Genomic and phylogenetic characterization of the 11 discovered viruses, which we coined plastroviruses (plant-associated astro-like viruses), showed that they are related to both astro- and potyviruses and allowed us to identify divergent Serine protease, RdRp and viral capsid domains encoded in the plastrovirus genome. Interestingly, some of the plastroviruses shared different features with potyviruses including the replacement of the catalytic Ser by a Cys residue in the protease active site. These results suggest that plastroviruses may have reached different points on an evolutionary trajectory from astro-like to poty-like genomes. A model how potyviruses might have emerged from (pl)astro-like ancestors in a multi-step process is discussed.

Preferential Homologous Chromosome Pairing in a Tetraploid Intergeneric Somatic Hybrid (Citrus reticulata + Poncirus trifoliata) Revealed by Molecular Marker Inheritance.

The creation of intergeneric somatic hybrids between Citrus and Poncirus is an efficient approach for citrus rootstock breeding, offering the possibility of combining beneficial traits from both genera into novel rootstock lineages. These somatic hybrids are also used as parents for further tetraploid sexual breeding. In order to optimize these latter breeding schemes, it is essential to develop knowledge on the mode of inheritance in the intergeneric tetraploid hybrids. We assessed the meiotic behavior of an intergeneric tetraploid somatic hybrid resulting from symmetric protoplast fusion of diploid Citrus reticulata and diploid Poncirus trifoliata. The analysis was based on the segregation patterns of 16 SSR markers and 9 newly developed centromeric/pericentromeric SNP markers, representing all nine linkage groups of the Citrus genetic map. We found strong but incomplete preferential pairing between homologues of the same ancestral genome. The proportion of gametes that can be explained by random meiotic chromosome associations (τ) varied significantly between chromosomes, from 0.09 ± 0.02 to 0.47 ± 0.09, respectively, in chromosome 2 and 1. This intermediate inheritance between strict disomy and tetrasomy, with global preferential disomic tendency, resulted in a high level of intergeneric heterozygosity of the diploid gametes. Although limited, intergeneric recombinations occurred, whose observed rates, ranging from 0.09 to 0.29, respectively, in chromosome 2 and 1, were significantly correlated with τ. Such inheritance is of particular interest for rootstock breeding because a large part of the multi-trait value selected at the teraploid parent level is transmitted to the progeny, while the potential for some intergeneric recombination offers opportunities for generating plants with novel allelic combinations that can be targeted by selection.

On the reconstruction of the ancestral bacterial genomes in genus Mycobacterium and Brucella

BackgroundTo reconstruct the evolution history of DNA sequences, novel models of increasing complexity regarding the number of free parameters taken into account in the sequence evolution, as well as faster and more accurate algorithms, and statistical and computational methods, are needed. More particularly, as the principal forces that have led to major structural changes are genome rearrangements (such as translocations, fusions, and so on), understanding their underlying mechanisms, among other things via the ancestral genome reconstruction, are essential. In this problem, since finding the ancestral genomes that minimize the number of rearrangements in a phylogenetic tree is known to be NP-hard for three or more genomes, heuristics are commonly chosen to obtain approximations of the exact solution. The aim of this work is to show that another path is possible.ResultsVarious algorithms and software already deal with the difficult nature of the problem of reconstruction of the ancestral genome, but they do not function with precision, in particular when indels or single nucleotide polymorphisms fall into repeated sequences. In this article, and despite the theoretical NP-hardness of the ancestral reconstruction problem, we show that an exact solution can be found in practice in various cases, encompassing organelles and some bacteria. A practical example proves that an accurate reconstruction, which also allows to highlight homoplasic events, can be obtained. This is illustrated by the reconstruction of ancestral genomes of two bacterial pathogens, belonging in Mycobacterium and Brucella genera.ConclusionsBy putting together automatically reconstructed ancestral regions with handmade ones for problematic cases, we show that an accurate reconstruction of ancestors of the Brucella genus and of the Mycobacterium tuberculosis complex is possible. By doing so, we are able to investigate the evolutionary history of each pathogen by computing their common ancestors. They can be investigated extensively, by studying the gene content evolution over time, the resistance acquisition, and the impacts of mobile elements on genome plasticity.

Scaffolding of Ancient Contigs and Ancestral Reconstruction in a Phylogenetic Framework

Ancestral genome reconstruction is an important task to analyze the evolution of genomes. Recent progress in sequencing ancient DNA led to the publication of so-called paleogenomes and allows the integration of this sequencing data in genome evolution analysis. However, the de novo assembly of ancient genomes is usually fragmented due to DNA degradation over time among others. Integrated phylogenetic assembly addresses the issue of genome fragmentation in the ancient DNA assembly while aiming to improve the reconstruction of all ancient genomes in the phylogeny simultaneously. The fragmented assembly of the ancient genome can be represented as an assembly graph, indicating contradicting ordering information of contigs. In this setting, our approach is to compare the ancient data with extant finished genomes. We generalize a reconstruction approach minimizing the Single-Cut-or-Join rearrangement distance towards multifurcating trees and include edge lengths to improve the reconstruction in practice. This results in a polynomial time algorithm that includes additional ancient DNA data at one node in the tree, resulting in consistent reconstructions of ancestral genomes.