Evolutionary Dynamics of Chromatin Structure and Duplicate Gene Expression in Diploid and Allopolyploid Cotton.

Polyploidy is a prominent mechanism of plant speciation and adaptation, yet the mechanistic understandings of duplicated gene regulation remain elusive. Chromatin structure dynamics are suggested to govern gene regulatory control. Here, we characterized genome-wide nucleosome organization and chromatin accessibility in allotetraploid cotton, Gossypium hirsutum (AADD, 2n = 4X = 52), relative to its two diploid parents (AA or DD genome) and their synthetic diploid hybrid (AD), using DNS-seq. The larger A-genome exhibited wider average nucleosome spacing in diploids, and this intergenomic difference diminished in the allopolyploid but not hybrid. Allopolyploidization also exhibited increased accessibility at promoters genome-wide and synchronized cis-regulatory motifs between subgenomes. A prominent cis-acting control was inferred for chromatin dynamics and demonstrated by transposable element removal from promoters. Linking accessibility to gene expression patterns, we found distinct regulatory effects for hybridization and later allopolyploid stages, including nuanced establishment of homoeolog expression bias and expression level dominance. Histone gene expression and nucleosome organization are coordinated through chromatin accessibility. Our study demonstrates the capability to track high-resolution chromatin structure dynamics and reveals their role in the evolution of cis-regulatory landscapes and duplicate gene expression in polyploids, illuminating regulatory ties to subgenomic asymmetry and dominance.

Integrating SNVs and CNAs on a phylogenetic tree from single-cell DNA sequencing data.

Single-cell DNA sequencing enables the construction of evolutionary trees that can reveal how tumors gain mutations and grow. Different whole-genome amplification procedures render genomic materials of different characteristics, often suitable for the detection of either single-nucleotide variation or copy number aberration, but not ideally for both. Consequently, this hinders the inference of a comprehensive phylogenetic tree and limits opportunities to investigate the interplay of SNVs and CNAs. Existing methods such as SCARLET and COMPASS require that the SNVs and CNAs are detected from the same sets of cells, which is technically challenging. Here we present a novel computational tool, SCsnvcna, that places SNVs on a tree inferred from CNA signals, whereas the sets of cells rendering the SNVs and CNAs are independent, offering a more practical solution in terms of the technical challenges. SCsnvcna is a Bayesian probabilistic model using both the genotype constraints on the tree and the cellular prevalence to search the optimal solution. Comprehensive simulations and comparison with seven state-of-the-art methods show that SCsnvcna is robust and accurate in a variety of circumstances. Particularly, SCsnvcna most frequently produces the lowest error rates, with ability to scale to a wide range of numerical values for leaf nodes in the tree, SNVs, and SNV cells. The application of SCsnvcna to two published colorectal cancer data sets shows highly consistent placement of SNV cells and SNVs with the original study while also supporting a refined placement of ATP7B, illustrating SCsnvcna's value in analyzing complex multitumor samples.

Maize (Zea mays L.) Nucleoskeletal Proteins Regulate Nuclear Envelope Remodeling and Function in Stomatal Complex Development and Pollen Viability.

In eukaryotes, the nuclear envelope (NE) encloses chromatin and separates it from the rest of the cell. The Linker of Nucleoskeleton and Cytoskeleton (LINC) complex physically bridges across the NE, linking nuclear and cytoplasmic components. In plants, these LINC complexes are beginning to be ascribed roles in cellular and nuclear functions, including chromatin organization, regulation of nuclei shape and movement, and cell division. Homologs of core LINC components, KASH and SUN proteins, have previously been identified in maize. Here, we characterized the presumed LINC-associated maize nucleoskeletal proteins NCH1 and NCH2, homologous to members of the plant NMCP/CRWN family, and MKAKU41, homologous to AtKAKU4. All three proteins localized to the nuclear periphery when transiently and heterologously expressed as fluorescent protein fusions in Nicotiana benthamiana. Overexpression of MKAKU41 caused dramatic changes in the organization of the nuclear periphery, including nuclear invaginations that stained positive for non-nucleoplasmic markers of the inner and outer NE membranes, and the ER. The severity of these invaginations was altered by changes in LINC connections and the actin cytoskeleton. In maize, MKAKU41 appeared to share genetic functions with other LINC components, including control of nuclei shape, stomatal complex development, and pollen viability. Overall, our data show that NCH1, NCH2, and MKAKU41 have characteristic properties of LINC-associated plant nucleoskeletal proteins, including interactions with NE components suggestive of functions at the nuclear periphery that impact the overall nuclear architecture.

