First genetic maps development and QTL mining in Ranunculus asiaticus L. through ddRADseq.

Persian Buttercup (Ranunculus asiaticus L.; 2x=2n=16; estimated genome size: 7.6Gb) is an ornamental and perennial crop native of Asia Minor and Mediterranean basin, marketed both as cut flower or potted plant. Currently new varieties are developed by selecting plants carrying desirable traits in segregating progenies obtained by controlled mating, which are propagated through rhizomes or micro-propagated in vitro. In order to escalate selection efficiency and respond to market requests, more knowledge of buttercup genetics would facilitate the identification of markers associated with loci and genes controlling key ornamental traits, opening the way for molecular assisted breeding programs. Reduced-representation sequencing (RRS) represents a powerful tool for plant genotyping, especially in case of large genomes such as the one of buttercup, and have been applied for the development of high-density genetic maps in several species. We report on the development of the first molecular-genetic maps in R. asiaticus based on of a two-way pseudo-testcross strategy. A double digest restriction-site associated DNA (ddRAD) approach was applied for genotyping two F1 mapping populations, whose female parents were a genotype of a so called 'ponpon' and of a 'double flower' varieties, while the common male parental ('Cipro') was a genotype producing a simple flower. The ddRAD generated a total of ~2Gb demultiplexed reads, resulting in an average of 8,3M reads per line. The sstacks pipeline was applied for the construction of a mock reference genome based on sequencing data, and SNP markers segregating in only one of the parents were retained for map construction by treating the F1 population as a backcross. The four parental maps (two of the female parents and two of the common male parent) were aligned with 106 common markers and 8 linkage groups were identified, corresponding to the haploid chromosome number of the species. An average of 586 markers were associated with each parental map, with a marker density ranging from 1 marker/cM to 4.4 markers/cM. The developed maps were used for QTL analysis for flower color, leading to the identification of major QTLs for purple pigmentation. These results contribute to dissect on the genetics of Persian buttercup, enabling the development of new approaches for future varietal development.

Genotyping-by-Sequencing and QTL Mapping of Biomass Yield in Two Switchgrass F1 Populations (Lowland x Coastal and Coastal x Upland).

The prevalence of genetic diversity in switchgrass germplasm can be exploited to capture favorable alleles that increase its range of adaptation and biomass yield. The objectives of the study were to analyze the extent of polymorphism and patterns of segregation distortion in two F1 populations and use the linkage maps to locate QTL for biomass yield. We conducted genotyping-by-sequencing on two populations derived from crosses between the allotetraploid lowland genotype AP13 (a selection from "Alamo") and coastal genotype B6 (a selection from PI 422001) with 285 progeny (AB population) and between B6 and the allotetraploid upland VS16 (a selection from "Summer") with 227 progeny (BV population). As predictable from the Euclidean distance between the parents, a higher number of raw variants was discovered in the coastal × upland BV cross (6 M) compared to the lowland × coastal AB cross (2.5 M). The final number of mapped markers was 3,107 on the BV map and 2,410 on the AB map. More segregation distortion of alleles was seen in the AB population, with 75% distorted loci compared to 11% distorted loci in the BV population. The distortion in the AB population was seen across all chromosomes in both the AP13 and B6 maps and likely resulted from zygotic or post-zygotic selection for increased levels of heterozygosity. Our results suggest lower genetic compatibility between the lowland AP13 and the coastal B6 ecotype than between B6 and the upland ecotype VS16. Four biomass QTLs were mapped in the AB population (LG 2N, 6K, 6N, and 8N) and six QTLs in the BV population [LG 1N (2), 8N (2), 9K, and 9N]. The QTL, with the largest and most consistent effect across years, explaining between 8.4 and 11.5% of the variation, was identified on 6N in the AP13 map. The cumulative effect of all the QTLs explained a sizeable portion of the phenotypic variation in both AB and BV populations and the markers associated with them may potentially be used for the marker-assisted improvement of biomass yield. Since switchgrass improvement is based on increasing favorable allele frequencies through recurrent selection, the transmission bias within individuals and loci needs to be considered as this may affect the genetic gain if the favorable alleles are distorted.

