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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08281-y.

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.

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.

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3765-8) contains supplementary material, which is available to authorized users.

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.

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.

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.

Thirteen nuclear microsatellite loci for butternut (Juglans cinerea L.)

Butternut (Juglans cinerea L.) is an eastern North American forest tree severely threatened by an exotic fungal pathogen, Sirococcus clavigignenti-juglandacearum. We report here 13 nuclear microsatellites for genetic evaluation of the remaining natural populations. Summary statistics are reported for individuals from a population of butternuts in central Kentucky (N = 63). All markers were polymorphic, with an average of 13.7 alleles per locus observed. Four loci exhibited significantly fewer heterozygotes than expected under Hardy-Weinberg equilibrium (P < 0.05).

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.

RNAi knock-downs support roles for the mucin-like (AeIMUC1) gene and short-chain dehydrogenase/reductase (SDR) gene in Aedes aegypti susceptibility to Plasmodium gallinaceum

The mosquito midgut represents the first barrier encountered by the Plasmodium parasite (Haemosporida: Plasmodiidae) when it is ingested in blood from an infected vertebrate. Previous studies identified the Aedes aegypti (L.) (Diptera: Culicidae) mucin-like (AeIMUC1) and short-chain dehydrogenase/reductase (SDR) genes as midgut-expressed candidate genes influencing susceptibility to infection by Plasmodium gallinaceum (Brumpt). We used RNA inference (RNAi) by double-stranded RNA (dsRNA) injections to examine ookinete survival to the oocyst stage following individual gene knock-downs. Double-stranded RNA gene knock-downs were performed 3 days prior to P. gallinaceum infection and oocyst development was evaluated at 7 days post-infection. Mean numbers of parasites developing to the oocyst stage were significantly reduced by 52.3% in dsAeIMUC1-injected females and by 36.5% in dsSDR-injected females compared with females injected with a dsβ-gal control. The prevalence of infection was significantly reduced in dsAeIMUC1- and dsSDR-injected females compared with females injected with dsβ-gal; these reductions resulted in a two- and three-fold increase in the number of uninfected individuals, respectively. Overall, these results suggest that both AeIMUC1 and SDR play a role in Ae. aegypti vector competence to P. gallinaceum.

Complete sequences of mitochondria genomes of Aedes aegypti and Culex quinquefasciatus and comparative analysis of mitochondrial DNA fragments inserted in the nuclear genomes

We present complete sequences of the mitochondrial genomes for two important mosquitoes, Aedes aegypti and Culex quinquefasciatus, that are major vectors of dengue virus and lymphatic filariasis, respectively. The A. aegypti mitochondrial genome is 16,655 bp in length and that of C. quinquefasciatus is 15,587 bp, yet both contain 13 protein coding genes, 22 transfer RNA (tRNA) genes, one 12S ribosomal RNA (rRNA) gene, one 16S rRNA gene and a control region (CR) in the same order. The difference in the genome size is due to the difference in the length of the control region. We also analyzed insertions of nuclear copies of mtDNA-like sequences (NUMTs) in a comparative manner between the two mosquitoes. The NUMT sequences occupy ~0.008% of the A. aegypti genome and ~0.001% of the C. quinquefasciatus genome. Several NUMTs were found localized in the introns of predicted protein coding genes in both genomes (32 genes in A. aegypti but only four in C. quinquefasciatus). None of these NUMT-containing genes had an ortholog between the two species or had paralogous copies within a genome that was also NUMT-containing. It was further observed that the NUMT-containing genes were relatively longer but had lower GC content compared to the NUMT-less paralogous copies. Moreover, stretches of homologies are present among the genic and non-genic NUMTs that may play important roles in genomic rearrangement of NUMTs in these genomes. Our study provides new insights on understanding the roles of nuclear mtDNA sequences in genome complexities of these mosquitoes.

Heritability and adaptive phenotypic plasticity of adult body size in the mosquito Aedes aegypti with implications for dengue vector competence

