Intragenomic Variation Research Articles

Genetics: Feedforward loop for diversity.

DNA-sequence analysis suggests that genetic mutations arise at elevated rates in genomes harbouring high levels of heterozygosity — the state in which the two copies of a genetic region contain sequence differences. See Letter p.463 To discover why mutations rates vary within genomes, Laurence Hurst and colleagues examined intragenomic variation in mutation rate directly in Arabidopsis, rice and the honey bee using a parent–offspring sequencing strategy. They find that mutation rates are higher in heterozygotes and in proximity to crossover events. Mutations occur disproportionately more often in heterozygous than in homozygous domains and gene clusters under purifying selection (commonly homozygous) and under balancing selection (mainly heterozygous) have low and high mutation rates, respectively. The authors suggest that extremely weak selection on the mutation rate may therefore not be necessary to explain why mutational hot and cold spots might correspond to regions under positive/balancing and purifying selection, respectively.

A Male-Specific Genetic Map of the Microcrustacean Daphnia pulex Based on Single-Sperm Whole-Genome Sequencing.

Genetic linkage maps are critical for assembling draft genomes to a meaningful chromosome level and for deciphering the genomic underpinnings of biological traits. The estimates of recombination rates derived from genetic maps also play an important role in understanding multiple aspects of genomic evolution such as nucleotide substitution patterns and accumulation of deleterious mutations. In this study, we developed a high-throughput experimental approach that combines fluorescence-activated cell sorting, whole-genome amplification, and short-read sequencing to construct a genetic map using single-sperm cells. Furthermore, a computational algorithm was developed to analyze single-sperm whole-genome sequencing data for map construction. These methods allowed us to rapidly build a male-specific genetic map for the freshwater microcrustacean Daphnia pulex, which shows significant improvements compared to a previous map. With a total of mapped 1672 haplotype blocks and an average intermarker distance of 0.87 cM, this map spans a total genetic distance of 1451 Kosambi cM and comprises 90% of the resolved regions in the current Daphnia reference assembly. The map also reveals the mistaken mapping of seven scaffolds in the reference assembly onto chromosome II by a previous microsatellite map based on F2 crosses. Our approach can be easily applied to many other organisms and holds great promise for unveiling the intragenomic and intraspecific variation in the recombination rates.

New Ceratocystis species from Eucalyptus and Cunninghamia in South China.

During routine surveys for possible fungal pathogens in the rapidly expanding plantations of Eucalyptus and Cunninghamia lanceolata in China, numerous isolates of unknown species in the genus Ceratocystis (Microascales) were obtained from tree wounds. In this study we identified the Ceratocystis isolates from Eucalyptus and Cunninghamia in the GuangDong, GuangXi, FuJian and HaiNan Provinces of South China based on morphology and through comparisons of DNA sequence data for the ITS, partial β-tubulin and TEF-1α gene regions. Morphological and DNA sequence comparisons revealed two previously unknown species residing in the Indo-Pacific Clade. These are described here as Ceratocystis cercfabiensis sp. nov. and Ceratocystis collisensis sp. nov. Isolates of Ceratocystis cercfabiensis showed intragenomic variation in their ITS sequences and four strains were selected for cloning of the ITS gene region. Twelve ITS haplotypes were obtained from 17 clones selected for sequencing, differing in up to seven base positions and representing two separate phylogenetic groups. This is the first evidence of multiple ITS types in isolates of Ceratocystis residing in the Indo-Pacific Clade. Caution should thus be exercised when using the ITS gene region as a barcoding marker for Ceratocystis species in this clade. This study also represents the first record of a species of Ceratocystis from Cunninghamia.

Intrapatient diversity of Achromobacter spp. involved in chronic colonization of Cystic Fibrosis airways