HGGA: hierarchical guided genome assembler.

BACKGROUND: De novo genome assembly typically produces a set of contigs instead of the complete genome. Thus additional data such as genetic linkage maps, optical maps, or Hi-C data is needed to resolve the complete structure of the genome. Most of the previous work uses the additional data to order and orient contigs. RESULTS: Here we introduce a framework to guide genome assembly with additional data. Our approach is based on clustering the reads, such that each read in each cluster originates from nearby positions in the genome according to the additional data. These sets are then assembled independently and the resulting contigs are further assembled in a hierarchical manner. We implemented our approach for genetic linkage maps in a tool called HGGA. CONCLUSIONS: Our experiments on simulated and real Pacific Biosciences long reads and genetic linkage maps show that HGGA produces a more contiguous assembly with less contigs and from 1.2 to 9.8 times higher NGA50 or N50 than a plain assembly of the reads and 1.03 to 6.5 times higher NGA50 or N50 than a previous approach integrating genetic linkage maps with contig assembly. Furthermore, also the correctness of the assembly remains similar or improves as compared to an assembly using only the read data.

Genetic mapping of loci involved in oil tocopherol composition control in Russian sunflower (Helianthus annuus L.) lines.

Tocopherols are antioxidants that preserve oil lipids against oxidation and serve as a natural source of vitamin E in the human diet. Compared with other major oilseeds like rapeseed and soybean, sunflower (Helianthus annuus L.) exhibits low phenotypic diversity of tocopherol composition, both in wild and cultivated accessions from germplasm collections. Two major mutations that alter tocopherol composition were identified in genetic collections, and several studies suggested additional loci controlling tocopherol composition, with their expression possibly depending on the genetic background. In the present study, we performed QTL mapping of tocopherol composition in two independent F2 crosses between lines with contrasting tocopherol composition from the Pustovoit All-Russia Research Institute of Oil Crops (VNIIMK) collection. We used genotyping-bysequencing (GBS) to construct single nucleotide polymorphism-based genetic maps, and performed QTL mapping using quantitative and qualitative encoding for phenotypic traits. Our results support the notion that the tocopherol composition in the assessed crosses is controlled by two loci. We additionally selected and validated two single nucleotide polymorphism markers for each cross which could be used for marker-assisted selection.

Large-scale advances in SSR markers with high-throughput sequencing in Euphorbia fischeriana Steud

BACKGROUND: Euphorbia fischeriana Steud is a very important medicinal herb and has significant medical value for healing cancer, edema and tuberculosis in China. The lack of molecular markers for Euphorbia fischeriana Steud is a dominant barrier to genetic research. For the purpose of developing many simple sequence repeat (SSR) molecular markers, we completed transcriptome analysis with the Illumina HiSeq 2000 platform. RESULTS: Approximately 9.1 million clean reads were acquired and then assembled into approximately 186.3 thousand nonredundant unigenes, 53,146 of which were SSR-containing unigenes. A total of 76,193 SSR loci were identified. Of these SSR loci, 28,491 were detected at the terminal position of ESTs, which made it difficult to design SSR primers for these SSR-containing sequences, and the residual SSRs were thus used to design primer pairs. Analyzing the results of these markers revealed that the mononucleotide motif A/T (44,067, 57.83% of all SSRs) was the most abundant, followed by the dinucleotide type AG/CT (9430, 12.38%). Using 100 randomly selected primer pairs, 77 primers were successfully amplified in Euphorbia fischeriana Steud, and 79 were successfully amplified in three other related species. The markers developed displayed relatively high quality and cross-species transferability. CONCLUSIONS: The large number of EST-SSRs exploited successfully in Euphorbia fischeriana Steud for the first time could provide genetic information for research on linkage maps, variety identification, genetic diversity analysis, and molecular marker-assisted breeding.