Adaptive phenotypic plasticity is particularly important to organisms with developmental cycles that undergo ontogenetic niche shifts that differentially subject individual life stages to heterogeneous and often stressful environmental conditions. The yellow fever and dengue fever vector mosquito, Aedes aegypti, typically breeds in small water-filled containers that expose the developing aquatic larvae to competition for resources with conspecifics and high probabilities for habitat drying. Here we investigated the heritability (h(2)) and phenotypic plasticity among A. aegypti laboratory populations and field populations from Trinidad, West Indies. Heritability for body size was moderate or completely eroded among the laboratory populations, while field populations contained high genetic variation among both males and females. Norms of reactions based on optimum vs. deficient larval conditions for artificial sibling families representing Trinidad field populations suggested significant gene × environment interactions influence body size and that there may be sex specific differences in allocation of resources. Individuals reared under optimum laboratory conditions were significantly larger and showed much less variability in body size plasticity than their field reared cohorts, suggesting that exposure to environmental stress may be common for A. aegypti larval development and would undoubtedly impact other traits, including arbovirus vector competence among adult females, in a similar fashion. Broad genetic variance in body size and other characters is likely maintained by balancing selection. Our results also suggest the need for caution in translating conclusions from experiments with laboratory colonies to natural populations. These would likely be more informative to expected phenotypes under natural conditions if conducted over a range of conditions that simulate environmental stress.

Development of bacterial artificial chromosome library resources for parasitoid Hymenoptera (Nasonia vitripennis and Nasonia giraulti: Pteromalidae)

The species of the genus Nasonia possess qualities that make them excellent candidates for genetic and genomic studies. To increase the wealth of genomic resources for the genus we constructed publicly available bacterial artificial chromosome (BAC) libraries for Nasonia vitripennis and Nasonia giraulti. Libraries have 36 864 clones each, empty-vector contents of approximately 2% and average insert sizes of 113.1 and 97.7 Kb, respectively, representing 12 and 11 genome equivalents. The N. vitripennis library was used for genome sequence assembly and in efforts at positional cloning of a developmental gene. The genome assembly of N. vitripennis is currently composed on 6181 un-joined scaffolds. These BAC libraries can be used to identify and close regions between scaffolds of the genome assemblies of both species.

Functional and evolutionary insights from the genomes of three parasitoid Nasonia species

We report here genome sequences and comparative analyses of three closely related parasitoid wasps: Nasonia vitripennis, N. giraulti, and N. longicornis. Parasitoids are important regulators of arthropod populations, including major agricultural pests and disease vectors, and Nasonia is an emerging genetic model, particularly for evolutionary and developmental genetics. Key findings include the identification of a functional DNA methylation tool kit; hymenopteran-specific genes including diverse venoms; lateral gene transfers among Pox viruses, Wolbachia, and Nasonia; and the rapid evolution of genes involved in nuclear-mitochondrial interactions that are implicated in speciation. Newly developed genome resources advance Nasonia for genetic research, accelerate mapping and cloning of quantitative trait loci, and will ultimately provide tools and knowledge for further increasing the utility of parasitoids as pest insect-control agents.

Geographically extensive hybridization between the forest trees American butternut and Japanese walnut

We investigate the question of naturally occurring interspecific hybrids between two forest trees: the native North American butternut (Juglans cinerea L.) and the introduced Japanese walnut (Juglans ailantifolia Carrière). Using nuclear and chloroplast DNA markers, we provide evidence for 29 F(1) and 22 advanced generation hybrids in seven locations across the eastern and southern range of the native species. Two locations show extensive admixture (95% J. ailantifolia and hybrids) while other locations show limited admixture. Hybridization appears to be asymmetrical with 90.9 per cent of hybrids having J. ailantifolia as the maternal parent. This is, to our knowledge, the first genetic data supporting natural hybridization between these species. The long-term outcome of introgression could include loss of native diversity, but could also include transfer of useful traits from the introduced species.

The CHD3 remodeler PICKLE promotes trimethylation of histone H3 lysine 27

