Celine, a long interspersed nuclear element retrotransposon, colonizes in the centromeres of poplar chromosomes.

Centromeres in most multicellular eukaryotes are composed of long arrays of repetitive DNA sequences. Interestingly, several transposable elements, including the well-known long terminal repeat (LTR) retrotransposon CRM (centromeric retrotransposon of maize), were found to be enriched in functional centromeres marked by the centromeric histone H3 (CENH3). Here we report a centromeric long interspersed nuclear element (LINE), Celine, in Populus species. Celine has colonized preferentially in the CENH3-associated chromatin of every poplar chromosome, with 84% of the Celine elements localized in the CENH3-binding domains. By contrast, only 51% of the CRM elements were bound to CENH3 domains in Populus trichocarpa. These results suggest different centromere targeting mechanisms employed by Celine and CRM elements. Nevertheless, the high target specificity seems to be detrimental to further amplification of the Celine elements, leading to a shorter life span and patchy distribution among plant species compared to the CRM elements. Using a phylogenetically guided approach we were able to identify Celine-like LINE elements in tea plant (Camellia sinensis) and green ash tree (Fraxinus pennsylvanica). The centromeric localization of these Celine-like LINEs was confirmed in both species. We demonstrate that the centromere targeting property of Celine-like LINEs is of primitive origin and has been conserved among distantly related plant species.

Chromosome ends initiate homologous chromosome pairing during rice meiosis.

During meiotic prophase I, chromosomes undergo large-scale dynamics to allow homologous chromosome pairing, prior to which chromosome ends attach to the inner nuclear envelope and form a chromosomal bouquet. Chromosome pairing is crucial for homologous recombination and accurate chromosome segregation during meiosis. However, the specific mechanism by which homologous chromosomes recognize each other is poorly understood. Here, we investigated the process of homologous chromosome pairing during early prophase I of meiosis in rice (Oryza sativa) using pooled oligo probes specific to an entire chromosome or chromosome arm. We revealed that chromosome pairing begins from both ends and extends towards the center from early zygotene through late zygotene. Genetic analysis of both trisomy and autotetraploidy also showed that pairing initiation is induced by both ends of a chromosome. However, healed ends that lack the original terminal regions on telocentric and acrocentric chromosomes cannot initiate homologous chromosome pairing, even though they may still enter the telomere clustering region at the bouquet stage. Furthermore, a chromosome that lacks the distal parts on both sides loses the ability to pair with other intact chromosomes. Thus, the native ends of chromosomes play a crucial role in initiating homologous chromosome pairing during meiosis and likely have a substantial impact on genome differentiation.

Dynamics of accessible chromatin regions and subgenome dominance in octoploid strawberry.

Subgenome dominance has been reported in diverse allopolyploid species, where genes from one subgenome are preferentially retained and are more highly expressed than those from other subgenome(s). However, the molecular mechanisms responsible for subgenome dominance remain poorly understood. Here, we develop genome-wide map of accessible chromatin regions (ACRs) in cultivated strawberry (2n = 8x = 56, with A, B, C, D subgenomes). Each ACR is identified as an MNase hypersensitive site (MHS). We discover that the dominant subgenome A contains a greater number of total MHSs and MHS per gene than the submissive B/C/D subgenomes. Subgenome A suffers fewer losses of MHS-related DNA sequences and fewer MHS fragmentations caused by insertions of transposable elements. We also discover that genes and MHSs related to stress response have been preferentially retained in subgenome A. We conclude that preservation of genes and their cognate ACRs, especially those related to stress responses, play a major role in the establishment of subgenome dominance in octoploid strawberry.

Molecular dissection of an intronic enhancer governing cold-induced expression of the vacuolar invertase gene in potato.

