Concerted evolution: molecular mechanism and biological implications

Genetic and epigenetic characteristics of FSHD-associated 4q and 10q D4Z4 that are distinct from non-4q/10q D4Z4 homologs

Facioscapulohumeral dystrophy (FSHD) is one of the most prevalent muscular dystrophies. The majority of FSHD cases are linked to a decreased copy number of D4Z4 macrosatellite repeats on chromosome 4q (FSHD1). Less than 5% of FSHD cases have no repeat contraction (FSHD2), most of which are associated with mutations of SMCHD1. FSHD is associated with the transcriptional derepression of DUX4 encoded within the D4Z4 repeat, and SMCHD1 contributes to its regulation. We previously found that the loss of heterochromatin mark (i.e., histone H3 lysine 9 tri-methylation (H3K9me3)) at D4Z4 is a hallmark of both FSHD1 and FSHD2. However, whether this loss contributes to DUX4 expression was unknown. Furthermore, additional D4Z4 homologs exist on multiple chromosomes, but they are largely uncharacterized and their relationship to 4q/10q D4Z4 was undetermined. We found that the suppression of H3K9me3 results in displacement of SMCHD1 at D4Z4 and increases DUX4 expression in myoblasts. The DUX4 open reading frame (ORF) is disrupted in D4Z4 homologs and their heterochromatin is unchanged in FSHD. The results indicate the significance of D4Z4 heterochromatin in DUX4 gene regulation and reveal the genetic and epigenetic distinction between 4q/10q D4Z4 and the non-4q/10q homologs, highlighting the special role of the 4q/10q D4Z4 chromatin and the DUX4 ORF in FSHD.

Is homologous recombination really an error-free process?

Homologous recombination (HR) is an evolutionarily conserved process that plays a pivotal role in the equilibrium between genetic stability and diversity. HR is commonly considered to be error-free, but several studies have shown that HR can be error-prone. Here, we discuss the actual accuracy of HR. First, we present the product of genetic exchanges (gene conversion, GC, and crossing over, CO) and the mechanisms of HR during double strand break repair and replication restart. We discuss the intrinsic capacities of HR to generate genome rearrangements by GC or CO, either during DSB repair or replication restart. During this process, abortive HR intermediates generate genetic instability and cell toxicity. In addition to genome rearrangements, HR also primes error-prone DNA synthesis and favors mutagenesis on single stranded DNA, a key DNA intermediate during the HR process. The fact that cells have developed several mechanisms protecting against HR excess emphasize its potential risks. Consistent with this duality, several pro-oncogenic situations have been consistently associated with either decreased or increased HR levels. Nevertheless, this versatility also has advantages that we outline here. We conclude that HR is a double-edged sword, which on one hand controls the equilibrium between genome stability and diversity but, on the other hand, can jeopardize the maintenance of genomic integrity. Therefore, whether non-homologous end joining (which, in contrast with HR, is not intrinsically mutagenic) or HR is the more mutagenic process is a question that should be re-evaluated. Both processes can be "Dr. Jekyll" in maintaining genome stability/variability and "Mr. Hyde" in jeopardizing genome integrity.

Glossina palpalis palpalis populations from Equatorial Guinea belong to distinct allopatric clades

Background

Luba is one of the four historical foci of Human African Trypanosomiasis (HAT) on Bioko Island, in Equatorial Guinea. Although no human cases have been detected since 1995, T. b. gambiense was recently observed in the vector Glossina palpalis palpalis. The existence of cryptic species within this vector taxon has been previously suggested, although no data are available regarding the evolutionary history of tsetse flies populations in Bioko.

Methods

A phylogenetic analysis of 60 G. p. palpalis from Luba was performed sequencing three mitochondrial (COI, ND2 and 16S) and one nuclear (rDNA-ITS1) DNA markers. Phylogeny reconstruction was performed by Distance Based, Maximum Likelihood and Bayesian Inference methods.

Results

The COI and ND2 mitochondrial genes were concatenated and revealed 10 closely related haplotypes with a dominant one found in 61.1% of the flies. The sequence homology of the other 9 haplotypes compared to the former ranged from 99.6 to 99.9%. Phylogenetic analysis clearly clustered all island samples with flies coming from the Western African Clade (WAC), and separated from the flies belonging to the Central Africa Clade (CAC), including samples from Mbini and Kogo, two foci of mainland Equatorial Guinea. Consistent with mitochondrial data, analysis of the microsatellite motif present in the ITS1 sequence exhibited two closely related genotypes, clearly divergent from the genotypes previously identified in Mbini and Kogo.

Conclusions

We report herein that tsetse flies populations circulating in Equatorial Guinea are composed of two allopatric subspecies, one insular and the other continental. The presence of these two G. p. palpalis cryptic taxa in Equatorial Guinea should be taken into account to accurately manage vector control strategy, in a country where trypanosomiasis transmission is controlled but not definitively eliminated yet.

Disparate molecular evolution of two types of repetitive DNAs in the genome of the grasshopper Eyprepocnemis plorans

Wide arrays of repetitive DNA sequences form an important part of eukaryotic genomes. These repeats appear to evolve as coherent families, where repeats within a family are more similar to each other than to other orthologous representatives in related species. The continuous homogenization of repeats, through selective and non-selective processes, is termed concerted evolution. Ascertaining the level of variation between repeats is crucial to determining which evolutionary model best explains the homogenization observed for these sequences. Here, for the grasshopper Eyprepocnemis plorans, we present the analysis of intragenomic diversity for two repetitive DNA sequences (a satellite DNA (satDNA) and the 45S rDNA) resulting from the independent microdissection of several chromosomes. Our results show different homogenization patterns for these two kinds of paralogous DNA sequences, with a high between-chromosome structure for rDNA but no structure at all for the satDNA. This difference is puzzling, considering the adjacent localization of the two repetitive DNAs on paracentromeric regions in most chromosomes. The disparate homogenization patterns detected for these two repetitive DNA sequences suggest that several processes participate in the concerted evolution in E. plorans, and that these mechanisms might not work as genome-wide processes but rather as sequence-specific ones.

Organization and biology of the porcine serum amyloid A (SAA) gene cluster: isoform specific responses to bacterial infection

