Spontaneous mutations in maize pollen are frequent in some lines and arise mainly from retrotranspositions and deletions.

While studying spontaneous mutations at the maize bronze (bz) locus, we made the unexpected discovery that specific low-copy number retrotransposons are mobile in the pollen of some maize lines, but not of others. We conducted large-scale genetic experiments to isolate new bz mutations from several Bz stocks and recovered spontaneous stable mutations only in the pollen parent in reciprocal crosses. Most of the new stable bz mutations resulted from either insertions of low-copy number long terminal repeat (LTR) retrotransposons or deletions, the same two classes of mutations that predominated in a collection of spontaneous wx mutations [Wessler S (1997) The Mutants of Maize, pp 385-386]. Similar mutations were recovered at the closely linked sh locus. These events occurred with a frequency of 2-4 × 10-5 in two lines derived from W22 and in 4Co63, but not at all in B73 or Mo17, two inbreds widely represented in Corn Belt hybrids. Surprisingly, the mutagenic LTR retrotransposons differed in the active lines, suggesting differences in the autonomous element make-up of the lines studied. Some active retrotransposons, like Hopscotch, Magellan, and Bs2, a Bs1 variant, were described previously; others, like Foto and Focou in 4Co63, were not. By high-throughput sequencing of retrotransposon junctions, we established that retrotranposition of Hopscotch, Magellan, and Bs2 occurs genome-wide in the pollen of active lines, but not in the female germline or in somatic tissues. We discuss here the implications of these results, which shed light on the source, frequency, and nature of spontaneous mutations in maize.

Large-scale reconstruction of chromatin structures of maize temperate and tropical inbred lines.

Maize is a model plant species often used for genetics and genomics research because of its genetic diversity. There are prominent morphological, genetic, and epigenetic variations between tropical and temperate maize lines. However, the genome-wide chromatin conformations of these two maize types remain unexplored. We applied a Hi-C approach to compare the genome-wide chromatin interactions between temperate inbred line D132 and tropical line CML288. A reconstructed maize three-dimensional genome model revealed the spatial segregation of the global A and B compartments. The A compartments contain enriched genes and active epigenome marks, whereas the B compartments are gene-poor, transcriptionally silent chromatin regions. Whole-genome analyses indicated that the global A compartment content of CML288 was 3.12% lower than that of D132. Additionally, global and A/B sub-compartments were associated with differential gene expression and epigenetic changes between two inbred lines. About 25.3% of topologically associating domains (TADs) were determined to be associated with complex domain-level modifications that induced transcriptional changes, indicative of a large-scale reorganization of chromatin structures between the inbred maize lines. Furthermore, differences in chromatin interactions between the two lines correlated with epigenetic changes. These findings provide a solid foundation for the wider plant community to further investigate the genome-wide chromatin structures in other plant species.

Identification of ZmNF-YC2 and its regulatory network for maize flowering time.

Flowering time is an important agronomic trait that determines the distribution and adaptation of plants. The accurate prediction of flowering time in elite germplasm is critical for maize breeding. However, the molecular mechanisms underlying the photoperiod response remain elusive in maize. Here we cloned the flowering time-controlling gene, ZmNF-YC2, by map-based cloning and confirmed that ZmNF-YC2 is the nuclear transcription factor Y subunit C-2 protein and a positive regulator of flowering time in maize under long-day conditions. Our results show that ZmNF-YC2 promotes the expression of ZmNF-YA3. ZmNF-YA3 negatively regulates the transcription of ZmAP2. ZmAP2 suppresses the expression of ZMM4 to delay flowering time. We then developed a gene regulatory model of flowering time in maize using ZmNF-YC2, ZmNF-YA3, ZmAP2, ZMM4, and other key genes. The cascading regulation by ZmNF-YC2 of maize flowering time has not been reported in other species.

The wild sweetpotato (Ipomoea trifida) genome provides insights into storage root development.