Achromobacter spp. are increasingly identified in Cystic Fibrosis (CF) patients and their ability to persistently colonize the CF respiratory tract (CFRT) suggests that Achromobacter species possess adaptive characteristics. We studied genome dynamics in 118 isolates recovered from 13 patients with Achromobacter chronic colonization (5–26 isolates per patient recovered over 13–61months). Isolates were identified to species level by nrdA gene sequencing, subjected to Pulsed-Field Gel Electrophoresis (PFGE) and multiplex rep-PCR (MR-PCR), and rrs intragenomic diversity was studied by PCR-Temporal Temperature Gel Electrophoresis (TTGE). Intrapatient diversity was assessed: (i) from dynamics of XbaI and/or SpeI-based pulsotypes, (ii) from comparison of MR-PCR profiles, and (iii) by longitudinal analysis of rrs intragenomic diversity. Patients were chronically colonized by Achromobacter xylosoxidans (n=10), Achromobacter dolens (n=1) or Achromobacter insuavis (n=2). All strains displayed genomic diversification over time but A. insuavis showed higher pulsotype diversity compared to other species. Intragenomic rrs heterogeneity was found in strains from 6 of 13 patients and may be persistently observed. Achromobacter genome evolution observed during chronic colonization of the CFRT warrants further investigation of the adaptation features of the different species, as well as of the selective forces driving this adaptation in the CFRT.

Intra-genomic variation in symbiotic dinoflagellates: recent divergence or recombination between lineages?

BackgroundThe symbiosis between corals and the dinoflagellate alga Symbiodinium is essential for the development and survival of coral reefs. Yet this fragile association is highly vulnerable to environmental disturbance. A coral’s ability to tolerate temperature stress depends on the fitness of its resident symbionts, whose thermal optima vary extensively between lineages. However, the in hospite population genetic structure of Symbiodinium is poorly understood and mostly based on analysis of bulk DNA extracted from thousands to millions of cells. Using quantitative single-cell PCR, we enumerated DNA polymorphisms in the symbionts of the reef-building coral Pocillopora damicornis, and applied a model selection approach to explore the potential for recombination between coexisting Symbiodinium populations.ResultsTwo distinct Symbiodinium ITS2 sequences (denoted C100 and C109) were retrieved from all P. damicornis colonies analysed. However, the symbiont assemblage consisted of three distinct Symbiodinium populations: cells featuring pure arrays of ITS2 type C109, near-homogeneous cells of type C100 (with trace ITS2 copies of type C109), and those with co-dominant C100 and C109 ITS2 repeats. The symbiont consortia of some colonies consisted almost entirely of these putative C100 × C109 recombinants.ConclusionsOur results are consistent with the occurrence of sexual recombination between Symbiodinium types C100 and C109. While the multiple-copy nature of the ITS2 dictates that the observed pattern of intra-genomic co-dominance may be a result of incomplete concerted evolution of intra-genomic polymorphisms, this is a less likely explanation given the occurrence of homogeneous cells of the C109 type. Conclusive evidence for inter-lineage recombination and introgression in this genus will require either direct observational evidence or a single-cell genotyping approach targeting multiple, single-copy loci.Electronic supplementary materialThe online version of this article (doi:10.1186/s12862-015-0325-1) contains supplementary material, which is available to authorized users.

Non-random expression of ribosomal DNA units in a grasshopper showing high intragenomic variation for the ITS2 region

We analyse intragenomic variation of the ITS2 internal transcribed spacer of ribosomal DNA (rDNA) in the grasshopper Eyprepocnemis plorans, by means of tagged PCR 454 amplicon sequencing performed on both genomic DNA (gDNA) and RNA-derived complementary DNA (cDNA), using part of the ITS2 flanking coding regions (5.8S and 28S rDNA) as an internal control for sequencing errors. Six different ITS2 haplotypes (i.e. variants for at least one nucleotide in the complete ITS2 sequence) were found in a single population, one of them (Hap4) being specific to a supernumerary (B) chromosome. The analysis of both gDNA and cDNA from the same individuals provided an estimate of the expression efficiency of the different haplotypes. We found random expression (i.e. about similar recovery in gDNA and cDNA) for three haplotypes (Hap1, Hap2 and Hap5), but significant underexpression for three others (Hap3, Hap4 and Hap6). Hap4 was the most extremely underexpressed and, remarkably, it showed the lowest sequence conservation for the flanking 5.8-28S coding regions in the gDNA reads but the highest conservation (100%) in the cDNA ones, suggesting the preferential expression of mutation-free rDNA units carrying this ITS2 haplotype. These results indicate that the ITS2 region of rDNA is far from complete homogenization in this species, and that the different rDNA units are not expressed at random, with some of them being severely downregulated.

