Cyclical regression covariates remove the major confounding effect of cyclical developmental gene expression with strain-specific drug response in the malaria parasite Plasmodium falciparum.

BACKGROUND: The cyclical nature of gene expression in the intraerythrocytic development cycle (IDC) of the malaria parasite, Plasmodium falciparum, confounds the accurate detection of specific transcriptional differences, e.g. as provoked by the development of drug resistance. In lab-based studies, P. falciparum cultures are synchronized to remove this confounding factor, but the rapid detection of emerging resistance to artemisinin therapies requires rapid analysis of transcriptomes extracted directly from clinical samples. Here we propose the use of cyclical regression covariates (CRC) to eliminate the major confounding effect of developmentally driven transcriptional changes in clinical samples. We show that elimination of this confounding factor reduces both Type I and Type II errors and demonstrate the effectiveness of this approach using a published dataset of 1043 transcriptomes extracted directly from patient blood samples with different patient clearance times after treatment with artemisinin. RESULTS: We apply this method to two publicly available datasets and demonstrate its ability to reduce the confounding of differences in transcript levels due to misaligned intraerythrocytic development time. Adjusting the clinical 1043 transcriptomes dataset with CRC results in detection of fewer functional categories than previously reported from the same data set adjusted using other methods. We also detect mostly the same functional categories, but observe fewer genes within these categories. Finally, the CRC method identifies genes in a functional category that was absent from the results when the dataset was adjusted using other methods. Analysis of differential gene expression in the clinical data samples that vary broadly for developmental stage resulted in the detection of far fewer transcripts in fewer functional categories while, at the same time, identifying genes in two functional categories not present in the unadjusted data analysis. These differences are consistent with the expectation that CRC reduces both false positives and false negatives with the largest effect on datasets from samples with greater variance in developmental stage. CONCLUSIONS: Cyclical regression covariates have immediate application to parasite transcriptome sequencing directly from clinical blood samples and to cost-constrained in vitro experiments.

Transcriptome profiles of Quercus rubra responding to increased O3 stress.

BACKGROUND: Climate plays an essential role in forest health, and climate change may increase forest productivity losses due to abiotic and biotic stress. Increased temperature leads to the increased formation of ozone (O3). Ozone is formed by the interaction of sunlight, molecular oxygen and by the reactions of chemicals commonly found in industrial and automobile emissions such as nitrogen oxides and volatile organic compounds. Although it is well known that productivity of Northern red oak (Quercus rubra) (NRO), an ecologically and economically important species in the forests of eastern North America, is reduced by exposure to O3, limited information is available on its responses to exogenous stimuli at the level of gene expression. RESULTS: RNA sequencing yielded more than 323 million high-quality raw sequence reads. De novo assembly generated 52,662 unigenes, of which more than 42,000 sequences could be annotated through homology-based searches. A total of 4140 differential expressed genes (DEGs) were detected in response to O3 stress, as compared to their respective controls. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of the O3-response DEGs revealed perturbation of several biological pathways including energy, lipid, amino acid, carbohydrate and terpenoid metabolism as well as plant-pathogen interaction. CONCLUSION: This study provides the first reference transcriptome for NRO and initial insights into the genomic responses of NRO to O3. Gene expression profiling reveals altered primary and secondary metabolism of NRO seedlings, including known defense responses such as terpenoid biosynthesis.

High-quality genetic mapping with ddRADseq in the non-model tree Quercus rubra.

