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

Postzygotic barriers persist despite ongoing introgression in hybridizing Mimulus species

The evolution of postzygotic isolation is thought to be a key step in maintaining species boundaries upon secondary contact, yet the dynamics and persistence of hybrid incompatibilities in sympatric species are not well understood.Here, we explore these issues using genetic mapping in three populations of recombinant inbred lines between naturally hybridizing monkeyflowers Mimulus guttatus and M. nasutus from the sympatric Catherine Creek population.The three M. guttatus founders differ dramatically in admixture history. Comparative genetic mapping also reveals three putative inversions segregating among the M. guttatus founders, two due to admixture. We observe strong, genome-wide transmission ratio distortion, but patterns are highly variable among populations. Some distortion is explained by epistatic selection favoring parental genotypes, but tests of inter-chromosomal linkage disequilibrium also reveal multiple candidate Dobzhansky-Muller incompatibilities. We also map several genetic loci for hybrid fertility, including two interacting pairs coinciding with peaks of distortion.Remarkably, in this limited sample of M. guttatus, we discover abundant segregating variation for hybrid incompatibilities with M. nasutus, suggesting this population harbors diverse contributors to postzygotic isolation. Moreover, even with substantial admixture, hybrid incompatibilities between Mimulus species persist, suggesting postzygotic isolation might be a potent force in maintaining species barriers in this system.

Unravelling transmission ratio distortion across the bovine genome: identification of candidate regions for reproduction defects

BACKGROUND

Biological mechanisms affecting gametogenesis, embryo development and postnatal viability have the potential to alter Mendelian inheritance expectations resulting in observable transmission ratio distortion (TRD). Although the discovery of TRD cases have been around for a long time, the current widespread and growing use of DNA technologies in the livestock industry provides a valuable resource of large genomic data with parent-offspring genotyped trios, enabling the implementation of TRD approach. In this research, the objective is to investigate TRD using SNP-by-SNP and sliding windows approaches on 441,802 genotyped Holstein cattle and 132,991 (or 47,910 phased) autosomal SNPs.

RESULTS

The TRD was characterized using allelic and genotypic parameterizations. Across the whole genome a total of 604 chromosomal regions showed strong significant TRD. Most (85%) of the regions presented an allelic TRD pattern with an under-representation (reduced viability) of carrier (heterozygous) offspring or with the complete or quasi-complete absence (lethality) for homozygous individuals. On the other hand, the remaining regions with genotypic TRD patterns exhibited the classical recessive inheritance or either an excess or deficiency of heterozygote offspring. Among them, the number of most relevant novel regions with strong allelic and recessive TRD patterns were 10 and 5, respectively. In addition, functional analyses revealed candidate genes regulating key biological processes associated with embryonic development and survival, DNA repair and meiotic processes, among others, providing additional biological evidence of TRD findings.

CONCLUSIONS

Our results revealed the importance of implementing different TRD parameterizations to capture all types of distortions and to determine the corresponding inheritance pattern. Novel candidate genomic regions containing lethal alleles and genes with functional and biological consequences on fertility and pre- and post-natal viability were also identified, providing opportunities for improving breeding success in cattle.

Diversity and determinants of recombination landscapes in flowering plants

During meiosis, crossover rates are not randomly distributed along the chromosome and their location may have a strong impact on the functioning and evolution of the genome. To date, the broad diversity of recombination landscapes among plants has rarely been investigated and a formal comparative genomic approach is still needed to characterize and assess the determinants of recombination landscapes among species and chromosomes. We gathered genetic maps and genomes for 57 flowering plant species, corresponding to 665 chromosomes, for which we estimated large-scale recombination landscapes. We found that the number of crossover per chromosome spans a limited range (between one to five/six) whatever the genome size, and that there is no single relationship across species between genetic map length and chromosome size. Instead, we found a general relationship between the relative size of chromosomes and recombination rate, while the absolute length constrains the basal recombination rate for each species. At the chromosome level, we identified two main patterns (with a few exceptions) and we proposed a conceptual model explaining the broad-scale distribution of crossovers where both telomeres and centromeres play a role. These patterns correspond globally to the underlying gene distribution, which affects how efficiently genes are shuffled at meiosis. These results raised new questions not only on the evolution of recombination rates but also on their distribution along chromosomes.

Meiotic recombination is a universal feature of sexually reproducing species. During meiosis, crossovers play a fundamental role for the proper segregation of chromosomes during meiosis and reshuffles alleles among chromosomes. How much variation in recombination is expected within a genome and among different species remains a central question for understanding the evolution of recombination. We characterized and compared recombination landscapes in a large set of plant species with a wide range of genome size. We found that the number of crossovers varied little among species, from one mandatory to no more than five or six crossovers per chromosomes, whatever the genome size. However, we identified two main patterns of variation along chromosomes (with a few exceptions) that can be explained by a new conceptual model where chromosome length, chromosome structure and gene density play a role. The strong association between gene density and recombination was already known, but raised new questions not only about the evolution of recombination rates but also on their distribution along chromosomes.

OUP accepted manuscript

Substantial morphological variation in land plants remains inaccessible to genetic analysis because current models lack variation in important ecological and agronomic traits. The genus Gilia was historically a model for biosystematics studies and includes variation in morphological traits that are poorly understood at the genetic level. We assembled a chromosome-scale reference genome of G. yorkii and used it to investigate genome evolution in the Polemoniaceae. We performed QTL (quantitative trait loci) mapping in a G. yorkii×G. capitata interspecific population for traits related to inflorescence architecture and flower color. The genome assembly spans 2.75 Gb of the estimated 2.80-Gb genome, with 96.7% of the sequence contained in the nine largest chromosome-scale scaffolds matching the haploid chromosome number. Gilia yorkii experienced at least one round of whole-genome duplication shared with other Polemoniaceae after the eudicot paleohexaploidization event. We identified QTL linked to variation in inflorescence architecture and petal color, including a candidate for the major flower color QTL—a tandem duplication of flavanol 3′,5′-hydroxylase. Our results demonstrate the utility of Gilia as a forward genetic model for dissecting the evolution of development in plants including the causal loci underlying inflorescence architecture transitions.

Genomic Resources for Darters (Percidae: Etheostominae) Provide Insight into Postzygotic Barriers Implicated in Speciation

Comparative genomic approaches are increasingly being used to study the evolution of reproductive barriers in nonmodel species. Although numerous studies have examined prezygotic isolation in darters (Percidae), investigations into postzygotic barriers have remained rare due to long generation times and a lack of genomic resources. Orangethroat and rainbow darters naturally hybridize and provide a remarkable example of male-driven speciation via character displacement. Backcross hybrids suffer from high mortality, which appears to promote behavioral isolation in sympatry. To investigate the genomic architecture of postzygotic isolation, we used Illumina and PacBio sequencing to generate a chromosome-level, annotated assembly of the orangethroat darter genome and high-density linkage maps for orangethroat and rainbow darters. We also analyzed genome-wide RADseq data from wild-caught adults of both species and laboratory-generated backcrosses to identify genomic regions associated with hybrid incompatibles. Several putative chromosomal translocations and inversions were observed between orangethroat and rainbow darters, suggesting structural rearrangements may underlie postzygotic isolation. We also found evidence of selection against recombinant haplotypes and transmission ratio distortion in backcross hybrid genomes, providing further insight into the genomic architecture of genetic incompatibilities. Notably, regions with high levels of genetic divergence between species were enriched for genes associated with developmental and meiotic processes, providing strong candidates for postzygotic isolating barriers. These findings mark significant contributions to our understanding of the genetic basis of reproductive isolation between species undergoing character displacement. Furthermore, the genomic resources presented here will be instrumental for studying speciation in darters, the most diverse vertebrate group in North America.

Patterns of Hybrid Seed Inviability in the Mimulus guttatus sp. Complex Reveal a Potential Role of Parental Conflict in Reproductive Isolation