Causes and consequences of crossing-over evidenced via a high-resolution recombinational landscape of the honey bee.

BackgroundSocial hymenoptera, the honey bee (Apis mellifera) in particular, have ultra-high crossover rates and a large degree of intra-genomic variation in crossover rates. Aligned with haploid genomics of males, this makes them a potential model for examining the causes and consequences of crossing over. To address why social insects have such high crossing-over rates and the consequences of this, we constructed a high-resolution recombination atlas by sequencing 55 individuals from three colonies with an average marker density of 314 bp/marker.ResultsWe find crossing over to be especially high in proximity to genes upregulated in worker brains, but see no evidence for a coupling with immune-related functioning. We detect only a low rate of non-crossover gene conversion, contrary to current evidence. This is in striking contrast to the ultrahigh crossing-over rate, almost double that previously estimated from lower resolution data. We robustly recover the predicted intragenomic correlations between crossing over and both population level diversity and GC content, which could be best explained as indirect and direct consequences of crossing over, respectively.ConclusionsOur data are consistent with the view that diversification of worker behavior, but not immune function, is a driver of the high crossing-over rate in bees. While we see both high diversity and high GC content associated with high crossing-over rates, our estimate of the low non-crossover rate demonstrates that high non-crossover rates are not a necessary consequence of high recombination rates.Electronic supplementary materialThe online version of this article (doi:10.1186/s13059-014-0566-0) contains supplementary material, which is available to authorized users.

Characterization of retrotransposon Bov-B LINE reverse transcriptase gene sequences in parthenogenetic lizards Darevskia unisexualis and bisexual species D. nairensis and D. valentini

Cloning and sequencing of the partial reverse transcriptase gene (750 bp) of the Bov-B LINE retrotransposon have been held in parthenogenetic lizards Darevskia unisexualis and its assumed parental bisexual species D. nairensis and D.valentini. It was shown that the percentage of transcriptionally active copies of this gene, which does not contain a stop codon, is almost the same in the three species and is about 75%. The intragenomic divergence level of these sequences is low and was found to be 2.6% in D. unisexualis, 1.9% in D. nairensis, and 1.6% in D. valentini. The phylogenetic analysis shows the distribution of copies of D. unisexualis in each of the two clusters of RT sequences characteristic of D. nairensis and D. valentini. This result supports the view of the hybrid origin of D. unisexualis and does not exclude intraspecific hybridization between D. nairensis and D. valentini.

Assessing Symbiodinium diversity in scleractinian corals via next-generation sequencing-based genotyping of the ITS2 rDNA region.

The persistence of coral reef ecosystems relies on the symbiotic relationship between scleractinian corals and intracellular, photosynthetic dinoflagellates in the genus Symbiodinium. Genetic evidence indicates that these symbionts are biologically diverse and exhibit discrete patterns of environmental and host distribution. This makes the assessment of Symbiodinium diversity critical to understanding the symbiosis ecology of corals. Here, we applied pyrosequencing to the elucidation of Symbiodinium diversity via analysis of the internal transcribed spacer 2 (ITS2) region, a multicopy genetic marker commonly used to analyse Symbiodinium diversity. Replicated data generated from isoclonal Symbiodinium cultures showed that all genomes contained numerous, yet mostly rare, ITS2 sequence variants. Pyrosequencing data were consistent with more traditional denaturing gradient gel electrophoresis (DGGE) approaches to the screening of ITS2 PCR amplifications, where the most common sequences appeared as the most intense bands. Further, we developed an operational taxonomic unit (OTU)-based pipeline for Symbiodinium ITS2 diversity typing to provisionally resolve ecologically discrete entities from intragenomic variation. A genetic distance cut-off of 0.03 collapsed intragenomic ITS2 variants of isoclonal cultures into single OTUs. When applied to the analysis of field-collected coral samples, our analyses confirm that much of the commonly observed SymbiodiniumITS2 diversity can be attributed to intragenomic variation. We conclude that by analysing Symbiodinium populations in an OTU-based framework, we can improve objectivity, comparability and simplicity when assessing ITS2 diversity in field-based studies.

Wide Occurrence of SSU rDNA Intragenomic Polymorphism in Foraminifera and its Implications for Molecular Species Identification