BACKGROUND: Restriction site associated DNA sequencing (RADseq) has the potential to be a broadly applicable, low-cost approach for high-quality genetic linkage mapping in forest trees lacking a reference genome. The statistical inference of linear order must be as accurate as possible for the correct ordering of sequence scaffolds and contigs to chromosomal locations. Accurate maps also facilitate the discovery of chromosome segments containing allelic variants conferring resistance to the biotic and abiotic stresses that threaten forest trees worldwide. We used ddRADseq for genetic mapping in the tree Quercus rubra, with an approach optimized to produce a high-quality map. Our study design also enabled us to model the results we would have obtained with less depth of coverage. RESULTS: Our sequencing design produced a high sequencing depth in the parents (248×) and a moderate sequencing depth (15×) in the progeny. The digital normalization method of generating a de novo reference and the SAMtools SNP variant caller yielded the most SNP calls (78,725). The major drivers of map inflation were multiple SNPs located within the same sequence (77% of SNPs called). The highest quality map was generated with a low level of missing data (5%) and a genome-wide threshold of 0.025 for deviation from Mendelian expectation. The final map included 849 SNP markers (1.8% of the 78,725 SNPs called). Downsampling the individual FASTQ files to model lower depth of coverage revealed that sequencing the progeny using 96 samples per lane would have yielded too few SNP markers to generate a map, even if we had sequenced the parents at depth 248×. CONCLUSIONS: The ddRADseq technology produced enough high-quality SNP markers to make a moderately dense, high-quality map. The success of this project was due to high depth of coverage of the parents, moderate depth of coverage of the progeny, a good framework map, an optimized bioinformatics pipeline, and rigorous premapping filters. The ddRADseq approach is useful for the construction of high-quality genetic maps in organisms lacking a reference genome if the parents and progeny are sequenced at sufficient depth. Technical improvements in reduced representation sequencing (RRS) approaches are needed to reduce the amount of missing data.

Balancing selection maintains diversity in a cold tolerance gene in broadly distributed live oaks.

Cold poses major physiological challenges to plants, especially long-lived trees. In trees occurring along variable temperature clines, the expected direction and consequences of selection on cold acclimation ability and freezing tolerance are not straightforward. Here we estimated selection in cold acclimation genes at two evolutionary timescales in all seven species of the American live oaks (Quercus subsection Virentes). Two cold response candidate genes were chosen: ICE1, a key gene in the cold acclimation pathway, and HOS1, which modulates cold response by negatively regulating ICE1. Two housekeeping genes, GAPDB and CHR11, were also analyzed. At the shallow evolutionary timescale, we demonstrate that HOS1 experienced recent balancing selection in the two most broadly distributed species, Q. virginiana and Q. oleoides. At a deeper evolutionary scale, a codon-based model of evolution revealed the signature of negative selection in ICE1. In contrast, three positively selected codons have been identified in HOS1, possibly a signature of the diversification of Virentes into warmer climates from a freezing adapted lineage of oaks. Our findings indicate that evolution has favored diversity in cold tolerance modulation through balancing selection in HOS1 while maintaining core cold acclimation ability, as evidenced by purifying selection in ICE1.

Back to America: tracking the origin of European introduced populations of Quercus rubra L.

Quercus rubra has been introduced in Europe since the end of the 17th century. It is widely distributed today across this continent and considered invasive in some countries. Here, we investigated the distribution of genetic diversity of both native and introduced populations with the aim of tracing the origin of introduced populations. A large sampling of 883 individuals from 73 native and 38 European locations were genotyped at 69 SNPs. In the natural range, we found a continuous geographic gradient of variation with a predominant latitudinal component. We explored the existence of ancestral populations by performing Bayesian clustering analysis and found support for two or three ancestral genetic clusters. Approximate Bayesian Computations analyses based on these two or three clusters support recent extensive secondary contacts between them, suggesting that present-day continuous genetic variation resulted from recent admixture. In the introduced range, one main genetic cluster was not recovered in Europe, suggesting that source populations were preferentially located in the northern part of the natural distribution. However, our results cannot refute the introduction of populations from the southern states that did not survive in Europe.

Influence of mosquito genotype on transcriptional response to dengue virus infection.