Genomic conflicts may play a central role in the evolution of reproductive barriers. Theory predicts that early-onset hybrid inviability may stem from conflict between parents for resource allocation to offspring. Here, we describe M. decorus: a group of cryptic species within the M. guttatus species complex that are largely reproductively isolated by hybrid seed inviability (HSI). HSI between M. guttatus and M. decorus is common and strong, but populations of M. decorus vary in the magnitude and directionality of HSI with M. guttatus. Patterns of HSI between M. guttatus and M. decorus, as well as within M. decorus, conform to the predictions of parental conflict: first, reciprocal F1s exhibit size differences and parent-of-origin-specific endosperm defects; second, the extent of asymmetry between reciprocal F1 seed size is correlated with asymmetry in HSI; and third, inferred differences in the extent of conflict predict the extent of HSI between populations. We also find that HSI is rapidly evolving, as populations that exhibit the most HSI are each others' closest relative. Lastly, although all populations appear largely outcrossing, we find that the differences in the inferred strength of conflict scale positively with π, suggesting that demographic or life history factors other than transitions to self-fertilization may influence the rate of parental-conflict-driven evolution. Overall, these patterns suggest the rapid evolution of parent-of-origin-specific resource allocation alleles coincident with HSI within and between M. guttatus and M. decorus. Parental conflict may therefore be an important evolutionary driver of reproductive isolation.

Implementation of Bayesian methods to identify SNP and haplotype regions with transmission ratio distortion across the whole genome: TRDscan v.1.0

Realized deviations from the expected Mendelian inheritance of alleles from heterozygous parents have been previously reported in a broad range of organisms (i.e., transmission ratio distortion; TRD). Various biological mechanisms affecting gametes, embryos, fetuses, or even postnatal offspring can produce patterns of TRD. However, knowledge about its prevalence and potential causes in livestock species is still scarce. Specific Bayesian models have been recently developed for the analyses of TRD for biallelic loci, which accommodated a wide range of population structures, enabling TRD investigation in livestock populations. The parameterization of these models is flexible and allows the study of overall (parent-unspecific) TRD and sire- and dam-specific TRD. This research aimed at deriving Bayesian models for fitting TRD on the basis of haplotypes, testing the models for both haplotype- and SNP-based methods in simulated data and actual Holstein genotypes, and developing a specific software for TRD analyses. Results obtained on simulated data sets showed that the statistical power of the analysis increased with sample size of trios (n), proportion of heterozygous parents, and the magnitude of the TRD. On the other hand, the statistical power to detect TRD decreased with the number of alleles at each loci. Bayesian analyses showed a strong Pearson correlation coefficient (≥0.97) between simulated and estimated TRD that reached the significance level of Bayes factor ≥10 for both single-marker and haplotype analyses when n ≥ 25. Moreover, the accuracy in terms of the mean absolute error decreased with the increase of the sample size and increased with the number of alleles at each loci. Using real data (55,732 genotypes of Holstein trios), SNP- and haplotype-based distortions were detected with overall TRD, sire-TRD, or dam-TRD, showing different magnitudes of TRD and statistical relevance. Additionally, the haplotype-based method showed more ability to capture TRD compared with individual SNP. To discard possible random TRD in real data, an approximate empirical null distribution of TRD was developed. The program TRDscan v.1.0 was written in Fortran 2008 language and provides a powerful statistical tool to scan for TRD regions across the whole genome. This developed program is freely available at http://www.casellas.info/files/TRDscan.zip.

Completing the hybridization triangle: the inheritance of genetic incompatibilities during homoploid hybrid speciation in ragworts (Senecio)

A new homoploid hybrid lineage needs to establish a degree of reproductive isolation from its parent species if it is to persist as an independent entity, but the role hybridization plays in this process is known in only a handful of cases. The homoploid hybrid ragwort species, Senecio squalidus (Oxford ragwort), originated following the introduction of hybrid plants to the UK approximately 320 years ago. The source of the hybrid plants was from a naturally occurring hybrid zone between S. aethnensis and S. chrysanthemifolius on Mount Etna, Sicily. Previous studies of the parent species found evidence for multiple incompatibility loci causing transmission ratio distortion of genetic markers in their hybrid progeny. This study closes the hybridization triangle by reporting a genetic mapping analysis of the remaining two paired cross combinations between S. squalidus and its parents. Genetic maps produced from F₂ mapping families were generally collinear but with half of the linkage groups showing evidence of genomic reorganization between genetic maps. The new maps produced from crosses between S. squalidus and each parent showed multiple incompatibility loci distributed across the genome, some of which co-locate with previously reported incompatibility loci between the parents. These findings suggest that this young homoploid hybrid species has inherited a unique combination of genomic rearrangements and incompatibilities from its parents that contribute to its reproductive isolation.

Patterns of Transmission Ratio Distortion in Interspecific Lettuce Hybrids Reveal a Sex-Independent Gametophytic Barrier

Interspecific crosses can result in progeny with reduced vitality or fertility due to genetic incompatibilities between species, a phenomenon known as hybrid incompatibility (HI). HI is often caused by a bias against deleterious allele combinations, which results in transmission ratio distortion (TRD). Here, we determined the genome-wide distribution of HI between wild lettuce, Lactuca saligna, and cultivated lettuce, L. sativa, in a set of backcross inbred lines (BILs) with single introgression segments from L. saligna introgressed into a L. sativa genetic background. Almost all BILs contained an introgression segment in a homozygous state except a few BILs, for which we were able to obtain only a single heterozygous introgression. Their inbred progenies displayed severe TRD with a bias toward the L. sativa allele and complete nontransmission of the homozygous L. saligna introgression, i.e., absolute HI. These HI might be caused by deleterious heterospecific allele combinations at two loci. We used an multilocus segregating interspecific F2 population to identify candidate conspecific loci that can nullify the HI in BILs. Segregation analysis of developed double-introgression progenies showed nullification of three HI and proved that these HI are explained by nuclear pairwise incompatibilities. One of these digenic HI showed 29% reduced seed set and its pattern of TRD pointed to a sex-independent gametophytic barrier. Namely, this HI was caused by complete nontransmission of one heterospecific allele combination at the haploid stage, surprisingly in both male and female gametophytes. Our study shows that two-locus incompatibility systems contribute to reproductive barriers among Lactuca species.

Interspecific genetic maps in Miscanthus floridulus and M. sacchariflorus accelerate detection of QTLs associated with plant height and inflorescence

Miscanthus is recognized as a promising lignocellulosic crop for the production of bioethanol and bioproducts worldwide. To facilitate the identification of agronomical important traits and establish genetics knowledge, two genetic maps were developed from a controlled interspecific cross between M. floridulus and M. sacchariflorus. A total of 650 SSR markers were mapped in M. floridulus, spanning 19 linkage groups and 2053.31 cM with an average interval of 3.25 cM. The map of M. sacchariflorus comprised 495 SSR markers in 19 linkage groups covering 1684.86 cM with an average interval of 3.54 cM. The estimation on genome length indicated that the genome coverage of parental genetic maps were 93.87% and 89.91%, respectively. Eighty-eight bi-parental common markers were allowed to connect the two maps, and six pairs of syntenic linkage groups were recognized. Furthermore, quantitative trait loci (QTL) mapping of three agronomic traits, namely, plant height (PH), heading time (HT), and flowering time (FT), demonstrated that a total of 66 QTLs were identified in four consecutive years using interval mapping and multiple-QTL model. The LOD value of these QTLs ranged from 2.51 to 10.60, and the phenotypic variation explained varied from 9.50 to 37.10%. QTL cluster in syntenic groups MF19/MS7 contained six stable QTLs associated with PH, HT, and FT. In conclusion, we report for the first time the genetic mapping of biomass traits in M. floridulus and M. sacchariflorus. These results will be a valuable genetic resource, facilitating the discovery of essential genes and breeding of Miscanthus.

Linkage Map of Lissotriton Newts Provides Insight into the Genetic Basis of Reproductive Isolation