Potato (Solanum tuberosum) is the third most important food crop in the world. Potato tubers must be stored at cold temperatures to minimize sprouting and losses due to disease. However, cold temperatures strongly induce the expression of the potato vacuolar invertase gene (VInv) and cause reducing sugar accumulation. This process, referred to as "cold-induced sweetening," is a major postharvest problem for the potato industry. We discovered that the cold-induced expression of VInv is controlled by a 200 bp enhancer, VInvIn2En, located in its second intron. We identified several DNA motifs in VInvIn2En that bind transcription factors involved in the plant cold stress response. Mutation of these DNA motifs abolished VInvIn2En function as a transcriptional enhancer. We developed VInvIn2En deletion lines in both diploid and tetraploid potato using clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated nuclease 9 (Cas9)-mediated gene editing. VInv transcription in cold-stored tubers was significantly reduced in the deletion lines. Interestingly, the VInvIn2En sequence is highly conserved among distantly related Solanum species, including tomato (Solanum lycopersicum) and other non-tuber-bearing species. We conclude that the VInv gene and the VInvIn2En enhancer have adopted distinct roles in the cold stress response in tubers of tuber-bearing Solanum species.

Translocations and inversions: major chromosomal rearrangements during Vigna (Leguminosae) evolution.

KEY MESSAGE: Inversions and translocations are the major chromosomal rearrangements involved in Vigna subgenera evolution, being Vigna vexillata the most divergent species. Centromeric repositioning seems to be frequent within the genus. Oligonucleotide-based fluorescence in situ hybridization (Oligo-FISH) provides a powerful chromosome identification system for inferring plant chromosomal evolution. Aiming to understand macrosynteny, chromosomal diversity, and the evolution of bean species from five Vigna subgenera, we constructed cytogenetic maps for eight taxa using oligo-FISH-based chromosome identification. We used oligopainting probes from chromosomes 2 and 3 of Phaseolus vulgaris L. and two barcode probes designed from V. unguiculata (L.) Walp. genome. Additionally, we analyzed genomic blocks among the Ancestral Phaseoleae Karyotype (APK), two V. unguiculata subspecies (V. subg. Vigna), and V. angularis (Willd.) Ohwi & Ohashi (V. subg. Ceratotropis). We observed macrosynteny for chromosomes 2, 3, 4, 6, 7, 8, 9, and 10 in all investigated taxa except for V. vexillata (L.) A. Rich (V. subg. Plectrotropis), in which only chromosomes 4, 7, and 9 were unambiguously identified. Collinearity breaks involved with chromosomes 2 and 3 were revealed. We identified minor differences in the painting pattern among the subgenera, in addition to multiple intra- and interblock inversions and intrachromosomal translocations. Other rearrangements included a pericentric inversion in chromosome 4 (V. subg. Vigna), a reciprocal translocation between chromosomes 1 and 5 (V. subg. Ceratotropis), a potential deletion in chromosome 11 of V. radiata (L.) Wilczek, as well as multiple intrablock inversions and centromere repositioning via genomic blocks. Our study allowed the visualization of karyotypic patterns in each subgenus, revealing important information for understanding intrageneric karyotypic evolution, and suggesting V. vexillata as the most karyotypically divergent species.

Classified functional mixed effects model prediction.

In nowadays biomedical research, there has been a growing demand for making accurate prediction at subject levels. In many of these situations, data are collected as longitudinal curves and display distinct individual characteristics. Thus, prediction mechanisms accommodated with functional mixed effects models (FMEM) are useful. In this paper, we developed a classified functional mixed model prediction (CFMMP) method, which adapts classified mixed model prediction (CMMP) to the framework of FMEM. Performance of CFMMP against functional regression prediction based on simulation studies and the consistency property of CFMMP estimators are explored. Real-world applications of CFMMP are illustrated using real world examples including data from the hormone research menstrual cycles and the diffusion tensor imaging.

Cold-induced deposition of bivalent H3K4me3-H3K27me3 modification and nucleosome depletion in Arabidopsis.

Epigenetic regulation of gene expression plays a crucial role in plant development and environmental adaptation. The H3K4me3 and H3K27me3 have not only been discovered in the regulation of gene expression in multiple biological processes but also in responses to abiotic stresses in plants. However, evidence for the presence of both H3K4me3 and H3K27me3 on the same nucleosome is sporadic. Cold-induced deposition of bivalent H3K4me3-H3K27me3 modifications and nucleosome depletion over a considerable number of active genes is documented in potato tubers and provides clues on an additional role of the bivalent modifications. Limited by the available information of genes encoding PcG/TrxG proteins as well as their corresponding mutants in potatoes, the molecular mechanism underlying the cold-induced deposition of the bivalent mark remains elusive. In this study, we found a similar deposition of the bivalent H3K4me3-H3K27me3 mark over 2129 active genes in cold-treated Arabidopsis Col-0 seedlings. The expression levels of the bivalent mark-associated genes tend to be independent of bivalent modification levels. However, these genes were associated with greater chromatin accessibility, presumably to provide a distinct chromatin environment for gene expression. In mutants clf28 and lhp1, failure to deposit H3K27me3 in active genes upon cold treatment implies that the CLF is potentially involved in cold-induced deposition of H3K27me3, with assistance from LHP1. Failure to deposit H3K4me3 during cold treatment in atx1-2 suggests a regulatory role of ATX1 in the deposition of H3K4me3. In addition, we observed a cold-induced global reduction in nucleosome occupancy, which is potentially mediated by LHP1 in an H3K27me3-dependent manner.