BACKGROUND: Sweetpotato (Ipomoea batatas (L.) Lam.) is the seventh most important crop in the world and is mainly cultivated for its underground storage root (SR). The genetic studies of this species have been hindered by a lack of high-quality reference sequence due to its complex genome structure. Diploid Ipomoea trifida is the closest relative and putative progenitor of sweetpotato, which is considered a model species for sweetpotato, including genetic, cytological, and physiological analyses. RESULTS: Here, we generated the chromosome-scale genome sequence of SR-forming diploid I. trifida var. Y22 with high heterozygosity (2.20%). Although the chromosome-based synteny analysis revealed that the I. trifida shared conserved karyotype with Ipomoea nil after the separation, I. trifida had a much smaller genome than I. nil due to more efficient eliminations of LTR-retrotransposons and lack of species-specific amplification bursts of LTR-RTs. A comparison with four non-SR-forming species showed that the evolution of the beta-amylase gene family may be related to SR formation. We further investigated the relationship of the key gene BMY11 (with identity 47.12% to beta-amylase 1) with this important agronomic trait by both gene expression profiling and quantitative trait locus (QTL) mapping. And combining SR morphology and structure, gene expression profiling and qPCR results, we deduced that the products of the activity of BMY11 in splitting starch granules and be recycled to synthesize larger granules, contributing to starch accumulation and SR swelling. Moreover, we found the expression pattern of BMY11, sporamin proteins and the key genes involved in carbohydrate metabolism and stele lignification were similar to that of sweetpotato during the SR development. CONCLUSIONS: We constructed the high-quality genome reference of the highly heterozygous I. trifida through a combined approach and this genome enables a better resolution of the genomics feature and genome evolutions of this species. Sweetpotato SR development genes can be identified in I. trifida and these genes perform similar functions and patterns, showed that the diploid I. trifida var. Y22 with typical SR could be considered an ideal model for the studies of sweetpotato SR development.

Multiple EIT and EIA in optical microresonators.

Multiple-path interference plays a fundamental role in classical and quantum physics. In this work, we propose two general schemes to realize multiple electromagnetically induced transparency (EIT) and electromagnetically induced absorption (EIA) in coupled microresonators and optomechanical systems. We give explicit physical descriptions and find out that these two schemes are essentially equivalent to each other. More importantly, we experimentally demonstrate both multiple EIT and EIA by coupling a microtoroid to a microsphere that supports multiple high Q optical modes with dense modes distributions. The theory fits well with the experimental results. We believe that our study and experimental results lay a foundation for realizing arbitrary multiple pathways interference in applications.

Highly interwoven communities of a gene regulatory network unveil topologically important genes for maize seed development.

The complex interactions between transcription factors (TFs) and their target genes in a spatially and temporally specific manner are crucial to all cellular processes. Reconstruction of gene regulatory networks (GRNs) from gene expression profiles can help to decipher TF-gene regulations in a variety of contexts; however, the inevitable prediction errors of GRNs hinder optimal data mining of RNA-Seq transcriptome profiles. Here we perform an integrative study of Zea mays (maize) seed development in order to identify key genes in a complex developmental process. First, we reverse engineered a GRN from 78 maize seed transcriptome profiles. Then, we studied collective gene interaction patterns and uncovered highly interwoven network communities as the building blocks of the GRN. One community, composed of mostly unknown genes interacting with opaque2, brittle endosperm1 and shrunken2, contributes to seed phenotypes. Another community, composed mostly of genes expressed in the basal endosperm transfer layer, is responsible for nutrient transport. We further integrated our inferred GRN with gene expression patterns in different seed compartments and at various developmental stages and pathways. The integration facilitated a biological interpretation of the GRN. Our yeast one-hybrid assays verified six out of eight TF-promoter bindings in the reconstructed GRN. This study identified topologically important genes in interwoven network communities that may be crucial to maize seed development.

Rolling-circle amplification of centromeric Helitrons in plant genomes.

The unusual eukaryotic Helitron transposons can readily capture host sequences and are, thus, evolutionarily important. They are presumed to amplify by rolling-circle replication (RCR) because some elements encode predicted proteins homologous to RCR prokaryotic transposases. In support of this replication mechanism, it was recently shown that transposition of a bat Helitron generates covalently closed circular intermediates. Another strong prediction is that RCR should generate tandem Helitron concatemers, yet almost all Helitrons identified to date occur as solo elements in the genome. To investigate alternative modes of Helitron organization in present-day genomes, we have applied the novel computational tool HelitronScanner to 27 plant genomes and have uncovered numerous tandem arrays of partially decayed, truncated Helitrons in all of them. Strikingly, most of these Helitron tandem arrays are interspersed with other repeats in centromeres. Many of these arrays have multiple Helitron 5' ends, but a single 3' end. The number of repeats in any one array can range from a handful to several hundreds. We propose here an RCR model that conforms to the present Helitron landscape of plant genomes. Our study provides strong evidence that plant Helitrons amplify by RCR and that the tandemly arrayed replication products accumulate mostly in centromeres.