An Improved Oil Palm Genome Assembly as a Valuable Resource for Crop Improvement and Comparative Genomics in the Arecoideae Subfamily.

Oil palm (Elaeis guineensis Jacq.) is the most traded crop among the economically important palm species. Here, we report an extended version genome of E. guineensis that is 1.2 Gb in length, an improvement of the physical genome coverage to 79% from the previous 43%. The improvement was made by assigning an additional 1968 originally unplaced scaffolds that were available publicly into the physical genome. By integrating three ultra-dense linkage maps and using them to place genomic scaffolds, the 16 pseudomolecules were extended. As we show, the improved genome has enhanced the mapping resolution for genome-wide association studies (GWAS) and permitted further identification of candidate genes/protein-coding regions (CDSs) and any non-coding RNA that may be associated with them for further studies. We then employed the new physical map in a comparative genomics study against two other agriculturally and economically important palm species-date palm (Phoenix dactylifera L.) and coconut palm (Cocos nucifera L.)-confirming the high level of conserved synteny among these palm species. We also used the improved oil palm genome assembly version as a palm genome reference to extend the date palm physical map. The improved genome of oil palm will enable molecular breeding approaches to expedite crop improvement, especially in the largest subfamily of Arecoideae, which consists of 107 species belonging to Arecaceae.

Identification of a candidate gene associated with isoflavone content in soybean seeds using genome-wide association and linkage mapping.

Isoflavone, a secondary metabolite produced by Glycine max (L.) Merr. (soybean), is valuable for human and plant health. The genetic architecture of soybean isoflavone content remains unclear, however, despite several mapping studies. We generated genomic data for 200 soybean cultivars and 150 recombinant inbred lines (RILs) to localize putative loci associated with soybean seed isoflavone content. Using a genome-wide association study (GWAS), we identified 87 single-nucleotide polymorphisms (SNPs) that were significantly associated with isoflavone concentration. Using linkage mapping, we identified 37 quantitative trait loci (QTLs) underlying the content of four isoflavones found in the RILs. A major locus on chromosome 8 (qISO8-1) was co-located by both the GWAS and linkage mapping. qISO8-1 was fine mapped to a 99.5-kb region, flanked by SSR_08_1651 and SSR_08_1656, in a BC2 F5 population. GmMPK1, encoding a mitogen-activated protein kinase, was identified as the causal gene in qISO8-1, and two natural GmMPK1 polymorphisms were significantly associated with isoflavone content. Overexpression of GmMPK1 in soybean hairy roots resulted in increased isoflavone concentrations. Overexpressing GmMPK1 in transgenic soybeans had greater resistance to Phytophthora root rot, suggesting that GmMPK1 might increase soybean resistance to biotic stress by influencing isoflavone content. Our results not only increase our understanding of the genetic architecture of soybean seed isoflavone content, but also provide a framework for the future marker-assisted breeding of high isoflavone content in soybean cultivars.

Major QTLs for Trunk Height and Correlated Agronomic Traits Provide Insights into Multiple Trait Integration in Oil Palm Breeding.

Superior oil yield is always the top priority of the oil palm industry. Short trunk height (THT) and compactness traits have become increasingly important to improve harvesting efficiency since the industry started to suffer yield losses due to labor shortages. Breeding populations with low THT and short frond length (FL) are actually available, such as Dumpy AVROS pisifera (DAV) and Gunung Melayu dura (GM). However, multiple trait stacking still remains a challenge for oil palm breeding, which usually requires 12-20 years to complete a breeding cycle. In this study, yield and height increment in the GM × GM (GM-3341) and the GM × DAV (GM-DAV-3461) crossing programs were evaluated and palms with good yield and smaller height increment were identified. In the GM-3341 family, non-linear THT growth between THT_2008 (seven years old) and THT_2014 (13 years old) was revealed by a moderate correlation, suggesting that inter-palm competition becomes increasingly important. In total, 19 quantitative trait loci (QTLs) for THT_2008 (8), oil per palm (O/P) (7) and FL (4) were localized on the GM-3341 linkage map, with an average mapping interval of 2.01 cM. Three major QTLs for THT_2008, O/P and FL are co-located on chromosome 11 and reflect the correlation of THT_2008 with O/P and FL. Multiple trait selection for high O/P and low THT (based on the cumulative effects of positive alleles per trait) identified one palm from 100 palms, but with a large starting population of 1000-1500 seedling per cross, this low frequency could be easily compensated for during breeding selection.