Ribosomal DNA is commonly used as a marker for protist phylogeny and taxonomy because of its ubiquity and its expected species specificity thanks to the mechanism of concerted evolution. However, numerous studies reported the occurrence of intragenomic (intra-individual) polymorphism in various protists and particularly in Foraminifera. To infer to what extent the SSU rDNA intragenomic variability occurs in Foraminifera, we studied 16 foraminiferal species belonging to single-chambered monothalamids and multi-chambered Globothalamea, with one to six individuals per species. We performed single-cell DNA extractions and PCRs of a 600bp fragment of SSU rDNA, and sequenced 9 to 23 clones per individual for a total of 818 sequences. We found intragenomic variability in almost all species, even after excluding singleton mutations. Intra-individual sequence divergence ranged from 0 to 5.15% and was higher than 1% in 11 species. Variability was usually located at the end of stem-loop structures and included compensatory single nucleotide polymorphisms and expansion segments polymorphisms. However, the polymorphisms did not change the secondary structure of the rRNA. Our results suggest a non-concerted evolution of rRNA genes in Foraminifera. The origin of this variability and its implications for species identification in environmental DNA studies are discussed.

Cryptic species revealed by molecular phylogenetic analysis of sequences obtained from basidiomata of Tulasnella

Delimitation of species and the search for a proper threshold for defining phylogenetic species in fungi are under discussion. In this study, morphological and molecular data are correlated to delimit species of Tulasnella, the most important mycobionts of Orchidaceae, which suffer from poor taxonomy. Resupinate basidiomata of Tulasnella species were collected in Ecuador and Germany, and 11 specimens (seven from Ecuador, four from Germany) were assigned to traditional species concepts by use of morphological keys. The specimens were compared by micro-anatomical examination with 75 specimens of Tulasnella borrowed from fungaria to obtain better insights on variation of characters. Sequences of the ITS region (127) were obtained after cloning from the fresh basidiomata and from pure cultures. Proportional variability of ITS sequences was analyzed within and among the cultures and the specimens designated to different morphospecies. Results suggested an intragenomic variation of less than 2%, an intraspecific variation of up to 4% and an interspecific divergence of more than 9% in Tulasnella. Cryptic species in Tulasnella, mostly from Ecuador, were revealed by phylogenetic analyses with 4% intraspecific divergence as a minimum threshold for delimiting species. Conventional diagnostic morphological characters appeared insufficient for species characterization. Arguments are presented for molecular delimitation of the established species Tulasnella albida, T. asymmetrica, T. eichleriana, T. cf. pinicola, T. tomaculum and T. violea.

Ribosomal DNA Sequence Heterogeneity Reflects Intraspecies Phylogenies and Predicts Genome Structure in Two Contrasting Yeast Species

The ribosomal RNA encapsulates a wealth of evolutionary information, including genetic variation that can be used to discriminate between organisms at a wide range of taxonomic levels. For example, the prokaryotic 16S rDNA sequence is very widely used both in phylogenetic studies and as a marker in metagenomic surveys and the internal transcribed spacer region, frequently used in plant phylogenetics, is now recognized as a fungal DNA barcode. However, this widespread use does not escape criticism, principally due to issues such as difficulties in classification of paralogous versus orthologous rDNA units and intragenomic variation, both of which may be significant barriers to accurate phylogenetic inference. We recently analyzed data sets from the Saccharomyces Genome Resequencing Project, characterizing rDNA sequence variation within multiple strains of the baker's yeast Saccharomyces cerevisiae and its nearest wild relative Saccharomyces paradoxus in unprecedented detail. Notably, both species possess single locus rDNA systems. Here, we use these new variation datasets to assess whether a more detailed characterization of the rDNA locus can alleviate the second of these phylogenetic issues, sequence heterogeneity, while controlling for the first. We demonstrate that a strong phylogenetic signal exists within both datasets and illustrate how they can be used, with existing methodology, to estimate intraspecies phylogenies of yeast strains consistent with those derived from whole-genome approaches. We also describe the use of partial Single Nucleotide Polymorphisms, a type of sequence variation found only in repetitive genomic regions, in identifying key evolutionary features such as genome hybridization events and show their consistency with whole-genome Structure analyses. We conclude that our approach can transform rDNA sequence heterogeneity from a problem to a useful source of evolutionary information, enabling the estimation of highly accurate phylogenies of closely related organisms, and discuss how it could be extended to future studies of multilocus rDNA systems. [concerted evolution; genome hydridisation; phylogenetic analysis; ribosomal DNA; whole genome sequencing; yeast]