HelitronScanner uncovers a large overlooked cache of Helitron transposons in many plant genomes.

Transposons make up the bulk of eukaryotic genomes, but are difficult to annotate because they evolve rapidly. Most of the unannotated portion of sequenced genomes is probably made up of various divergent transposons that have yet to be categorized. Helitrons are unusual rolling circle eukaryotic transposons that often capture gene sequences, making them of considerable evolutionary importance. Unlike other DNA transposons, Helitrons do not end in inverted repeats or create target site duplications, so they are particularly challenging to identify. Here we present HelitronScanner, a two-layered local combinational variable (LCV) tool for generalized Helitron identification that represents a major improvement over previous identification programs based on DNA sequence or structure. HelitronScanner identified 64,654 Helitrons from a wide range of plant genomes in a highly automated way. We tested HelitronScanner's predictive ability in maize, a species with highly heterogeneous Helitron elements. LCV scores for the 5' and 3' termini of the predicted Helitrons provide a primary confidence level and element copy number provides a secondary one. Newly identified Helitrons were validated by PCR assays or by in silico comparative analysis of insertion site polymorphism among multiple accessions. Many new Helitrons were identified in model species, such as maize, rice, and Arabidopsis, and in a variety of organisms where Helitrons had not been reported previously to our knowledge, leading to a major upward reassessment of their abundance in plant genomes. HelitronScanner promises to be a valuable tool in future comparative and evolutionary studies of this major transposon superfamily.

iTRAQ-Based Proteomics Analysis and Network Integration for Kernel Tissue Development in Maize.

Grain weight is one of the most important yield components and a developmentally complex structure comprised of two major compartments (endosperm and pericarp) in maize (Zea mays L.), however, very little is known concerning the coordinated accumulation of the numerous proteins involved. Herein, we used isobaric tags for relative and absolute quantitation (iTRAQ)-based comparative proteomic method to analyze the characteristics of dynamic proteomics for endosperm and pericarp during grain development. Totally, 9539 proteins were identified for both components at four development stages, among which 1401 proteins were non-redundant, 232 proteins were specific in pericarp and 153 proteins were specific in endosperm. A functional annotation of the identified proteins revealed the importance of metabolic and cellular processes, and binding and catalytic activities for the tissue development. Three and 76 proteins involved in 49 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were integrated for the specific endosperm and pericarp proteins, respectively, reflecting their complex metabolic interactions. In addition, four proteins with important functions and different expression levels were chosen for gene cloning and expression analysis. Different concordance between mRNA level and the protein abundance was observed across different proteins, stages, and tissues as in previous research. These results could provide useful message for understanding the developmental mechanisms in grain development in maize.

Three-pathway electromagnetically induced transparency in coupled-cavity optomechanical system.

Recently Qu and Agarwal [Phys. Rev. A 22, 031802 (2013)] found a three-pathway electromagnetically induced absorption (TEIA) phenomenon within a mechanically coupled two-cavity system, where there exist a sharp EIA dip in the broad electromagnetically induced transparency peak in the transmission spectrum. In this work, we study the response of a probe light in a pair of directly coupled microcavities with one mechanical mode. We find that in addition to the sharp TEIA dip within a broad EIT window as found by Qu and Agarwal, three-pathway electromagnetically induced transparency (TEIT) within the broad EIT window could also exist under certain conditions. We give explicit physical explanations and detailed calculations. Our results provide a method for controlling transition between TEIA and TEIT in coupled optomechanical systems, and reveal the multiple pathways interference is versatile for controlling light.

InsertionMapper: a pipeline tool for the identification of targeted sequences from multidimensional high throughput sequencing data.

