Wayne Xu

Parent-of-Origin Effects on Gene Expression and DNA Methylation in the Maize Endosperm

This work uses deep sequencing of RNAs from maize endosperm and embryo to identify 54 maternally expressed genes and 46 paternally expressed genes and then examines genome-wide DNA methylation and gene expression, finding hypomethylation of the maternal allele and endosperm-specific expression for many of the imprinted genes. Imprinting describes the differential expression of alleles based on their parent of origin. Deep sequencing of RNAs from maize (Zea mays) endosperm and embryo tissue 14 d after pollination was used to identify imprinted genes among a set of ~12,000 genes that were expressed and contained sequence polymorphisms between the B73 and Mo17 genotypes. The analysis of parent-of-origin patterns of expression resulted in the identification of 100 putative imprinted genes in maize endosperm, including 54 maternally expressed genes (MEGs) and 46 paternally expressed genes (PEGs). Three of these genes have been previously identified as imprinted, while the remaining 97 genes represent novel imprinted maize genes. A genome-wide analysis of DNA methylation identified regions with reduced endosperm DNA methylation in, or near, 19 of the 100 imprinted genes. The reduced levels of DNA methylation in endosperm are caused by hypomethylation of the maternal allele for both MEGs and PEGs in all cases tested. Many of the imprinted genes with reduced DNA methylation levels also show endosperm-specific expression patterns. The imprinted maize genes were compared with imprinted genes identified in genome-wide screens of rice (Oryza sativa) and Arabidopsis thaliana, and at least 10 examples of conserved imprinting between maize and each of the other species were identified.

Structural variants in the soybean genome localize to clusters of biotic stress response genes

Genome-wide structural and gene content variations are hypothesized to drive important phenotypic variation within a species. Structural and gene content variations were assessed among four soybean (Glycine max) genotypes using array hybridization and targeted resequencing. Many chromosomes exhibited relatively low rates of structural variation (SV) among genotypes. However, several regions exhibited both copy number and presence-absence variation, the most prominent found on chromosomes 3, 6, 7, 16, and 18. Interestingly, the regions most enriched for SV were specifically localized to gene-rich regions that harbor clustered multigene families. The most abundant classes of gene families associated with these regions were the nucleotide-binding and receptor-like protein classes, both of which are important for plant biotic defense. The colocalization of SV with plant defense response signal transduction pathways provides insight into the mechanisms of soybean resistance gene evolution and may inform the development of new approaches to resistance gene cloning.

Phenotypic and genomic analyses of a fast neutron mutant population resource in soybean

Mutagenized populations have become indispensable resources for introducing variation and studying gene function in plant genomics research. In this study, fast neutron (FN) radiation was used to induce deletion mutations in the soybean (Glycine max) genome. Approximately 120,000 soybean seeds were exposed to FN radiation doses of up to 32 Gray units to develop over 23,000 independent M2 lines. Here, we demonstrate the utility of this population for phenotypic screening and associated genomic characterization of striking and agronomically important traits. Plant variation was cataloged for seed composition, maturity, morphology, pigmentation, and nodulation traits. Mutants that showed significant increases or decreases in seed protein and oil content across multiple generations and environments were identified. The application of comparative genomic hybridization (CGH) to lesion-induced mutants for deletion mapping was validated on a midoleate x-ray mutant, M23, with a known FAD2-1A (for fatty acid desaturase) gene deletion. Using CGH, a subset of mutants was characterized, revealing deletion regions and candidate genes associated with phenotypes of interest. Exome resequencing and sequencing of PCR products confirmed FN-induced deletions detected by CGH. Beyond characterization of soybean FN mutants, this study demonstrates the utility of CGH, exome sequence capture, and next-generation sequencing approaches for analyses of mutant plant genomes. We present this FN mutant soybean population as a valuable public resource for future genetic screens and functional genomics research.

The composition and origins of genomic variation among individuals of the soybean reference cultivar Williams 82