CHD3 proteins are ATP-dependent chromatin remodelers that contribute to repression of developmentally regulated genes in both animal and plant systems. In animals, this repression has been linked to a multiple subunit complex, Mi-2/NuRD, whose constituents include a CHD3 protein, a histone deacetylase, and a methyl-CpG-binding domain protein. In Arabidopsis, PICKLE (PKL) codes for a CHD3 protein that acts during germination to repress expression of seed-associated genes. Repression of seed-associated traits is promoted in pkl seedlings by the plant growth regulator gibberellin (GA). We undertook a microarray analysis to determine how PKL and GA act to promote the transition from seed to seedling. We found that PKL and GA act in separate pathways to repress expression of seed-specific genes. Comparison of genomic datasets revealed that PKL-dependent genes are enriched for trimethylation of histone H3 lysine 27 (H3K27me3), a repressive epigenetic mark. Chromatin immunoprecipitation studies demonstrate that PKL promotes H3K27me3 in both germinating seedlings and in adult plants but do not identify a connection between PKL-dependent expression and acetylation levels. Taken together, our analyses illuminate a new pathway by which CHD3 remodelers contribute to repression in eukaryotes.

Quantitative trait loci determining autogeny and body size in the Asian tiger mosquito (Aedes albopictus)

The majority of mosquito species require a blood meal to stimulate vitellogenesis and subsequent oviposition (anautogeny), but some autogenous individuals complete their first ovarian cycle without a blood meal. Autogeny may be facultative or obligatory. In this study, we selected for an autogenous strain in the Asian tiger mosquito Aedes albopictus and examined an F(1) intercross population for quantitative trait loci (QTL) determining the autogeny trait as well as wing length as a proxy for body size. Using composite interval mapping, we identified four QTL for each trait and observed considerable overlap in genome positions between each QTL for autogeny (follicle size) and wing length. Most QTL were minor in magnitude, individually explaining <10% of the phenotypic variation. Alleles from the autogenous parent generally showed a dominance or overdominance effect on both phenotypes. Strong genetic and phenotypic correlations indicate that autogeny and wing length are determined by up to four clusters of tightly linked genes or the potential pleiotropic effects of single genes. Although females from the autogenous strain produced approximately fivefold more eggs following a blood meal than through autogeny, we suggest that the maintenance of alleles for autogeny in natural populations is likely due to balancing selection. Autogeny should be favored under conditions of limited host availability for blood feeding or increased defensive behavior by the host and adequate larval nutrition. Correlation between autogeny and body size may reflect an increased ability for larger females to accumulate sufficient nutrient reserves to support oogenesis without the requirement for a blood meal.

Investigations of dengue-2 susceptibility and body size among Aedes aegypti populations

The mosquito Aedes aegypti (L.) (Diptera: Culicidae) is the primary global vector for dengue virus (DENV), yet considerable genetic variation exists among populations in terms of its competence to vector DENV. Variability in adult body size has also been observed among various mosquito populations and several studies have reported a relationship between body size and arbovirus dissemination, although most of these relied on artificially derived variation in body size. Here we examine the relationship between body size and disseminated DENV infection among 10 Ae. aegypti populations reared under optimum laboratory conditions. Body size variability was inferred from wing length measurements and DENV competence was evaluated as the proportion of individuals with disseminated infections following exposure to the dengue-2 JAM1409 strain. There were significant differences in mean wing lengths among populations (anova, F(9,22)= 7.10, P < 0.0001), ranging from 2.16 mm (Bangkok population) to 2.79 mm (MOYO-S [susceptible] population). We also observed significant differences among some populations in mean DENV infection rates (Waller-Duncan K-ratio t-test), ranging from 19.54% (MOYO-R [refractory] population) to 56.60% (MOYO-S population). However, we did not observe evidence for significant interactions between body size and DENV dissemination. We suggest that either the two traits are genetically independent or that our ability to detect interactions between them was limited by their respective inheritances as quantitative traits.

Genetic basis for reproductive diapause is correlated with life history traits within the Culex pipiens complex