Construction of three half-sib SSR linkage maps derived from overwintering cultivated rice and segregation distortion loci mapping.

Segregation distortion is a common phenomenon that has been observed in genetics and plant breeding; however, the mechanism of segregation distortion is unknown. In the present study, three half-sib F2 populations derived from three japonica overwinter (perennial) rice varieties (W1, W2, and W3) crossed to the indica rice variety Minghui725 (MH725) were developed to construct three half-sib linkage maps. We established linkage map lengths of 2032.8, 2317.4, and 2108.7 cM with average intervals of 20.1, 20.5, and 19.7 cM using 101, 113, and 107 SSR markers in W1/MH725, W2/MH725, and W3/MH725, respectively. Discrepancies in marker order and genetic linkage distance occurred in the three half-sib linkage maps due to segregation distortion. A total of 88 markers exhibited segregation distortion across the three linkage maps at P < 0.01 level, 42 segregation distortion loci (SDLs) were detected across the three half-sib populations and exhibited variable LOD value that ranged from 3.2 (SDL2f) to 30.1 (SDL5d), and 13 of the 42 SDLs were repeatedly located at the same chromosomal regions of the previously published hybrid sterility quantitative trait loci. Data from this study provide an extensive archive for investigating the genetic characteristic of overwintering cultivated rice and the future exploration and innovation of overwintering rice breeding.

Updates on Genomic Resources in Chickpea for Crop Improvement.

In recent years, rapid advancement has been done in generation of genomic resources for the important legume crop chickpea. Here, we provide an update on important advancements made on availability of genomic resources for this crop. The availability of reference genome and transcriptome sequences, and resequencing of several accessions have enabled the discovery of gene space and molecular markers in chickpea. These resources have helped in elucidating evolutionary relationship and identification of quantitative trait loci for important agronomic traits. Gene expression in different tissues/organs during development and under abiotic/biotic stresses has been interrogated. In addition, single-base resolution DNA methylation patterns in different organs have been analyzed to understand gene regulation. Overall, we provide a consolidated overview of available genomic resources of chickpea that may help in fulfilling the promises for improvement of this important crop.

De Novo Genome Assembly of Populus simonii Further Supports That Populus simonii and Populus trichocarpa Belong to Different Sections.

Populus simonii is an important tree in the genus Populus, widely distributed in the Northern Hemisphere and having a long cultivation history. Although this species has ecologically and economically important values, its genome sequence is currently not available, hindering the development of new varieties with wider adaptive and commercial traits. Here, we report a chromosome-level genome assembly of P. simonii using PacBio long-read sequencing data aided by Illumina paired-end reads and related genetic linkage maps. The assembly is 441.38 Mb in length and contain 686 contigs with a contig N50 of 1.94 Mb. With the linkage maps, 336 contigs were successfully anchored into 19 pseudochromosomes, accounting for 90.2% of the assembled genome size. Genomic integrity assessment showed that 1,347 (97.9%) of the 1,375 genes conserved among all embryophytes can be found in the P. simonii assembly. Genomic repeat analysis revealed that 41.47% of the P. simonii genome is composed of repetitive elements, of which 40.17% contained interspersed repeats. A total of 45,459 genes were predicted from the P. simonii genome sequence and 39,833 (87.6%) of the genes were annotated with one or more related functions. Phylogenetic analysis indicated that P. simonii and Populus trichocarpa should be placed in different sections, contrary to the previous classification according to morphology. The genome assembly not only provides an important genetic resource for the comparative and functional genomics of different Populus species, but also furnishes one of the closest reference sequences for identifying genomic variants in an F1 hybrid population derived by crossing P. simonii with other Populus species.