Soybean (Glycine max) is a self-pollinating species that has relatively low nucleotide polymorphism rates compared with other crop species. Despite the low rate of nucleotide polymorphisms, a wide range of heritable phenotypic variation exists. There is even evidence for heritable phenotypic variation among individuals within some cultivars. Williams 82, the soybean cultivar used to produce the reference genome sequence, was derived from backcrossing a Phytophthora root rot resistance locus from the donor parent Kingwa into the recurrent parent Williams. To explore the genetic basis of intracultivar variation, we investigated the nucleotide, structural, and gene content variation of different Williams 82 individuals. Williams 82 individuals exhibited variation in the number and size of introgressed Kingwa loci. In these regions of genomic heterogeneity, the reference Williams 82 genome sequence consists of a mosaic of Williams and Kingwa haplotypes. Genomic structural variation between Williams and Kingwa was maintained between the Williams 82 individuals within the regions of heterogeneity. Additionally, the regions of heterogeneity exhibited gene content differences between Williams 82 individuals. These findings show that genetic heterogeneity in Williams 82 primarily originated from the differential segregation of polymorphic chromosomal regions following the backcross and single-seed descent generations of the breeding process. We conclude that soybean haplotypes can possess a high rate of structural and gene content variation, and the impact of intracultivar genetic heterogeneity may be significant. This detailed characterization will be useful for interpreting soybean genomic data sets and highlights important considerations for research communities that are developing or utilizing a reference genome sequence.

Complementary genetic and genomic approaches help characterize the linkage group I seed protein QTL in soybean

BackgroundThe nutritional and economic value of many crops is effectively a function of seed protein and oil content. Insight into the genetic and molecular control mechanisms involved in the deposition of these constituents in the developing seed is needed to guide crop improvement. A quantitative trait locus (QTL) on Linkage Group I (LG I) of soybean (Glycine max (L.) Merrill) has a striking effect on seed protein content.ResultsA soybean near-isogenic line (NIL) pair contrasting in seed protein and differing in an introgressed genomic segment containing the LG I protein QTL was used as a resource to demarcate the QTL region and to study variation in transcript abundance in developing seed. The LG I QTL region was delineated to less than 8.4 Mbp of genomic sequence on chromosome 20. Using Affymetrix® Soy GeneChip and high-throughput Illumina® whole transcriptome sequencing platforms, 13 genes displaying significant seed transcript accumulation differences between NILs were identified that mapped to the 8.4 Mbp LG I protein QTL region.ConclusionsThis study identifies gene candidates at the LG I protein QTL for potential involvement in the regulation of protein content in the soybean seed. The results demonstrate the power of complementary approaches to characterize contrasting NILs and provide genome-wide transcriptome insight towards understanding seed biology and the soybean genome.

Uncovering the Arabidopsis thaliana nectary transcriptome: investigation of differential gene expression in floral nectariferous tissues

BackgroundMany flowering plants attract pollinators by offering a reward of floral nectar. Remarkably, the molecular events involved in the development of nectaries, the organs that produce nectar, as well as the synthesis and secretion of nectar itself, are poorly understood. Indeed, to date, no genes have been shown to directly affect the de novo production or quality of floral nectar. To address this gap in knowledge, the ATH1 Affymetrix® GeneChip array was used to systematically investigate the Arabidopsis nectary transcriptome to identify genes and pathways potentially involved in nectar production.ResultsIn this study, we identified a large number of genes differentially expressed between secretory lateral nectaries and non-secretory median nectary tissues, as well as between mature lateral nectaries (post-anthesis) and immature lateral nectaries (pre-anthesis). Expression within nectaries was also compared to thirteen non-nectary reference tissues, from which 270 genes were identified as being significantly upregulated in nectaries. The expression patterns of 14 nectary-enriched genes were also confirmed via RT PCR. Upon looking into functional groups of upregulated genes, pathways involved in gene regulation, carbohydrate metabolism, and lipid metabolism were particularly enriched in nectaries versus reference tissues.ConclusionA large number of genes preferentially expressed in nectaries, as well as between nectary types and developmental stages, were identified. Several hypotheses relating to mechanisms of nectar production and regulation thereof are proposed, and provide a starting point for reverse genetics approaches to determine molecular mechanisms underlying nectar synthesis and secretion.

An Integrative Approach to Genomic Introgression Mapping