Linkage maps are widely used to investigate structure, function, and evolution of genomes. In speciation research, maps facilitate the study of the genetic architecture of reproductive isolation by allowing identification of genomic regions underlying reduced fitness of hybrids. Here we present a linkage map for European newts of the Lissotriton vulgaris species complex, constructed using two families of F2 L. montandoni × L. vulgaris hybrids. The map consists of 1146 protein-coding genes on 12 linkage groups, equal to the haploid chromosome number, with a total length of 1484 cM (1.29 cM per marker). It is notably shorter than two other maps available for salamanders, but the differences in map length are consistent with cytogenetic estimates of the number of chiasmata per chromosomal arm. Thus, large salamander genomes do not necessarily translate into long linkage maps, as previously suggested. Consequently, salamanders are an excellent model to study evolutionary consequences of recombination rate variation in taxa with large genomes and a similar number of chromosomes. A complex pattern of transmission ratio distortion (TRD) was detected: TRD occurred mostly in one family, in one breeding season, and was clustered in two genomic segments. This is consistent with environment-dependent mortality of individuals carrying L. montandoni alleles in these two segments and suggests a role of TRD blocks in reproductive isolation. The reported linkage map will empower studies on the genomic architecture of divergence and interactions between the genomes of hybridizing newts.

High-density genetic map construction and QTLs identification for plant height in white jute (Corchorus capsularis L.) using specific locus amplified fragment (SLAF) sequencing

Background

Genetic mapping and quantitative trait locus (QTL) detection are powerful methodologies in plant improvement and breeding. White jute (Corchorus capsularis L.) is an important industrial raw material fiber crop because of its elite characteristics. However, construction of a high-density genetic map and identification of QTLs has been limited in white jute due to a lack of sufficient molecular markers. The specific locus amplified fragment sequencing (SLAF-seq) strategy combines locus-specific amplification and high-throughput sequencing to carry out de novo single nuclear polymorphism (SNP) discovery and large-scale genotyping. In this study, SLAF-seq was employed to obtain sufficient markers to construct a high-density genetic map for white jute. Moreover, with the development of abundant markers, genetic dissection of fiber yield traits such as plant height was also possible. Here, we present QTLs associated with plant height that were identified using our newly constructed genetic linkage groups.

Results

An F₈ population consisting of 100 lines was developed. In total, 69,446 high-quality SLAFs were detected of which 5,074 SLAFs were polymorphic; 913 polymorphic markers were used for the construction of a genetic map. The average coverage for each SLAF marker was 43-fold in the parents, and 9.8-fold in each F₈ individual. A linkage map was constructed that contained 913 SLAFs on 11 linkage groups (LGs) covering 1621.4 cM with an average density of 1.61 cM per locus. Among the 11 LGs, LG1 was the largest with 210 markers, a length of 406.34 cM, and an average distance of 1.93 cM between adjacent markers. LG11 was the smallest with only 25 markers, a length of 29.66 cM, and an average distance of 1.19 cM between adjacent markers. ‘SNP_only’ markers accounted for 85.54% and were the predominant markers on the map. QTL mapping based on the F₈ phenotypes detected 11 plant height QTLs including one major effect QTL across two cultivation locations, with each QTL accounting for 4.14–15.63% of the phenotypic variance.

Conclusions

To our knowledge, the linkage map constructed here is the densest one available to date for white jute. This analysis also identified the first QTL in white jute. The results will provide an important platform for gene/QTL mapping, sequence assembly, genome comparisons, and marker-assisted selection breeding for white jute.

Bayesian analysis of parent-specific transmission ratio distortion in seven Spanish beef cattle breeds

Transmission ratio distortion (TRD) is the departure from the expected Mendelian ratio in offspring, a poorly investigated biological phenomenon in livestock species. Given the current availability of specific parametric methods for the analysis of segregation data, this study focused on the screening of TRD in 602 402 single nucleotide polymorphisms covering all autosomal chromosomes in seven Spanish beef cattle breeds. On average, 0.13% (n = 786) and 0.01% (n = 29) of genetic markers evidenced sire- or dam-specific TRD respectively. There were no single nucleotide polymorphisms accounting for both sire- and dam-specific TRD at the same time, and only one marker (rs43147474) accounted for (sire-specific) TRD in all seven breeds. It must be noted that rs43147474 is located in the fourth intronic region of the GTP-binding protein 10 gene, and this locus has been previously linked to the maintenance of mitochondria and nucleolar architectures. Alternatively, other candidate genes surround this hot-spot for sire-specific TRD in the cattle genome, and they are related to embryonic and postnatal lethality as well as prostate cancer, among others. This research characterized the distribution of TRD in the bovine genome, highlighting heterogeneous results when comparing across breeds.

Comparative Analysis of Regions with Distorted Segregation in Three Diploid Populations of Potato

Genes associated with gametic and zygotic selection could underlie segregation distortion, observed as alterations of expected Mendelian genotypic frequencies in mapping populations. We studied highly dense genetic maps based on single nucleotide polymorphisms to elucidate the genetic nature of distorted segregation in potato. Three intra- and interspecific diploid segregating populations were used. DRH and D84 are crosses between the sequenced doubled monoploid DM 1-3 516 R44 Solanum tuberosum Group Phureja and either RH89-039-16 S. tuberosum or 84SD22, a S. tuberosum × S. chacoense hybrid. MSX902 is an interspecific cross between 84SD22 and Ber83 S. berthaultii × 2 × species mosaic. At the 0.05 significance level, 21%, 57%, and 51% of the total markers mapped in DRH, D84, and MSX902 exhibited distorted segregation, respectively. Segregation distortion regions for DRH were located on chromosomes 9 and 12; for D84 on chromosomes 2, 3, 4, 6, 7, and 8; and on chromosomes 1, 2, 7, 9, and 12 for MSX902. In general, each population had unique segregation distortion regions and directions of distortion. Interspecific crosses showed greater levels of distorted segregation and lower recombination rates as determined from the male parents. The different genomic regions where the segregation distortion regions occurred in the three populations likely reflect unique genetic combinations producing distorted segregation.

Genetic loci with parent-of-origin effects cause hybrid seed lethality in crosses between Mimulus species

In flowering plants, F1 hybrid seed lethality is a common outcome of crosses between closely related diploid species, but the genetic basis of this early-acting and potentially widespread form of postzygotic reproductive isolation is largely unknown. We intercrossed two closely related species of monkeyflower, Mimulus guttatus and Mimulus tilingii, to characterize the mechanisms and strength of postzygotic reproductive isolation. Then, using a reciprocal backcross design, we performed high-resolution genetic mapping to determine the genetic architecture of hybrid seed lethality and directly test for loci with parent-of-origin effects. We found that F1 hybrid seed lethality is an exceptionally strong isolating barrier between Mimulus species, with reciprocal crosses producing < 1% viable seeds. This form of postzygotic reproductive isolation appears to be highly polygenic, indicating that multiple incompatibility loci have accumulated rapidly between these closely related Mimulus species. It is also primarily caused by genetic loci with parent-of-origin effects, suggesting a possible role for imprinted genes in the evolution of Mimulus hybrid seed lethality. Our findings suggest that divergence in loci with parent-of-origin effects, which is probably driven by genomic coevolution within lineages, might be an important source of hybrid incompatibilities between flowering plant species.

The relationship of recombination rate, genome structure, and patterns of molecular evolution across angiosperms

Background

Although homologous recombination affects the efficacy of selection in populations, the pattern of recombination rate evolution and its effects on genome evolution across plants are largely unknown. Recombination can reduce genome size by enabling the removal of LTR retrotransposons, alter codon usage by GC biased gene conversion, contribute to complex histories of gene duplication and loss through tandem duplication, and enhance purifying selection on genes. Therefore, variation in recombination rate across species may explain some of the variation in genomic architecture as well as rates of molecular evolution. We used phylogenetic comparative methods to investigate the evolution of global meiotic recombination rate in angiosperms and its effects on genome architecture and selection at the molecular level using genetic maps and genome sequences from thirty angiosperm species.

Results

Recombination rate is negatively correlated with genome size, which is likely caused by the removal of LTR retrotransposons. After correcting recombination rates for euchromatin content, we also found an association between global recombination rate and average gene family size. This suggests a role for recombination in the preservation of duplicate genes or expansion of gene families. An analysis of the correlation between the ratio of nonsynonymous to synonymous substitution rates (dN/dS) and recombination rate in 3748 genes indicates that higher recombination rates are associated with an increased efficacy of purifying selection, suggesting that global recombination rates affect variation in rates of molecular evolution across distantly related angiosperm species, not just between populations. We also identified shifts in dN/dS for recombination proteins that are associated with shifts in global recombination rate across our sample of angiosperms.

