Genetics of adaptive radiation in Hawaiian and Cook Islands species of Tetramolopium (Asteraceae). II. Genetic linkage map and its implications for interspecific breeding barriers

High-density interspecific genetic linkage mapping provides insights into genomic incompatibility between channel catfish and blue catfish

Catfish is the leading aquaculture species in the United States. The interspecific hybrid catfish produced by mating female channel catfish with male blue catfish outperform both of their parent species in a number of traits. However, mass production of the hybrids has been difficult because of reproductive isolation. Investigations of genome structure and organization of the hybrids provide insights into the genetic basis for maintenance of species divergence in the face of gene flow, thereby helping develop strategies for introgression and efficient production of the hybrids for aquaculture. In this study, we constructed a high-density genetic linkage map using the hybrid catfish system with the catfish 250K SNP array. A total of 26,238 SNPs were mapped to 29 linkage groups, with 12,776 unique marker positions. The linkage map spans approximately 3240 cM with an average intermarker distance of 0.25 cM. A fraction of markers (986 of 12,776) exhibited significant deviation from the expected Mendelian ratio of segregation, and they were clustered in major genomic blocks across 15 LGs, most notably LG9 and LG15. The distorted markers exhibited significant bias for maternal alleles among the backcross progenies, suggesting strong selection against the blue catfish alleles. The clustering of distorted markers within genomic blocks should lend insights into speciation as marked by incompatibilities between the two species. Such findings should also have profound implications for understanding the genomic evolution of closely related species as well as the introgression of hybrid production programs in aquaculture.

A hybrid genetic linkage map of two ecologically and morphologically divergent Midas cichlid fishes (Amphilophus spp.) obtained by massively parallel DNA sequencing (ddRADSeq)

Cichlid fishes are an excellent model system for studying speciation and the formation of adaptive radiations because of their tremendous species richness and astonishing phenotypic diversity. Most research has focused on African rift lake fishes, although Neotropical cichlid species display much variability as well. Almost one dozen species of the Midas cichlid species complex (Amphilophus spp.) have been described so far and have formed repeated adaptive radiations in several Nicaraguan crater lakes. Here we apply double-digest restriction-site associated DNA sequencing to obtain a high-density linkage map of an interspecific cross between the benthic Amphilophus astorquii and the limnetic Amphilophus zaliosus, which are sympatric species endemic to Crater Lake Apoyo, Nicaragua. A total of 755 RAD markers were genotyped in 343 F₂ hybrids. The map resolved 25 linkage groups and spans a total distance of 1427 cM with an average marker spacing distance of 1.95 cM, almost matching the total number of chromosomes (n = 24) in these species. Regions of segregation distortion were identified in five linkage groups. Based on the pedigree of parents to F₂ offspring, we calculated a genome-wide mutation rate of 6.6 × 10⁻⁸ mutations per nucleotide per generation. This genetic map will facilitate the mapping of ecomorphologically relevant adaptive traits in the repeated phenotypes that evolved within the Midas cichlid lineage and, as the first linkage map of a Neotropical cichlid, facilitate comparative genomic analyses between African cichlids, Neotropical cichlids and other teleost fishes.

Mapping unexplored genomes II: genetic architecture of species differences in the woody Sonchus alliance (Asteraceae) in the Macaronesian Islands

Despite numerous, well-documented evolutionary histories of plant groups which underwent rapid radiation in various oceanic archipelagos, very little is known about the genetic basis of species differences and adaptive radiation. This paper represents the first such study in the Macaronesian Islands using non-model endemic plants, the woody Sonchus alliance. Here I inferred the genetic basis of species differences between two Canary Island endemics, the herbaceous perennial, shade tolerant Lactucosonchus webbii and the woody, coastal desert perennial Sonchus radicatus by quantitative trait locus (QTL) mapping using AFLP markers. A total of 23 QTL (7.3-23.8% PVE; phenotypic variance explained) for 11 morphological traits were found, one for flowering time (31% PVE), and five QTL (7-10.7% PVE) for two physiological traits (intrinsic water use efficiency and stomatal conductance). Interpreted cautiously, these results suggest that major morphological and some physiological differences between the two species are controlled by numerous genes with small to moderate effect. This implies that major morphological changes in island plants can be more complex than suggested by other studies, such as in Tetramolopium in the Hawaiian Islands. The genetic basis of arborescence on islands, one of the most spectacular convergent features of plants across different lineages and archipelagos, is also discussed.