Chromosome-scale genome assembly of the 'Munstead' cultivar of Lavandula angustifolia.

OBJECTIVES: Lavandula angustifolia (English lavender) is commercially important not only as an ornamental species but also as a major source of fragrances. To better understand the genomic basis of chemical diversity in lavender, we sequenced, assembled, and annotated the 'Munstead' cultivar of L. angustifolia. DATA DESCRIPTION: A total of 80 Gb of Oxford Nanopore Technologies reads was used to assemble the 'Munstead' genome using the Canu genome assembler software. Following multiple rounds of error correction and scaffolding using Hi-C data, the final chromosome-scale assembly represents 795,075,733 bp across 25 chromosomes with an N50 scaffold length of 31,371,815 bp. Benchmarking Universal Single Copy Orthologs analysis revealed 98.0% complete orthologs, indicative of a high-quality assembly representative of genic space. Annotation of protein-coding sequences revealed 58,702 high-confidence genes encoding 88,528 gene models. Access to the 'Munstead' genome will permit comparative analyses within and among lavender accessions and provides a pivotal species for comparative analyses within Lamiaceae.

Cold stress induces differential gene expression of retained homeologs in Camelina sativa cv Suneson.

Camelina sativa (L.) Crantz, a member of the Brassicaceae, has potential as a biofuel feedstock which is attributable to the production of fatty acids in its seeds, its fast growth cycle, and low input requirements. While a genome assembly is available for camelina, it was generated from short sequence reads and is thus highly fragmented in nature. Using long read sequences, we generated a chromosome-scale, highly contiguous genome assembly (644,491,969 bp) for the spring biotype cultivar 'Suneson' with an N50 contig length of 12,031,512 bp and a scaffold N50 length of 32,184,682 bp. Annotation of protein-coding genes revealed 91,877 genes that encode 133,355 gene models. We identified a total of 4,467 genes that were significantly up-regulated under cold stress which were enriched in gene ontology terms associated with "response to cold" and "response to abiotic stress". Coexpression analyses revealed multiple coexpression modules that were enriched in genes differentially expressed following cold stress that had putative functions involved in stress adaptation, specifically within the plastid. With access to a highly contiguous genome assembly, comparative analyses with Arabidopsis thaliana revealed 23,625 A. thaliana genes syntenic with 45,453 Suneson genes. Of these, 24,960 Suneson genes were syntenic to 8,320 A. thaliana genes reflecting a 3 camelina homeolog to 1 Arabidopsis gene relationship and retention of all three homeologs. Some of the retained triplicated homeologs showed conserved gene expression patterns under control and cold-stressed conditions whereas other triplicated homeologs displayed diverged expression patterns revealing sub- and neo-functionalization of the homeologs at the transcription level. Access to the chromosome-scale assembly of Suneson will enable both basic and applied research efforts in the improvement of camelina as a sustainable biofuel feedstock.

Dynamics of cis-regulatory sequences and transcriptional divergence of duplicated genes in soybean.

