Disentangling structural genomic and behavioural barriers in a sea of connectivity

Genetic divergence among populations arises through natural selection or drift and is counteracted by connectivity and gene flow. In sympatric populations, isolating mechanisms are thus needed to limit the homogenizing effects of gene flow to allow for adaptation and speciation. Chromosomal inversions act as an important mechanism maintaining isolating barriers, yet their role in sympatric populations and divergence with gene flow is not entirely understood. Here, we revisit the question of whether inversions play a role in the divergence of connected populations of the marine fish Atlantic cod (Gadus morhua), by exploring a unique data set combining whole‐genome sequencing data and behavioural data obtained with acoustic telemetry. Within a confined fjord environment, we find three genetically differentiated Atlantic cod types belonging to the oceanic North Sea population, the western Baltic population and a local fjord‐type cod. Continuous behavioural tracking over 4 year revealed temporally stable sympatry of these types within the fjord. Despite overall weak genetic differentiation consistent with high levels of gene flow, we detected significant frequency shifts of three previously identified inversions, indicating an adaptive barrier to gene flow. In addition, behavioural data indicated that North Sea cod and individuals homozygous for the LG12 inversion had lower fitness in the fjord environment. However, North Sea and fjord‐type cod also occupy different depths, possibly contributing to prezygotic reproductive isolation and representing a behavioural barrier to gene flow. Our results provide the first insights into a complex interplay of genomic and behavioural isolating barriers in Atlantic cod and establish a new model system towards an understanding of the role of genomic structural variants in adaptation and diversification.     

Genomic divergence between the migratory and stationary ecotypes of Atlantic cod

Atlantic cod displays a range of phenotypic and genotypic variations, which includes the differentiation into coastal stationary and offshore migratory types of cod that co‐occur in several parts of its distribution range and are often sympatric on the spawning grounds. Differentiation of these ecotypes may involve both historical separation and adaptation to ecologically distinct environments, the genetic basis of which is now beginning to be unravelled. Genomic analyses based on recent sequencing advances are able to document genomic divergence in more detail and may facilitate the exploration of causes and consequences of genome‐wide patterns. We examined genomic divergence between the stationary and migratory types of cod in the Northeast Atlantic, using next‐generation sequencing of pooled DNA from each of two population samples. Sequence data was mapped to the published cod genome sequence, arranged in more than 6000 scaffolds (611 Mb). We identified 25 divergent scaffolds (26 Mb) with a higher than average gene density, against a backdrop of overall moderate genomic differentiation. Previous findings of localized genomic divergence in three linkage groups were confirmed, including a large (15 Mb) genomic region, which seems to be uniquely involved in the divergence of migratory and stationary cod. The results of the pooled sequencing approach support and extend recent findings based on single‐nucleotide polymorphism markers and suggest a high degree of reproductive isolation between stationary and migratory cod in the North‐east Atlantic.     

Genomics of speciation and introgression in Princess cichlid fishes from Lake Tanganyika

How variation in the genome translates into biological diversity and new species originate has endured as the mystery of mysteries in evolutionary biology. African cichlid fishes are prime model systems to address speciation‐related questions for their remarkable taxonomic and phenotypic diversity, and the possible role of gene flow in this process. Here, we capitalize on genome sequencing and phylogenomic analyses to address the relative impacts of incomplete lineage sorting, introgression and hybrid speciation in the Neolamprologus savoryi‐complex (the ‘Princess cichlids’) from Lake Tanganyika. We present a time‐calibrated species tree based on whole‐genome sequences and provide strong evidence for incomplete lineage sorting in the early phases of diversification and multiple introgression events affecting different stages. Importantly, we find that the Neolamprologus chromosomes show centre‐to‐periphery biases in nucleotide diversity, sequence divergence, GC content, incomplete lineage sorting and rates of introgression, which are likely modulated by recombination density and linked selection. The detection of heterogeneous genomic landscapes has strong implications on the genomic mechanisms involved in speciation. Collinear chromosomal regions can be protected from gene flow and harbour incompatibility genes if they reside in lowly recombining regions, and coupling can evolve between nonphysically linked genomic regions (chromosome centres in particular). Simultaneously, higher recombination towards chromosome peripheries makes these more dynamic, evolvable regions where adaptation polymorphisms have a fertile ground. Hence, differences in genome architecture could explain the levels of taxonomic and phenotypic diversity seen in taxa with collinear genomes and might have contributed to the spectacular cichlid diversity observed today.     