Conclusions

Although our analyses only reveal correlations, not mechanisms, and do not include potential covariates of recombination rate, like effective population size, they suggest that global recombination rates may play an important role in shaping the macroevolutionary patterns of gene and genome evolution in plants. Interspecific recombination rate variation is tightly correlated with genome size as well as variation in overall LTR retrotransposon abundances. Recombination may shape gene-to-gene variation in dN/dS between species, which might impact the overall gene duplication and loss rates.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0473-3) contains supplementary material, which is available to authorized users.

Genetic Map Construction and Detection of Genetic Loci Underlying Segregation Distortion in an Intraspecific Cross of Populus deltoides

Based on a two-way pseudo-testcross strategy, high density and complete coverage linkage maps were constructed for the maternal and paternal parents of an intraspecific F2 pedigree of Populus deltoides. A total of 1,107 testcross markers were obtained, and the mapping population consisted of 376 progeny. Among these markers, 597 were from the mother, and were assigned into 19 linkage groups, spanning a total genetic distance of 1,940.3 cM. The remaining 519 markers were from the father, and were also were mapped into 19 linkage groups, covering 2,496.3 cM. The genome coverage of both maps was estimated as greater than 99.9% at 20 cM per marker, and the numbers of linkage groups of both maps were in accordance with the 19 haploid chromosomes in Populus. Marker segregation distortion was observed in large contiguous blocks on some of the linkage groups. Subsequently, we mapped the segregation distortion loci in this mapping pedigree. Altogether, eight segregation distortion loci with significant logarithm of odds supports were detected. Segregation distortion indicated the uneven transmission of the alternate alleles from the mapping parents. The corresponding genome regions might contain deleterious genes or be associated with hybridization incompatibility. In addition to the detection of segregation distortion loci, the established genetic maps will serve as a basic resource for mapping genetic loci controlling traits of interest in future studies.

Spodoptera frugiperda (Lepidoptera: Noctuidae) host-plant variants: two host strains or two distinct species?

The moth Spodoptera frugiperda is a well-known pest of crops throughout the Americas, which consists of two strains adapted to different host-plants: the first feeds preferentially on corn, cotton and sorghum whereas the second is more associated with rice and several pasture grasses. Though morphologically indistinguishable, they exhibit differences in their mating behavior, pheromone compositions, and show development variability according to the host-plant. Though the latter suggest that both strains are different species, this issue is still highly controversial because hybrids naturally occur in the wild, not to mention the discrepancies among published results concerning mating success between the two strains. In order to clarify the status of the two host-plant strains of S. frugiperda, we analyze features that possibly reflect the level of post-zygotic isolation: (1) first generation (F1) hybrid lethality and sterility; (2) patterns of meiotic segregation of hybrids in reciprocal second generation (F2), as compared to the meiosis of the two parental strains. We found a significant reduction of mating success in F1 in one direction of the cross and a high level of microsatellite markers showing transmission ratio distortion in the F2 progeny. Our results support the existence of post-zygotic reproductive isolation between the two laboratory strains and are in accordance with the marked level of genetic differentiation that was recovered between individuals of the two strains collected from the field. Altogether these results provide additional evidence in favor of a sibling species status for the two strains.

Sequencing, assembling, and correcting draft genomes using recombinant populations

Current de novo whole-genome sequencing approaches often are inadequate for organisms lacking substantial preexisting genetic data. Problems with these methods are manifest as: large numbers of scaffolds that are not ordered within chromosomes or assigned to individual chromosomes, misassembly of allelic sequences as separate loci when the individual(s) being sequenced are heterozygous, and the collapse of recently duplicated sequences into a single locus, regardless of levels of heterozygosity. Here we propose a new approach for producing de novo whole-genome sequences—which we call recombinant population genome construction—that solves many of the problems encountered in standard genome assembly and that can be applied in model and nonmodel organisms. Our approach takes advantage of next-generation sequencing technologies to simultaneously barcode and sequence a large number of individuals from a recombinant population. The sequences of all recombinants can be combined to create an initial de novo assembly, followed by the use of individual recombinant genotypes to correct assembly splitting/collapsing and to order and orient scaffolds within linkage groups. Recombinant population genome construction can rapidly accelerate the transformation of nonmodel species into genome-enabled systems by simultaneously producing a high-quality genome assembly and providing genomic tools (e.g., high-confidence single-nucleotide polymorphisms) for immediate applications. In populations segregating for important functional traits, this approach also enables simultaneous mapping of quantitative trait loci. We demonstrate our method using simulated Illumina data from a recombinant population of Caenorhabditis elegans and show that the method can produce a high-fidelity, high-quality genome assembly for both parents of the cross.

When virulence originates from non-agricultural hosts: new insights into plant breeding

Monogenic plant resistance breakdown is a model for testing evolution in action in pathogens. As a rule, plant pathologists argue that virulence - the allele that allows pathogens to overcome resistance - is due to a new mutation at the avirulence locus within the native/endemic population that infects susceptible crops. In this article, we develop an alternative and neglected scenario where a given virulence pre-exists in a non-agricultural host and might be accidentally released or introduced on the matching resistant cultivar in the field. The main difference between the two scenarios is the divergence time expected between the avirulent and the virulent populations. As a consequence, population genetic approaches such as genome scans and Approximate Bayesian Computation methods allow explicit testing of the two scenarios by timing the divergence. This review then explores the fundamental implications of this alternative scenario for plant breeding, including the invasion of virulence or the evolution of more aggressive hybrids, and proposes concrete solutions to achieve a sustainable resistance.

Selection on plant male function genes identifies candidates for reproductive isolation of yellow monkeyflowers

Understanding the genetic basis of reproductive isolation promises insight into speciation and the origins of biological diversity. While progress has been made in identifying genes underlying barriers to reproduction that function after fertilization (post-zygotic isolation), we know much less about earlier acting pre-zygotic barriers. Of particular interest are barriers involved in mating and fertilization that can evolve extremely rapidly under sexual selection, suggesting they may play a prominent role in the initial stages of reproductive isolation. A significant challenge to the field of speciation genetics is developing new approaches for identification of candidate genes underlying these barriers, particularly among non-traditional model systems. We employ powerful proteomic and genomic strategies to study the genetic basis of conspecific pollen precedence, an important component of pre-zygotic reproductive isolation among yellow monkeyflowers (Mimulus spp.) resulting from male pollen competition. We use isotopic labeling in combination with shotgun proteomics to identify more than 2,000 male function (pollen tube) proteins within maternal reproductive structures (styles) of M. guttatus flowers where pollen competition occurs. We then sequence array-captured pollen tube exomes from a large outcrossing population of M. guttatus, and identify those genes with evidence of selective sweeps or balancing selection consistent with their role in pollen competition. We also test for evidence of positive selection on these genes more broadly across yellow monkeyflowers, because a signal of adaptive divergence is a common feature of genes causing reproductive isolation. Together the molecular evolution studies identify 159 pollen tube proteins that are candidate genes for conspecific pollen precedence. Our work demonstrates how powerful proteomic and genomic tools can be readily adapted to non-traditional model systems, allowing for genome-wide screens towards the goal of identifying the molecular basis of genetically complex traits.

Genetic mapping of two components of reproductive isolation between two sibling species of moths, Ostrinia nubilalis and O. scapulalis

We report the quantitative trait loci (QTL) mapping of reproductive isolation traits between Ostrinia nubilalis (the European corn borer) and its sibling species O. scapulalis (the Adzuki bean borer), focusing on two traits: mating isolation (mi) and pheromone production (Pher). Four genetic maps were generated from two backcross families, with two maps (one chromosomal map and one linkage map) per backcross. We located 165-323 AFLP markers on these four maps, resulting in the identification of 27-31 linkage groups, depending on the map considered. No-choice mating experiments with the offspring of each backcross led to the detection of at least two QTLs for mi in different linkage groups. QTLs underlying Pher were located in a third linkage group. The Z heterochromosome was identified by a specific marker (Tpi) and did not carry any of these QTLs. Finally, we considered the global divergence between the two sibling species, distortions of segregation throughout the genome, and the location and effect of mi and Pher QTLs in light of the known candidate genes for reproductive isolation within the genus Ostrinia and, more broadly, in phytophagous insects.