Transcriptional divergence of duplicated genes after whole genome duplication (WGD) has been described in many plant lineages and is often associated with subgenome dominance, a genome-wide mechanism. However, it is unknown what underlies the transcriptional divergence of duplicated genes in polyploid species that lack subgenome dominance. Soybean is a paleotetraploid with a WGD that occurred 5 to 13 Mya. Approximately 50% of the duplicated genes retained from this WGD exhibit transcriptional divergence. We developed accessible chromatin region (ACR) datasets from leaf, flower, and seed tissues using MNase-hypersensitivity sequencing. We validated enhancer function of several ACRs associated with known genes using CRISPR/Cas9-mediated genome editing. The ACR datasets were used to examine and correlate the transcriptional patterns of 17,111 pairs of duplicated genes in different tissues. We demonstrate that ACR dynamics are correlated with divergence of both expression level and tissue specificity of individual gene pairs. Gain or loss of flanking ACRs and mutation of cis-regulatory elements (CREs) within the ACRs can change the balance of the expression level and/or tissue specificity of the duplicated genes. Analysis of DNA sequences associated with ACRs revealed that the extensive sequence rearrangement after the WGD reshaped the CRE landscape, which appears to play a key role in the transcriptional divergence of duplicated genes in soybean. This may represent a general mechanism for transcriptional divergence of duplicated genes in polyploids that lack subgenome dominance.

Genome-wide chromatin accessibility landscape and dynamics of transcription factor networks during ovule and fiber development in cotton.

BACKGROUND: The development of cotton fiber is regulated by the orchestrated binding of regulatory proteins to cis-regulatory elements associated with developmental genes. The cis-trans regulatory dynamics occurred throughout the course of cotton fiber development are elusive. Here we generated genome-wide high-resolution DNase I hypersensitive sites (DHSs) maps to understand the regulatory mechanisms of cotton ovule and fiber development. RESULTS: We generated DNase I hypersensitive site (DHS) profiles from cotton ovules at 0 and 3 days post anthesis (DPA) and fibers at 8, 12, 15, and 18 DPA. We obtained a total of 1185 million reads and identified a total of 199,351 DHSs through ~ 30% unique mapping reads. It should be noted that more than half of DNase-seq reads mapped multiple genome locations and were not analyzed in order to achieve a high specificity of peak profile and to avoid bias from repetitive genomic regions. Distinct chromatin accessibilities were observed in the ovules (0 and 3 DPA) compared to the fiber elongation stages (8, 12, 15, and 18 DPA). Besides, the chromatin accessibility during ovules was particularly elevated in genomic regions enriched with transposable elements (TEs) and genes in TE-enriched regions were involved in ovule cell division. We analyzed cis-regulatory modules and revealed the influence of hormones on fiber development from the regulatory divergence of transcription factor (TF) motifs. Finally, we constructed a reliable regulatory network of TFs related to ovule and fiber development based on chromatin accessibility and gene co-expression network. From this network, we discovered a novel TF, WRKY46, which may shape fiber development by regulating the lignin content. CONCLUSIONS: Our results not only reveal the contribution of TEs in fiber development, but also predict and validate the TFs related to fiber development, which will benefit the research of cotton fiber molecular breeding.

Genome-Wide Identification of Open Chromatin in Plants Using MH-Seq.

Functional cis-regulatory elements (CREs) act as precise transcriptional switches for fine-tuning gene transcription. Identification of CREs is critical for understanding regulatory mechanisms of gene expression associated with various biological processes in eukaryotes. It is well known that CREs reside in open chromatin that exhibits hypersensitivity to enzyme cleavage and physical shearing. Currently, high-throughput methodologies, such as DNase-seq, ATAC-seq, and FAIRE-seq, have been widely applied in mapping open chromatin in various eukaryotic genomes. More recently, differential MNase (micrococcal nuclease) treatment has been successfully employed to map open chromatin in addition to profiling nucleosome landscape in both mammalian and plant species. We have developed a MNase hypersensitivity sequencing (MH-seq) technique in plants. The MH-seq procedure includes plant nuclei fixation and purification, differential treatments of purified nuclei with MNase, specific recovery of MNase-trimmed small DNA fragments within 20~100 bp in length, and MH-seq library construction followed by Illumina sequencing and data analysis. MH-seq has been successfully applied for global identification of open chromatin in both Arabidopsis thaliana and maize. It has been proven to be an attractive alternative for profiling open chromatin. Thus, MH-seq is expected to be valuable in probing chromatin accessibility on a genome-wide scale for other plants with sequenced genomes. Moreover, MHS data allow to implement footprinting assays to unveil binding sites of transcription factors.

Identification and functional validation of super-enhancers in Arabidopsis thaliana.