The evolution of late season reproductive arrest (diapause) among female Culex pipiens mosquitoes allows them to overwinter in temperate climates, while females of the sibling species Culex quinquefasciatus do not exhibit the diapause phenotype. We present results for quantitative trait loci (QTL) analyses of two independent segregating populations derived from crosses between C. pipiens and C. quinquefasciatus. QTL for diapause and three life history traits were identified and compared for genome positions and gene effects. Using a combination of composite interval mapping, single-marker analysis and all possible subsets regression, we identified multiple QTL for each trait, totalling 14 and 17 QTL for each population, respectively. Individual QTL across traits often mapped to similar genome locations, suggesting these traits may be controlled in part by genes with pleiotropic effects or multiple tightly linked genes. The majority of QTL were intermediate in magnitude in that they explained 10-25% of the phenotypic variation. The majority of QTL showed overdominance effects. We suggest that this could impact natural populations, as increased heterosis across hybrid zones may allow populations to adapt to environmental conditions via stabilizing selection, and yet maintain species identity outside these regions because of strong morphological integration, wherein related traits evolve as an integrated unit.

Sequencing of a new target genome: the Pediculus humanus humanus (Phthiraptera: Pediculidae) genome project

The human body louse, Pediculus humanus humanus (L.), and the human head louse, Pediculus humanus capitis, belong to the hemimetabolous order Phthiraptera. The body louse is the primary vector that transmits the bacterial agents of louse-borne relapsing fever, trench fever, and epidemic typhus. The genomes of the bacterial causative agents of several of these aforementioned diseases have been sequenced. Thus, determining the body louse genome will enhance studies of host-vector-pathogen interactions. Although not important as a major disease vector, head lice are of major social concern. Resistance to traditional pesticides used to control head and body lice have developed. It is imperative that new molecular targets be discovered for the development of novel compounds to control these insects. No complete genome sequence exists for a hemimetabolous insect species primarily because hemimetabolous insects often have large (2000 Mb) to very large (up to 16,300 Mb) genomes. Fortuitously, we determined that the human body louse has one of the smallest genome sizes known in insects, suggesting it may be a suitable choice as a minimal hemimetabolous genome in which many genes have been eliminated during its adaptation to human parasitism. Because many louse species infest birds and mammals, the body louse genome-sequencing project will facilitate studies of their comparative genomics. A 6-8X coverage of the body louse genome, plus sequenced expressed sequence tags, should provide the entomological, evolutionary biology, medical, and public health communities with useful genetic information.

Application of One Sided t-tests and a Generalized Experiment Wise Error Rate to High-Density Oligonucleotide Microarray Experiments: An Example Using Arabidopsis

MOTIVATION: A formidable challenge in the analysis of microarray data is the identification of those genes that exhibit differential expression. The objectives of this research were to examine the utility of simple ANOVA, one sided t tests, natural log transformation, and a generalized experiment wise error rate methodology for analysis of such experiments. As a test case, we analyzed a Affymetrix GeneChip microarray experiment designed to test for the effect of a CHD3 chromatin remodeling factor, PICKLE, and an inhibitor of the plant hormone gibberellin (GA), on the expression of 8256 Arabidopsis thaliana genes. RESULTS: The GFWER(k) is defined as the probability of rejecting k or more true null hypothesis at a given p level. Computing probabilities by GFWER(k) was shown to be simple to apply and, depending on the value of k, can greatly increase power. A k value as small as 2 or 3 was concluded to be adequate for large or small experiments respectively. A one sided t-test along with GFWER(2)=.05 identified 43 genes as exhibiting PICKLE-dependent expression. Expression of all 43 genes was re-examined by qRT-PCR, of which 36 (83.7%) were confirmed to exhibit PICKLE-dependent expression.

PICKLE acts throughout the plant to repress expression of embryonic traits and may play a role in gibberellin-dependent responses