A second-generation anchored genetic linkage map of the tammar wallaby (Macropus eugenii)

Background

The tammar wallaby, Macropus eugenii, a small kangaroo used for decades for studies of reproduction and metabolism, is the model Australian marsupial for genome sequencing and genetic investigations. The production of a more comprehensive cytogenetically-anchored genetic linkage map will significantly contribute to the deciphering of the tammar wallaby genome. It has great value as a resource to identify novel genes and for comparative studies, and is vital for the ongoing genome sequence assembly and gene ordering in this species.

Results

A second-generation anchored tammar wallaby genetic linkage map has been constructed based on a total of 148 loci. The linkage map contains the original 64 loci included in the first-generation map, plus an additional 84 microsatellite loci that were chosen specifically to increase coverage and assist with the anchoring and orientation of linkage groups to chromosomes. These additional loci were derived from (a) sequenced BAC clones that had been previously mapped to tammar wallaby chromosomes by fluorescence in situ hybridization (FISH), (b) End sequence from BACs subsequently FISH-mapped to tammar wallaby chromosomes, and (c) tammar wallaby genes orthologous to opossum genes predicted to fill gaps in the tammar wallaby linkage map as well as three X-linked markers from a published study. Based on these 148 loci, eight linkage groups were formed. These linkage groups were assigned (via FISH-mapped markers) to all seven autosomes and the X chromosome. The sex-pooled map size is 1402.4 cM, which is estimated to provide 82.6% total coverage of the genome, with an average interval distance of 10.9 cM between adjacent markers. The overall ratio of female/male map length is 0.84, which is comparable to the ratio of 0.78 obtained for the first-generation map.

Conclusions

Construction of this second-generation genetic linkage map is a significant step towards complete coverage of the tammar wallaby genome and considerably extends that of the first-generation map. It will be a valuable resource for ongoing tammar wallaby genetic research and assembling the genome sequence. The sex-pooled map is available online at http://compldb.angis.org.au/.

Modeling segregation distortion for viability selection. I. Reconstruction of linkage maps with distorted markers

Molecular markers have been widely used to map quantitative trait loci (QTL). The QTL mapping partly relies on accurate linkage maps. The non-Mendelian segregation of markers, which affects not only the estimation of genetic distance between two markers but also the order of markers on a same linkage group, is usually observed in QTL analysis. However, these distorted markers are often ignored in the real data analysis of QTL mapping so that some important information may be lost. In this paper, we developed a multipoint approach via Hidden Markov chain model to reconstruct the linkage maps given a specified gene order while simultaneously making use of distorted, dominant and missing markers in an F(2) population. The new method was compared with the methods in the MapManager and Mapmaker programs, respectively, and verified by a series of Monte Carlo simulation experiments along with a working example. Results showed that the adjusted linkage maps can be used for further QTL or segregation distortion locus (SDL) analysis unless there are strong evidences to prove that all markers show normal Mendelian segregation.

Linkage maps of the dwarf and Normal lake whitefish (Coregonus clupeaformis) species complex and their hybrids reveal the genetic architecture of population divergence

Elucidating the genetic architecture of population divergence may reveal the evolution of reproductive barriers and the genomic regions implicated in the process. We assembled genetic linkage maps for the dwarf and Normal lake whitefish species complex and their hybrids. A total of 877 AFLP loci and 30 microsatellites were positioned. The homology of mapped loci between families supported the existence of 34 linkage groups (of 40n expected) exhibiting 83% colinearity among linked loci between these two families. Classes of AFLP markers were not randomly distributed among linkage groups. Both AFLP and microsatellites exhibited deviations from Mendelian expectations, with 30.4% exhibiting significant segregation distortion across 28 linkage groups of the four linkage maps in both families (P < 0.00001). Eight loci distributed over seven homologous linkage groups were significantly distorted in both families and the level of distortion, when comparing homologous loci of the same phase between families, was correlated (Spearman R = 0.378, P = 0.0021). These results suggest that substantial divergence incurred during allopatric glacial separation and subsequent sympatric ecological specialization has resulted in several genomic regions that are no longer complementary between dwarf and Normal populations issued from different evolutionary glacial lineages.