Intraspecific and intragenomic variability of ITS rDNA sequences reveals taxonomic problems in Ceratocystis fimbriata sensu stricto.

Fourteen new species in the Latin American Clade (LAC) of the Ceratocystis fimbriata complex recently were distinguished from C. fimbriata sensu stricto largely based on variation in ITS rDNA sequences. Among the 116 isolates representing the LAC, there were 41 ITS haplotypes. Maximum parsimony (MP) analysis of ITS sequences produced poorly resolved trees. In contrast, analyses of mating-type genes (MAT1-1-2 and MAT1-2-1) resolved a single MP tree with branches of high bootstrap and posterior probability support. Four isolates showed intragenomic variation in ITS sequences. Cloning and sequencing of PCR products from a single haploid strain identified two or more ITS sequences differing at up to 16 base positions and representing two described species. Isolates from introduced populations that appeared to be clonal based on microsatellite markers varied at up to 14 bp in ITS sequence. Strains of seven Brazilian ITS haplotypes and an isolate from Ipomoea batatas (on which the species name C. fimbriata was based) were fully interfertile in sexual crosses. These analyses support three phylogenetic species that differ in pathogenicity: C. platani, C. cacaofunesta and C. colombiana. Five ITS species (C. manginecans, C. mangicola, C. mangivora, C. acaciivora, C. eucalypticola) appear to be ITS haplotypes that have been moved from or within Brazil on nursery stock. The taxonomic status of other species delineated primarily by ITS sequences (C. diversiconidia, C. papillata, C. neglecta, C. ecuadoriana, C. fimbriatomima, C. curvata) needs further study, but they are considered doubtful species.

Intraspecific and intragenomic variability of ITS rDNA sequences reveals taxonomic problems in Ceratocystis fimbriata sensu stricto

Fourteen new species in the Latin American Clade (LAC) of the Ceratocystis fimbriata complex recently were distinguished from C. fimbriata sensu stricto largely based on variation in ITS rDNA sequences. Among the 116 isolates representing the LAC, there were 41 ITS haplotypes. Maximum parsimony (MP) analysis of ITS sequences produced poorly resolved trees. In contrast, analyses of mating-type genes (MAT1-1-2 and MAT1-2-1) resolved a single MP tree with branches of high bootstrap and posterior probability support. Four isolates showed intragenomic variation in ITS sequences. Cloning and sequencing of PCR products from a single haploid strain identified two or more ITS sequences differing at up to 16 base positions and representing two described species. Isolates from introduced populations that appeared to be clonal based on microsatellite markers varied at up to 14 bp in ITS sequence. Strains of seven Brazilian ITS haplotypes and an isolate from Ipomoea batatas (on which the species name C. fimbriata was based) were fully interfertile in sexual crosses. These analyses support three phylogenetic species that differ in pathogenicity: C. platani, C. cacaofunesta and C. colombiana. Five ITS species (C. manginecans, C. mangicola, C. mangivora, C. acaciivora, C. eucalypticola) appear to be ITS haplotypes that have been moved from or within Brazil on nursery stock. The taxonomic status of other species delineated primarily by ITS sequences (C. diversiconidia, C. papillata, C. neglecta, C. ecuadoriana, C. fimbriatomima, C. curvata) needs further study, but they are considered doubtful species.

Extreme phenotypic polymorphism in the coral genus <I>Pocillopora;</I> micro-morphology corresponds to mitochondrial groups, while colony morphology does not