Serum amyloid A (SAA) is a prominent acute phase protein. Although its biological functions are debated, the wide species distribution of highly homologous SAA proteins and their uniform behavior in response to injury or inflammation in itself suggests a significant role for this protein. The pig is increasingly being used as a model for the study of inflammatory reactions, yet only little is known about how specific SAA genes are regulated in the pig during acute phase responses and other responses induced by pro-inflammatory host mediators. We designed SAA gene specific primers and quantified the gene expression of porcine SAA1, SAA2, SAA3, and SAA4 by reverse transcriptase quantitative polymerase chain reaction (RT-qPCR) in liver, spleen, and lung tissue from pigs experimentally infected with the Gram-negative swine specific bacterium Actinobacillus pleuropneumoniae, as well as from pigs experimentally infected with the Gram-positive bacterium Staphylococcus aureus. Our results show that: 1) SAA1 may be a pseudogene in pigs; 2) we were able to detect two previously uncharacterized SAA transcripts, namely SAA2 and SAA4, of which the SAA2 transcript is primarily induced in the liver during acute infection and presumably contributes to circulating SAA in pigs; 3) Porcine SAA3 transcription is induced both hepatically and extrahepatically during acute infection, and may be correlated to local organ affection; 4) Hepatic transcription of SAA4 is markedly induced in pigs infected with A. pleuropneumoniae, but only weakly in pigs infected with S. aureus. These results for the first time establish the infection response patterns of the four porcine SAA genes which will be of importance for the use of the pig as a model for human inflammatory responses, e.g. within sepsis, cancer, and obesity research.

Twenty years of physical mapping of major ribosomal RNA genes across the teleosts: A review of research

Molecular cytogenetic data on the number and position of 45S ribosomal RNA genes (rDNA; located in nucleolus organizing regions, NORs) detected by FISH in 330 species of 77 families and 22 orders of bony fishes (Teleostei) and, additionally, 11 species of basal ray-finned fishes are compiled and analyzed. The portion of species with single rDNA sites in the sample amounts to 72%. The percentage of species with multiple NORs decreases with increasing numbers of rDNA loci per genome, i.e. scarcely 3% of species carry 4 or more rDNA-bearing chromosome pairs. 43% of all rDNA sites analyzed occur terminally on the short arms of chromosomes or constitute them. In general, terminal rDNA sites account for 87% of all examined cases. Interspecific variation in the location of single rDNA sites among related taxa, polymorphisms of multiple NORs in some groups of teleosts and analytical outcomes on the subject are reviewed.

Elongation factor-1α, a putative single-copy nuclear gene, has divergent sets of paralogs in an arachnid

Identification of paralogy in candidate nuclear loci is an important prerequisite in phylogenetics and statistical phylogeography, but one that is often overlooked. One marker commonly assumed to be a single-copy gene and claimed to harbor great utility for inferring recent divergences is elongation factor-1α (EF-1α). To test this hypothesis, we systematically cloned EF-1α in three disjunct populations of the harvestman Metasiro americanus. Here we show that EF-1α has a large number of paralogs in this species. The paralogs do not evolve in a concerted manner, and the paralogs diverged prior to the population divergence. Moreover, the paralogs of M. americanus are not comparable to the highly divergent EF-1α paralogs found in bees and spiders, which are easily recognized and separated through the use of specific primers. We demonstrate statistically that our detection of paralogs cannot be attributed to amplification error. The presence of EF-1α paralogs in M. americanus prevents its use in statistical phylogeography, and the presence of out-paralogs argues against its use in phylogenetic inference among recently diverged clades. These data contradict the common assumption that EF-1α is for most or all taxa a single-copy gene, or that it has a small number of paralogs that are homogenized through gene conversion, unequal crossing over, or other processes.

Employing 454 amplicon pyrosequencing to reveal intragenomic divergence in the internal transcribed spacer rDNA region in fungi

The rDNA internal transcribed spacer (ITS) region has been accepted as a DNA barcoding marker for fungi and is widely used in phylogenetic studies; however, intragenomic ITS variability has been observed in a broad range of taxa, including prokaryotes, plants, animals, and fungi, and this variability has the potential to inflate species richness estimates in molecular investigations of environmental samples. In this study 454 amplicon pyrosequencing of the ITS1 region was applied to 99 phylogenetically diverse axenic single-spore cultures of fungi (Dikarya: Ascomycota and Basidiomycota) to investigate levels of intragenomic variation. Three species (one Basidiomycota and two Ascomycota), in addition to a positive control species known to contain ITS paralogs, displayed levels of molecular variation indicative of intragenomic variation; taxon inflation due to presumed intragenomic variation was ≈9%. Intragenomic variability in the ITS region appears to be widespread but relatively rare in fungi (≈3–5% of species investigated in this study), suggesting this problem may have minor impacts on species richness estimates relative to PCR and/or pyrosequencing errors. Our results indicate that 454 amplicon pyrosequencing represents a powerful tool for investigating levels of ITS intragenomic variability across taxa, which may be valuable for better understanding the fundamental mechanisms underlying concerted evolution of repetitive DNA regions.

Non-concerted ITS evolution in fungi, as revealed from the important medicinal fungus Ophiocordyceps sinensis

The internal transcribed spacer (ITS) of nuclear ribosomal DNA (nrDNA) has been widely used as a molecular marker in phylogenetic studies and has been selected as a DNA barcode for fungi. It is generally believed that nrDNA conforms to concerted evolution in most eukaryotes; however, intraindividual-intraspecific polymorphisms of this region were reported in various organisms, suggesting a non-concerted evolutionary process. In Ophiocordyceps sinensis, one of the most valuable medicinal fungi, a remarkable variation of the ITS region has been revealed. Some highly divergent sequences were thought to represent cryptic species, different species or genotypes in previous studies. To clarify the unusual ITS polymorphisms observed in O. sinensis, specific primers were designed to amplify ITS paralogs from pure cultures of both single-ascospore and tissue isolates in this study. All of the available ITS sequences, including those generated by this group and those in GenBank, were analyzed. Several AT-biased ITS paralogs were classified as pseudogenes based on their nucleotide compositions, secondary structures and minimum free energies of their 5.8S rRNAs, substitution rates, phylogenetic positions and gene expression analyses. Furthermore, ITS pseudogenes were amplified with specific primers from 10 of the 28 strains tested, including eight single-ascospore and two tissue isolates. Divergent ITS paralogs were proved to coexist in individual genomes, suggesting a non-concerted mechanism of evolution in the ITS region of O. sinensis. The hypotheses that divergent ITS paralogs represent cryptic or other species or different genotypes were thus rejected.

Concerted evolution in the ribosomal RNA cistron

BACKGROUND

Gene conversion is the mechanism proposed to be responsible for the homogenization of multigene families such as the nuclear ribosomal gene clusters. This concerted evolutionary process prevents individual genes in gene clusters from accumulating mutations. The mechanism responsible for concerted evolution is not well understood but recombination during meiosis has been hypothesized to play a significant role in this homogenization. In this study we tested the hypothesis of unequal crossing over playing a significant role in gene conversion events within the ribosomal RNA cistron during meiosis, mitosis or both life stages in the fungal tree pathogen Ceratocystis manginecans.