A seed marks the transition between two developmental states; a plant is an embryo during seed formation, whereas it is a seedling after emergence from the seed. Two factors have been identified in Arabidopsis that play a role in establishment of repression of the embryonic state: PKL (PICKLE), which codes for a putative CHD3 chromatin remodeling factor, and gibberellin (GA), a plant growth regulator. Previous observations have also suggested that PKL mediates some aspects of GA responsiveness in the adult plant. To investigate possible mechanisms by which PKL and GA might act to repress the embryonic state, we further characterized the ability of PKL and GA to repress embryonic traits and reexamined the role of PKL in mediating GA-dependent responses. We found that PKL acts throughout the seedling to repress expression of embryonic traits. Although the ability of pkl seedlings to express embryonic traits is strongly induced by inhibiting GA biosynthesis, it is only marginally responsive to abscisic acid and SPY (SPINDLY), factors that have previously been demonstrated to inhibit GA-dependent responses during germination. We also observed that pkl plants exhibit the phenotypic hallmarks of a mutation in a positive regulator of a GA response pathway including reduced GA responsiveness and increased synthesis of bioactive GAs. These observations indicate that PKL may mediate a subset of GA-dependent responses during shoot development.

Genomic analysis in the sting-2 quantitative trait locus for defensive behavior in the honey bee, Apis mellifera

We have sequenced an 81-kb genomic region from the honey bee, Apis mellifera, associated with a quantitative trait locus (QTL) sting-2 for aggressive behavior. This sequence represents the first extensive study of the honey-bee genome structure encompassing putative genes in a QTL for a behavioral trait. Expression of 13 putative genes, as well as two transcripts that were present in a honey-bee EST database, was confirmed through reverse transcription analysis of mRNA from the honey-bee head. Whereas most transcripts exhibited little or no variation between European and Africanized honey-bee alleles, one transcript demonstrated significant nonsynonymous substitutions, deletions, and insertions. All 13 putative genes lacked similarity to known invertebrate or vertebrate proteins or transcripts. This observation may be reflective of the processes that determine the genomic evolution of an insect with social behavior and/or haplo-diploidy and are an indication of the unique nature of the honey-bee genome. These results make this sequence an invaluable research tool for the ongoing honey-bee whole-genome sequencing effort.

Transcriptome identification of putative genes involved in protein catabolism and innate immune response in human body louse (Pediculicidae: Pediculus humanus)

Genomics information relating to human body lice is surprisingly scarce, and this has constrained studies of their physiology, immunology and vector biology. To identify novel body louse genes, we used engorged adult lice to generate a cDNA library. Initially, 1152 clones were screened for inserts, edited for removal of vector sequences and base pairs of poor quality, and viewed for splicing variations, gene families and polymorphism. Computational methods identified 506 inferred open reading frames including the first predicted louse defensin. The inferred defensin aligns well with other insect defensins and has highly conserved cysteine residues, as are known for other defensin sequences. Two cysteine and five serine proteinases were categorized according to their inferred catalytic sites. We also discovered seven putative ubiquitin-pathway genes and four iron metabolizing deduced enzymes. Finally, glutathione-S-transferases and cytochrome P450 genes were among the detoxification enzymes found. Results from this first systematic effort to discover human body louse genes should promote further studies in Phthiraptera and lice.

Transcriptome analysis of the cowpea weevil bruchid: identification of putative proteinases and alpha-amylases associated with food breakdown

We describe here the first systematic work to discover insect genes involved in food breakdown using a cDNA library enriched for gut-expressed transcripts from Callosobruchus maculatus. A total of 1056 clones were screened for cDNA insert-containing plasmids, and 503 nonredundant open reading frames were discovered. Twenty-three inferred genes potentially involved in digestive processes in cowpea weevil were identified, including proteinases and amylases. The predicted catalytic sites were identified in the inferred cysteine and aspartic acid proteinases, and in alpha-amylases. Transcriptome analysis of the cowpea bruchid will potentially permit gene discovery in other beetles, an insect order of major economic and ecological importance that is poorly represented in genomic databases.

Genetic basis of heterosis explored by simple sequence repeat markers in a random-mated maize population