The mosquito Aedes aegypti is the principal vector that transmits dengue virus (DENV) to humans. The primary factors that trigger a susceptible or refractory interaction of A. aegypti with DENV are not well understood. In this study, our aim is to characterize the influence of vector genotype on differential gene expression of susceptible vs. refractory A. aegypti strains to DENV infection. To accomplish that, we identified differential expression of a set of complementary DNAs (cDNAs; n = 9,504) of the D2S3 (susceptible) and Moyo-D (refractory) strains of A. aegypti to DENV serotype 2 (JAM1409) and compared these results to the differential expression of cDNAs in a different susceptible vector genotype (Moyo-S) relative to the same refractory genotype (Moyo-D) identified from our previous study. We observed that, although the number of differentially expressed transcripts (DETs) was similar in both the studies, about ~95% of the DETs were distinct between Moyo-D/D2S3 vs. Moyo-D/Moyo-S. This suggested that A. aegypti response, to infection of a given genotype of dengue, is largely dependent upon the vector genotype. However, we observed a set of common DETs among the vector strains that were associated with predicted functions such as endocytosis, regulation of autophagy, peroxisome, and lipid metabolism that may be relatively universal in conferring mosquito response to DENV infection.

Gene flow from Fraxinus cultivars into natural stands of Fraxinus pennsylvanica occurs range-wide, is regionally extensive, and is associated with a loss of allele richness.

In North America, a comparatively small number of Fraxinus (ash) cultivars were planted in large numbers in both urban and rural environments across the entire range of Fraxinus pennsylvanica Marsh (green ash) over the last 80 years. Undetected cultivar gene flow, if extensive, could significantly lower genetic diversity within populations, suppress differentiation between populations, generate interspecific admixture not driven by long-standing natural processes, and affect the impact of abiotic and biotic threats. In this investigation we generated the first range-wide genetic assessment of F. pennsylvanica to detect the extent of cultivar gene flow into natural stands. We used 16 EST-SSR markers to genotype 48 naturally regenerated populations of F. pennsylvanica distributed across the native range (1291 trees), 19 F. pennsylvanica cultivars, and one F. americana L. (white ash) cultivar to detect cultivar propagule dispersal into these populations. We detected first generation cultivar parentage with high confidence in 171 individuals in 34 of the 48 populations and extensive cultivar parentage (23-50%) in eight populations. The incidence of cultivar parentage was negatively associated with allele richness (R2 = 0.151, p = 0.006). The evidence for a locally high frequency of cultivar propagule dispersal and the interspecific admixture in eastern populations will inform Fraxinus gene pool conservation strategies and guide the selection of individuals for breeding programs focused on increasing resistance to the emerald ash borer (Agrilus planipennis Fairmaire), an existential threat to the Fraxinus species of North America.

Conserving a threatened North American walnut: a chromosome-scale reference genome for butternut (Juglans cinerea).

With the advent of affordable and more accurate third-generation sequencing technologies, and the associated bioinformatic tools, it is now possible to sequence, assemble, and annotate more species of conservation concern than ever before. Juglans cinerea, commonly known as butternut or white walnut, is a member of the walnut family, native to the Eastern United States and Southeastern Canada. The species is currently listed as Endangered on the IUCN Red List due to decline from an invasive fungus known as Ophiognomonia clavigignenti-juglandacearum (Oc-j) that causes butternut canker. Oc-j creates visible sores on the trunks of the tree which essentially starves and slowly kills the tree. Natural resistance to this pathogen is rare. Conserving butternut is of utmost priority due to its critical ecosystem role and cultural significance. As part of an integrated undergraduate and graduate student training program in biodiversity and conservation genomics, the first reference genome for Juglans cinerea is described here. This chromosome-scale 539 Mb assembly was generated from over 100 × coverage of Oxford Nanopore long reads and scaffolded with the Juglans mandshurica genome. Scaffolding with a closely related species oriented and ordered the sequences in a manner more representative of the structure of the genome without altering the sequence. Comparisons with sequenced Juglandaceae revealed high levels of synteny and further supported J. cinerea's recent phylogenetic placement. Comparative assessment of gene family evolution revealed a significant number of contracting families, including several associated with biotic stress response.

Stage-specific transcription during development of Aedes aegypti.