METHODS

Ceratocystis manginecans, a haploid ascomycete, reproduces homothallically and was found to have two distinct sequences within the internally transcribed spacer (ITS) region of the ribosomal RNA cistron. The different ITS types were scored using PCR-RFLP assays and chi-square analyses to determine the level of significance of the changes in the ratios of the ITS types.

RESULTS

The relative ratios of the two ITS sequence types changed when the fungal isolates were cultured vegetatively or allowed to produced sexual structures and spores. These active changes were shown to occur more frequently during meiosis than mitosis.

CONCLUSION

The evidence presented provides concrete support for homogenization in the rRNA gene clusters found in this fungus and that the most reasonable explanation for this process is unequal crossing over.

Duplication and concerted evolution in a master sex determiner under balancing selection

The transformer (tra) gene is a key regulator in the signalling hierarchy controlling all aspects of somatic sexual differentiation in Drosophila and other insects. Here, we show that six of the seven sequenced ants have two copies of tra. Surprisingly, the two paralogues are always more similar within species than among species. Comparative sequence analyses indicate that this pattern is owing to the ongoing concerted evolution after an ancestral duplication rather than independent duplications in each of the six species. In particular, there was strong support for inter-locus recombination between the paralogues of the ant Atta cephalotes. In the five species where the location of paralogues is known, they are adjacent to each other in four cases and separated by only few genes in the fifth case. Because there have been extensive genomic rearrangements in these lineages, this suggests selection acting to conserve their synteny. In three species, we also find a signature of positive selection in one of the paralogues. In three bee species where information is available, the tra gene is also duplicated, the copies are adjacent and in at least one species there was recombination between paralogues. These results suggest that concerted evolution plays an adaptive role in the evolution of this gene family.

ITS2 characterization and Anopheles species identification of the subgenus Cellia

In Mizoram, the origin and molecular nature of Anopheles species is poorly understood, despite the region having high malarial incidence and rich biodiversity. A diagnostic PCR assay for distinguishing the Cellia subgenera members of Anopheles species was developed based on the interspecific ITS2 variation. No intraspecific variation was found and the size (362-604bp) and GC content (48.8-58.9%) of the ITS2 were highly variable among Anophelines. The ITS2 of A. vagus is significantly longer than those of other Anopheles species. Significant relationship was observed among repeats, minimum free energy and RNA secondary structures. Different types of microsatellites were identified and among them dinucleotide, pentanucleotide and polynucleotide microsatellites were predominant. Variation in the length of the ITS2 between species was due to indels in simple repeats. Four domain types of RNA secondary structures were identified and the lowest free energy values were predicted using the computer software, RNAfold. Types I and II were observed only in Neocellia and Myzomyia series and Types III and IV were common in Neocellia and Pyretophorus series. ITS2-based PCR protocol provides a means for vector ecologists, malaria epidemiologists and control personnel to accurately identify members of the subgenera Cellia and a better understanding of their genomic status in Mizoram.

The fester locus in Botryllus schlosseri experiences selection

Background

Allorecognition, the ability of an organism to distinguish self from non-self, occurs throughout the entire tree of life. Despite the prevalence and importance of allorecognition systems, the genetic basis of allorecognition has rarely been characterized outside the well-known MHC (Major Histocompatibility Complex) in vertebrates and SI (Self-Incompatibility) in plants. Where loci have been identified, their evolutionary history is an open question. We have previously identified the genes involved in self/non-self recognition in the colonial ascidian Botryllus schlosseri, and we can now begin to investigate their evolution. In B. schlosseri, colonies sharing 1 or more alleles of a gene called FuHC (Fusion Histocompatibility) will fuse. Protein products of a locus called fester, located ~300 kb from FuHC, have been shown to play multiple roles in the histocompatibility reaction, as activating and/or inhibitory receptors. We test whether the proteins encoded by this locus are evolving neutrally or are experiencing balancing, directional, or purifying selection.

Results

Nearly all of the variation in the fester locus resides within populations. The 13 housekeeping genes (12 nuclear genes and mitochondrial cytochrome oxidase I) have substantially more structure among populations within groups and among groups than fester. All polymorphism statistics (Tajima's D, Fu and Li's D* and F*) are significantly negative for the East Coast A-type alleles, and Fu and Li's F* statistic is significantly negative for the West Coast A-type alleles. These results are likely due to selection rather than demography, given that 10 of the housekeeping loci have no populations with significant values for any of the polymorphism statistics. The majority of codons in the fester proteins have ω values < 1, but 15–27 codons have > 95% posterior probability of ω values > 1.

Conclusion

Fester proteins are evolving non-neutrally. The polymorphism statistics are consistent with either purifying selection or directional selection. The ω statistics show that the majority of the protein is experiencing purifying selection (ω < 1), but that 15–27 codons are undergoing either balancing or directional selection: ω > 1 is compatible with either scenario. The distribution of variation within and among populations points towards balancing selection and away from directional selection. While these data do not provide unambiguous support for a specific type of selection, they contribute to our evolutionary understanding of a critical biological process by determining the forces that affect loci involved in allorecognition.

A nuclear single-nucleotide polymorphism (SNP) potentially useful for the separation of Rhodnius prolixus from members of the Rhodnius robustus cryptic species complex (Hemiptera: Reduviidae)

The design and application of rational strategies that rely on accurate species identification are pivotal for effective vector control. When morphological identification of the target vector species is impractical, the use of molecular markers is required. Here we describe a non-coding, single-copy nuclear DNA fragment that contains a single-nucleotide polymorphism (SNP) with the potential to distinguish the important domestic Chagas disease vector, Rhodnius prolixus, from members of the four sylvatic Rhodnius robustus cryptic species complex. A total of 96 primer pairs obtained from whole genome shotgun sequencing of the R. prolixus genome (12,626 random reads) were tested on 43 R. prolixus and R. robustus s.l. samples. One of the seven amplicons selected (AmpG) presented a SNP, potentially diagnostic for R. prolixus, on the 280th site. The diagnostic nature of this SNP was then confirmed based on the analysis of 154 R. prolixus and R. robustus s.l. samples representing the widest possible geographic coverage. The results of a 60% majority-rule Bayesian consensus tree and a median-joining network constructed based on the genetic variability observed reveal the paraphyletic nature of the R. robustus species complex, with respect to R. prolixus. The AmpG region is located in the fourth intron of the Transmembrane protein 165 gene, which seems to be in the R. prolixus X chromosome. Other possible chromosomal locations of the AmpG region in the R. prolixus genome are also presented and discussed.