Cytoplasmic male sterility contributes to hybrid incompatibility between subspecies of Arabidopsis lyrata

In crosses between evolutionarily diverged populations, genomic incompatibilities may result in sterile hybrids, indicating evolution of reproductive isolation. In several plant families, crosses within a population can also lead to male sterile progeny because of conflict between the maternally and biparentally inherited genomes. We examined hybrid fertility between subspecies of the perennial outcrossing self-incompatible Lyrate rockcress (Arabidopsis lyrata) in large reciprocal F2 progenies and three generations of backcrosses. In one of the reciprocal F2 progenies, almost one-fourth of the plants were male-sterile. Correspondingly, almost one-half of the plants in one of the four reciprocal backcross progenies expressed male sterility. In an additional four independent F2 and backcross families, three segregated male sterility. The observed asymmetrical hybrid incompatibility is attributable to male sterility factors in one cytoplasm, for which the other population lacks effective fertility restorers. Genotyping of 96 molecular markers and quantitative trait locus mapping revealed that only 60% of the plants having the male sterile cytoplasm and lacking the corresponding restorers were phenotypically male-sterile. Genotyping data showed that there is only one restorer locus, which mapped to a 600-kb interval at the top of chromosome 2 in a region containing a cluster of pentatricopeptide repeat genes. Male fertility showed no trade-off with seed production. We discuss the role of cytoplasm and genomic conflict in incipient speciation and conclude that cytoplasmic male sterility–lowering hybrid fitness is a transient effect with limited potential to form permanent reproductive barriers between diverged populations of hermaphrodite self-incompatible species.

Chromosome rearrangements, recombination suppression, and limited segregation distortion in hybrids between Yellowstone cutthroat trout (Oncorhynchus clarkii bouvieri) and rainbow trout (O. mykiss)

Background

Introgressive hybridization is an important evolutionary process that can lead to the creation of novel genome structures and thus potentially new genetic variation for selection to act upon. On the other hand, hybridization with introduced species can threaten native species, such as cutthroat trout (Oncorhynchus clarkii) following the introduction of rainbow trout (O. mykiss). Neither the evolutionary consequences nor conservation implications of rainbow trout introgression in cutthroat trout is well understood. Therefore, we generated a genetic linkage map for rainbow-Yellowstone cutthroat trout (O. clarkii bouvieri) hybrids to evaluate genome processes that may help explain how introgression affects hybrid genome evolution.

Results

The hybrid map closely aligned with the rainbow trout map (a cutthroat trout map does not exist), sharing all but one linkage group. This linkage group (RYHyb20) represented a fusion between an acrocentric (Omy28) and a metacentric chromosome (Omy20) in rainbow trout. Additional mapping in Yellowstone cutthroat trout indicated the two rainbow trout homologues were fused in the Yellowstone genome. Variation in the number of hybrid linkage groups (28 or 29) likely depended on a Robertsonian rearrangement polymorphism within the rainbow trout stock. Comparison between the female-merged F₁ map and a female consensus rainbow trout map revealed that introgression suppressed recombination across large genomic regions in 5 hybrid linkage groups. Two of these linkage groups (RYHyb20 and RYHyb25_29) contained confirmed chromosome rearrangements between rainbow and Yellowstone cutthroat trout indicating that rearrangements may suppress recombination. The frequency of allelic and genotypic segregation distortion varied among parents and families, suggesting few incompatibilities exist between rainbow and Yellowstone cutthroat trout genomes.

Conclusions

Chromosome rearrangements suppressed recombination in the hybrids. This result supports several previous findings demonstrating that recombination suppression restricts gene flow between chromosomes that differ by arrangement. Conservation of synteny and map order between the hybrid and rainbow trout maps and minimal segregation distortion in the hybrids suggest rainbow and Yellowstone cutthroat trout genomes freely introgress across chromosomes with similar arrangement. Taken together, these results suggest that rearrangements impede introgression. Recombination suppression across rearrangements could enable large portions of non-recombined chromosomes to persist within admixed populations.

Transmission ratio distortion of molecular markers in a doubled haploid population originated from a natural hybrid between Osmunda japonica and O. lancea

In ferns, intra-gametophytic selfing occurs as a mode of reproduction where two gametes from the same gametophyte form a completely homozygous sporophyte. Intra-gametophytic selfing is considered to be prevented by lethal or deleterious recessive genes in several diploid species. In order to investigate the modes and tempo of selection acting different developmental stages, doubled haploids obtained from intra-gametophytic selfing within isolated gametophytes of a putative F1 hybrid between Osmunda japonica and O. lancea were analyzed with EST_derived molecular markers, and the distribution pattern of transmission ratio distortion (TRD) along linkage map was clarified. As the results, the markers with skewness were clustered in two linkage groups. For the two highly distorted regions, gametophytes and F2 population were also examined. The markers skewed towards O. japonica on a linkage group (LG_2) showed skewness also in gametophytes, and the TRD was generated in the process of spore formation or growth of gametophytes. Also, selection appeared to be operating in the gametophytic stage. The markers on other linkage group (LG_11) showed highest skewness towards O. lancea in doubled haploids, and it was suggested that the segregation of LG_11 were influenced by zygotic lethality or genotypic evaluation and that some deleterious recessive genes exist in LG_11 and reduce the viability of homozygotes with O. japonica alleles. It is very likely that a region of LG_11were responsible for the low frequencies of intra-gametophytic selfing in O. japonica.

Investigating incipient speciation in Arabidopsis lyrata from patterns of transmission ratio distortion

Our understanding of the development of intrinsic reproductive isolation is still largely based on theoretical models and thorough empirical studies on a small number of species. Theory suggests that reproductive isolation develops through accumulation of epistatic genic incompatibilities, also known as Bateson-Dobzhansky-Muller (BDM) incompatibilities. We can detect these from marker transmission ratio distortion (TRD) in hybrid progenies of crosses between species or populations, where TRD is expected to result from selection against heterospecific allele combinations in hybrids. TRD may also manifest itself because of intragenomic conflicts or competition between gametes or zygotes. We studied early stage speciation in Arabidopsis lyrata by investigating patterns of TRD across the genome in F2 progenies of three reciprocal crosses between four natural populations. We found that the degree of TRD increases with genetic distance between crossed populations, but also that reciprocal progenies may differ substantially in their degree of TRD. Chromosomes AL6 and especially AL1 appear to be involved in many single- and two-locus distortions, but the location and source of TRD vary between crosses and between reciprocal progenies. We also found that the majority of single- and two-locus TRD appears to have a gametic, as opposed to zygotic, origin. Thus, while theory on BDM incompatibilities is typically illustrated with derived nuclear alleles proving incompatible in hybrid zygotes, our results suggest a prominent role for distortions emerging before zygote formation.

Indirect evolution of hybrid lethality due to linkage with selected locus in Mimulus guttatus

Most species are superbly and intricately adapted to the environments in which they live. Adaptive evolution by natural selection is the primary force shaping biological diversity. Differences between closely related species in ecologically selected characters such as habitat preference, reproductive timing, courtship behavior, or pollinator attraction may prevent interbreeding in nature, causing reproductive isolation. But does ecological adaptation cause reproductive incompatibilities such as hybrid sterility or lethality? Although several genes causing hybrid incompatibilities have been identified, there is intense debate over whether the genes that contribute to ecological adaptations also cause hybrid incompatibilities. Thirty years ago, a genetic study of local adaptation to copper mine soils in the wildflower Mimulus guttatus identified a locus that appeared to cause copper tolerance and hybrid lethality in crosses to other populations. But do copper tolerance and hybrid lethality have the same molecular genetic basis? Here we show, using high-resolution genome mapping, that copper tolerance and hybrid lethality are not caused by the same gene but are in fact separately controlled by two tightly linked loci. We further show that selection on the copper tolerance locus indirectly caused the hybrid incompatibility allele to go to high frequency in the copper mine population because of hitchhiking. Our results provide a new twist on Darwin's original supposition that hybrid incompatibilities evolve as an incidental by-product of ordinary adaptation to the environment.

