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.

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.

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.

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 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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Development of genomic microsatellites in Gleditsia triacanthos (Fabaceae) using Illumina sequencing.

PREMISE OF THE STUDY: Fourteen genomic microsatellite markers were developed and characterized in honey locust, Gleditsia triacanthos, using Illumina sequencing. Due to their high variability, these markers can be applied in analyses of genetic diversity and structure, and in mating system and gene flow studies. • METHODS AND RESULTS: Thirty-six individuals from across the species range were included in a genetic diversity analysis and yielded three to 20 alleles per locus. Observed heterozygosity and expected heterozygosity ranged from 0.214 to 0.944 and from 0.400 to 0.934, respectively, with minimal occurrence of null alleles. Regular segregation of maternal alleles was observed at seven loci and moderate segregation distortion at four of 11 loci that were heterozygous in the seed parent. • CONCLUSIONS: Honey locust is an important agroforestry tree capable of very fast growth and tolerance of poor site conditions. This is the first report of genomic microsatellites for this species.

Human-impacted landscapes facilitate hybridization between a native and an introduced tree.

Spatial and temporal dynamics of hybridization, in particular the influence of local environmental conditions, are well studied for sympatric species but less is known for native-introduced systems, especially for long-lived species. We used microsatellite and chloroplast DNA markers to characterize the influence of anthropogenic landscapes on the extent, direction, and spatial distribution of hybridization between a native North American tree Juglans cinerea (butternut) and an introduced tree Juglans ailantifolia (Japanese walnut) for 1363 trees at 48 locations across the native range of butternut. Remarkably, admixture in anthropogenic sites reached nearly 70%, while fragmented and continuous forests showed minimal admixture (<8%). Furthermore, more hybrids in anthropogenic sites had J. ailantifolia seed parents (95%) than hybrids in fragmented and continuous forests (69% and 59%, respectively). Our results show a strong influence of landscape type on rate and direction of realized gene flow. While hybrids are common in anthropogenic landscapes, our results suggest that even small forested landscapes serve as substantial barriers to hybrid establishment, a key consideration for butternut conservation planning, a species already exhibiting severe decline, and for other North American forest trees that hybridize with introduced congeners.

Comparative expression profiles of midgut genes in dengue virus refractory and susceptible Aedes aegypti across critical period for virus infection.

BACKGROUND: Aedes aegypti is the primary mosquito vector for dengue virus (DENV) worldwide. Infectivity of dengue virus varies among natural populations of this mosquito. How A. aegypti responds to DENV infection relative to which genes and associated pathways contribute to its differential susceptibility as a vector is not well defined. METHODS/PRINCIPAL FINDINGS: Here, we used custom cDNA microarrays to identify groups of genes that were differentially expressed in midgut tissues between susceptible and refractory strains in a highly time specific manner. While genes involved in protein processing in the endoplasmic reticulum, mRNA surveillance, and the proteasome were significantly up-regulated in the susceptible strain, several metabolic processes including glycolysis, glycan biosynthesis and Wnt pathway were active in the refractory strain. In addition, several key signaling genes were expressed as common responsive genes in both susceptible and refractory mosquitoes that may be necessary for signal transduction to trigger the appropriate host response to the viral infection. These are coordinately expressed in the form of tight gene networks and expression clusters that may be necessary to differentially contribute to the progression of dengue infection between the two strains. CONCLUSIONS: Our data show that highly correlated differential expression of responsive genes throughout the post infection period in A. aegypti midgut tissues is necessary for a coordinated transcriptional response of the mosquito genes to host or defend the viral infection.