Near-isogenic lines (NILs) are valuable genetic resources for many crop species, including soybean (Glycine max). The development of new molecular platforms promises to accelerate the mapping of genetic introgressions in these materials. Here, we compare some existing and emerging methodologies for genetic introgression mapping: single-feature polymorphism analysis, Illumina GoldenGate single nucleotide polymorphism (SNP) genotyping, and de novo SNP discovery via RNA-Seq analysis of next-generation sequence data. We used these methods to map the introgressed regions in an iron-inefficient soybean NIL and found that the three mapping approaches are complementary when utilized in combination. The comparative RNA-Seq approach offers several additional advantages, including the greatest mapping resolution, marker depth, and de novo marker utility for downstream fine-mapping analysis. We applied the comparative RNA-Seq method to map genetic introgressions in an additional pair of NILs exhibiting differential seed protein content. Furthermore, we attempted to optimize the comparative RNA-Seq approach by assessing the impact of sequence depth, SNP identification methodology, and post hoc analyses on SNP discovery rates. We conclude that the comparative RNA-Seq approach can be optimized with sufficient sampling and by utilizing a post hoc correction accounting for gene density variation that controls for false discoveries.

A new recombinant mosquito salivary allergen, rAed a 2: allergenicity, clinical relevance, and cross‐reactivity

Background:  Mosquito salivary proteins cause allergic reactions in humans. The allergenicity, clinical relevance, and species cross‐reactivity of a new 37‐kDa recombinant mosquito (Aedes aegypti) salivary allergen, rAed a 2, were evaluated.

RNA-dependent dynamic histone acetylation regulates MCL1 alternative splicing

Histone deacetylases (HDACs) and lysine acetyltransferases (KATs) catalyze dynamic histone acetylation at regulatory and coding regions of transcribed genes. Highly phosphorylated HDAC2 is recruited within corepressor complexes to regulatory regions, while the nonphosphorylated form is associated with the gene body. In this study, we characterized the nonphosphorylated HDAC2 complexes recruited to the transcribed gene body and explored the function of HDAC-complex-mediated dynamic histone acetylation. HDAC1 and 2 were coimmunoprecipitated with several splicing factors, including serine/arginine-rich splicing factor 1 (SRSF1) which has roles in alternative splicing. The co-chromatin immunoprecipitation of HDAC1/2 and SRSF1 to the gene body was RNA-dependent. Inhibition of HDAC activity and knockdown of HDAC1, HDAC2 or SRSF1 showed that these proteins were involved in alternative splicing of MCL1. HDAC1/2 and KAT2B were associated with nascent pre-mRNA in general and with MCL1 pre-mRNA specifically. Inhibition of HDAC activity increased the occupancy of KAT2B and acetylation of H3 and H4 of the H3K4 methylated alternative MCL1 exon 2 nucleosome. Thus, nonphosphorylated HDAC1/2 is recruited to pre-mRNA by splicing factors to act at the RNA level with KAT2B and other KATs to catalyze dynamic histone acetylation of the MCL1 alternative exon and alter the splicing of MCL1 pre-mRNA.

Transcriptome mapping of pAR060302, a blaCMY-2-positive broad-host-range IncA/C plasmid

ABSTRACT The multidrug resistance-encoding plasmids belonging to the IncA/C incompatibility group have recently emerged among Escherichia coli and Salmonella enterica strains in the United States. These plasmids have a unique genetic structure compared to other enterobacterial plasmid types, a broad host range, and a propensity to acquire large numbers of antimicrobial resistance genes via their accessory regions. Using E. coli strain DH5α harboring the prototype IncA/C plasmid pAR060302, we sought to define the baseline transcriptome of IncA/C plasmids under laboratory growth and in the face of selective pressure. The effects of ampicillin, florfenicol, or streptomycin exposure were compared to those on cells left untreated at logarithmic phase using Illumina platform-based RNA sequencing (RNA-Seq). Under growth in Luria-Bertani broth lacking antibiotics, much of the backbone of pAR060302 was transcriptionally inactive, including its putative transfer regions. A few plasmid backbone genes of interest were highly transcribed, including genes of a putative toxin-antitoxin system and an H-NS-like transcriptional regulator. In contrast, numerous genes within the accessory regions of pAR060302 were highly transcribed, including the resistance genes floR, bla CMY-2, aadA, and aacA. Treatment with ampicillin or streptomycin resulted in no genes being differentially expressed compared to controls lacking antibiotics, suggesting that many of the resistance-associated genes are not differentially expressed due to exposure to these antibiotics. In contrast, florfenicol treatment resulted in the upregulation of floR and numerous chromosomal genes. Overall, the transcriptome mapping of pAR060302 suggests that it mitigates the fitness costs of carrying resistance-associated genes through global regulation with its transcriptional regulators.