Improving Illumina assemblies with Hi‐C and long reads: An example with the North African dromedary

Researchers have assembled thousands of eukaryotic genomes using Illumina reads, but traditional mate‐pair libraries cannot span all repetitive elements, resulting in highly fragmented assemblies. However, both chromosome conformation capture techniques, such as Hi‐C and Dovetail Genomics Chicago libraries and long‐read sequencing, such as Pacific Biosciences and Oxford Nanopore, help span and resolve repetitive regions and therefore improve genome assemblies. One important livestock species of arid regions that does not have a high‐quality contiguous reference genome is the dromedary (Camelus dromedarius). Draft genomes exist but are highly fragmented, and a high‐quality reference genome is needed to understand adaptation to desert environments and artificial selection during domestication. Dromedaries are among the last livestock species to have been domesticated, and together with wild and domestic Bactrian camels, they are the only representatives of the Camelini tribe, which highlights their evolutionary significance. Here we describe our efforts to improve the North African dromedary genome. We used Chicago and Hi‐C sequencing libraries from Dovetail Genomics to resolve the order of previously assembled contigs, producing almost chromosome‐level scaffolds. Remaining gaps were filled with Pacific Biosciences long reads, and then scaffolds were comparatively mapped to chromosomes. Long reads added 99.32 Mbp to the total length of the new assembly. Dovetail Chicago and Hi‐C libraries increased the longest scaffold over 12‐fold, from 9.71 Mbp to 124.99 Mbp and the scaffold N50 over 50‐fold, from 1.48 Mbp to 75.02 Mbp. We demonstrate that Illumina de novo assemblies can be substantially upgraded by combining chromosome conformation capture and long‐read sequencing.     

Targeted multiplex next‐generation sequencing: advances in techniques of mitochondrial and nuclear DNA sequencing for population genomics

Next‐generation sequencing (NGS) is emerging as an efficient and cost‐effective tool in population genomic analyses of nonmodel organisms, allowing simultaneous resequencing of many regions of multi‐genomic DNA from multiplexed samples. Here, we detail our synthesis of protocols for targeted resequencing of mitochondrial and nuclear loci by generating indexed genomic libraries for multiplexing up to 100 individuals in a single sequencing pool, and then enriching the pooled library using custom DNA capture arrays. Our use of DNA sequence from one species to capture and enrich the sequencing libraries of another species (i.e. cross‐species DNA capture) indicates that efficient enrichment occurs when sequences are up to about 12% divergent, allowing us to take advantage of genomic information in one species to sequence orthologous regions in related species. In addition to a complete mitochondrial genome on each array, we have included between 43 and 118 nuclear loci for low‐coverage sequencing of between 18 kb and 87 kb of DNA sequence per individual for single nucleotide polymorphisms discovery from 50 to 100 individuals in a single sequencing lane. Using this method, we have generated a total of over 500 whole mitochondrial genomes from seven cetacean species and green sea turtles. The greater variation detected in mitogenomes relative to short mtDNA sequences is helping to resolve genetic structure ranging from geographic to species‐level differences. These NGS and analysis techniques have allowed for simultaneous population genomic studies of mtDNA and nDNA with greater genomic coverage and phylogeographic resolution than has previously been possible in marine mammals and turtles.     

High‐contiguity genome assembly of the chemosynthetic gammaproteobacterial endosymbiont of the cold seep tubeworm Lamellibrachia barhami