BACKGROUND: The advent of next-generation high-throughput technologies has revolutionized whole genome sequencing, yet some experiments require sequencing only of targeted regions of the genome from a very large number of samples. These regions can be amplified by PCR and sequenced by next-generation methods using a multidimensional pooling strategy. However, there is at present no available generalized tool for the computational analysis of target-enriched NGS data from multidimensional pools. RESULTS: Here we present InsertionMapper, a pipeline tool for the identification of targeted sequences from multidimensional high throughput sequencing data. InsertionMapper consists of four independently working modules: Data Preprocessing, Database Modeling, Dimension Deconvolution and Element Mapping. We illustrate InsertionMapper with an example from our project 'New reverse genetics resources for maize', which aims to sequence-index a collection of 15,000 independent insertion sites of the transposon Ds in maize. Identified sequences are validated by PCR assays. This pipeline tool is applicable to similar scenarios requiring analysis of the tremendous output of short reads produced in NGS sequencing experiments of targeted genome sequences. CONCLUSIONS: InsertionMapper is proven efficacious for the identification of target-enriched sequences from multidimensional high throughput sequencing data. With adjustable parameters and experiment configurations, this tool can save great computational effort to biologists interested in identifying their sequences of interest within the huge output of modern DNA sequencers. InsertionMapper is freely accessible at https://sourceforge.net/p/insertionmapper and http://bo.csam.montclair.edu/du/insertionmapper.

Genome resources for the elite bread wheat cultivar Aikang 58 and mining of elite homeologous haplotypes for accelerating wheat improvement.

Despite recent progress in crop genomics studies, the genomic changes brought about by modern breeding selection are still poorly understood, thus hampering genomics-assisted breeding, especially in polyploid crops with compound genomes such as common wheat (Triticum aestivum). In this work, we constructed genome resources for the modern elite common wheat variety Aikang 58 (AK58). Comparative genomics between AK58 and the landrace cultivar Chinese Spring (CS) shed light on genomic changes that occurred through recent varietal improvement. We also explored subgenome diploidization and divergence in common wheat and developed a homoeologous locus-based genome-wide association study (HGWAS) approach, which was more effective than single homoeolog-based GWAS in unraveling agronomic trait-associated loci. A total of 123 major HGWAS loci were detected using a genetic population derived from AK58 and CS. Elite homoeologous haplotypes (HHs), formed by combinations of subgenomic homoeologs of the associated loci, were found in both parents and progeny, and many could substantially improve wheat yield and related traits. We built a website where users can download genome assembly sequence and annotation data for AK58, perform blast analysis, and run JBrowse. Our work enriches genome resources for wheat, provides new insights into genomic changes during modern wheat improvement, and suggests that efficient mining of elite HHs can make a substantial contribution to genomics-assisted breeding in common wheat and other polyploid crops.

The polychromatic Helitron landscape of the maize genome.

Maize Helitron transposons are intriguing because of their notable ability to capture gene fragments and move them around the genome. To document more extensively their variability and their contribution to the remarkable genome structure variation of present-day maize, we have analyzed their composition, copy number, timing of insertion, and chromosomal distribution. First, we searched 2.4 Gb of sequences generated by the Maize Genome Sequencing Project with our HelitronFinder program. We identified 2,791 putative nonautonomous Helitrons and manually curated a subset of 272. The predicted Helitrons measure 11.9 kb on average and carry from zero to nine gene fragments, captured from 376 different genes. Although the diversity of Helitron gene fragments in maize is greater than in other species, more than one-third of annotated Helitrons carry fragments derived from just one of two genes. Most members in these two subfamilies inserted in the genome less than one million years ago. Second, we conducted a BLASTN search of the maize sequence database with queries from two previously described agenic Helitrons not detected by HelitronFinder. Two large subfamilies of Helitrons or Helitron-related transposons were identified. One subfamily, termed Cornucopious, consists of thousands of copies of an approximately 1.0-kb agenic Helitron that may be the most abundant transposon in maize. The second subfamily consists of >150 copies of a transposon-like sequence, termed Heltir, that has terminal inverted repeats resembling Helitron 3' termini. Nonautonomous Helitrons make up at least 2% of the maize genome and most of those tested show +/- polymorphisms among modern inbred lines.

Student Attitudes Contribute to the Effectiveness of a Genomics CURE.