Coral species are difficult to identify due to few distinct morphological features and high levels of phenotypic variation; furthermore, hybridization is thought to play an important role in coral evolution. The genus Pocillopora is one of the most abundant reef-building corals in the Pacific, yet species boundaries within this important group remain enigmatic, with conflicting interpretations of molecular data. To further examine the relationship between morphology and genetic variation, we examined in situ photographs, microscopic images, and molecular markers (ORF, CR, ITS-2) from 59 Pocillopora colonies representing a range of atypical morphologies thought to be rare or endemic to remote locations (the Hawaiian Archipelago, Johnston Atoll, Palmyra Atoll, and Howland Island). Additionally, we combined our data with previously published genetic data from the western Indian Ocean, Taiwan, the Great Barrier Reef, and the tropical east Pacific (n = 298) to examine phylogeographic relationships more broadly. Mitochondrial markers revealed six distinct clades, while nuclear data (ITS-2) resolved no clear groups due in part to intragenomic variation. Although colony-level morphology is highly variable within each mitochondrial clade, microscopic corallite-level morphology showed less variation and stronger concordance with mitochondrial groups, and one clade was found only in the Hawaiian Archipelago. These findings will assist with efforts to delineate species boundaries in the genus Pocillopora and contribute to understanding the scale of species-level variation in reef building corals. Accurate determination of species boundaries will improve conservation efforts and further our understanding of the evolution and biodiversity of reef building corals.

The utility of PacBio circular consensus sequencing for characterizing complex gene families in non-model organisms.

BackgroundMolecular characterization of highly diverse gene families can be time consuming, expensive, and difficult, especially when considering the potential for relatively large numbers of paralogs and/or pseudogenes. Here we investigate the utility of Pacific Biosciences single molecule real-time (SMRT) circular consensus sequencing (CCS) as an alternative to traditional cloning and Sanger sequencing PCR amplicons for gene family characterization. We target vomeronasal gene receptors, one of the most diverse gene families in mammals, with the goal of better understanding intra-specific V1R diversity of the gray mouse lemur (Microcebus murinus). Our study compares intragenomic variation for two V1R subfamilies found in the mouse lemur. Specifically, we compare gene copy variation within and between two individuals of M. murinus as characterized by different methods for nucleotide sequencing. By including the same individual animal from which the M. murinus draft genome was derived, we are able to cross-validate gene copy estimates from Sanger sequencing versus CCS methods.ResultsWe generated 34,088 high quality circular consensus sequences of two diverse V1R subfamilies (here referred to as V1RI and V1RIX) from two individuals of Microcebus murinus. Using a minimum threshold of 7× coverage, we recovered approximately 90% of V1RI sequences previously identified in the draft M. murinus genome (59% being identical at all nucleotide positions). When low coverage sequences were considered (i.e. < 7× coverage) 100% of V1RI sequences identified in the draft genome were recovered. At least 13 putatively novel V1R loci were also identified using CCS technology.ConclusionsRecent upgrades to the Pacific Biosciences RS instrument have improved the CCS technology and offer an alternative to traditional sequencing approaches. Our results suggest that the Microcebus murinus V1R repertoire has been underestimated in the draft genome. In addition to providing an improved understanding of V1R diversity in the mouse lemur, this study demonstrates the utility of CCS technology for characterizing complex regions of the genome. We anticipate that long-read sequencing technologies such as PacBio SMRT will allow for the assembly of multigene family clusters and serve to more accurately characterize patterns of gene copy variation in large gene families, thus revealing novel micro-evolutionary patterns within non-model organisms.Electronic supplementary materialThe online version of this article (doi:10.1186/1471-2164-15-720) contains supplementary material, which is available to authorized users.

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.

Multiple ITS Haplotypes in the Genome of the Lichenized Basidiomycete Cora inversa (Hygrophoraceae): Fact or Artifact?