Symbiotic relationships between vestimentiferan tubeworms and chemosynthetic Gammaproteobacteria build the foundations of many hydrothermal vent and hydrocarbon seep ecosystems in the deep sea. The association between the vent tubeworm Riftia pachyptila and its endosymbiont Candidatus Endoriftia persephone has become a model system for symbiosis research in deep‐sea vestimentiferans, while markedly fewer studies have investigated symbiotic relationships in other tubeworm species, especially at cold seeps. Here we sequenced the endosymbiont genome of the tubeworm Lamellibrachia barhami from a cold seep in the Gulf of California, using short‐ and long‐read sequencing technologies in combination with Hi‐C and Dovetail Chicago libraries. Our final assembly had a size of ~4.17 MB, a GC content of 54.54%, 137X coverage, 4153 coding sequences, and a CheckM completeness score of 97.19%. A single scaffold contained 99.51% of the genome. Comparative genomic analyses indicated that the L. barhami symbiont shares a set of core genes and many metabolic pathways with other vestimentiferan symbionts, while containing 433 unique gene clusters that comprised a variety of transposases, defence‐related genes and a lineage‐specific CRISPR/Cas3 system. This assembly represents the most contiguous tubeworm symbiont genome resource to date and will be particularly valuable for future comparative genomic studies investigating structural genome evolution, physiological adaptations and host‐symbiont communication in chemosynthetic animal‐microbe symbioses.     

Reverse genetics in complex multigene operons by co‐transformation of the plastid genome and its application to the open reading frame previously designated psbN

Reverse genetics approaches have contributed enormously to the elucidation of gene functions in plastid genomes and the determination of structure–function relationships in chloroplast multiprotein complexes. Gene knock‐outs are usually performed by disrupting the reading frame of interest with a selectable marker cassette. Site‐directed mutagenesis is done by placing the marker into the adjacent intergenic spacer and relying on co‐integration of the desired mutation by homologous recombination. These strategies are not applicable to genes residing in large multigene operons or other gene‐dense genomic regions, because insertion of the marker cassette into an operon‐internal gene or into the nearest intergenic spacer is likely to interfere with expression of adjacent genes in the operon or disrupt cis‐elements for the expression of neighboring genes and operons. Here we have explored the possibility of using a co‐transformation strategy to mutate a small gene of unknown function (psbN) that is embedded in a complex multigene operon. Although inactivation of psbN resulted in strong impairment of photosynthesis, homoplasmic knock‐out lines were readily recovered by co‐transformation with a selectable marker integrating >38 kb away from the targeted psbN. Our results suggest co‐transformation as a suitable strategy for the functional analysis of plastid genes and operons, which allows the recovery of unselected homoplasmic mutants even if the introduced mutations entail a significant selective disadvantage. Moreover, our data provide evidence for involvement of the psbN gene product in the biogenesis of both photosystem I and photosystem II. We therefore propose to rename the gene product ‘photosystem biogenesis factor 1′ and the gene pbf1.     

Karyotyping and cytogenetic mapping of Atlantic cod (Gadus morhua Linnaeus, 1758)

The Atlantic cod (Gadus morhua) is an important natural resource for northern societies and is now also considered to be a promising candidate for aquaculture. In recent years, much effort has been directed towards the development of genomic tools, and genome initiatives for Atlantic cod have been established. Despite the growing attention devoted to the Atlantic cod genomics, basic aspects of its genome structure and organization remain unknown. Thus, the present work aims to study cytogenetic features of the Atlantic cod as a contribution to the knowledge of this species' genome. The Atlantic cod displays a diploid number of 46 chromosomes, with a karyotypic formula 16 m/sm + 30 st/t. Conventional karyotyping was improved by chromosomal mapping of two classes of repetitive sequences. 18S rDNA clusters were assigned to pairs 2 and 4; small amounts of 18S rDNA clusters were occasionally detected on pair 5. These findings could not be related to the geographical origin of the specimens, but were consistent with the variability of these repeated genes in fish in general. 5S ribosomal gene clusters, apparently corresponding to a single 5S rDNA class, were detected on twelve chromosomes (pairs 11, 12, 14, 17, 20 and 21). The present update of the existing but meagre information on the karyotype of Atlantic cod, plus the first physical mapping of repetitive genes in this species herein, opens the way for an integrated approach that combines genetic and physical mapping with the assembly of the genome of this commercially important species.     