The genetic basis of heterosis in crop plants has not been completely resolved. Our objective in this study was to determine the level of dominance for quantitative trait loci (QTLs) that underlie heterosis in maize (Zea mays L.). An F2 population of an elite maize single cross, LH200 x LH216, was random mated for three generations in an attempt to break up repulsion linkages that might lead to pseudo-overdominance. The population was analyzed with 160 simple-sequence repeat markers. Phenotypic data analyses indicated overdominance for grain yield and partial dominance for plant height, grain moisture and stalk lodging. A total of 28 QTLs were identified for grain yield, 16 for grain moisture, 8 for stalk lodging, and 11 for plant height. For grain yield, 24 QTLs (86%) showed overdominance. In contrast, most of the QTLs for plant height, grain moisture and stalk lodging showed partial to complete dominance. Little epistasis was detected among the QTLs for any of the traits. Our results can be interpreted in one of two ways, or a combination of both: (1) QTLs for grain yield in maize exhibit true overdominance, or (2) QTLs for grain yield show partial to complete dominance, but they are so tightly linked such that three generations of random mating failed to separate their individual effects.

Coordinate repression of regulators of embryonic identity by PICKLE during germination in Arabidopsis

In angiosperms, germination represents an important developmental transition during which embryonic identity is repressed and vegetative identity emerges. PICKLE (PKL) encodes a CHD3-chromatin-remodeling factor necessary for the repression of expression of LEAFY COTYLEDON1 (LEC1), a central regulator of embryogenesis. A candidate gene approach and microarray analysis identified nine additional genes that exhibit PKL-dependent repression of expression during germination. Transcripts for all three LEAFY COTYLEDON genes, LEC1, LEC2, and FUS3, exhibit PKL-dependent repression, and all three transcripts are elevated more than 100-fold in pkl primary roots that inappropriately express embryonic traits (pickle roots). Three other genes that exhibit PKL-dependent regulation have expression patterns correlated with zygotic or somatic embryogenesis, and one gene encodes a putative Lin-11, Isl-1, MEC-3 (LIM) domain transcriptional regulator that is preferentially expressed in siliques. Genes that exhibit PKL-dependent repression during germination are not necessarily regulated by PKL at other points in development. Our data suggest that PKL selectively regulates a suite of genes during germination to repress embryonic identity. In particular, we propose that PKL acts as a master regulator of the LEAFY COTYLEDON genes, and that joint derepression of these genes is likely to contribute substantially to expression of embryonic identity in pkl seedlings.

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.

Inversions and gene order shuffling in Anopheles gambiae and A. funestus

In tropical Africa, Anopheles funestus is one of the three most important malaria vectors. We physically mapped 157 A. funestus complementary DNAs (cDNAs) to the polytene chromosomes of this species. Sequences of the cDNAs were mapped in silico to the A. gambiae genome as part of a comparative genomic study of synteny, gene order, and sequence conservation between A. funestus and A. gambiae. These species are in the same subgenus and diverged about as recently as humans and chimpanzees. Despite nearly perfect preservation of synteny, we found substantial shuffling of gene order along corresponding chromosome arms. Since the divergence of these species, at least 70 chromosomal inversions have been fixed, the highest rate of rearrangement of any eukaryote studied to date. The high incidence of paracentric inversions and limited colinearity suggests that locating genes in one anopheline species based on gene order in another may be limited to closely related taxa.

Chromosomal regions associated with segregation distortion in maize