The genetic basis of intrinsic and extrinsic post-zygotic reproductive isolation jointly promoting speciation in the lake whitefish species complex (Coregonus clupeaformis)

Understanding the genetic architecture of reproductive barriers and the evolutionary forces that drove their divergence represents a considerable challenge towards understanding speciation. The objective of this study was to determine the genetic basis of intrinsic and extrinsic post-zygotic isolation in diverging populations of dwarf and normal lake whitefish with allopatric glacial origins. We found that the rate of embryonic mortality was 5.3-6.5 times higher in dwarf-normal hybrid backcrosses during development than in F1 dwarf and normal crosses. When comparing embryos that died during development against larvae that successfully hatched, patterns of Mendelian segregation at 101 loci whose linkage is known identified 13 loci distributed over seven linkage groups that exhibited significant shifts in segregation ratios leading to significant segregation distortion at these loci in the surviving progeny. Controlled crosses and quantitative trait loci analysis revealed a significant genetic basis for developmental time until emergence, a trait critical to fish larval survival in nature. Hatching backcross progeny exhibited asynchronous emergence and transgressive segregation, suggesting that extrinsic post-zygotic isolation may select against hybridization in specific environmental contexts. Evidence of a genetic basis for increased embryonic mortality followed by asynchronous emergence indicated that intrinsic and extrinsic mechanisms are not mutually exclusive in the formation and maintenance of reproductive isolation, but may be jointly promoting population divergence and ultimately speciation.

A genetic linkage map ofSilene vulgarisbased on AFLP markers

A genetic linkage map of an intraspecific cross between 2 Silene vulgaris s.l. ecotypes is presented. Three-hundred AFLP markers from 2 different restriction enzyme combinations were used to genotype an F2mapping population. Maternal and paternal pure-coupling phase maps with 114 and 186 markers on 12 and 13 linkage groups, respectively, were constructed. Total map length of the paternal and maternal maps are 547 and 446 Kosambi cM, respectively. Nearly half of the markers (49%) exhibited significant transmission ratio distortion. Genome coverage and potential causes of the observed segregation ratio distortions are discussed. The maps represent a first step towards the identification of quantitative trait loci associated with habitat adaptation in the non-model species Silene vulgaris.Key words: AFLP, genome mapping, segregation distortion, Silene vulgaris.

Genetic mapping of species boundaries in Louisiana irises using IRRE retrotransposon display markers

Genetic mapping studies provide insight into the pattern and extent of genetic incompatibilities affecting hybridization between closely related species. Genetic maps of two species of Louisiana Irises, Iris fulva and I. brevicaulis, were constructed from transposon-based molecular markers segregating in reciprocal backcross (BC1) interspecific hybrids and used to investigate genomic patterns of species barriers inhibiting introgression. Linkage mapping analyses indicated very little genetic incompatibility between I. fulva and I. brevicaulis in the form of map regions exhibiting transmission ratio distortion, and this was confirmed using a Bayesian multipoint mapping analysis. These results demonstrate the utility of transposon-based marker systems for genetic mapping studies of wild plant species and indicate that the genomes of I. fulva and I. brevicaulis are highly permeable to gene flow and introgression from one another via backcrossing.

Transmission ratio distortion in intraspecific hybrids of Mimulus guttatus: implications for genomic divergence

We constructed a genetic linkage map between two divergent populations of Mimulus guttatus. We genotyped an F(2) mapping population (N = 539) at 154 AFLP, microsatellite, and gene-based markers. A framework map was constructed consisting of 112 marker loci on 14 linkage groups with a total map length of 1518 cM Kosambi. Nearly half of all markers (48%) exhibited significant transmission ratio distortion (alpha = 0.05). By using a Bayesian multipoint mapping method and visual inspection of significantly distorted markers, we detected 12 transmission ratio distorting loci (TRDL) throughout the genome. The high degree of segregation distortion detected in this intraspecific map indicates substantial genomic divergence that perhaps suggests genomic incompatibilities between these two populations. We compare the pattern of transmission ratio distortion in this map to an interspecific map constructed between M. guttatus and M. nasutus. A similar level of segregation distortion is detected in both maps. Collinear regions between maps are compared to determine if there are shared genetic patterns of non-Mendelian segregation distortion within and among Mimulus species.