A RAD-based linkage map and comparative genomics in the gudgeons (genus Gnathopogon, Cyprinidae)

Background

The construction of linkage maps is a first step in exploring the genetic basis for adaptive phenotypic divergence in closely related species by quantitative trait locus (QTL) analysis. Linkage maps are also useful for comparative genomics in non-model organisms. Advances in genomics technologies make it more feasible than ever to study the genetics of adaptation in natural populations. Restriction-site associated DNA (RAD) sequencing in next-generation sequencers facilitates the development of many genetic markers and genotyping. We aimed to construct a linkage map of the gudgeons of the genus Gnathopogon (Cyprinidae) for comparative genomics with the zebrafish Danio rerio (a member of the same family as gudgeons) and for the future QTL analysis of the genetic architecture underlying adaptive phenotypic evolution of Gnathopogon.

Results

We constructed the first genetic linkage map of Gnathopogon using a 198 F₂ interspecific cross between two closely related species in Japan: river-dwelling Gnathopogon elongatus and lake-dwelling Gnathopogon caerulescens. Based on 1,622 RAD-tag markers, a linkage map spanning 1,390.9 cM with 25 linkage groups and an average marker interval of 0.87 cM was constructed. We also identified a region involving female-specific transmission ratio distortion (TRD). Synteny and collinearity were extensively conserved between Gnathopogon and zebrafish.

Conclusions

The dense SNP-based linkage map presented here provides a basis for future QTL analysis. It will also be useful for transferring genomic information from a “traditional” model fish species, zebrafish, to screen candidate genes underlying ecologically important traits of the gudgeons.

Genomic collinearity and the genetic architecture of floral differences between the homoploid hybrid species Iris nelsonii and one of its progenitors, Iris hexagona

Hybrid speciation represents a relatively rapid form of diversification. Early models of homoploid hybrid speciation suggested that reproductive isolation between the hybrid species and progenitors primarily resulted from karyotypic differences between the species. However, genic incompatibilities and ecological divergence may also be responsible for isolation. Iris nelsonii is an example of a homoploid hybrid species that is likely isolated from its progenitors primarily by strong prezygotic isolation, including habitat divergence, floral isolation and post-pollination prezygotic barriers. Here, we used linkage mapping and quantitative trait locus (QTL) mapping approaches to investigate genomic collinearity and the genetic architecture of floral differences between I. nelsonii and one of its progenitor species I. hexagona. The linkage map produced from this cross is highly collinear with another linkage map produced between I. fulva and I. brevicaulis (the two other species shown to have contributed to the genomic makeup of I. nelsonii), suggesting that karyotypic differences do not contribute substantially to isolation in this homoploid hybrid species. Similar to other studies of the genetic architecture of floral characteristics, at least one QTL was found that explained >20% variance in each color trait, while minor QTLs were detected for each morphological trait. These QTLs will serve as hypotheses for regions under selection by pollinators.

Patterns of shoot architecture in locally adapted populations are linked to intraspecific differences in gene regulation

• Shoot architecture, including the number and location of branches, is a crucial aspect of plant function, morphological diversification, life history evolution and crop domestication. • Genes controlling shoot architecture are well characterized in, and largely conserved across, model flowering plant species. The role of these genes in the evolution of morphological diversity in natural populations, however, has not been explored. • We identify axillary meristem outgrowth as a primary driver of divergent branch number and life histories in two locally adapted populations of the monkeyflower, Mimulus guttatus. • Furthermore, we show that MORE AXILLARY GROWTH (MAX) gene expression strongly correlates with natural variation in branch outgrowth in this species, linking modification of the MAX-dependent pathway to the evolutionary diversification of shoot architecture.

Development of SSR markers and construction of a linkage map in jute

Jute is an important natural fibre crop, which is only second to cotton in its importance at the global level. It is mostly grown in Indian subcontinent and has been recently used for the development of genomics resources.We recently initiated a programme to develop simple sequence repeat markers and reported a set of 2469 SSR that were developed using four SSR-enriched libraries (Mir et al. 2009). In this communication, we report an additional set of 607 novel SSR in 393 SSR containing sequences. However, primers could be designed for only 417 potentially useful SSR. Polymorphism survey was carried out for 374 primer pairs using two parental genotypes (JRO 524 and PPO4) of a mapping population developed for fibre fineness; only 66 SSR were polymorphic. Owing to a low level of polymorphism between the parental genotypes and a high degree of segregation distortion in recombinant inbred lines, genotypic data of only 53 polymorphic SSR on the mapping population consisting of 120 RIL could be used for the construction of a linkage map; 36 SSR loci were mapped on six linkage groups that covered a total genetic distance of 784.3 cM. Hopefully, this map will be enriched with more SSR loci in future and will prove useful for identification of quantitative trait loci/genes for molecular breeding involving improvement of fibre fineness and other related traits in jute.

Genome scans for transmission ratio distortion regions in mice

Transmission ratio distortion (TRD) is the departure from the expected genotypic frequencies under Mendelian inheritance. This departure can be due to multiple physiological mechanisms during gametogenesis, fertilization, fetal and embryonic development, and early neonatal life. Although a few TRD loci have been reported in mouse, inheritance patterns have never been evaluated for TRD. In this article, we developed a Bayesian binomial model accounting for additive and dominant deviation TRD mechanisms. Moreover, this model was used to perform genome-wide scans for TRD quantitative trait loci (QTL) on six F2 mouse crosses involving between 296 and 541 mice and between 72 and 1854 genetic markers. Statistical significance of each model was checked at each genetic marker with Bayes factors. Genome scans revealed overdominance TRD QTL located in mouse chromosomes 1, 2, 12, 13, and 14 and additive TRD QTL in mouse chromosomes 2, 3, and 15, although these results did not replicate across mouse crosses. This research contributes new statistical tools for the analysis of specific genetic patterns involved in TRD in F2 populations, our results suggesting a relevant incidence of TRD phenomena in mouse with important implications for both statistical analyses and biological research.

Rapid evolution caused by pollinator loss in Mimulus guttatus

Anthropogenic perturbations including habitat loss and emerging disease are changing pollinator communities and generating novel selection pressures on plant populations. Disruption of plant-pollinator relationships is predicted to cause plant mating system evolution, although this process has not been directly observed. This study demonstrates the immediate evolutionary effects of pollinator loss within experimental populations of a predominately outcrossing wildflower. Initially equivalent populations evolved for five generations within two pollination treatments: abundant bumblebee pollinators versus no pollinators. The populations without pollinators suffered greatly reduced fitness in early generations but rebounded as they evolved an improved ability to self-fertilize. All populations diverged in floral, developmental, and life-history traits, but only a subset of characters showed clear association with pollination treatment. Pronounced treatment effects were noted for anther-stigma separation and autogamous seed set. Dramatic allele frequency changes at two chromosomal polymorphisms occurred in the no pollinator populations, explaining a large fraction of divergence in pollen viability. The pattern of phenotypic and genetic changes in this experiment favors a sequential model for the evolution of the multitrait "selfing syndrome" observed throughout angiosperms.

Node-specific branching and heterochronic changes underlie population-level differences in Mimulus guttatus (Phrymaceae) shoot architecture

PREMISE OF THE STUDY

Shoot architecture is a fundamentally developmental aspect of plant biology with implications for plant form, function, reproduction, and life history evolution. Mimulus guttatus is morphologically diverse and becoming a model for evolutionary biology. Shoot architecture, however, has never been studied from a developmental perspective in M. guttatus.

METHODS

We examined the development of branches and flowers in plants from two locally adapted populations of M. guttatus with contrasting flowering times, life histories, and branch numbers. We planted second-generation seed in growth chambers to control for maternal and environmental effects.

KEY RESULTS

Most branches occurred at nodes one and two of the main axis. Onset of branching occurred earlier and at a greater frequency in perennials than in annuals. In perennials, almost all flowers occurred at the fifth or more distal nodes. In annuals, most flowers occurred at the third and more distal nodes. Accessory axillary meristems and higher-order branching did not influence shoot architecture.