Kermit: linkage map guided long read assembly.

BACKGROUND: With long reads getting even longer and cheaper, large scale sequencing projects can be accomplished without short reads at an affordable cost. Due to the high error rates and less mature tools, de novo assembly of long reads is still challenging and often results in a large collection of contigs. Dense linkage maps are collections of markers whose location on the genome is approximately known. Therefore they provide long range information that has the potential to greatly aid in de novo assembly. Previously linkage maps have been used to detect misassemblies and to manually order contigs. However, no fully automated tools exist to incorporate linkage maps in assembly but instead large amounts of manual labour is needed to order the contigs into chromosomes. RESULTS: We formulate the genome assembly problem in the presence of linkage maps and present the first method for guided genome assembly using linkage maps. Our method is based on an additional cleaning step added to the assembly. We show that it can simplify the underlying assembly graph, resulting in more contiguous assemblies and reducing the amount of misassemblies when compared to de novo assembly. CONCLUSIONS: We present the first method to integrate linkage maps directly into genome assembly. With a modest increase in runtime, our method improves contiguity and correctness of genome assembly.

A European Whitefish Linkage Map and Its Implications for Understanding Genome-Wide Synteny Between Salmonids Following Whole Genome Duplication.

Genomic datasets continue to increase in number due to the ease of production for a wider selection of species including non-model organisms. For many of these species, especially those with large or polyploid genomes, highly contiguous and well-annotated genomes are still rare due to the complexity and cost involved in their assembly. As a result, a common starting point for genomic work in non-model species is the production of a linkage map. Dense linkage maps facilitate the analysis of genomic data in a variety of ways, from broad scale observations regarding genome structure e.g., chromosome number and type or sex-related structural differences, to fine scale patterns e.g., recombination rate variation and co-localization of differentiated regions. Here we present both sex-averaged and sex-specific linkage maps for Coregonus sp. "Albock", a member of the European whitefish lineage (C. lavaretus spp. complex), containing 5395 single nucleotide polymorphism (SNP) loci across 40 linkage groups to facilitate future investigation into the genomic basis of whitefish adaptation and speciation. The map was produced using restriction-site associated digestion (RAD) sequencing data from two wild-caught parents and 156 F1 offspring. We discuss the differences between our sex-averaged and sex-specific maps and identify genome-wide synteny between C. sp. "Albock" and Atlantic Salmon (Salmo salar), which have diverged following the salmonid-specific whole genome duplication. Our analysis confirms that many patterns of synteny observed between Atlantic Salmon and Oncorhynchus and Salvelinus species are also shared by members of the Coregoninae subfamily. We also show that regions known for their species-specific rediploidization history can pose challenges for synteny identification since these regions have diverged independently in each salmonid species following the salmonid-specific whole genome duplication. The European whitefish map provided here will enable future studies to understand the distribution of loci of interest, e.g., FST outliers, along the whitefish genome as well as assisting with the de novo assembly of a whitefish reference genome.

Accounting for Errors in Low Coverage High-Throughput Sequencing Data When Constructing Genetic Maps Using Biparental Outcrossed Populations.

Next-generation sequencing is an efficient method that allows for substantially more markers than previous technologies, providing opportunities for building high-density genetic linkage maps, which facilitate the development of nonmodel species' genomic assemblies and the investigation of their genes. However, constructing genetic maps using data generated via high-throughput sequencing technology (e.g., genotyping-by-sequencing) is complicated by the presence of sequencing errors and genotyping errors resulting from missing parental alleles due to low sequencing depth. If unaccounted for, these errors lead to inflated genetic maps. In addition, map construction in many species is performed using full-sibling family populations derived from the outcrossing of two individuals, where unknown parental phase and varying segregation types further complicate construction. We present a new methodology for modeling low coverage sequencing data in the construction of genetic linkage maps using full-sibling populations of diploid species, implemented in a package called GUSMap. Our model is based on the Lander-Green hidden Markov model but extended to account for errors present in sequencing data. We were able to obtain accurate estimates of the recombination fractions and overall map distance using GUSMap, while most existing mapping packages produced inflated genetic maps in the presence of errors. Our results demonstrate the feasibility of using low coverage sequencing data to produce genetic maps without requiring extensive filtering of potentially erroneous genotypes, provided that the associated errors are correctly accounted for in the model.