Whole‐exome targeted sequencing of the uncharacterized pine genome

The large genome size of many species hinders the development and application of genomic tools to study them. For instance, loblolly pine (Pinus taeda L.), an ecologically and economically important conifer, has a large and yet uncharacterized genome of 21.7 Gbp. To characterize the pine genome, we performed exome capture and sequencing of 14 729 genes derived from an assembly of expressed sequence tags. Efficiency of sequence capture was evaluated and shown to be similar across samples with increasing levels of complexity, including haploid cDNA, haploid genomic DNA and diploid genomic DNA. However, this efficiency was severely reduced for probes that overlapped multiple exons, presumably because intron sequences hindered probe:exon hybridizations. Such regions could not be entirely avoided during probe design, because of the lack of a reference sequence. To improve the throughput and reduce the cost of sequence capture, a method to multiplex the analysis of up to eight samples was developed. Sequence data showed that multiplexed capture was reproducible among 24 haploid samples, and can be applied for high‐throughput analysis of targeted genes in large populations. Captured sequences were de novo assembled, resulting in 11 396 expanded and annotated gene models, significantly improving the knowledge about the pine gene space. Interspecific capture was also evaluated with over 98% of all probes designed from P. taeda that were efficient in sequence capture, were also suitable for analysis of the related species Pinus elliottii Engelm.     

Long-Range Genomic Enrichment,Sequencing, and Assembly to Determine Unknown Sequences Flanking a Known microRNA

Conserved plant microRNAs (miRNAs) modulate important biological processes but little is known about conserved cis-regulatory elements (CREs) surrounding MIRNA genes. We developed a solution-based targeted genomic enrichment methodology to capture, enrich, and sequence flanking genomic regions surrounding conserved MIRNA genes with a locked-nucleic acid (LNA)-modified, biotinylated probe complementary to the mature miRNA sequence. Genomic DNA bound by the probe is captured by streptavidin-coated magnetic beads, amplified, sequenced and assembled de novo to obtain genomic DNA sequences flanking MIRNA locus of interest. We demonstrate the sensitivity and specificity of this enrichment methodology in Arabidopsis thaliana to enrich targeted regions spanning 10–20 kb surrounding known MIR166 and MIR165 loci. Assembly of the sequencing reads successfully recovered all targeted loci. While further optimization for larger, more complex genomes is needed, this method may enable determination of flanking genomic DNA sequence surrounding a known core (like a conserved mature miRNA) from multiple species that currently don''t have a full genome assembly available.     

Targeted Sequencing of Large Genomic Regions with CATCH-Seq

Current target enrichment systems for large-scale next-generation sequencing typically require synthetic oligonucleotides used as capture reagents to isolate sequences of interest. The majority of target enrichment reagents are focused on gene coding regions or promoters en masse. Here we introduce development of a customizable targeted capture system using biotinylated RNA probe baits transcribed from sheared bacterial artificial chromosome clone templates that enables capture of large, contiguous blocks of the genome for sequencing applications. This clone adapted template capture hybridization sequencing (CATCH-Seq) procedure can be used to capture both coding and non-coding regions of a gene, and resolve the boundaries of copy number variations within a genomic target site. Furthermore, libraries constructed with methylated adapters prior to solution hybridization also enable targeted bisulfite sequencing. We applied CATCH-Seq to diverse targets ranging in size from 125 kb to 3.5 Mb. Our approach provides a simple and cost effective alternative to other capture platforms because of template-based, enzymatic probe synthesis and the lack of oligonucleotide design costs. Given its similarity in procedure, CATCH-Seq can also be performed in parallel with commercial systems.     