BACKGROUND: Aedes aegypti is the most important global vector of dengue virus infection in humans. Availability of the draft genome sequence of this mosquito provides unique opportunities to study different aspects of its biology, including identification of genes and pathways relevant to the developmental processes associated with transition across individual life stages. However, detailed knowledge of gene expression patterns pertaining to developmental stages of A. aegypti is largely lacking. RESULTS: We performed custom cDNA microarray analyses to examine the expression patterns among six developmental stages: early larvae, late larvae, early pupae, late pupae, and adult male and female mosquitoes. Results revealed 1,551 differentially expressed transcripts (DETs) showing significant differences in levels of expression between these life stages. The data suggests that most of the differential expression occurs in a stage specific manner in A. aegypti. Based on hierarchical clustering of expression levels, correlated expression patterns of DETs were also observed among developmental stages. Weighted gene correlation network analysis revealed modular patterns of expression among the DETs. We observed that hydrolase activity, membrane, integral to membrane, DNA binding, translation, ribosome, nucleoside-triphosphatase activity, structural constituent of ribosome, ribonucleoprotein complex and receptor activity were among the top ten ranked GO (Gene Ontology) terms associated with DETs. Significant associations of DETs were also observed with specific KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway modules. Finally, comparisons with the previously reported developmental transcriptome of the malaria vector, Anopheles gambiae, indicated that gene expression patterns during developmental processes reflect both species-specific as well as common components of the two mosquito species. CONCLUSIONS: Our study shows that genes involved in the developmental life cycle of A. aegypti are expressed in a highly stage-specific manner. This suggests that transcriptional events associated with transition through larval, pupal and adult stages are largely discrete.

Genome sequences of the human body louse and its primary endosymbiont provide insights into the permanent parasitic lifestyle.

As an obligatory parasite of humans, the body louse (Pediculus humanus humanus) is an important vector for human diseases, including epidemic typhus, relapsing fever, and trench fever. Here, we present genome sequences of the body louse and its primary bacterial endosymbiont Candidatus Riesia pediculicola. The body louse has the smallest known insect genome, spanning 108 Mb. Despite its status as an obligate parasite, it retains a remarkably complete basal insect repertoire of 10,773 protein-coding genes and 57 microRNAs. Representing hemimetabolous insects, the genome of the body louse thus provides a reference for studies of holometabolous insects. Compared with other insect genomes, the body louse genome contains significantly fewer genes associated with environmental sensing and response, including odorant and gustatory receptors and detoxifying enzymes. The unique architecture of the 18 minicircular mitochondrial chromosomes of the body louse may be linked to the loss of the gene encoding the mitochondrial single-stranded DNA binding protein. The genome of the obligatory louse endosymbiont Candidatus Riesia pediculicola encodes less than 600 genes on a short, linear chromosome and a circular plasmid. The plasmid harbors a unique arrangement of genes required for the synthesis of pantothenate, an essential vitamin deficient in the louse diet. The human body louse, its primary endosymbiont, and the bacterial pathogens that it vectors all possess genomes reduced in size compared with their free-living close relatives. Thus, the body louse genome project offers unique information and tools to use in advancing understanding of coevolution among vectors, symbionts, and pathogens.

A haplotype-resolved chromosome-scale genome for Quercus rubra L. provides insights into the genetics of adaptive traits for red oak species.

Northern red oak (Quercus rubra L.) is an ecologically and economically important forest tree native to North America. We present a chromosome-scale genome of Q. rubra generated by the combination of PacBio sequences and chromatin conformation capture (Hi-C) scaffolding. This is the first reference genome from the red oak clade (section Lobatae). The Q. rubra assembly spans 739 Mb with 95.27% of the genome in 12 chromosomes and 33,333 protein-coding genes. Comparisons to the genomes of Quercus lobata and Quercus mongolica revealed high collinearity, with intrachromosomal structural variants present. Orthologous gene family analysis with other tree species revealed that gene families associated with defense response were expanding and contracting simultaneously across the Q. rubra genome. Quercus rubra had the most CC-NBS-LRR and TIR-NBS-LRR resistance genes out of the 9 species analyzed. Terpene synthase gene family comparisons further reveal tandem gene duplications in TPS-b subfamily, similar to Quercus robur. Phylogenetic analysis also identified 4 subfamilies of the IGT/LAZY gene family in Q. rubra important for plant structure. Single major QTL regions were identified for vegetative bud break and marcescence, which contain candidate genes for further research, including a putative ortholog of the circadian clock constituent cryptochrome (CRY2) and 8 tandemly duplicated genes for serine protease inhibitors, respectively. Genome-environment associations across natural populations identified candidate abiotic stress tolerance genes and predicted performance in a common garden. This high-quality red oak genome represents an essential resource to the oak genomic community, which will expedite comparative genomics and biological studies in Quercus species.