The Genomics Education Partnership (GEP) engages students in a course-based undergraduate research experience (CURE). To better understand the student attributes that support success in this CURE, we asked students about their attitudes using previously published scales that measure epistemic beliefs about work and science, interest in science, and grit. We found, in general, that the attitudes students bring with them into the classroom contribute to two outcome measures, namely, learning as assessed by a pre- and postquiz and perceived self-reported benefits. While the GEP CURE produces positive outcomes overall, the students with more positive attitudes toward science, particularly with respect to epistemic beliefs, showed greater gains. The findings indicate the importance of a student's epistemic beliefs to achieving positive learning outcomes.

The complete Ac/Ds transposon family of maize.

BACKGROUND: The nonautonomous maize Ds transposons can only move in the presence of the autonomous element Ac. They comprise a heterogeneous group that share 11-bp terminal inverted repeats (TIRs) and some subterminal repeats, but vary greatly in size and composition. Three classes of Ds elements can cause mutations: Ds-del, internal deletions of the 4.6-kb Ac element; Ds1, ~400-bp in size and sharing little homology with Ac, and Ds2, variably-sized elements containing about 0.5 kb from the Ac termini and unrelated internal sequences. Here, we analyze the entire complement of Ds-related sequences in the genome of the inbred B73 and ask whether additional classes of Ds-like (Ds-l) elements, not uncovered genetically, are mobilized by Ac. We also compare the makeup of Ds-related sequences in two maize inbreds of different origin. RESULTS: We found 903 elements with 11-bp Ac/Ds TIRs flanked by 8-bp target site duplications. Three resemble Ac, but carry small rearrangements. The others are much shorter, once extraneous insertions are removed. There are 331 Ds1 and 39 Ds2 elements, many of which are likely mobilized by Ac, and two novel classes of Ds-l elements. Ds-l3 elements lack subterminal homology with Ac, but carry transposase gene fragments, and represent decaying Ac elements. There are 44 such elements in B73. Ds-l4 elements share little similarity with Ac outside of the 11-bp TIR, have a modal length of ~1 kb, and carry filler DNA which, in a few cases, could be matched to gene fragments. Most Ds-related elements in B73 (486/903) fall in this class. None of the Ds-l elements tested responded to Ac. Only half of Ds insertion sites examined are shared between the inbreds B73 and W22. CONCLUSIONS: The majority of Ds-related sequences in maize correspond to Ds-l elements that do not transpose in the presence of Ac. Unlike actively transposing elements, many Ds-l elements are inserted in repetitive DNA, where they probably become methylated and begin to decay. The filler DNA present in most elements is occasionally captured from genes, a rare feature in transposons of the hAT superfamily to which Ds belongs. Maize inbreds of different origin are highly polymorphic in their DNA transposon makeup.

Long-read sequencing reveals genomic structural variations that underlie creation of quality protein maize.

Mutation of o2 doubles maize endosperm lysine content, but it causes an inferior kernel phenotype. Developing quality protein maize (QPM) by introgressing o2 modifiers (Mo2s) into the o2 mutant benefits millions of people in developing countries where maize is a primary protein source. Here, we report genome sequence and annotation of a South African QPM line K0326Y, which is assembled from single-molecule, real-time shotgun sequencing reads collinear with an optical map. We achieve a N50 contig length of 7.7 million bases (Mb) directly from long-read assembly, compared to those of 1.04 Mb for B73 and 1.48 Mb for Mo17. To characterize Mo2s, we map QTLs to chromosomes 1, 6, 7, and 9 using an F2 population derived from crossing K0326Y and W64Ao2. RNA-seq analysis of QPM and o2 endosperms reveals a group of differentially expressed genes that coincide with Mo2 QTLs, suggesting a potential role in vitreous endosperm formation.

Facilitating Growth through Frustration: Using Genomics Research in a Course-Based Undergraduate Research Experience.