iSeg: an efficient algorithm for segmentation of genomic and epigenomic data.

BACKGROUND: Identification of functional elements of a genome often requires dividing a sequence of measurements along a genome into segments where adjacent segments have different properties, such as different mean values. Despite dozens of algorithms developed to address this problem in genomics research, methods with improved accuracy and speed are still needed to effectively tackle both existing and emerging genomic and epigenomic segmentation problems. RESULTS: We designed an efficient algorithm, called iSeg, for segmentation of genomic and epigenomic profiles. iSeg first utilizes dynamic programming to identify candidate segments and test for significance. It then uses a novel data structure based on two coupled balanced binary trees to detect overlapping significant segments and update them simultaneously during searching and refinement stages. Refinement and merging of significant segments are performed at the end to generate the final set of segments. By using an objective function based on the p-values of the segments, the algorithm can serve as a general computational framework to be combined with different assumptions on the distributions of the data. As a general segmentation method, it can segment different types of genomic and epigenomic data, such as DNA copy number variation, nucleosome occupancy, nuclease sensitivity, and differential nuclease sensitivity data. Using simple t-tests to compute p-values across multiple datasets of different types, we evaluate iSeg using both simulated and experimental datasets and show that it performs satisfactorily when compared with some other popular methods, which often employ more sophisticated statistical models. Implemented in C++, iSeg is also very computationally efficient, well suited for large numbers of input profiles and data with very long sequences. CONCLUSIONS: We have developed an efficient general-purpose segmentation tool and showed that it had comparable or more accurate results than many of the most popular segment-calling algorithms used in contemporary genomic data analysis. iSeg is capable of analyzing datasets that have both positive and negative values. Tunable parameters allow users to readily adjust the statistical stringency to best match the biological nature of individual datasets, including widely or sparsely mapped genomic datasets or those with non-normal distributions.

The maize (Zea mays L.) nucleoside diphosphate kinase1 (ZmNDPK1) gene encodes a human NM23-H2 homologue that binds and stabilizes G-quadruplex DNA.

Noncanonical forms of DNA like the guanine quadruplex (G4) play important roles in regulating transcription and translation through interactions with their protein partners. Although potential G4 elements have been identified in or near genes from species diverse as bacteria, mammals, and plants, little is known about how they might function as cis-regulatory elements or as binding sites for trans-acting protein partners. In fact, until now no G4 binding partners have been identified in the plant kingdom. Here, we report on the cloning and characterization of the first plant-kingdom gene known to encode a G4-binding protein, maize (Zea mays L.) nucleoside diphosphate kinase1 (ZmNDPK1). Structural characterization by X-ray crystallography reveals that it is a homohexamer, akin to other known NDPKs like the human homologue NM23-H2. Further probing into the G4-binding properties of both NDPK homologues suggests that ZmNDPK1 possesses properties distinct from that of NM23-H2, which is known to interact with a G-rich sequence element upstream of the c-myc gene and, in doing so, modulate its expression. Indeed, ZmNDPK1 binds the folded G4 with low nanomolar affinity but corresponding unfolded G-rich DNA more weakly, whereas NM23-H2 binds both folded and unfolded G4 with low nanomolar affinities; nonetheless, both homologues appear to stabilize folded DNAs whether they were prefolded or not. We also demonstrate that the G4-binding activity of ZmNDPK1 is independent of nucleotide binding and kinase activity, suggesting that the G4-binding region and the enzyme active sites are separate. Together, these findings establish a broad evolutionary conservation of some NDPKs as G4-DNA binding enzymes, but with potentially distinct biochemical properties that may reflect divergent evolution or species-specific deployment of these elements in gene regulatory processes.

QTL Mapping and Candidate Gene Analysis of Telomere Length Control Factors in Maize (Zea mays L.).