The taxon-specific paralogs of grapevine PRLIP genes are highly induced upon powdery mildew infection

PRLIP (pathogenesis-related lipase) is a gene family encoding class 3 lipase-like proteins originally described and first characterized in Arabidopsis thaliana. Nine paralog genes of Arabidopsis can be separated into two groups based on expression characteristics and pathogen responses. Genes of Group 1 are clustered on chromosome 5 and show either high inducibility to different stress hormones and in response to pathogen attack or are undetectable at the transcript level. Group 2 contains the remaining genes, spread over the genome and are expressed constitutively in all the tissues tested. The aim of the present study was to determine the distribution of these two groups among plants, and to verify their differential expression. Orthologs of constitutively active members (Group 2) were found in all angiosperms, with available genome sequences. They are referred to as "core PRLIPs". In contrast, the gene cluster containing the pathogen-inducible PRLIPs (Group 1) was unique for Arabidopsis. Among other angiosperms, grapevine also possesses such a unique genome-specific group of PRLIP genes. To investigate whether these genes are also counterparts in pathogen responses, their expression pattern was tested under stress conditions. Two of the specific Vitis PRLIPs were highly induced in response to both powdery mildew infection and benzothiadiazole (BTH) treatment. Core Vitis PRLIPs, however, were not responsive to either pathogen attack or the chemical inducer. Our data provide insights into the distribution of a pathogenesis-related gene family in different plant lineages, and might reveal common characteristics with other inducible defense-related gene families.

Evolutionary patterns of non-coding RNAs

A plethora of new functions of non-coding RNAs (ncRNAs) have been discovered in past few years. In fact, RNA is emerging as the central player in cellular regulation, taking on active roles in multiple regulatory layers from transcription, RNA maturation, and RNA modification to translational regulation. Nevertheless, very little is known about the evolution of this "Modern RNA World" and its components. In this contribution, we attempt to provide at least a cursory overview of the diversity of ncRNAs and functional RNA motifs in non-translated regions of regular messenger RNAs (mRNAs) with an emphasis on evolutionary questions. This survey is complemented by an in-depth analysis of examples from different classes of RNAs focusing mostly on their evolution in the vertebrate lineage. We present a survey of Y RNA genes in vertebrates and study the molecular evolution of the U7 snRNA, the snoRNAs E1/U17, E2, and E3, the Y RNA family, the let-7 microRNA (miRNA) family, and the mRNA-like evf-1 gene. We furthermore discuss the statistical distribution of miRNAs in metazoans, which suggests an explosive increase in the miRNA repertoire in vertebrates. The analysis of the transcription of ncRNAs suggests that small RNAs in general are genetically mobile in the sense that their association with a hostgene (e.g. when transcribed from introns of a mRNA) can change on evolutionary time scales. The let-7 family demonstrates, that even the mode of transcription (as intron or as exon) can change among paralogous ncRNA.

Genetic variants in REC8, RNF212, and PRDM9 influence male recombination in cattle

We use >250,000 cross-over events identified in >10,000 bovine sperm cells to perform an extensive characterization of meiotic recombination in male cattle. We map Quantitative Trait Loci (QTL) influencing genome-wide recombination rate, genome-wide hotspot usage, and locus-specific recombination rate. We fine-map three QTL and present strong evidence that genetic variants in REC8 and RNF212 influence genome-wide recombination rate, while genetic variants in PRDM9 influence genome-wide hotspot usage.

Homologous recombination is an essential cellular process that determines proper chromosome segregation during meiosis, affects fertility, and influences evolvability. Nevertheless, the components of the recombination apparatus remain incompletely characterized in mammals. One approach to identify such components is to identify the genes that underlie inherited variation in recombination phenotypes. In addition to providing mechanistic insights, this would allow the study of the evolutionary forces that shape the recombination process. In this paper, we take advantage of genotypes for 50,000 genome-wide SNP markers to measure four recombination phenotypes (genome-wide recombination rate, genome-wide hotspot usage, locus-specific recombination rate, genome-wide cross-over interference) for >750 bulls on the basis of >250,000 cross-overs detected in sperm cells transmitted to >10,000 sons. We quantify the heritability and scan the genome for Quantitative Trait Loci (QTL) influencing each one of these recombination phenotypes. We perform a detailed genetic analysis of three such QTL, thereby providing evidence that genetic variants in REC8 and RNF212 influence genome-wide recombination rate, while genetic variants in an X-linked PRDM9 paralogue influence genome-wide hotspot usage.

Tracing back the nascence of a new sex-determination pathway to the ancestor of bees and ants

In several Hymenoptera, sexual fate is determined by the allelic composition at the complementary sex-determiner locus, a sex-determination mechanism that can strongly affect population dynamics. To date, the molecular identification of complementary sex determiner has only been achieved in the honeybee, where the complementary sex-determiner gene was reported to have arisen from duplication of the feminizer gene. Strikingly, the complementary sex-determiner gene was also proposed to be unique to the honeybee lineage. Here we identify feminizer and complementary sex-determiner orthologues in bumble bees and ants. We further demonstrate that the duplication of feminizer that produced complementary sex determiner occurred before the divergence of Aculeata species (~120 Myr ago). Finally, we provide evidence that the two genes evolved concertedly through gene conversion, complementary sex-determiner evolution being additionally shaped by mosaic patterns of selection. Thus, the complementary sex-determiner gene likely represents the molecular basis for single locus-complementary sex determination in the Aculeata infra-order, and possibly, in the entire Hymenoptera order.

In several Hymenoptera species - ants, bees and wasps - sexual fate is determined by the allelic composition at the complementary sex - determiner locus. This study identifies the honeybee complementary sex - determiner in bumble bee and ant orthologues, previously thought to be unique to the honeybee lineage.

Ancient gene paralogy may mislead inference of plastid phylogeny

Because of its ancient origin more than 1 billion years ago, the highly reduced plastid genomes of Plantae (e.g., plant chloroplasts) provide limited insights into the initial stages of endosymbiont genome reduction. The photosynthetic amoeba Paulinella provides a more useful model to study this process because its alpha-cyanobacterium-derived plastid originated ∼60 Ma and the genome still contains ∼1,000 genes. Here, we compared and contrasted features associated with genome reduction due to primary endosymbiosis in Paulinella plastids and in marine, free-living strains of the picocyanobacterium, Prochlorococcus. Both types of genomes show gene inactivation, concerted evolution, and contraction of gene families that impact highly conserved single-copy phylogenetic markers in the plastid such as psbA, psbC, and psbD. Our data suggest that these photosystem II genes may provide misleading phylogenetic signal because each of the constituent Plantae lineages has likely undergone a different, independent series of events that led to their reduction to a single copy. This issue is most problematic for resolving basal Plantae relationships when differential plastid gene loss was presumably ongoing, as we observe in Paulinella species. Our work uncovers a key, previously unappreciated aspect of organelle genome reduction and demonstrates "work-in-progress" models such as Paulinella to be critical to gain a fuller understanding of algal and plant genome evolution.