Construction of Agropyron Gaertn. genetic linkage maps using a wheat 660K SNP array reveals a homoeologous relationship with the wheat genome.

Agropyron Gaertn. (P genome) is a wild relative of wheat that harbours many genetic variations that could be used to increase the genetic diversity of wheat. To agronomically transfer important genes from the P genome to a wheat chromosome by induced homoeologous pairing and recombination, it is necessary to determine the chromosomal relationships between Agropyron and wheat. Here, we report using the wheat 660K single nucleotide polymorphism (SNP) array to genotype a segregating Agropyron F1 population derived from an interspecific cross between two cross-pollinated diploid collections 'Z1842' [A. cristatum (L.) Beauv.] (male parent) and 'Z2098' [A. mongolicum Keng] (female parent) and 35 wheat-A. cristatum addition/substitution lines. Genetic linkage maps were constructed using 913 SNP markers distributed among seven linkage groups spanning 839.7 cM. The average distance between adjacent markers was 1.8 cM. The maps identified the homoeologous relationship between the P genome and wheat and revealed that the P and wheat genomes are collinear and relatively conserved. In addition, obvious rearrangements and introgression spread were observed throughout the P genome compared with the wheat genome. Combined with genotyping data, the complete set of wheat-A. cristatum addition/substitution lines was characterized according to their homoeologous relationships. In this study, the homoeologous relationship between the P genome and wheat was identified using genetic linkage maps, and the detection mean for wheat-A. cristatum introgressions might significantly accelerate the introgression of genetic variation from Agropyron into wheat for exploitation in wheat improvement programmes.

A high-resolution reference genetic map positioning 8.8 K genes for the conifer white spruce: structural genomics implications and correspondence with physical distance.

Over the last decade, extensive genetic and genomic resources have been developed for the conifer white spruce (Picea glauca, Pinaceae), which has one of the largest plant genomes (20 Gbp). Draft genome sequences of white spruce and other conifers have recently been produced, but dense genetic maps are needed to comprehend genome macrostructure, delineate regions involved in quantitative traits, complement functional genomic investigations, and assist the assembly of fragmented genomic sequences. A greatly expanded P. glauca composite linkage map was generated from a set of 1976 full-sib progeny, with the positioning of 8793 expressed genes. Regions with significant low or high gene density were identified. Gene family members tended to be mapped on the same chromosomes, with tandemly arrayed genes significantly biased towards specific functional classes. The map was integrated with transcriptome data surveyed across eight tissues. In total, 69 clusters of co-expressed and co-localising genes were identified. A high level of synteny was found with pine genetic maps, which should facilitate the transfer of structural information in the Pinaceae. Although the current white spruce genome sequence remains highly fragmented, dozens of scaffolds encompassing more than one mapped gene were identified. From these, the relationship between genetic and physical distances was examined and the genome-wide recombination rate was found to be much smaller than most estimates reported for angiosperm genomes. This gene linkage map shall assist the large-scale assembly of the next-generation white spruce genome sequence and provide a reference resource for the conifer genomics community.

QTL variations for growth-related traits in eight distinct families of common carp (Cyprinus carpio).