Genetics of postzygotic isolation in Eucalyptus: whole-genome analysis of barriers to introgression in a wide interspecific cross of Eucalyptus grandis and E. globulus

The genetic architecture of hybrid fitness characters can provide valuable insights into the nature and evolution of postzygotic reproductive barriers in diverged species. We determined the genome-wide distribution of barriers to introgression in an F(1) hybrid of two Eucalyptus tree species, Eucalyptus grandis (W. Hill ex Maiden.) and E. globulus (Labill.). Two interspecific backcross families (N = 186) were used to construct comparative, single-tree, genetic linkage maps of an F(1) hybrid individual and two backcross parents. A total of 1354 testcross AFLP marker loci were evaluated in the three parental maps and a substantial proportion (27.7% average) exhibited transmission ratio distortion (alpha = 0.05). The distorted markers were located in distinct regions of the parental maps and marker alleles within each region were all biased toward either of the two parental species. We used a Bayesian approach to estimate the position and effect of transmission ratio distorting loci (TRDLs) in the distorted regions of each parental linkage map. The relative viability of TRDL alleles ranged from 0.20 to 0.72. Contrary to expectation, heterospecific (donor) alleles of TRDLs were favored as often as recurrent alleles in both backcrosses, suggesting that positive and negative heterospecific interactions affect introgression rates in this wide interspecific pedigree.

Comparing the linkage maps of the close relatives Arabidopsis lyrata and A. thaliana

We have constructed a genetic map of Arabidopsis lyrata, a self-incompatible relative of the plant model species A. thaliana. A. lyrata is a diploid (n = 8) species that diverged from A. thaliana (n = 5) approximately 5 MYA. Mapping was conducted in a full-sib progeny of two unrelated F(1) hybrids between two European populations of A. lyrata ssp. petraea. We used the least-squares method of the Joinmap program for map construction. The gross chromosomal differences between the two species were most parsimoniously explained with three fusions, two reciprocal translocations, and one inversion. The total map length was 515 cM, and the distances were 12% larger than those between corresponding markers in the linkage map of A. thaliana. The 72 markers, consisting of microsatellites and gene-based markers, were spaced on average every 8 cM. Transmission ratio distortion was extensive, and most distortions were specific to each reciprocal cross, suggesting cytoplasmic interactions. We estimate locations and most probable genotype frequencies of transmission ratio distorting loci (TRDL) with a Bayesian method and discuss the possible reasons for the observed distortions.

Genetic mapping in (Populus tomentosa x Populus bolleana) and P. tomentosa Carr. using AFLP markers

The AFLP genetic linkage maps for two poplar cultivars were constructed with the pseudo-test-cross mapping strategy. The hybrids were derived from an interspecific backcross between the female hybrid clone "TB01" (Populus tomentosa x Populus bolleana) and the male clone "LM50" ( P. tomentosa). A total of 782 polymorphic fragments were obtained with a PCR-based strategy using 49 enzyme-nested ( EcoRI/ MseI) primer combinations. Six hundred and thirty two of these fragments segregated in a 1:1 ratio ( P<0.01), indicating that these DNA polymorphisms are heterozygous in one parent and null in the other. The linkage analysis was performed using Mapmaker version 3.0 with LOD 5.0 and a maximum recombination fraction (theta) of 0.3. Map distances were estimated using the Kosambi mapping function. In the framework map for "LM50" (P. tomentosa), 218 markers were aligned in 19 major linkage groups. The linked loci spanned approximately 2,683 cM of the poplar genome, with an average distance of 12.3 cM between adjacent markers. For "TB01" (P. tomentosa x P. bolleana), the analysis revealed 144 loci, which were mapped to 19 major linkage groups and covered about 1,956 cM, with an average distance of 13.6 cM between adjacent markers. These maps covered about 87% and 77% of the estimated genome size of parents "LM50" and "TB01", respectively. The maps developed in this study lay an important foundation for future genomics research in poplar, providing a means for localizing genes controlling economically important traits in P. tomentosa.