A novel method for simultaneous Enterococcus species identification/typing and van genotyping by high resolution melt analysis

In order to develop a typing and identification method for van gene containing Enterococcus faecium, two multiplex PCR reactions were developed for use in HRM-PCR (High Resolution Melt-PCR): (i) vanA, vanB, vanC, vanC23 to detect van genes from different Enterococcus species; (ii) ISR (intergenic spacer region between the 16S and 23S rRNA genes) to detect all Enterococcus species and obtain species and isolate specific HRM curves. To test and validate the method three groups of isolates were tested: (i) 1672 Enterococcus species isolates from January 2009 to December 2009; (ii) 71 isolates previously identified and typed by PFGE (pulsed-field gel electrophoresis) and MLST (multi-locus sequence typing); and (iii) 18 of the isolates from (i) for which ISR sequencing was done. As well as successfully identifying 2 common genotypes by HRM from the Austin Hospital clinical isolates, this study analysed the sequences of all the vanB genes deposited in GenBank and developed a numerical classification scheme for the standardised naming of these vanB genotypes. The identification of Enterococcus faecalis from E. faecium was reliable and stable using ISR PCR. The typing of E. faecium by ISR PCR: (i) detected two variable peaks corresponding to different copy numbers of insertion sequences I and II corresponding to peak I and II respectively; (ii) produced 7 melt profiles for E. faecium with variable copy numbers of sequences I and II; (iii) demonstrated stability and instability of peak heights with equal frequency within the patient sample (36.4±4.5 days and 38.6±5.8 days respectively for 192 patients); (iv) detected ISR-HRM types with as much discrimination as PFGE and more than MLST; and (v) detected ISR-HRM types that differentiated some isolates that were identical by PFGE and MLST. In conjunction with the rapid and accurate van genotyping method described here, this ISR-HRM typing and identification method can be used as a stable identification and typing method with predictable instability based on recombination and concerted evolution of the rrn operon that will complement existing typing methods.

Effect of individual suppression of RBCS multigene family on Rubisco contents in rice leaves

In higher plants, a small subunit of Rubisco is encoded for by an RBCS multigene family in the nuclear genome. However, it is unknown how each multigene member contributes to the accumulation of Rubisco holoenzyme. Here, four RBCS genes that are highly expressed in leaf blaedes of rice (Oryza sativa L.) were individually suppressed by RNAi, and the effects on leaf Rubisco content were examined at seedling, vegetative and reproductive stages. Rubisco contents in each transgenic line declined irrespective of growth stage, and the ratios of Rubisco-N to total N were 66-96% of wild-type levels. The mRNA levels of the suppressed RBCS genes declined significantly, whereas those of the unsuppressed ones did not change drastically. These results indicate that four RBCS genes all contribute to accumulation of Rubisco holoenzyme irrespective of growth stage and that suppression of one RBCS gene is not fully compensated by other RBCS genes. Additionally, the mRNA levels of the large subunit of Rubisco showed a change similar to that of total RBCS mRNA level irrespective of genotype and growth stage. These results suggest that gene expression of RBCS and RBCL is regulated in a coordinated manner at the transcript level in rice.

The dispanins: a novel gene family of ancient origin that contains 14 human members

The Interferon induced transmembrane proteins (IFITM) are a family of transmembrane proteins that is known to inhibit cell invasion of viruses such as HIV-1 and influenza. We show that the IFITM genes are a subfamily in a larger family of transmembrane (TM) proteins that we call Dispanins, which refers to a common 2TM structure. We mined the Dispanins in 36 eukaryotic species, covering all major eukaryotic groups, and investigated their evolutionary history using Bayesian and maximum likelihood approaches to infer a phylogenetic tree. We identified ten human genes that together with the known IFITM genes form the Dispanin family. We show that the Dispanins first emerged in eukaryotes in a common ancestor of choanoflagellates and metazoa, and that the family later expanded in vertebrates where it forms four subfamilies (A-D). Interestingly, we also find that the family is found in several different phyla of bacteria and propose that it was horizontally transferred to eukaryotes from bacteria in the common ancestor of choanoflagellates and metazoa. The bacterial and eukaryotic sequences have a considerably conserved protein structure. In conclusion, we introduce a novel family, the Dispanins, together with a nomenclature based on the evolutionary origin.

Nucleotide and phylogenetic analyses of the Chlamydia trachomatis ompA gene indicates it is a hotspot for mutation

Background

Serovars of the human pathogen Chlamydia trachomatis occupy one of three specific tissue niches. Genomic analyses indicate that the serovars have a phylogeny congruent with their pathobiology and have an average substitution rate of less than one nucleotide per kilobase. In contrast, the gene that determines serovar specificity, ompA, has a phylogenetic association that is not congruent with tissue tropism and has a degree of nucleotide variability much higher than other genomic loci. The ompA gene encodes the major surface-exposed antigenic determinant, and the observed nucleotide diversity at the ompA locus is thought to be due to recombination and host immune selection pressure. The possible contribution of a localized increase in mutation rate, however, has not been investigated.

Results

Nucleotide diversity and phylogenetic relationships of the five constant and four variable domains of the ompA gene, as well as several loci surrounding ompA, were examined for each serovar. The loci flanking the ompA gene demonstrated that nucleotide diversity increased monotonically as ompA is approached and that their gene trees are not congruent with either ompA or tissue tropism. The variable domains of the ompA gene had a very high level of non-synonymous change, which is expected as these regions encode the surface-exposed epitopes and are under positive selection. However, the synonymous changes are clustered in the variable regions compared to the constant domains; if hitchhiking were to account for the increase in synonymous changes, these substitutions should be more evenly distributed across the gene. Recombination also cannot entirely account for this increase as the phylogenetic relationships of the constant and variable domains are congruent with each other.

Conclusions

The high number of synonymous substitutions observed within the variable domains of ompA appears to be due to an increased mutation rate within this region of the genome, whereas the increase in nucleotide substitution rate and the lack of phylogenetic congruence in the regions flanking ompA are characteristic motifs of gene conversion. Together, the increased mutation rate in the ompA gene, in conjunction with gene conversion and positive selection, results in a high degree of variability that promotes host immune evasion.