A hallmark of the research experience is encountering difficulty and working through those challenges to achieve success. This ability is essential to being a successful scientist, but replicating such challenges in a teaching setting can be difficult. The Genomics Education Partnership (GEP) is a consortium of faculty who engage their students in a genomics Course-Based Undergraduate Research Experience (CURE). Students participate in genome annotation, generating gene models using multiple lines of experimental evidence. Our observations suggested that the students' learning experience is continuous and recursive, frequently beginning with frustration but eventually leading to success as they come up with defendable gene models. In order to explore our "formative frustration" hypothesis, we gathered data from faculty via a survey, and from students via both a general survey and a set of student focus groups. Upon analyzing these data, we found that all three datasets mentioned frustration and struggle, as well as learning and better understanding of the scientific process. Bioinformatics projects are particularly well suited to the process of iteration and refinement because iterations can be performed quickly and are inexpensive in both time and money. Based on these findings, we suggest that a dynamic of "formative frustration" is an important aspect for a successful CURE.

Computational prediction and molecular confirmation of Helitron transposons in the maize genome.

BACKGROUND: Helitrons represent a new class of transposable elements recently uncovered in plants and animals. One remarkable feature of Helitrons is their ability to capture gene sequences, which makes them of considerable potential evolutionary importance. However, because Helitrons lack the typical structural features of other DNA transposable elements, identifying them is a challenge. Currently, most researchers identify Helitrons manually by comparing sequences. With the maize whole genome sequencing project underway, an automated computational Helitron searching tool is needed. The characterization of Helitron activities in maize needs to be addressed in order to better understand the impact of Helitrons on the organization of the genome. RESULTS: We developed and implemented a heuristic searching algorithm in PERL for identifying Helitrons. Our HelitronFinder program will (i) take FASTA-formatted DNA sequences as input and identify the hairpin looping patterns, and (ii) exploit the consensus 5' and 3' end sequences of known Helitrons to identify putative ends. We randomly selected five predicted Helitrons from the program's high quality output for molecular verification. Four out of the five predicted Helitrons were confirmed by PCR assays and DNA sequencing in different maize inbred lines. The HelitronFinder program identified two head-to-head dissimilar Helitrons in a maize BAC sequence. CONCLUSION: We have identified 140 new Helitron candidates in maize with our computational tool HelitronFinder by searching maize DNA sequences currently available in GenBank. Four out of five candidates were confirmed to be real by empirical methods, thus validating the predictions of HelitronFinder. Additional points to emerge from our study are that Helitrons do not always insert at an AT dinucleotide in the host sequences, that they can insert immediately adjacent to an existing Helitron, and that their movement may cause changes in the flanking region, such as deletions.

The maize W22 genome provides a foundation for functional genomics and transposon biology.

The maize W22 inbred has served as a platform for maize genetics since the mid twentieth century. To streamline maize genome analyses, we have sequenced and de novo assembled a W22 reference genome using short-read sequencing technologies. We show that significant structural heterogeneity exists in comparison to the B73 reference genome at multiple scales, from transposon composition and copy number variation to single-nucleotide polymorphisms. The generation of this reference genome enables accurate placement of thousands of Mutator (Mu) and Dissociation (Ds) transposable element insertions for reverse and forward genetics studies. Annotation of the genome has been achieved using RNA-seq analysis, differential nuclease sensitivity profiling and bisulfite sequencing to map open reading frames, open chromatin sites and DNA methylation profiles, respectively. Collectively, the resources developed here integrate W22 as a community reference genome for functional genomics and provide a foundation for the maize pan-genome.

Integrated regulatory network reveals the early salt tolerance mechanism of Populus euphratica.

Salinization is one of the major factors that threaten the existence of plants worldwide. Populus euphratica has been deemed to be a promising candidate for stress response research because of its high capacity to tolerate extreme salt stress. We carried out a genome-wide transcriptome analysis to identify the differentially expressed genes (DEGs) response to salt shock and elucidate the early salt tolerance mechanisms in P. euphratica. Both hierarchical clustering and DEG analysis demonstrated a predominant variation from time-course rather than NaCl intensity within 24 hours salt shock. Among the identified 1,678 salt-responsive DEGs, 74.1% (1,244) have not been reported before. We further created an integrated regulatory gene network of the salt response in P. euphratica by combining DEGs, transcription factors (TFs), Helitrons, miRNAs and their targets. The prominent pathways in this network are plant hormone transduction, starch and sucrose metabolism, RNA transport, protein processing in endoplasmic reticulum, etc. In addition, the network indicates calcium-related genes play key roles in P. euphratica response to salt shock. These results illustrated an overview of the systematic molecular response in P. euphratica under different intensities of salt shock and revealed the complex regulatory mechanism.