A high-quality reference genome for Fraxinus pennsylvanica for ash species restoration and research.

Green ash (Fraxinus pennsylvanica) is the most widely distributed ash tree in North America. Once common, it has experienced high mortality from the non-native invasive emerald ash borer (EAB; Agrilus planipennis). A small percentage of native green ash trees that remain healthy in long-infested areas, termed "lingering ash," display partial resistance to the insect, indicating that breeding and propagating populations with higher resistance to EAB may be possible. To assist in ash breeding, ecology and evolution studies, we report the first chromosome-level assembly from the genus Fraxinus for F. pennsylvanica with over 99% of bases anchored to 23 haploid chromosomes, spanning 757 Mb in total, composed of 49.43% repetitive DNA, and containing 35,470 high-confidence gene models assigned to 22,976 Asterid orthogroups. We also present results of range-wide genetic variation studies, the identification of candidate genes for important traits including potential EAB-resistance genes, and an investigation of comparative genome organization among Asterids based on this reference genome platform. Residual duplicated regions within the genome probably resulting from a recent whole genome duplication event in Oleaceae were visualized in relation to wild olive (Olea europaea var. sylvestris). We used our F. pennsylvanica chromosome assembly to construct reference-guided assemblies of 27 previously sequenced Fraxinus taxa, including F. excelsior. Thus, we present a significant step forward in genomic resources for research and protection of Fraxinus species.

A Roadmap for Participatory Chestnut Breeding for Nut Production in the Eastern United States.

Chestnut cultivation for nut production is increasing in the eastern half of the United States. Chinese chestnuts (Castanea mollissima Blume), or Chinese hybrids with European (C. sativa Mill.) and Japanese chestnuts (C. crenata Sieb. & Zucc.), are cultivated due to their high kernel quality, climatic adaptation, and disease resistance. Several hundred thousand pounds of high-quality fresh nuts are taken to market every fall, and several hundred additional orchards are entering bearing years. Grower-led on-farm improvement has largely facilitated this growth. A lack of significant investments in chestnut breeding in the region, paired with issues of graft incompatibility, has led many growers to cultivate seedlings of cultivars rather than grafted cultivars. After decades of evaluation, selection, and sharing of plant materials, growers have reached a threshold of improvement where commercial seedling orchards can be reliably established by planting offspring from elite selected parents. Growers recognize that if cooperation persists and university expertise and resources are enlisted, improvement can continue and accelerate. To this end, the University of Missouri Center for Agroforestry (UMCA) and chestnut growers throughout the eastern United States are partnering to formalize a participatory breeding program - the Chestnut Improvement Network. This partnership entails the UMCA providing an organizational structure and leadership to coordinate on-farm improvement, implement strategic crossing schemes, and integrate genetic tools. Chestnut growers offer structural capacity by cultivating seedling production orchards that provide financial support for the grower but also house segregating populations with improved individuals, in situ repositories, and selection trials, creating great value for the industry.

A Seed-Endophytic Bacillus safensis Strain With Antimicrobial Activity Has Genes for Novel Bacteriocin-Like Antimicrobial Peptides.