Evaluation of the internal transcribed spacer 2 (ITS2) as a molecular marker for phylogenetic inference using sequence and secondary structure information in blow flies (Diptera: Calliphoridae)

The internal transcribed spacer 2 (ITS2) is a small non-coding region located inside the nuclear ribosomal DNA cluster. ITS2 sequence variability is thought to be appropriate to differentiate species and for phylogenetic reconstructions analyses, which can be further improved if structural information is considered. We evaluated the potential of ITS2 as a molecular marker for phylogenetic inference in Calliphoridae (Diptera: Brachycera) using a broad range of inference methods and different substitution models, accounting or not for structural information. Sequence analyses revealed a hierarchically organized pattern of sequence variation and a small level of nucleotide substitution saturation. Intragenomic variation due to small sequence repeats was found mainly in the most variable domain (IV), but it has no significant impact on the phylogenetic signal at the species level. Inferred secondary structures revealed that GC pairs are more frequently found flanking bulges and loops regions in more conserved domains, thus ensuring structure stability. In the phylogenetic analyses, the use of substitution models accounting for structural information significantly improves phylogenetic inference in both neighbour-joining and Bayesian analyses, although the former provides limited resolution for dealing with highly divergent sequences. For Bayesian analyses, a significant improvement in likelihood was observed when considering structure information, although with small changes in topology and overall support, probably reflecting better evolutionary rates estimates. Based on these findings, ITS2 is a suitable molecular marker for phylogenetic analyses in Calliphoridae, at both species and generic level.

Polyphasic analysis of intraspecific diversity in Epicoccum nigrum warrants reclassification into separate species

BACKGROUND

Epicoccum nigrum Link (syn. E. purpurascens Ehrenb. ex Schlecht) is a saprophytic ascomycete distributed worldwide which colonizes a myriad of substrates. This fungus has been known as a biological control agent for plant pathogens and produces a variety of secondary metabolites with important biological activities as well as biotechnological application. E. nigrum produces darkly pigmented muriform conidia on short conidiophores on sporodochia and is a genotypically and phenotypically highly variable species. Since different isolates identified as E. nigrum have been evaluated as biological control agents and used for biocompound production, it is highly desirable that this species name refers to only one lineage. However, according to morphological and genetic variation, E. nigrum present two genotypes that may comprise more than one species.

METHODOLOGY/PRINCIPAL FINDINGS

We report the application of combined molecular (ITS and β-tubulin gene sequence analysis, PCR-RFLP and AFLP techniques), morphometric, physiological, genetic compatibility and recombination analysis to study the taxonomic relationships within an endophytic population that has been identified as E. nigrum. This combined analysis established two genotypes showing morphological, physiological and genetic divergence as well as genetic incompatibility characterized by colony inhibition, strongly indicating that these genotypes correspond to different species. Genotype 1 corresponds to E. nigrum while genotype 2 represents a new species, referred to in this study as Epicoccum sp.

CONCLUSIONS/SIGNIFICANCE

This research contributes to the knowledge of the Epicoccum genus and asserts that the classification of E. nigrum as a single variable species should be reassessed. In fact, based on the polyphasic approach we suggest the occurrence of cryptic species within E. nigrum and also that many of the sequences deposited as E. nigrum in GenBank and culture collection of microbial strains should be reclassified, including the reference strain CBS 161.73 sequenced in this work. In addition, this study provides valuable tools for differentiation of Epicoccum species.

Genome-level homology and phylogeny of Shewanella (Gammaproteobacteria: lteromonadales: Shewanellaceae)

BACKGROUND

The explosion in availability of whole genome data provides the opportunity to build phylogenetic hypotheses based on these data as well as the ability to learn more about the genomes themselves. The biological history of genes and genomes can be investigated based on the taxomonic history provided by the phylogeny. A phylogenetic hypothesis based on complete genome data is presented for the genus Shewanella (Gammaproteobacteria: Alteromonadales: Shewanellaceae). Nineteen taxa from Shewanella (16 species and 3 additional strains of one species) as well as three outgroup species representing the genera Aeromonas (Gammaproteobacteria: Aeromonadales: Aeromonadaceae), Alteromonas (Gammaproteobacteria: Alteromonadales: Alteromonadaceae) and Colwellia (Gammaproteobacteria: Alteromonadales: Colwelliaceae) are included for a total of 22 taxa.

RESULTS

Putatively homologous regions were found across unannotated genomes and tested with a phylogenetic analysis. Two genome-wide data-sets are considered, one including only those genomic regions for which all taxa are represented, which included 3,361,015 aligned nucleotide base-pairs (bp) and a second that additionally includes those regions present in only subsets of taxa, which totaled 12,456,624 aligned bp. Alignment columns in these large data-sets were then randomly sampled to create smaller data-sets. After the phylogenetic hypothesis was generated, genome annotations were projected onto the DNA sequence alignment to compare the historical hypothesis generated by the phylogeny with the functional hypothesis posited by annotation.

CONCLUSIONS

Individual phylogenetic analyses of the 243 locally co-linear genome regions all failed to recover the genome topology, but the smaller data-sets that were random samplings of the large concatenated alignments all produced the genome topology. It is shown that there is not a single orthologous copy of 16S rRNA across the taxon sampling included in this study and that the relationships among the multiple copies are consistent with 16S rRNA undergoing concerted evolution. Unannotated whole genome data can provide excellent raw material for generating hypotheses of historical homology, which can be tested with phylogenetic analysis and compared with hypotheses of gene function.

Intragenomic variation in the ITS rDNA region obscures phylogenetic relationships and inflates estimates of operational taxonomic units in genus Laetiporus

Regions of rDNA are commonly used to infer phylogenetic relationships among fungal species and as DNA barcodes for identification. These regions occur in large tandem arrays, and concerted evolution is believed to reduce intragenomic variation among copies within these arrays, although some variation still might exist. Phylogenetic studies typically use consensus sequencing, which effectively conceals most intragenomic variation, but cloned sequences containing intragenomic variation are becoming prevalent in DNA databases. To understand effects of using cloned rDNA sequences in phylogenetic analyses we amplified and cloned the ITS region from pure cultures of six Laetiporus species and one Wolfiporia species (Basidiomycota, Polyporales). An average of 66 clones were selected randomly and sequenced from 21 cultures, producing a total of 1399 interpretable sequences. Significant variation (≥ 5% variation in sequence similarity) was observed among ITS copies within six cultures from three species clades (L. cincinnatus, L. sp. clade J, and Wolfiporia dilatohypha) and phylogenetic analyses with the cloned sequences produced different trees relative to analyses with consensus sequences. Cloned sequences from L. cincinnatus fell into more than one species clade and numerous cloned L. cincinnatus sequences fell into entirely new clades, which if analyzed on their own most likely would be recognized as "undescribed" or "novel" taxa. The use of a 95% cut off for defining operational taxonomic units (OTUs) produced seven Laetiporus OTUs with consensus ITS sequences and 20 OTUs with cloned ITS sequences. The use of cloned rDNA sequences might be problematic in fungal phylogenetic analyses, as well as in fungal bar-coding initiatives and efforts to detect fungal pathogens in environmental samples.