Segregation distortion skews the genotypic frequencies from their Mendelian expectations. Our objectives in this study were to assess the frequency of occurrence of segregation distortion in maize, identify chromosomal regions consistently associated with segregation distortion, and examine the effects of gametophytic factors on linkage mapping. We constructed a simple sequence repeat (SSR) linkage map for a LH200/LH216 F(2)Syn3 (i.e., random-mated three times) population, and compared the segregation distortion in this map with the segregation distortion in three published linkage maps. Among 1,820 codominant markers across the four mapping populations, 301 (17%) showed segregation distortion ( P < 0.05). The frequency of markers showing segregation distortion ranged from 19% in the Tx303/CO159 mapping population to 36% in the B73/Mo17 mapping population. A positive relationship was found between the number of meioses and the frequency of segregation distortion detected in a population. On a given chromosome, nearly all of the markers showing segregation distortion favored the allele from the same parent. A total of 18 chromosomal regions on the ten maize chromosomes were associated with segregation distortion. The consistent location of these chromosomal regions in four populations suggested the presence of segregation distortion regions (SDRs). Three known gametophytic factors are possible genetic causes of these SDRs. As shown in previous research, segregation distortion does not affect the estimate of map distance when only one gametophytic factor is present in an SDR.

Population genetics of the yellow fever mosquito in Trinidad: comparisons of amplified fragment length polymorphism (AFLP) and restriction fragment length polymorphism (RFLP) markers

Recent development of DNA markers provides powerful tools for population genetic analyses. Amplified fragment length polymorphism (AFLP) markers result from a polymerase chain reaction (PCR)-based DNA fingerprinting technique that can detect multiple restriction fragments in a single polyacrylamide gel, and thus are potentially useful for population genetic studies. Because AFLP markers have to be analysed as dominant loci in order to estimate population genetic diversity and genetic structure parameters, one must assume that dominant (amplified) alleles are identical in state, recessive (unamplified) alleles are identical in state, AFLP fragments segregate according to Mendelian expectations and that the genotypes of an AFLP locus are in Hardy-Weinberg equilibrium (HWE). The HWE assumption is untestable for natural populations using dominant markers. Restriction fragment length polymorphism (RFLP) markers segregate as codominant alleles, and can therefore be used to test the HWE assumption that is critical for analysing AFLP data. This study examined whether the dominant AFLP markers could provide accurate estimates of genetic variability for the Aedes aegypti mosquito populations of Trinidad, West Indies, by comparing genetic structure parameters using AFLP and RFLP markers. For AFLP markers, we tested a total of five primer combinations and scored 137 putative loci. For RFLP, we examined a total of eight mapped markers that provide a broad coverage of mosquito genome. The estimated average heterozygosity with AFLP markers was similar among the populations (0.39), and the observed average heterozygosity with RFLP markers varied from 0.44 to 0.58. The average FST (standardized among-population genetic variance) estimates were 0.033 for AFLP and 0.063 for RFLP markers. The genotypes at several RFLP loci were not in HWE, suggesting that the assumption critical for analysing AFLP data was invalid for some loci of the mosquito populations in Trinidad. Therefore, the results suggest that, compared with dominant molecular markers, codominant DNA markers provide better estimates of population genetic variability, and offer more statistical power for detecting population genetic structure.

Amplification of DNA markers from evolutionarily diverse genomes using single primers of simple-sequence repeats

The abundance and scattered distribution of simple-sequence repeats (SSR) in eukaryotic genomes prompted us to explore the use of SSR-based oligonucleotide primers in single primer amplification reactions. In a pilot experiment, 23 primers were used across a panel of evolutionarily diverse eukaryotic genomes, including grapes, lettuce, tomato, pine, maize, salmon, chicken, Holstein cows and humans. The primers were 16-20 bases in length and represented SSRs of di-, tri-, tetra-, and pentanucleotide repeats. The results showed that tetranucleotide repeat primers were most effective in amplifying polymorphic patterns. Of 11 such primers tested, 70% produced polymorphic patterns from the DNA of one or more species. Primers representing a combination of two tetranucleotide repeats, or compound microsatellites, were equally effective. The polymorphisms contained in such fingerprints were able to identify individuals of vertebrate species as well as lines or varieties of plants. Inheritance of the polymorphic bands was studied in a maize recombinant inbred population, DE811 x B73. Thirty-two polymorphic bands, derived from two amplification patterns, were mapped as dominant markers on an existing RFLP map of the same population. The bands were distributed across nine of the ten chromosomes.

Pinning down loose ends: mapping telomeres and factors affecting their length