Transfer of the methicillin resistance genomic island among staphylococci by conjugation

Methicillin resistance creates a major obstacle for treatment of Staphylococcus aureus infections. The resistance gene, mecA, is carried on a large (20 kb to > 60 kb) genomic island, staphylococcal cassette chromosome mec (SCCmec), that excises from and inserts site‐specifically into the staphylococcal chromosome. However, although SCCmec has been designated a mobile genetic element, a mechanism for its transfer has not been defined. Here we demonstrate the capture and conjugative transfer of excised SCCmec. SCCmec was captured on pGO400, a mupirocin‐resistant derivative of the pGO1/pSK41 staphylococcal conjugative plasmid lineage, and pGO400::SCCmec (pRM27) was transferred by filter‐mating into both homologous and heterologous S. aureus recipients representing a range of clonal complexes as well as S. epidermidis. The DNA sequence of pRM27 showed that SCCmec had been transferred in its entirety and that its capture had occurred by recombination between IS257/431 elements present on all SCCmec types and pGO1/pSK41 conjugative plasmids. The captured SCCmec excised from the plasmid and inserted site‐specifically into the chromosomal att site of both an isogenic S. aureus and a S. epidermidis recipient. These studies describe a means by which methicillin resistance can be environmentally disseminated and a novel mechanism, IS‐mediated recombination, for the capture and conjugative transfer of genomic islands.     

Phylogenetic relationships among North Atlantic redfish (genus Sebastes) as revealed by mitochondrial DNA sequence analyses

Species structure and phylogenetic relationships among the four currently recognised North Atlantic redfish species, Sebastes mentella (‘deep‐sea’ and ‘oceanic’ phenotype), S. marinus (including ‘giant’ phenotype), S. fasciatus and S. viviparus, were investigated by sequencing the complete mitochondrial NADH dehydrogenase subunit 3 (ND3) gene for 313 individuals (333 including outliers). Low levels of sequence divergence of 0.3–2.3% between the North Atlantic Sebastes haplotypes, low resolution in phylogenetic trees and a star‐like statistical parsimony network indicate a population bottleneck and/or founder event with a subsequent population expansion and a rapid speciation in North Atlantic Sebastes. Diagnostic mutations within the ND3 gene were recognised for some of the species. The results lead to the hypothesis that the ancestral haplotype lineage is now very rare or has suffered extinction and that the most basal lineage is found in S. viviparus. The results further indicate direct descendants in the S. fasciatus and S. marinus lineages from the S. mentella lineage, whereas the S. viviparus lineage originates from an earlier division. The observation of two haplotype lineages in S. marinus with a relatively high level of intraspecific genetic divergence, points towards a possible cryptic speciation.     

BAC-based sequencing of behaviorally-relevant genes in the prairie vole

The prairie vole (Microtus ochrogaster) is an important model organism for the study of social behavior, yet our ability to correlate genes and behavior in this species has been limited due to a lack of genetic and genomic resources. Here we report the BAC-based targeted sequencing of behaviorally-relevant genes and flanking regions in the prairie vole. A total of 6.4 Mb of non-redundant or haplotype-specific sequence assemblies were generated that span the partial or complete sequence of 21 behaviorally-relevant genes as well as an additional 55 flanking genes. Estimates of nucleotide diversity from 13 loci based on alignments of 1.7 Mb of haplotype-specific assemblies revealed an average pair-wise heterozygosity (8.4×10(-3)). Comparative analyses of the prairie vole proteins encoded by the behaviorally-relevant genes identified >100 substitutions specific to the prairie vole lineage. Finally, our sequencing data indicate that a duplication of the prairie vole AVPR1A locus likely originated from a recent segmental duplication spanning a minimum of 105 kb. In summary, the results of our study provide the genomic resources necessary for the molecular and genetic characterization of a high-priority set of candidate genes for regulating social behavior in the prairie vole.     