Incomplete concerted evolution and reproductive isolation at the rDNA locus uncovers nine cryptic species within Anopheles longirostris from Papua New Guinea

Background

Nuclear ribosomal DNA (rDNA) genes and transcribed spacers are highly utilized as taxonomic markers in metazoans despite the lack of a cohesive understanding of their evolution. Here we follow the evolution of the rDNA second internal transcribed spacer (ITS2) and the mitochondrial DNA cytochrome oxidase I subunit in the malaria mosquito Anopheles longirostris from Papua New Guinea (PNG). This morphospecies inhabits a variety of ecological environments indicating that it may comprise a complex of morphologically indistinguishable species. Using collections from over 70 sites in PNG, the mtDNA was assessed via direct DNA sequencing while the ITS2 was assessed at three levels - crude sequence variation through restriction digest, intragenomic copy variant organisation (homogenisation) through heteroduplex analysis and DNA sequencing via cloning.

Results

Genetic evaluation of over 300 individuals revealed that A. longirostris comprises eight ITS2 PCR-RFLP genotypes and nine ITS2 heteroduplex genotypes showing distinct copy variant organization profiles after PCR amplification. Seven of these nine genotypes were found to be sympatric with other genotypes. Phylogenetic analysis of cloned ITS2 PCR products and mtDNA COI confirmed all nine clades with evidence of reproductive isolation at the rDNA locus. Compensatory base changes in the ITS2 secondary structure or in pseudoknots were absent when closely related species were assessed. Individuals from each ITS2 genotype showed the same copy variant heteroduplex profile suggesting that the rDNA array is fixed within each genotype.

Conclusion

The centromere-proximal position of the rDNA array in Anopheles mosquitoes has probably reduced interchromosomal recombination leaving intrachromosomal events responsible for the observed pattern of concerted evolution we see in these mosquitoes. The stability of these intragenomic ITS2 copy variants within individuals and interbreeding populations suggests that rDNA is moving as a single evolutionary unit through natural populations to fixation and has provided a complementary diagnostic tool to the restriction digest for studying genetic discontinuities and species boundaries. In this, the utility of the ITS2 as a universal taxonomic marker is probably contingent on several factors pertaining to spacer dimensions and the genomic location of the rDNA array with respect to recombination and proximity to regions potentially under selection.

Molecular species identification of Central European ground beetles (Coleoptera: Carabidae) using nuclear rDNA expansion segments and DNA barcodes

BACKGROUND

The identification of vast numbers of unknown organisms using DNA sequences becomes more and more important in ecological and biodiversity studies. In this context, a fragment of the mitochondrial cytochrome c oxidase I (COI) gene has been proposed as standard DNA barcoding marker for the identification of organisms. Limitations of the COI barcoding approach can arise from its single-locus identification system, the effect of introgression events, incomplete lineage sorting, numts, heteroplasmy and maternal inheritance of intracellular endosymbionts. Consequently, the analysis of a supplementary nuclear marker system could be advantageous.

RESULTS

We tested the effectiveness of the COI barcoding region and of three nuclear ribosomal expansion segments in discriminating ground beetles of Central Europe, a diverse and well-studied invertebrate taxon. As nuclear markers we determined the 18S rDNA: V4, 18S rDNA: V7 and 28S rDNA: D3 expansion segments for 344 specimens of 75 species. Seventy-three species (97%) of the analysed species could be accurately identified using COI, while the combined approach of all three nuclear markers provided resolution among 71 (95%) of the studied Carabidae.

CONCLUSION

Our results confirm that the analysed nuclear ribosomal expansion segments in combination constitute a valuable and efficient supplement for classical DNA barcoding to avoid potential pitfalls when only mitochondrial data are being used. We also demonstrate the high potential of COI barcodes for the identification of even closely related carabid species.

Population genetic structure of the malaria vector Anopheles nili in sub-Saharan Africa

Background

Anopheles nili is a widespread efficient vector of human malaria parasites in the humid savannas and forested areas of sub-Saharan Africa. Understanding An. nili population structure and gene flow patterns could be useful for the development of locally-adapted vector control measures.

Methods

Polymorphism at eleven recently developed microsatelitte markers, and sequence variation in four genes within the 28s rDNA subunit (ITS2 and D3) and mtDNA (COII and ND4) were assessed to explore the level of genetic variability and differentiation among nine populations of An. nili from Senegal, Ivory Coast, Burkina Faso, Nigeria, Cameroon and the Democratic Republic of Congo (DRC).

Results

All microsatellite loci successfully amplified in all populations, showing high and very similar levels of genetic diversity in populations from West Africa and Cameroon (mean Rs = 8.10-8.88, mean He = 0.805-0.849) and much lower diversity in the Kenge population from DRC (mean Rs = 5.43, mean He = 0.594). Bayesian clustering analysis of microsatellite allelic frequencies revealed two main genetic clusters in the dataset. The first one included only the Kenge population and the second grouped together all other populations. High Fst estimates based on microsatellites (Fst > 0.118, P < 0.001) were observed in all comparisons between Kenge and all other populations. By contrast, low Fst estimates (Fst < 0.022, P < 0.05) were observed between populations within the second cluster. The correlation between genetic and geographic distances was weak and possibly obscured by demographic instability. Sequence variation in mtDNA genes matched these results, whereas low polymorphism in rDNA genes prevented detection of any population substructure at this geographical scale.

Conclusion

Overall, high genetic homogeneity of the An. nili gene pool was found across its distribution range in West and Central Africa, although demographic events probably resulted in a higher level of genetic isolation in the marginal population of Kenge (DRC). The role of the equatorial forest block as a barrier to gene flow and the implication of such findings for vector control are discussed.

Chromosomal conservation and sequence diversity of ribosomal RNA genes of two distant Oryza species