A degenerately repeated sequence, proximal to the telomere heptanucleotide repeat in maize, contains restriction enzyme sites that permit the separation of telomeres from the rest of the chromosomes. Probing with a telomere-specific oligonucleotide revealed genotype-dependent telomere lengths that vary more than 25-fold in maize among the 22 inbreds that have been surveyed. These lengths were found to segregate reproducibly in a recombinant inbred family where 50% of the variation can be accounted for by three loci. The dynamic control over telomere length in maize appears to act rapidly to achieve new genotypically determined telomere lengths in the F1. Clones of telomere proximal sequences were used to map restriction fragment length loci at the distal ends of eight of 20 chromosome arms.

Improved tissue culture response of an elite maize inbred through backcross breeding, and identification of chromosomal regions important for regeneration by RFLP analysis

The frequency of initiation of friable, embryogenic callus from immature embryos of the elite maize inbred line B73 was increased dramatically by introgression of chromosomal segments from the inbred line A188 through classical backcross breeding. Less than 0.2% of the immature B73 embryos tested (5 of 3,710) formed embryogenic callus. The breeding scheme consisted of six generations of backcrossing to B73 with selection at each generation for high frequency initiation of embryogenic cultures. BC6 individuals were selfed for four generations to select homozygous lines. The average embryogenic culture initiation frequency increased to 46% (256/561). Nearly all (91%) of the embryos from one BC6 S4 plant formed embryogenic cultures. RFLP analysis was used to determine the locations and effects of the introgressed A188 chromosomal segments. Five segments were retained through at least the fifth backcross generation. The hypothesis that one or more of these five regions contains genes controlling somatic embryogenesis in maize was tested using an F2 population of the cross A188 X Mo17. A set of five DNA markers (three of them linked) explained 82% of the observed phenotypic variance for percentage of immature embryos forming embryognic callus. Four of the five markers were located in or near introgressed A188 chromosome segments.The region marked by probe c595 on the long arm of chromosome 9 was highly associated with several measures of in vitro culture response (percent embryogenic embryos, plants per embryo, and plants per embryogenic embryo). We propose that there is a major gene (or genes) in this region in A188 that promotes embryogenic callus initiation and plant regeneration in B73, Mo17, and probably many other recalcitrant inbred lines of maize.

Effects on soybean seed yield by members of 123 serocluster and USDA 110 in the greenhouse and in soils low in indigenous Bradyrhizobium japonicum

USDA 123spc and USDA 110str were inoculated on three soybean cultivars in a pot experiment in the greenhouse. USDA 110str yielded 33% more seed than USDA 123spc (P < 0.01). In a low-Bradyrhizobium soil (10 cells/g) at Arlington, WI, three members of the 123 serocluster (USDA 123spc, NJ2-lastr, and I-123spc), USDA 110, and USDA 138 were inoculated at a rate of 1 × 108 cells/2.5-cm row on two soybean cultivars. The inoculum strains formed 100% of the nodules. There was no difference (P < 0.05) in seed yield between any of the inoculum strains on either cultivar. In a low-Bradyrhizobium (102 cells/g), low-N soil, in Georgia, USDA 110, USDA 110str, and two members of the 123 serocluster (USDA 123spc and NJ2-lastr) were inoculated at an inoculum rate of 1 × 108 cell/2.5-cm row on two soybean cultivars. Plots were fertilized with 156 kg of N ha−1, or no N was applied. Serogroup 110 formed 59–90% of the nodules and serocluster 123 formed 48–82% of the nodules on the cultivars. There were no differences in seed yield between the inoculum strains, except that on cv. Gordon with added N, USDA 110 resulted in higher seed yield than USDA 123spc. The uninoculated control plots with cv. Braxton (without added N) also had a better seed yield (P < 0.05) than plots inoculated with USDA 123spc. In summary, inoculation with USDA 110str resulted in higher soybean seed yield than inoculation with USDA 123spc under greenhouse conditions, but not always under field conditions. Key words: Bradyrhizobium japonicum, field trials, seed yield.