Super-enhancers (SEs) are exceptionally large enhancers and are recognized to play prominent roles in cell identity in mammalian species. We surveyed the genomic regions containing large clusters of accessible chromatin regions (ACRs) marked by deoxyribonuclease (DNase) I hypersensitivity in Arabidopsis thaliana. We identified a set of 749 putative SEs, which have a minimum length of 1.5 kilobases and represent the top 2.5% of the largest ACR clusters. We demonstrate that the genomic regions associating with these SEs were more sensitive to DNase I than other nonpromoter ACRs. The SEs were preferentially associated with topologically associating domains. Furthermore, the SEs and their predicted cognate genes were frequently associated with organ development and tissue identity in A. thaliana. Therefore, the A. thaliana SEs and their cognate genes mirror the functional characteristics of those reported in mammalian species. We developed CRISPR/Cas-mediated deletion lines of a 3,578-bp SE associated with the thalianol biosynthetic gene cluster (BGC). Small deletions (131-157 bp) within the SE resulted in distinct phenotypic changes and transcriptional repression of all five thalianol genes. In addition, T-DNA insertions in the SE region resulted in transcriptional alteration of all five thalianol genes. Thus, this SE appears to play a central role in coordinating the operon-like expression pattern of the thalianol BGC.

Chromosome-scale genome assemblies and annotations for Poales species Carex cristatella, Carex scoparia, Juncus effusus, and Juncus inflexus.

The majority of sequenced genomes in the monocots are from species belonging to Poaceae, which include many commercially important crops. Here, we expand the number of sequenced genomes from the monocots to include the genomes of 4 related cyperids: Carex cristatella and Carex scoparia from Cyperaceae and Juncus effusus and Juncus inflexus from Juncaceae. The high-quality, chromosome-scale genome sequences from these 4 cyperids were assembled by combining whole-genome shotgun sequencing of Nanopore long reads, Illumina short reads, and Hi-C sequencing data. Some members of the Cyperaceae and Juncaceae are known to possess holocentric chromosomes. We examined the repeat landscapes in our sequenced genomes to search for potential repeats associated with centromeres. Several large satellite repeat families, comprising 3.2-9.5% of our sequenced genomes, showed dispersed distribution of large satellite repeat clusters across all Carex chromosomes, with few instances of these repeats clustering in the same chromosomal regions. In contrast, most large Juncus satellite repeats were clustered in a single location on each chromosome, with sporadic instances of large satellite repeats throughout the Juncus genomes. Recognizable transposable elements account for about 20% of each of the 4 genome assemblies, with the Carex genomes containing more DNA transposons than retrotransposons while the converse is true for the Juncus genomes. These genome sequences and annotations will facilitate better comparative analysis within monocots.

Chromosome painting reveals inter-chromosomal rearrangements and evolution of subgenome D of wheat.

Aegilops species represent the most important gene pool for breeding bread wheat (Triticum aestivum). Thus, understanding the genome evolution, including chromosomal structural rearrangements and syntenic relationships among Aegilops species or between Aegilops and wheat, is important for both basic genome research and practical breeding applications. In the present study, we attempted to develop subgenome D-specific fluorescence in situ hybridization (FISH) probes by selecting D-specific oligonucleotides based on the reference genome of Chinese Spring. The oligo-based chromosome painting probes consisted of approximately 26 000 oligos per chromosome and their specificity was confirmed in both diploid and polyploid species containing the D subgenome. Two previously reported translocations involving two D chromosomes have been confirmed in wheat varieties and their derived lines. We demonstrate that the oligo painting probes can be used not only to identify the translocations involving D subgenome chromosomes, but also to determine the precise positions of chromosomal breakpoints. Chromosome painting of 56 accessions of Ae. tauschii from different origins led us to identify two novel translocations: a reciprocal 3D-7D translocation in two accessions and a complex 4D-5D-7D translocation in one accession. Painting probes were also used to analyze chromosomes from more diverse Aegilops species. These probes produced FISH signals in four different genomes. Chromosome rearrangements were identified in Aegilops umbellulata, Aegilops markgrafii, and Aegilops uniaristata, thus providing syntenic information that will be valuable for the application of these wild species in wheat breeding.

Inference about ratios of age-standardized rates with sampling errors in the population denominators for estimating both rates.