Contrary to the chromosomal polymorphism of 45S ribosomal genes (45S rDNA) loci in other Oryza species, each of Oryza australiensis and Oryza brachyantha has only one 45S rDNA locus at the most conserved position of 45S rDNAs in Oryza. O. australiensis and O. brachyantha are known phylogenetically distant and have extremely different genome sizes among diploid Oryza species. This study reveals that the sequences and organizations of intergenic spacer (IGS) for 45S rDNA of both O. australiensis and O. brachyantha are different from other Oryza species. The IGS of O. australiensis contains 13 tandem repeats and only one transcriptional initiation site, while there are four tandem repeats and three transcriptional initiation sites in the IGS of O. brachyantha. Our results suggest different evolution processes of orthologous rDNA loci in the genus Oryza. Here we also demonstrate an efficient strategy to study locus-specific IGS before whole genome sequences data are available.

Comparative genomics and evolution of the alpha-defensin multigene family in primates

Defensin genes encode small cationic antimicrobial peptides that form an important part of the innate immune system. They are divided into three families, alpha (α), beta (β), and theta (), according to arrangement of the disulfide bonding pattern between cysteine residues. Considering the functional importance of defensins, investigators have studied the evolution and the genomic organization of defensin genes. However, these studies have been restricted mainly to β-defensins. To understand the evolutionary dynamics of α-defensin genes among primates, we identified the α-defensin repertoires in human, chimpanzee, orangutan, macaque, and marmoset. The α-defensin genes in primates can be classified into three phylogenetic classes (class I, II, and III). The presence of all three classes in the marmoset indicates that their divergence occurred before the separation of New World and Old World monkeys. Comparative analysis of the α-defensin genomic clusters suggests that the makeup of the α-defensin gene repertoires between primates is quite different, as their genes have undergone dramatic birth-and-death evolution. Analysis of the encoded peptides of the α-defensin genes indicates that despite the overall high level of sequence divergence, certain amino acid residues or motifs are conserved within and between the three phylogenetic classes. The evolution of α-defensins in primates, therefore, appears to be governed by two opposing evolutionary forces. One force stabilizes specific amino acid residues and motifs to preserve the functional and structural integrity of the molecules and the other diversifies the sequences generating molecules with a wide range of activities against a large number of pathogens.

The evolutionary history of the Arabidopsis lyrata complex: a hybrid in the amphi-Beringian area closes a large distribution gap and builds up a genetic barrier

BACKGROUND

The genomes of higher plants are, on the majority, polyploid, and hybridisation is more frequent in plants than in animals. Both polyploidisation and hybridisation contribute to increased variability within species, and may transfer adaptations between species in a changing environment. Studying these aspects of evolution within a diversified species complex could help to clarify overall spatial and temporal patterns of plant speciation. The Arabidopsis lyrata complex, which is closely related to the model plant Arabidopsis thaliana, is a perennial, outcrossing, herbaceous species complex with a circumpolar distribution in the Northern Hemisphere as well as a disjunct Central European distribution in relictual habitats. This species complex comprises three species and four subspecies, mainly diploids but also several tetraploids, including one natural hybrid. The complex is ecologically, but not fully geographically, separated from members of the closely related species complex of Arabidopsis halleri, and the evolutionary histories of both species compexes have largely been influenced by Pleistocene climate oscillations.

RESULTS

Using DNA sequence data from the nuclear encoded cytosolic phosphoglucoisomerase and Internal Transcribed Spacers 1 and 2 of the ribosomal DNA, as well as the trnL/F region from the chloroplast genome, we unravelled the phylogeography of the various taxonomic units of the A. lyrata complex. We demonstrate the existence of two major gene pools in Central Europe and Northern America. These two major gene pools are constructed from different taxonomic units. We also confirmed that A. kamchatica is the allotetraploid hybrid between A. lyrata and A. halleri, occupying the amphi-Beringian area in Eastern Asia and Northern America. This species closes the large distribution gap of the various other A. lyrata segregates. Furthermore, we revealed a threefold independent allopolyploid origin of this hybrid species in Japan, China, and Kamchatka.

CONCLUSIONS

Unglaciated parts of the Eastern Austrian Alps and arctic Eurasia, including Beringia, served as major glacial refugia of the Eurasian A. lyrata lineage, whereas A. halleri and its various subspecies probably survived in refuges in Central Europe and Eastern Asia with a large distribution gap in between. The North American A. lyrata lineage probably survived the glaciation in the southeast of North America. The dramatic climatic changes during glaciation and deglaciation cycles promoted not only secondary contact and formation of the allopolyploid hybrid A. kamchatica, but also provided the environment that allowed this species to fill a large geographic gap separating the two genetically different A. lyrata lineages from Eurasia and North America. With our example focusing on the evolutionary history of the A. lyrata species complex, we add substantial information to a broad evolutionary framework for future investigations within this emerging model system in molecular and evolutionary biology.

Extrachromosomal circles of satellite repeats and 5S ribosomal DNA in human cells

BACKGROUND

Extrachomosomal circular DNA (eccDNA) is ubiquitous in eukaryotic organisms and was detected in every organism tested, including in humans. A two-dimensional gel electrophoresis facilitates the detection of eccDNA in preparations of genomic DNA. Using this technique we have previously demonstrated that most of eccDNA consists of exact multiples of chromosomal tandemly repeated DNA, including both coding genes and satellite DNA.

RESULTS

Here we report the occurrence of eccDNA in every tested human cell line. It has heterogeneous mass ranging from less than 2 kb to over 20 kb. We describe eccDNA homologous to human alpha satellite and the SstI mega satellite. Moreover, we show, for the first time, circular multimers of the human 5S ribosomal DNA (rDNA), similar to previous findings in Drosophila and plants. We further demonstrate structures that correspond to intermediates of rolling circle replication, which emerge from the circular multimers of 5S rDNA and SstI satellite.

CONCLUSIONS

These findings, and previous reports, support the general notion that every chromosomal tandem repeat is prone to generate eccDNA in eukryoric organisms including humans. They suggest the possible involvement of eccDNA in the length variability observed in arrays of tandem repeats. The implications of eccDNA on genome biology may include mechanisms of centromere evolution, concerted evolution and homogenization of tandem repeats and genomic plasticity.

The use and limits of ITS data in the analysis of intraspecific variation in Passiflora L. (Passifloraceae)

The discovery and characterization of informative intraspecific genetic markers is fundamental for evolutionary and conservation genetics studies. Here, we used nuclear ribosomal ITS sequences to access intraspecific genetic diversity in 23 species of the genus Passiflora L. Some degree of variation was detected in 21 of these. The Passiflora and Decaloba (DC.) Rchb. subgenera showed significant differences in the sizes of the two ITS regions and in GC content, which can be related to reproductive characteristics of species in these subgenera. Furthermore, clear geographical patterns in the spatial distribution of sequence types were identified in six species. The results indicate that ITS may be a useful tool for the evaluation of intraspecific genetic variation in Passiflora.