CONCLUSIONS

We found no evidence for trade-offs between flowers and branches because axillary meristem number was not limiting: a large number of meristems remained quiescent. If, however, quiescence is a component of meristem allocation strategy, then meristems may be limited despite presence of quiescent meristems. At the two basalmost nodes, branch number was determined by mechanisms governing either meristem initiation or outgrowth, rather than flowering vs. branching. At the third and more distal nodes, heterochronic processes contributed to flowering time and branch number differences between populations.

Relationship between transmission ratio distortion and genetic divergence in intraspecific rice crosses

The strength of reproductive isolation often correlates positively with parental divergence in both animals and plants. Here, we assess the relationship between transmission ratio distortion (TRD) of marker loci and parental divergence in intraspecific rice (Oryza sativa L.) crosses. We produced 10 diverse F(2) populations by crossing a temperate japonica reference accession with each of 10 donor accessions belonging to subpopulations different from the reference accession, and then genotyped the F(2) populations using molecular markers distributed across the entire genome. Significant TRDs (α = 0.05) were detected in 9 of the 10 F(2) populations. TRD regions on chromosomes 3 and 6 were common to several populations; in contrast, other TRD regions were unique to a single population, indicating the diversification of genomic location of TRDs among the populations. The level of TRD (estimated from the overall number and magnitude of TRDs) was significantly correlated with the genetic distance between the donor accessions and the reference accession. Our results suggest that in intraspecific rice crosses, parental divergence may result in diversification of the TRD pattern, followed by an increase in the level of TRD.

Genetic diversity of microsatellite loci in hierarchically structured populations

Microsatellite loci are widely used for investigating patterns of genetic variation within and among populations. Those patterns are in turn determined by population sizes, migration rates, and mutation rates. We provide exact expressions for the first two moments of the allele frequency distribution in a stochastic model appropriate for studying microsatellite evolution with migration, mutation, and drift under the assumption that the range of allele sizes is bounded. Using these results, we study the behavior of several measures related to Wright's F(ST), including Slatkin's R(ST). Our analytical approximations for F(ST) and R(ST) show that familiar relationships between N(e)m and F(ST) or R(ST) hold when the migration and mutation rates are small. Using the exact expressions for F(ST) and R(ST), our numerical results show that, when the migration and mutation rates are large, these relationships no longer hold. Our numerical results also show that the diversity measures most closely related to F(ST) depend on mutation rates, mutational models (stepwise versus two-phase), migration rates, and population sizes. Surprisingly, R(ST) is relatively insensitive to the mutation rates and mutational models. The differing behaviors of R(ST) and F(ST) suggest that properties of the among-population distribution of allele frequencies may allow the roles of mutation and migration in producing patterns of diversity to be distinguished, a topic of continuing investigation.

Prevalence of segregation distortion in diploid alfalfa and its implications for genetics and breeding applications

Segregation distortion (SD) is often observed in plant populations; its presence can affect mapping and breeding applications. To investigate the prevalence of SD in diploid alfalfa (Medicago sativa L.), we developed two unrelated segregating F(1) populations and one F(2) population. We genotyped all populations with SSR markers and assessed SD at each locus in each population. The three maps were syntenic and largely colinear with the Medicago truncatula genome sequence. We found genotypic SD for 24 and 34% of markers in the F(1) populations and 68% of markers in the F(2) population; distorted markers were identified on every linkage group. The smaller percentage of genotypic SD in the F(1) populations could be because they were non-inbred and/or due to non-fully informative markers. For the F(2) population, 60 of 90 mapped markers were distorted, and they clustered into eight segregation distortion regions (SDR). Most SDR identified in the F(1) populations were also identified in the F(2) population. Genotypic SD was primarily due to zygotic rather than allelic distortion, suggesting zygotic not gametic selection is the main cause of SD. On the F(2) linkage map, distorted markers in all SDR except two showed heterozygote excess. The severe SD in the F(2) population likely biased genetic distances among markers and possibly also marker ordering and could affect QTL mapping of agronomic traits. To reduce the effects of SD and non-fully informative markers, we suggest constructing linkage maps and conducting QTL mapping in advanced generation populations.

Caenorhabditis briggsae recombinant inbred line genotypes reveal inter-strain incompatibility and the evolution of recombination

The nematode Caenorhabditis briggsae is an emerging model organism that allows evolutionary comparisons with C. elegans and exploration of its own unique biological attributes. To produce a high-resolution C. briggsae recombination map, recombinant inbred lines were generated from reciprocal crosses between two strains and genotyped at over 1,000 loci. A second set of recombinant inbred lines involving a third strain was also genotyped at lower resolution. The resulting recombination maps exhibit discrete domains of high and low recombination, as in C. elegans, indicating these are a general feature of Caenorhabditis species. The proportion of a chromosome's physical size occupied by the central, low-recombination domain is highly correlated between species. However, the C. briggsae intra-species comparison reveals striking variation in the distribution of recombination between domains. Hybrid lines made with the more divergent pair of strains also exhibit pervasive marker transmission ratio distortion, evidence of selection acting on hybrid genotypes. The strongest effect, on chromosome III, is explained by a developmental delay phenotype exhibited by some hybrid F2 animals. In addition, on chromosomes IV and V, cross direction-specific biases towards one parental genotype suggest the existence of cytonuclear epistatic interactions. These interactions are discussed in relation to surprising mitochondrial genome polymorphism in C. briggsae, evidence that the two strains diverged in allopatry, the potential for local adaptation, and the evolution of Dobzhansky-Muller incompatibilities. The genetic and genomic resources resulting from this work will support future efforts to understand inter-strain divergence as well as facilitate studies of gene function, natural variation, and the evolution of recombination in Caenorhabditis nematodes.

Interpopulation hybridization results in widespread viability selection across the genome in Tigriopus californicus

Background

Genetic interactions within hybrids influence their overall fitness. Understanding the details of these interactions can improve our understanding of speciation. One experimental approach is to investigate deviations from Mendelian expectations (segregation distortion) in the inheritance of mapped genetic markers. In this study, we used the copepod Tigriopus californicus, a species which exhibits high genetic divergence between populations and a general pattern of reduced fitness in F2 interpopulation hybrids. Previous studies have implicated both nuclear-cytoplasmic and nuclear-nuclear interactions in causing this fitness reduction. We identified and mapped population-diagnostic single nucleotide polymorphisms (SNPs) and used these to examine segregation distortion across the genome within F2 hybrids.

Results

We generated a linkage map which included 45 newly elucidated SNPs and 8 population-diagnostic microsatellites used in previous studies. The map, the first available for the Copepoda, was estimated to cover 75% of the genome and included markers on all 12 T. californicus chromosomes. We observed little segregation distortion in newly hatched F2 hybrid larvae (fewer than 10% of markers at p < 0.05), but strikingly higher distortion in F2 hybrid adult males (45% of markers at p < 0.05). Hence, segregation distortion was primarily caused by selection against particular genetic combinations which acted between hatching and maturity. Distorted markers were not distributed randomly across the genome but clustered on particular chromosomes. In contrast to other studies in this species we found little evidence for cytonuclear coadaptation. Instead, different linkage groups exhibited markedly different patterns of distortion, which appear to have been influenced by nuclear-nuclear epistatic interactions and may also reflect genetic load carried within the parental lines.

Conclusion

Adult male F2 hybrids between two populations of T. californius exhibit dramatic segregation distortion across the genome. Distorted loci are clustered within specific linkage groups, and the direction of distortion differs between chromosomes. This segregation distortion is due to selection acting between hatching and adulthood.

A molecular recombination map of Antirrhinum majus

BACKGROUND

Genetic recombination maps provide important frameworks for comparative genomics, identifying gene functions, assembling genome sequences and for breeding. The molecular recombination map currently available for the model eudicot Antirrhinum majus is the result of a cross with Antirrhinum molle, limiting its usefulness within A. majus.

RESULTS