A rate ratio (RR) is an important metric for comparing cancer risks among different subpopulations. Inference for RR becomes complicated when populations used for calculating age-standardized cancer rates involve sampling errors, a situation that arises increasingly often when sample surveys must be used to obtain the population data. We compare a few strategies of estimating the standardized RR and propose bias-corrected ratio estimators as well as the corresponding variance estimators and confidence intervals that simultaneously consider the sampling error in estimating populations and the traditional Poisson error in the occurrence of cancer case or death. Performance of the proposed methods is evaluated empirically based on simulation studies. An application to immigration disparities in cancer mortality among Hispanic Americans is discussed. Our simulation studies show that a bias-corrected RR estimator performs the best in reducing the bias without increasing the coefficient of variation; the proposed variance estimators for the RR estimators and associated confidence intervals are fairly accurate. Finding of our application study are both interesting and consistent with the common sense as well as the results of our simulation studies.

Phased, chromosome-scale genome assemblies of tetraploid potato reveal a complex genome, transcriptome, and predicted proteome landscape underpinning genetic diversity.

Cultivated potato is a clonally propagated autotetraploid species with a highly heterogeneous genome. Phased assemblies of six cultivars including two chromosome-scale phased genome assemblies revealed extensive allelic diversity, including altered coding and transcript sequences, preferential allele expression, and structural variation that collectively result in a highly complex transcriptome and predicted proteome, which are distributed across the homologous chromosomes. Wild species contribute to the extensive allelic diversity in tetraploid cultivars, demonstrating ancestral introgressions predating modern breeding efforts. As a clonally propagated autotetraploid that undergoes limited meiosis, dysfunctional and deleterious alleles are not purged in tetraploid potato. Nearly a quarter of the loci bore mutations are predicted to have a high negative impact on protein function, complicating breeder's efforts to reduce genetic load. The StCDF1 locus controls maturity, and analysis of six tetraploid genomes revealed that 12 allelic variants of StCDF1 are correlated with maturity in a dosage-dependent manner. Knowledge of the complexity of the tetraploid potato genome with its rampant structural variation and embedded deleterious and dysfunctional alleles will be key not only to implementing precision breeding of tetraploid cultivars but also to the construction of homozygous, diploid potato germplasm containing favorable alleles to capitalize on heterosis in F1 hybrids.

Chromosome-specific painting unveils chromosomal fusions and distinct allopolyploid species in the Saccharum complex.

Karyotypes provide key cytogenetic information on the phylogenetic relationships and evolutionary origins in related eukaryotic species. Despite our knowledge of the chromosome numbers of sugarcane and its wild relatives, the chromosome composition and evolution among the species in the Saccharum complex have been elusive owing to the complex polyploidy and the large numbers of chromosomes of these species. Oligonucleotide-based chromosome painting has become a powerful tool of cytogenetic studies especially for plant species with large numbers of chromosomes. We developed oligo-based chromosome painting probes for all 10 chromosomes in Saccharum officinarum (2n = 8x = 80). The 10 painting probes generated robust fluorescence in situ hybridization signals in all plant species within the Saccharum complex, including species in the genera Saccharum, Miscanthus, Narenga and Erianthus. We conducted comparative chromosome analysis using the same set of probes among species from four different genera within the Saccharum complex. Excitingly, we discovered several novel cytotypes and chromosome rearrangements in these species. We discovered that fusion from two different chromosomes is a common type of chromosome rearrangement associated with the species in the Saccharum complex. Such fusion events changed the basic chromosome number and resulted in distinct allopolyploids in the Saccharum complex.

Variance Estimation and Confidence Intervals from Genome-wide Association Studies Through High-dimensional Misspecified Mixed Model Analysis.

We study variance estimation and associated confidence intervals for parameters characterizing genetic effects from genome-wide association studies (GWAS) in misspecified mixed model analysis. Previous studies have shown that, in spite of the model misspecification, certain quantities of genetic interests are consistently estimable, and consistent estimators of these quantities can be obtained using the restricted maximum likelihood (REML) method under a misspecified linear mixed model. However, the asymptotic variance of such a REML estimator is complicated and not ready to be implemented for practical use. In this paper, we develop practical and computationally convenient methods for estimating such asymptotic variances and constructing the associated confidence intervals. Performance of the proposed methods is evaluated empirically based on Monte-Carlo simulations and real-data application.