Integration of populations and differentiation of species

The framework for modern studies of speciation was established as part of the Neo-Darwinian synthesis of the early twentieth century. Here we evaluate this framework in the light of recent empirical and theoretical studies. Evidence from experimental studies of selection, quantitative genetic studies of species' differences, and the molecular evolution of 'isolation' genes, all agree that directional selection is the primary cause of speciation, as initially proposed by Darwin. Likewise, as suggested by Dobzhansky and Mayr, gene flow does hold species together, but probably more by facilitating the spread of beneficial mutants and associated hitchhiking events than by homogenizing neutral loci. Reproductive barriers are important as well in that they preserve adaptations, but as has been stressed by botanists for close to a century, they rarely protect the entire genome from gene flow in recently diverged species. Contrary to early views, it is now clear that speciation can occur in the presence of gene flow. However, recent theory does support the long-held view that population structure and small population size may increase speciation rates, but only under special conditions and not because of the increased efficacy of drift as suggested by earlier authors. Rather, low levels of migration among small populations facilitates the rapid accumulation of beneficial mutations that indirectly cause hybrid incompatibilities.

The potential of Prunus davidiana for introgression into peach [Prunus persica (L.) Batsch] assessed by comparative mapping

The potential for introgression of Prunus davidiana, a wild species related to peach, was evaluated with respect to problems of non-Mendelian segregation or suppressed recombination which often hamper breeding processes based on interspecific crosses. Three connected (F1, F2 and BC2) populations, derived from a cross between P. davidiana clone P1908 and the peach cultivar Summergrand were used. The intraspecific map of P. davidiana already established using the F1 progeny was complemented, and two interspecific maps, for the F2 and BC2 progenies, were built with a set of markers selected from the Prunus reference map. With the molecular data collected for the F2 map construction, regions with distorted marker segregation were detected on the genome; one third of all loci deviated significantly from the expected Mendelian ratios. However, some of these distorted segregations were probably not due to the interspecific cross. On linkage group 6, a skewed area under gametic selection was most likely influenced by the self-incompatibility gene of P. davidiana. Using anchor loci, a good colinearity between the three maps built and the Prunus reference map was demonstrated. Comparative mapping also revealed that homologous recombination occurred normally between P. davidiana and the Prunus persica genome. This confirmed the closeness of the two species. Higher recombination rates were generally observed between P. davidiana and P. persica than between Prunus amygdalus and P. persica. The consequences for plant breeding strategy are discussed. The three maps of the F1, F2 and BC2 progenies provide useful tools for QTL detection and marker-assisted selection, as well as for assessing the efficiency of the peach breeding scheme applied to introgress P. davidiana genes into peach cultivated varieties.

The First Comprehensive Genetic Linkage Map of a Marsupial: The Tammar Wallaby (Macropus eugenii)

The production of a marsupial genetic linkage map is perhaps one of the most important objectives in marsupial research. This study used a total of 353 informative meioses and 64 genetic markers to construct a framework genetic linkage map for the tammar wallaby (Macropus eugenii). Nearly all markers (93.8%) formed a significant linkage (LOD &gt; 3.0) with at least one other marker, indicating that the majority of the genome had been mapped. In fact, when compared with chiasmata data, &gt;70% (828 cM) of the genome has been covered. Nine linkage groups were identified, with all but one (LG7; X-linked) allocated to the autosomes. These groups ranged in size from 15.7 to 176.5 cM and have an average distance of 16.2 cM between adjacent markers. Of the autosomal linkage groups (LGs), LG2 and LG3 were assigned to chromosome 1 and LG4 localized to chromosome 3 on the basis of physical localization of genes. Significant sex-specific distortions toward reduced female recombination rates were revealed in 22% of comparisons. When comparing the X chromosome data to closely related species it is apparent that they are conserved in both synteny and gene order.

A genetic map in the Mimulus guttatus species complex reveals transmission ratio distortion due to heterospecific interactions