Bacteriocins are a highly diverse group of antimicrobial peptides that have been identified in a wide range of commensal and probiotic organisms, especially those resident in host microbiomes. Rising antibiotic resistance have fueled renewed research into new drug scaffolds such as antimicrobial peptides for use in therapeutics. In this investigation, we examined mung bean seeds for endophytes possessing activity against human and plant pathogens. We isolated a novel strain of Bacillus safensis, from the contents of surface-sterilized mung bean seed, which we termed B. safensis C3. Genome sequencing of C3 identified three distinct biosynthetic systems that produce bacteriocin-based peptides. C3 exhibited antibacterial activity against Escherichia coli, Xanthomonas axonopodis, and Pseudomonas syringae. Robust antimicrobial activity of B. safensis C3 was observed when C3 was co-cultured with Bacillus subtilis. Using the cell-free supernatant of C3 and cation exchange chromatography, we enriched a product that retained antimicrobial activity against B. subtilis. The peptide was found to be approximately 3.3 kDa in size by mass spectrometry, and resistant to proteolysis by Carboxypeptidase Y and Endoproteinase GluC, suggesting that it is a modified variant of an AS-48 like bacteriocin. Our findings open new avenues into further development of novel bacteriocin-based scaffolds for therapeutic development, as well as further investigations into how our discoveries of bacteriocin-producing plant commensal microorganisms may have the potential for an immediate impact on the safety of food supplies.

The power and promise of genetic mapping from Plasmodium falciparum crosses utilizing human liver-chimeric mice.

Genetic crosses are most powerful for linkage analysis when progeny numbers are high, parental alleles segregate evenly and numbers of inbred progeny are minimized. We previously developed a novel genetic crossing platform for the human malaria parasite Plasmodium falciparum, an obligately sexual, hermaphroditic protozoan, using mice carrying human hepatocytes (the human liver-chimeric FRG NOD huHep mouse) as the vertebrate host. We report on two genetic crosses-(1) an allopatric cross between a laboratory-adapted parasite (NF54) of African origin and a recently patient-derived Asian parasite, and (2) a sympatric cross between two recently patient-derived Asian parasites. We generated 144 unique recombinant clones from the two crosses, doubling the number of unique recombinant progeny generated in the previous 30 years. The allopatric African/Asian cross has minimal levels of inbreeding and extreme segregation distortion, while in the sympatric Asian cross, inbred progeny predominate and parental alleles segregate evenly. Using simulations, we demonstrate that these progeny provide the power to map small-effect mutations and epistatic interactions. The segregation distortion in the allopatric cross slightly erodes power to detect linkage in several genome regions. We greatly increase the power and the precision to map biomedically important traits with these new large progeny panels.

Range-wide distribution of genetic diversity in the North American tree Juglans cinerea: a product of range shifts, not ecological marginality or recent population decline.

The spatial distribution of genetic diversity is a product of recent and historical ecological processes, as well as anthropogenic activities. A current challenge in population and conservation genetics is to disentangle the relative effects of these processes, as a first step in predicting population response to future environmental change. In this investigation, we compare the influence of contemporary population decline, contemporary ecological marginality and postglacial range shifts. Using classical model comparison procedures and Bayesian methods, we have identified postglacial range shift as the clear determinant of genetic diversity, differentiation and bottlenecks in 29 populations of butternut, Juglans cinerea L., a North American outcrossing forest tree. Although butternut has experienced dramatic 20th century decline because of an introduced fungal pathogen, our analysis indicates that recent population decline has had less genetic impact than postglacial recolonization history. Location within the range edge vs. the range core also failed to account for the observed patterns of diversity and differentiation. Our results suggest that the genetic impact of large-scale recent population losses in forest trees should be considered in the light of Pleistocene-era large-scale range shifts that may have had long-term genetic consequences. The data also suggest that the population dynamics and life history of wind-pollinated forest trees may provide a buffer against steep population declines of short duration, a result having important implications for habitat management efforts, ex situ conservation sampling and population viability analysis.

Genomic Variation Among and Within Six Juglans Species.