Identifying the causes and consequences of assembly gaps using a multiplatform genome assembly of a bird‐of‐paradise

Genome assemblies are currently being produced at an impressive rate by consortia and individual laboratories. The low costs and increasing efficiency of sequencing technologies now enable assembling genomes at unprecedented quality and contiguity. However, the difficulty in assembling repeat‐rich and GC‐rich regions (genomic “dark matter”) limits insights into the evolution of genome structure and regulatory networks. Here, we compare the efficiency of currently available sequencing technologies (short/linked/long reads and proximity ligation maps) and combinations thereof in assembling genomic dark matter. By adopting different de novo assembly strategies, we compare individual draft assemblies to a curated multiplatform reference assembly and identify the genomic features that cause gaps within each assembly. We show that a multiplatform assembly implementing long‐read, linked‐read and proximity sequencing technologies performs best at recovering transposable elements, multicopy MHC genes, GC‐rich microchromosomes and the repeat‐rich W chromosome. Telomere‐to‐telomere assemblies are not a reality yet for most organisms, but by leveraging technology choice it is now possible to minimize genome assembly gaps for downstream analysis. We provide a roadmap to tailor sequencing projects for optimized completeness of both the coding and noncoding parts of nonmodel genomes.     

Low‐coverage whole‐genome sequencing reveals molecular markers for spawning season and sex identification in Gulf of Maine Atlantic cod (Gadus morhua,Linnaeus 1758)

Atlantic cod (Gadus morhua, Linnaeus 1758) in the western Gulf of Maine are managed as a single stock despite several lines of evidence supporting two spawning groups (spring and winter) that overlap spatially, while exhibiting seasonal spawning isolation. Low‐coverage whole‐genome sequencing was used to evaluate the genomic population structure of Atlantic cod spawning groups in the western Gulf of Maine and Georges Bank using 222 individuals collected over multiple years. Results indicated low total genomic differentiation, while also showing strong differentiation between spring and winter‐spawning groups at specific regions of the genome. Guided regularized random forest and ranked F _ST methods were used to select panels of single nucleotide polymorphisms (SNPs) that could reliably distinguish spring and winter‐spawning Atlantic cod (88.5% assignment rate), as well as males and females (95.0% assignment rate) collected in the western Gulf of Maine. These SNP panels represent a valuable tool for fisheries research and management of Atlantic cod in the western Gulf of Maine that will aid investigations of stock production and support accuracy of future assessments.     

Telomeric gene deletion and intrachromosomal amplification in antimony‐resistant Leishmania

Antimonials are still the mainstay of treatment against leishmaniasis but drug resistance is increasing. We carried out short read next‐generation sequencing (NGS) and comparative genomic hybridization (CGH) of three independent Leishmania major antimony‐resistant mutants. Copy number variations were consistently detected with both NGS and CGH. A major attribute of antimony resistance was a novel terminal deletion of variable length (67 kb to 204 kb) of the polyploid chromosome 31 in the three mutants. Terminal deletions in two mutants occurred at the level of inverted repeated sequences. The AQP1 gene coding for an aquaglyceroporin was part of the deleted region and its transfection into resistant mutants reverted resistance to SbIII. We also highlighted an intrachromosomal amplification of a subtelomeric locus on chromosome 34 in one mutant. This region encoded for ascorbate‐dependent peroxidase (APX) and glucose‐6‐phosphate dehydrogenase (G6PDH). Overexpression of these genes in revertant backgrounds demonstrated resistance to SbIII and protection from reactive oxygen species (ROS). Generation of a G6PDH null mutant in one revertant exhibited SbIII sensitivity and a decreased protection of ROS. Our genomic analyses and functional validation highlighted novel genomic rearrangements, functionally important resistant loci and the implication of new genes in antimony resistance in Leishmania.     