Telomere length is a quantitative trait important for many cellular functions. Failure to regulate telomere length contributes to genomic instability, cellular senescence, cancer, and apoptosis in humans, but the functional significance of telomere regulation in plants is much less well understood. To gain a better understanding of telomere biology in plants, we used quantitative trait locus (QTL) mapping to identify genetic elements that control telomere length variation in maize (Zea mays L.). For this purpose, we measured the median and mean telomere lengths from 178 recombinant inbred lines of the IBM mapping population and found multiple regions that collectively accounted for 33-38% of the variation in telomere length. Two-way analysis of variance revealed interaction between the quantitative trait loci at genetic bin positions 2.09 and 5.04. Candidate genes within these and other significant QTL intervals, along with select genes known a priori to regulate telomere length, were tested for correlations between expression levels and telomere length in the IBM population and diverse inbred lines by quantitative real-time PCR. A slight but significant positive correlation between expression levels and telomere length was observed for many of the candidate genes, but Ibp2 was a notable exception, showing instead a negative correlation. A rad51-like protein (TEL-MD_5.04) was strongly supported as a candidate gene by several lines of evidence. Our results highlight the value of QTL mapping plus candidate gene expression analysis in a genetically diverse model system for telomere research.

Maize ribosome-inactivating proteins (RIPs) with distinct expression patterns have similar requirements for proenzyme activation.

Ribosome-inactivating proteins (RIPs, EC 3.2.2.22) are potent naturally occurring toxins found in numerous and diverse plant species. The maize RIP is unusual among the plant RIPs because it is synthesized as an inactive precursor (also known as maize proRIP1 or b-32). The proenzyme undergoes proteolytic activation that results in the removal of the NH(2)-terminal, the COOH-terminal, and internal sequences to form a two-chain holoenzyme capable of irreversibly modifying the large rRNA. The characterization of a second maize RIP (RIP2), encoded by the gene designated Rip3:2 is described here. Low levels of Rip3:2 RNA were detected in roots, shoots, tassels, silks, and leaves, but the Rip3:2 gene, unlike the Rip3:1 gene, is not under the control of the transcriptional activator Opaque-2. Instead, its expression was up-regulated by drought. Rip3:2 encodes a 31.1 kDa polypeptide that is very similar to proRIP1 in regions corresponding to those found in the active protein and the NH(2)-terminal extension. A 19-amino-acid internal portion of proRIP2 has little similarity to the proRIP1 sequence except that both are very rich in acidic residues. RIP activity assays revealed that Rip3:2 encodes a polypeptide that acquires RNA-specific N-glycosidase activity after proteolytic cleavage. Accumulation as inactive proenzymes may therefore be a general feature of maize RIPs. Differential regulation of the two RIP genes suggests that the corresponding proteins may be involved in defence-related functions with one being regulated developmentally and the other being responsive to an environmental stimulus.

The Yin Yang-1 (YY1) protein undergoes a DNA-replication-associated switch in localization from the cytoplasm to the nucleus at the onset of S phase.

The essential Yin Yang-1 gene (YY1) encodes a ubiquitous, conserved, multifunctional zinc-finger transcription factor in animals. The YY1 protein regulates initiation, activation, or repression of transcription from a variety of genes required for cell growth, development, differentiation, or tumor suppression, as well as from genes in some retroviruses and DNA viruses. Among the specific functions attributed to YY1 is a role in cell-cycle-specific upregulation of the replication-dependent histone genes. The YY1 protein binds to the histone alpha element, a regulatory sequence found in all replication-dependent histone genes. We therefore examined the abundance, DNA-binding activity and localization of the YY1 protein throughout the cell cycle in unperturbed, shake-off-synchronized Chinese hamster ovary and HeLa cells. We found that, whereas the DNA-binding activity of YY1 increased dramatically early in S phase, the YY1 mRNA and protein levels did not. YY1 changed subcellular distribution patterns during the cell cycle, from mainly cytoplasmic at G1 to mainly nuclear at early and middle S phase, then back to primarily cytoplasmic later in S phase. Nuclear accumulation of YY1 near the G1/S boundary coincided with both an increase in YY1 DNA-binding activity and the coordinate up-regulation of the replication-dependent histone genes. The DNA synthesis inhibitor aphidicolin caused a nearly complete loss of nuclear YY1, whereas addition of caffeine or 2-aminopurine to aphidicolin-treated cells restored both DNA synthesis and YY1 localization in the nucleus. These findings reveal a mechanism by which YY1 localization is coupled to DNA synthesis and responsive to cell-cycle signaling pathways. Taken together, our results provide insight into how YY1 might participate in the cell-cycle control over a variety of nuclear events required for cell division and proliferation.