As part of a study of the genetics of floral adaptation and speciation in the Mimulus guttatus species complex, we constructed a genetic linkage map of an interspecific cross between M. guttatus and M. nasutus. We genotyped an F(2) mapping population (N = 526) at 255 AFLP, microsatellite, and gene-based markers and derived a framework map through repeated rounds of ordering and marker elimination. The final framework map consists of 174 marker loci on 14 linkage groups with a total map length of 1780 cM Kosambi. Genome length estimates (2011-2096 cM) indicate that this map provides thorough coverage of the hybrid genome, an important consideration for QTL mapping. Nearly half of the markers in the full data set (49%) and on the framework map (48%) exhibited significant transmission ratio distortion (alpha = 0.05). We localized a minimum of 11 transmission ratio distorting loci (TRDLs) throughout the genome, 9 of which generate an excess of M. guttatus alleles and a deficit of M. nasutus alleles. This pattern indicates that the transmission ratio distortion results from particular interactions between the heterospecific genomes and suggests that substantial genetic divergence has occurred between these Mimulus species. We discuss possible causes of the unequal representation of parental genomes in the F(2) generation.

Molecular phylogenetics of the sexually deceptive orchid genus Ophrys (Orchidaceae) based on nuclear and chloroplast DNA sequences

We present a phylogenetic analysis of the major lineages of the sexually deceptive orchid genus Ophrys based on nuclear ribosomal (nr) DNA (internal transcribed spacer region) and noncoding chloroplast (cp) DNA (trnL-trnF region) sequences. Sequence divergence within and among major Ophrys lineages was low for both nrDNA and cpDNA sequences. Separate analyses resulted in similar but poorly resolved trees. An incongruence length difference test revealed that nrDNA and cpDNA data sets were not incongruent. A combined analysis resulted in a better-resolved phylogenetic hypothesis of relationships among the major Ophrys lineages. Our data strongly support a division of Ophrys into two groups. These groups do not correspond to the earlier proposed sections Euophrys and Pseudophrys and are thus in conflict with traditional classifications. Our results support a well-resolved monophyletic group that contains the geographically widespread O. bombyliflora, O. speculum, O. tenthredinifera, and the O. fusca-lutea lineage. Relationships in the other group are poorly resolved. Based on our observations that taxa with identical sequences at presumably rapidly evolving loci clearly differ in floral morphology, we hypothesize that the diversity in the genus Ophrys is the result of a recent radiation in this orchid lineage.

Multipoint mapping of viability and segregation distorting loci using molecular markers

In line-crossing experiments, deviations from Mendelian segregation ratios are usually observed for some markers. We hypothesize that these deviations are caused by one or more segregation-distorting loci (SDL) linked to the markers. We develop both a maximum-likelihood (ML) method and a Bayesian method to map SDL using molecular markers. The ML mapping is implemented via an EM algorithm and the Bayesian method is performed via the Markov chain Monte Carlo (MCMC). The Bayesian mapping is computationally more intensive than the ML mapping but can handle more complicated models such as multiple SDL and variable number of SDL. Both methods are applied to a set of simulated data and real data from a cross of two Scots pine trees.

Genetic architecture of species differences in annual sunflowers: implications for adaptive trait introgression

Genetic architecture may profoundly influence the ability of adaptive traits to spread between species via introgressive hybridization. Here, we examine the genomic location of quantitative trait loci (QTL) associated with pollen sterility and morphological traits distinguishing two annual sunflowers, Helianthus annuus and H. debilis ssp. cucumerifolius. These species are of particular interest since they hybridize naturally, and the form of H. annuus in Texas (called ssp. texanus) is thought to have arisen through introgression. Analysis of 226 BC(1) progeny from a cross between H. annuus and H. debilis revealed 56 QTL for 15 morphological traits and 2 QTL for pollen sterility. Four morphological QTL are tightly linked (<10 cM) to one or more sterility factors and 7 are closely allied with underrepresented and presumably negatively selected chromosomal blocks. Although these 11 QTL seem unlikely to move between the species, no barrier to introgression was detected for the remaining 45 morphological QTL. In fact, due to widespread pleiotropy (or tight linkage), the introgression of just three small chromosomal blocks appears sufficient to largely recover the phenotype of ssp. texanus. Subsequent work will test for the occurrence and fitness consequences of the identified QTL in natural populations of ssp. texanus.