Genomic analysis in Juglans (walnuts) is expected to transform the breeding and agricultural production of both nuts and lumber. To that end, we report here the determination of reference sequences for six additional relatives of Juglans regia: Juglans sigillata (also from section Dioscaryon), Juglans nigra, Juglans microcarpa, Juglans hindsii (from section Rhysocaryon), Juglans cathayensis (from section Cardiocaryon), and the closely related Pterocarya stenoptera While these are 'draft' genomes, ranging in size between 640Mbp and 990Mbp, their contiguities and accuracies can support powerful annotations of genomic variation that are often the foundation of new avenues of research and breeding. We annotated nucleotide divergence and synteny by creating complete pairwise alignments of each reference genome to the remaining six. In addition, we have re-sequenced a sample of accessions from four Juglans species (including regia). The variation discovered in these surveys comprises a critical resource for experimentation and breeding, as well as a solid complementary annotation. To demonstrate the potential of these resources the structural and sequence variation in and around the polyphenol oxidase loci, PPO1 and PPO2 were investigated. As reported for other seed crops variation in this gene is implicated in the domestication of walnuts. The apparently Juglandaceae specific PPO1 duplicate shows accelerated divergence and an excess of amino acid replacement on the lineage leading to accessions of the domesticated nut crop species, Juglans regia and sigillata.

Evidence for inversion polymorphism related to sympatric host race formation in the apple maggot fly, Rhagoletis pomonella.

Evidence suggests that the apple maggot, Rhagoletis pomonella (Diptera: Tephritidae) is undergoing sympatric speciation (i.e., divergence without geographic isolation) in the process of shifting and adapting to a new host plant. Prior to the introduction of cultivated apples (Malus pumila) in North America, R. pomonella infested the fruit of native hawthorns (Crataegus spp.). However, sometime in the mid-1800s the fly formed a sympatric race on apple. The recently derived apple-infesting race shows consistent allele frequency differences from the hawthorn host race for six allozyme loci mapping to three different chromosomes. Alleles at all six of these allozymes correlate with the timing of adult eclosion, an event dependent on the duration of the overwintering pupal diapause. This timing difference differentially adapts the univoltine fly races to an approximately 3- to 4-week difference in the peak fruiting times of apple and hawthorn trees, partially reproductively isolating the host races. Here, we report finding substantial gametic disequilibrium among allozyme and complementary DNA (cDNA) markers encompassing the three chromosomal regions differentiating apple and hawthorn flies. The regions of disequilibrium extend well beyond the previously characterized six allozyme loci, covering substantial portions of chromosomes 1, 2, and 3 (haploid n = 6 in R. pomonella). Moreover, significant recombination heterogeneity and variation in gene order were observed among single-pair crosses for each of the three genomic regions, implying the existence of inversion polymorphism. We therefore have evidence that genes affecting diapause traits involved in host race formation reside within large complexes of rearranged genes. We explore whether these genomic regions (inversions) constitute coadapted gene complexes and discuss the implications of our findings for sympatric speciation in Rhagoletis.

Preliminary Genomic Characterization of Ten Hardwood Tree Species from Multiplexed Low Coverage Whole Genome Sequencing.

Forest health issues are on the rise in the United States, resulting from introduction of alien pests and diseases, coupled with abiotic stresses related to climate change. Increasingly, forest scientists are finding genetic/genomic resources valuable in addressing forest health issues. For a set of ten ecologically and economically important native hardwood tree species representing a broad phylogenetic spectrum, we used low coverage whole genome sequencing from multiplex Illumina paired ends to economically profile their genomic content. For six species, the genome content was further analyzed by flow cytometry in order to determine the nuclear genome size. Sequencing yielded a depth of 0.8X to 7.5X, from which in silico analysis yielded preliminary estimates of gene and repetitive sequence content in the genome for each species. Thousands of genomic SSRs were identified, with a clear predisposition toward dinucleotide repeats and AT-rich repeat motifs. Flanking primers were designed for SSR loci for all ten species, ranging from 891 loci in sugar maple to 18,167 in redbay. In summary, we have demonstrated that useful preliminary genome information including repeat content, gene content and useful SSR markers can be obtained at low cost and time input from a single lane of Illumina multiplex sequence.