BACKGROUND: Comparing QTL analyses of multiple pair-mating families can provide a better understanding of important allelic variations and distributions. However, most QTL mapping studies in common carp have been based on analyses of individual families. In order to improve our understanding of heredity and variation of QTLs in different families and identify important QTLs, we performed QTL analysis of growth-related traits in multiple segregating families. RESULTS: We completed a genome scan for QTLs that affect body weight (BW), total length (TL), and body thickness (BT) of 522 individuals from eight full-sib families using 250 microsatellites evenly distributed across 50 chromosomes. Sib-pair and half-sib model mapping identified 165 QTLs on 30 linkage groups. Among them, 10 (genome-wide P <0.01 or P < 0.05) and 28 (chromosome-wide P < 0.01) QTLs exhibited significant evidence of linkage, while the remaining 127 exhibited a suggestive effect on the above three traits at a chromosome-wide (P < 0.05) level. Multiple QTLs obtained from different families affect BW, TL, and BT and locate at close or identical positions. It suggests that same genetic factors may control variability in these traits. Furthermore, the results of the comparative QTL analysis of multiple families showed that one QTL was common in four of the eight families, nine QTLs were detected in three of the eight families, and 26 QTLs were found common to two of the eight families. These common QTLs are valuable candidates in marker-assisted selection. CONCLUSION: A large number of QTLs were detected in the common carp genome and associated with growth-related traits. Some of the QTLs of different growth-related traits were identified at similar chromosomal regions, suggesting a role for pleiotropy and/or tight linkage and demonstrating a common genetic basis of growth trait variations. The results have set up an example for comparing QTLs in common carp and provided insights into variations in the identified QTLs affecting body growth. Discovery of these common QTLs between families and growth-related traits represents an important step towards understanding of quantitative genetic variation in common carp.

High-density linkage maps fail to detect any genetic component to sex determination in a Rana temporaria family.

Sex chromosome differentiation in Rana temporaria varies strikingly among populations or families: whereas some males display well-differentiated Y haplotypes at microsatellite markers on linkage group 2 (LG2), others are genetically undistinguishable from females. We analysed with RADseq markers one family from a Swiss lowland population with no differentiated sex chromosomes, and where sibship analyses had failed to detect any association between the phenotypic sex of progeny and parental haplotypes. Offspring were reared in a common tank in outdoor conditions and sexed at the froglet stage. We could map a total of 2177 SNPs (1123 in the mother, 1054 in the father), recovering in both adults 13 linkage groups (= chromosome pairs) that were strongly syntenic to Xenopus tropicalis despite > 200 My divergence. Sexes differed strikingly in the localization of crossovers, which were uniformly distributed in the female but limited to chromosome ends in the male. None of the 2177 markers showed significant association with offspring sex. Considering the very high power of our analysis, we conclude that sex determination was not genetic in this family; which factors determined sex remain to be investigated.

CHROMOSOMAL REPATTERNING AND LINKAGE GROUP CONSERVATION IN MOSQUITO KARYOTYPIC EVOLUTION.

Chromosome number and morphology in mosquitoes is remarkably uniform: virtually all mosquitoes have a diploid chromosome number of six (2N = 6), and their chromosomes are invariably metacentric or submetacentric. Numerical changes obviously have not been important in mosquito chromosomal evolution, and because of the morphological similarity of their chromosomes, it appears that structural changes have played little or no role in mosquito karyotypic evolution. The goal of the present study was to identify the types and relative numbers of chromosomal changes in mosquito evolution and to extend the comparison where possible to the higher diptera. To do this, we compared the enzyme linkage maps of six species of Aedes to each other and to enzyme maps of seven other mosquito species and to Drosophila melanogaster. Our results indicate that Aedes chromosomes have been modified by inversions, most which were paracentric, and by translocations, most which were Robertsonian. Intrageneric comparison of Aedes enzyme maps also revealed groups of linked enzyme loci whose integrity has been maintained throughout Aedes evolution (conserved linkages/syntenies). Intergeneric comparisons of Aedes enzyme maps with those of species in the genera Culex, Anopheles, and Toxorhynchites disclosed conserved associations of enzyme loci between mosquito genera. These findings lead us to postulate that the ancestral mosquito karyotype consisted of six chromosomal elements which, other than being combined in different ways in various mosquito groups, have remained essentially intact during mosquito evolution. Furthermore, the identification of groups of linked enzyme loci common to mosquitoes and to D. melanogaster indicates that linkage group conservation may characterize the karyotypic evolution of all dipteran insects.