We created a molecular linkage map of A. majus based on segregation of markers in the F2 population of two inbred lab strains of A. majus. The resulting map consisted of over 300 markers in eight linkage groups, which could be aligned with a classical recombination map and the A. majus karyotype. The distribution of recombination frequencies and distorted transmission of parental alleles differed from those of a previous inter-species hybrid. The differences varied in magnitude and direction between chromosomes, suggesting that they had multiple causes. The map, which covered an estimated of 95% of the genome with an average interval of 2 cM, was used to analyze the distribution of a newly discovered family of MITE transposons and tested for its utility in positioning seven mutations that affect aspects of plant size.

CONCLUSIONS

The current map has an estimated interval of 1.28 Mb between markers. It shows a lower level of transmission ratio distortion and a longer length than the previous inter-species map, making it potentially more useful. The molecular recombination map further indicates that the IDLE MITE transposons are distributed throughout the genome and are relatively stable. The map proved effective in mapping classical morphological mutations of A. majus.

Extensive recombination rate variation in the house mouse species complex inferred from genetic linkage maps

The rate of recombination is a key genomic parameter that displays considerable variation among taxa. Species comparisons have demonstrated that the rate of evolution in recombination rate is strongly dependent on the physical scale of measurement. Individual recombination hotspots are poorly conserved among closely related taxa, whereas genomic-scale recombination rate variation bears a strong signature of phylogenetic history. In contrast, the mode and tempo of evolution in recombination rates measured on intermediate physical scales is poorly understood. Here, we conduct a detailed statistical comparison between two whole-genome F₂ genetic linkage maps constructed from experimental intercrosses between closely related house mouse subspecies (Mus musculus). Our two maps profile a common wild-derived inbred strain of M. m. domesticus crossed to distinct wild-derived inbred strains representative of two other house mouse subspecies, M. m. castaneus and M. m. musculus. We identify numerous orthologous genomic regions with significant map length differences between these two crosses. Because the genomes of these recently diverged house mice are highly collinear, observed differences in map length (centimorgans) are suggestive of variation in broadscale recombination rate (centimorgans per megabase) within M. musculus. Collectively, these divergent intervals span 19% of the house mouse genome, disproportionately aggregating on the X chromosome. In addition, we uncover strong statistical evidence for a large effect, sex-linked, site-specific modifier of recombination rate segregating within M. musculus. Our findings reveal considerable variation in the megabase-scale recombination landscape among recently diverged taxa and underscore the continued importance of genetic linkage maps in the post-genome era.

Progress and Promise in using Arabidopsis to Study Adaptation, Divergence, and Speciation

Fundamental questions remain to be answered on how lineages split and new species form. The Arabidopsis genus, with several increasingly well characterized species closely related to the model system A. thaliana, provides a rare opportunity to address key questions in speciation research. Arabidopsis species, and in some cases populations within a species, vary considerably in their habitat preferences, adaptations to local environments, mating system, life history strategy, genome structure and chromosome number. These differences provide numerous open doors for understanding the role these factors play in population divergence and how they may cause barriers to arise among nascent species. Molecular tools available in A. thaliana are widely applicable to its relatives, and together with modern comparative genomic approaches they will provide new and increasingly mechanistic insights into the processes underpinning lineage divergence and speciation. We will discuss recent progress in understanding the molecular basis of local adaptation, reproductive isolation and genetic incompatibility, focusing on work utilizing the Arabidopsis genus, and will highlight several areas in which additional research will provide meaningful insights into adaptation and speciation processes in this genus.

The first-generation Daphnia magna linkage map

Background

Daphnia magna is a well-established model species in ecotoxicology, ecology and evolution. Several new genomics tools are presently under development for this species; among them, a linkage map is a first requirement for estimating the genetic background of phenotypic traits in quantitative trait loci (QTL) studies and is also very useful in assembling the genome. It also enables comparative studies between D. magna and D. pulex, for which a linkage map already exists.

Results

Here we describe the first genetic linkage map of D. magna. We generated 214 F2 (intercross) clonal lines as the foundation of the linkage analysis. The linkage map itself is based on 109 microsatellite markers, which produced ten major linkage groups ranging in size from 31.1 cM to 288.5 cM. The total size of this linkage map extends to 1211.6 Kosambi cM, and the average interval for the markers within linkage groups is 15.1 cM. The F2 clones can be used to map QTLs for traits that differ between the parental clones. We successfully mapped the location of two loci with infertility alleles, one inherited from the paternal clone (Iinb1) and the other from the maternal clone (Xinb3).

Conclusions

The D. magna linkage map presented here provides extensive coverage of the genome and a given density of markers that enable us to detect QTLs of moderate to strong effects. It is similar in size to the linkage map of D. pulex.

A genetic model for the female sterility barrier between Asian and African cultivated rice species

S(1) is the most important locus acting as a reproductive barrier between Oryza sativa and O. glaberrima. It is a complex locus, with factors that may affect male and female fertility separately. Recently, the component causing the allelic elimination of pollen was fine mapped. However, the position and nature of the component causing female sterility remains unknown. To fine map the factor of the S(1) locus affecting female fertility, we developed a mapping approach based on the evaluation of the degree of female transmission ratio distortion (fTRD) of markers. Through implementing this methodology in four O. sativa x O. glaberrima crosses, the female component of the S(1) locus was mapped into a 27.8-kb (O. sativa) and 50.3-kb (O. glaberrima) region included within the interval bearing the male component of the locus. Moreover, evidence of additional factors interacting with S(1) was also found. In light of the available data, a model where incompatibilities in epistatic interactions between S(1) and the additional factors are the cause of the female sterility barrier between O. sativa and O. glaberrima was developed to explain the female sterility and the TRD mediated by S(1). According to our model, the recombination ratio and allelic combinations between these factors would determine the final allelic frequencies observed for a given cross.

Is local adaptation in Mimulus guttatus caused by trade-offs at individual loci?

Local adaptation is considered to be the result of fitness trade-offs for particular phenotypes across different habitats. However, it is unclear whether such phenotypic trade-offs exist at the level of individual genetic loci. Local adaptation could arise from trade-offs of alternative alleles at individual loci or by complementary sets of loci with different fitness effects of alleles in one habitat but selective neutrality in the alternative habitat. To evaluate the genome-wide basis of local adaptation, we performed a field-based quantitative trait locus (QTL) mapping experiment on recombinant inbred lines (RILs) created from coastal perennial and inland annual races of the yellow monkeyflower (Mimulus guttatus) grown reciprocally in native parental habitats. Overall, we detected 19 QTLs affecting one or more of 16 traits measured in two environments, most of small effect. We identified 15 additional QTL effects at two previously identified candidate QTLs [DIVERGENCE (DIV)]. Significant QTL by environment interactions were detected at the DIV loci, which was largely attributable to genotypic differences at a single field site. We found no detectable evidence for trade-offs for any one component of fitness, although DIV2 showed a trade-off involving different fitness traits between sites, suggesting that local adaptation is largely controlled by non-overlapping loci. This is surprising for an outcrosser, implying that reduced gene flow prevents the evolution of individuals adapted to multiple environments. We also determined that native genotypes were not uniformly adaptive, possibly reflecting fixed mutational load in one of the populations.

The genetics of phenotypic plasticity in plant defense: trichome production in Mimulus guttatus

Insect herbivory is a major driving force of plant evolution. Phenotypic plasticity and developmental variation provide a means for plants to cope with variable herbivory. We characterized the genetics of developmental variation and phenotypic plasticity in trichome density, a putative defensive trait of Mimulus guttatus (yellow monkeyflower). Our results are evaluated in relation to the optimal defense theory, which provides testable predictions for plastic and developmental patterns in defense traits. We found that both developmental stage and simulated insect damage affected trichome production, but in different ways. Plants were more likely to produce at least some trichomes on later leaves than on earlier leaves, regardless of damage. Damage did not affect the average probability of producing trichomes, but it did increase the density of hairs on trichome-positive plants. We mapped trichome quantitative trait loci (QTL) by selectively genotyping a large panel of recombinant inbred lines derived from two highly divergent populations. Several highly pleiotropic QTL influenced multiple aspects of the trichome phenotype (constitutive, developmental, and/or plastic responses). Only one of the QTL influenced trichome induction following damage. In a result that is consistent with a central prediction of optimal defense theory, the high allele at this location was from the ancestral population with low constitutive trichome production.