The internal transcribed spacer region (ITS) of the nuclear rDNA cistron represents the barcoding locus for Fungi. Intragenomic variation of this multicopy gene can interfere with accurate phylogenetic reconstruction of biological entities. We investigated the amount and nature of this variation for the lichenized fungus Cora inversa in the Hygrophoraceae (Basidiomycota: Agaricales), analyzing base call and length variation in ITS1 454 pyrosequencing data of three samples of the target mycobiont, for a total of 16,665 reads obtained from three separate repeats of the same samples under different conditions. Using multiple fixed alignment methods (PaPaRa) and maximum likelihood phylogenetic analysis (RAxML), we assessed phylogenetic relationships of the obtained reads, together with Sanger ITS sequences from the same samples. Phylogenetic analysis showed that all ITS1 reads belonged to a single species, C. inversa. Pyrosequencing data showed 266 insertion sites in addition to the 325 sites expected from Sanger sequences, for a total of 15,654 insertions (0.94 insertions per read). An additional 3,279 substitutions relative to the Sanger sequences were detected in the dataset, out of 5,461,125 bases to be called. Up to 99.3% of the observed indels in the dataset could be interpreted as 454 pyrosequencing errors, approximately 65% corresponding to incorrectly recovered homopolymer segments, and 35% to carry-forward-incomplete-extension errors. Comparison of automated clustering and alignment-based phylogenetic analysis demonstrated that clustering of these reads produced a 35-fold overestimation of biological diversity in the dataset at the 95% similarity threshold level, whereas phylogenetic analysis using a maximum likelihood approach accurately recovered a single biological entity. We conclude that variation detected in 454 pyrosequencing data must be interpreted with great care and that a combination of a sufficiently large number of reads per taxon, a set of Sanger references for the same taxon, and at least two runs under different emulsion PCR and sequencing conditions, are necessary to reliably separate biological variation from 454 sequencing errors. Our study shows that clustering methods are highly sensitive to artifactual sequence variation and inadequate to properly recover biological diversity in a dataset, if sequencing errors are substantial and not removed prior to clustering analysis.

Intra-genomic variation in the ribosomal repeats of nematodes.

Ribosomal loci represent a major tool for investigating environmental diversity and community structure via high-throughput marker gene studies of eukaryotes (e.g. 18S rRNA). Since the estimation of species’ abundance is a major goal of environmental studies (by counting numbers of sequences), understanding the patterns of rRNA copy number across species will be critical for informing such high-throughput approaches. Such knowledge is critical, given that ribosomal RNA genes exist within multi-copy repeated arrays in a genome. Here we measured the repeat copy number for six nematode species by mapping the sequences from whole genome shotgun libraries against reference sequences for their rRNA repeat. This revealed a 6-fold variation in repeat copy number amongst taxa investigated, with levels of intragenomic variation ranging from 56 to 323 copies of the rRNA array. By applying the same approach to four C. elegans mutation accumulation lines propagated by repeated bottlenecking for an average of ~400 generations, we find on average a 2-fold increase in repeat copy number (rate of increase in rRNA estimated at 0.0285-0.3414 copies per generation), suggesting that rRNA repeat copy number is subject to selection. Within each Caenorhabditis species, the majority of intragenomic variation found across the rRNA repeat was observed within gene regions (18S, 28S, 5.8S), suggesting that such intragenomic variation is not a product of selection for rRNA coding function. We find that the dramatic variation in repeat copy number among these six nematode genomes would limit the use of rRNA in estimates of organismal abundance. In addition, the unique pattern of variation within a single genome was uncorrelated with patterns of divergence between species, reflecting a strong signature of natural selection for rRNA function. A better understanding of the factors that control or affect copy number in these arrays, as well as their rates and patterns of evolution, will be critical for informing estimates of global biodiversity.

Extreme diversification of the mating type–high‐mobility group (MATA‐HMG) gene family in a plant‐associated arbuscular mycorrhizal fungus

Arbuscular mycorrhizal fungi (AMF) are important plant symbionts that have long been considered evolutionary anomalies because of their apparent long-term lack of sexuality, but recent explorations of available DNA sequence have challenged this notion by revealing the presence of homologues of fungal mating type-high-mobility group (MATA-HMG) and core meiotic genes in these organisms. To obtain more insights into the sexual potential of AMF, homologues of MATA-HMGs were sought in the transcriptome of three AMF isolates, and their functional and evolutionary trajectories were studied in genetically divergent strains of Rhizophagus irregularis using conventional and quantitative PCR procedures. Our analyses revealed the presence of at least 76 homologues of MATA-HMGs in R.irregularis isolates. None of these was found to be surrounded by genes generally found near other known fungal mating type loci, but here we report the presence of a 9-kb-long region in the AMF R.irregularis harbouring a total of four tandem-repeated MATA-HMGs; a feature that highlights a potentially elevated intragenomic diversity in this AMF species. The present study provides intriguing insights into the genome evolution of R.irregularis, and represents a stepping stone for understanding the potential of these fungi to undergo cryptic sex.