Characterizing restriction enzyme‐associated loci in historic ragweed (Ambrosia artemisiifolia) voucher specimens using custom‐designed RNA probes

Population genetic studies of nonmodel organisms frequently employ reduced representation library (RRL) methodologies, many of which rely on protocols in which genomic DNA is digested by one or more restriction enzymes. However, because high molecular weight DNA is recommended for these protocols, samples with degraded DNA are generally unsuitable for RRL methods. Given that ancient and historic specimens can provide key temporal perspectives to evolutionary questions, we explored how custom‐designed RNA probes could enrich for RRL loci (Restriction Enzyme‐Associated Loci baits, or REALbaits). Starting with genotyping‐by‐sequencing (GBS) data generated on modern common ragweed (Ambrosia artemisiifolia L.) specimens, we designed 20 000 RNA probes to target well‐characterized genomic loci in herbarium voucher specimens dating from 1835 to 1913. Compared to shotgun sequencing, we observed enrichment of the targeted loci at 19‐ to 151‐fold. Using our GBS capture pipeline on a data set of 38 herbarium samples, we discovered 22 813 SNPs, providing sufficient genomic resolution to distinguish geographic populations. For these samples, we found that dilution of REALbaits to 10% of their original concentration still yielded sufficient data for downstream analyses and that a sequencing depth of ~7m reads was sufficient to characterize most loci without wasting sequencing capacity. In addition, we observed that targeted loci had highly variable rates of success, which we primarily attribute to similarity between loci, a trait that ultimately interferes with unambiguous read mapping. Our findings can help researchers design capture experiments for RRL loci, thereby providing an efficient means to integrate samples with degraded DNA into existing RRL data sets.     

Molecular systematics of Brassica and allied genera in subtribes Brassicinae, Raphaninae, Moricandiinae, and Cakilinae (Brassicaceae, tribe Brassiceae); the organization and evolution of ribosomal gene families

DNA restriction endonuclease fragment analysis was used to obtain new information on the genomic organization of ribosomal DNA (rDNA) of Brassica and allied genera. The total genomic DNA of 95 accessions of 52 species representing 16 genera was restricted with six enzymes, and the restriction fragments were probed with three ribosomal clones (pTA71, Ver 18‐6, and Ver 6‐5). Eleven repeat unit length classes were recognized. The repeat unit size classes of 8.9 kb and 9.5 kb were observed most commonly, being represented in 17 and 14 species, respectively. The restriction enzyme SacI produced three to six (generally three) bands with detectable hybridization to the probe pTA71. This probe–enzyme combination indicated a remarkable uniformity amongst Brassica and allied genera in the coding region of repeat units. By contrast, an extensive size variation in the restriction fragments could be localized in the intergenic spacer (IGS) region. Eleven IGS‐containing length variants were detected. Complex hybridization patterns, resulting from extensive repeat unit heterogeneity and taxon‐specific methylation of one or more cleavage sites, were obtained with the EcoRI + pTA71 combination. The relative homologies between the coding regions were evident from the presence of 1.5 kb in all the taxa, and 0.4‐, 1.3‐, and 1.7‐kb fragments in 33, 27, and 24 species, respectively. The SacI + pTA71 and EcoRI + pTA71 combinations were generally able to distinguish taxa both within and between genera. Three restriction endonuclease digests probed with three ribosomal clones yielded essentially identical fragmentation patterns across all the accessions within the cultivated species Brassica campestris, B. oleracea, and B. juncea. In B. napus, three and seven accessions exhibited restriction profiles similar to one and both diploid progenitor species, respectively. Overall, rDNA repeat unit length polymorphism showed good correlation with the cytodeme‐based classification of Brassica and allied genera. © 2008 The Linnean Society of London, Botanical Journal of the Linnean Society, 2008, 157 , 545–557.