Genomic, regulatory and epigenetic mechanisms underlying duplicated gene evolution in the natural allotetraploid Oryza minuta

RNAcontacts: A Pipeline for Predicting Contacts from RNA Proximity Ligation Assays

High-throughput RNA proximity ligation assays are molecular methods that are used to simultaneously analyze the spatial proximity of many RNAs in living cells. Their principle is based on cross-linking, fragmentation, and subsequent religation of RNAs, followed by high-throughput sequencing. The generated fragments have two different types of splits, one resulting from pre-mRNA splicing and the other formed by the ligation of spatially close RNA strands. Here, we present RNAcontacts, a universal pipeline for detecting RNA-RNA contacts in high-throughput RNA proximity ligation assays. RNAcontacts circumvents the inherent problem of mapping sequences with two distinct types of splits using a two-pass alignment, in which splice junctions are inferred from a control RNA-seq experiment on the first pass and then provided to the aligner as bona fide introns on the second pass. Compared to previously developed methods, our approach allows for a more sensitive detection of RNA contacts and has a higher specificity with respect to splice junctions that are present in the biological sample. RNAcontacts automatically extracts contacts, clusters their ligation points, computes the read support, and generates tracks for visualizing through the UCSC Genome Browser. The pipeline is implemented in Snakemake, a reproducible and scalable workflow management system for rapid and uniform processing of multiple datasets. RNAcontacts is a generic pipeline for the detection of RNA contacts that can be used with any proximity ligation method as long as one of the interacting partners is RNA. RNAcontacts is available via the GitHub repository https://github.com/smargasyuk/ RNAcontacts/.

Introgression of Maize Diversity for Drought Tolerance: Subtropical Maize Landraces as Source of New Positive Variants

Current climate change models predict an increased frequency and intensity of drought for much of the developing world within the next 30 years. These events will negatively affect maize yields, potentially leading to economic and social instability in many smallholder farming communities. Knowledge about the genetic resources available for traits related to drought tolerance has great importance in developing breeding program strategies. The aim of this research was to study a maize landrace introgression panel to identify chromosomal regions associated with a drought tolerance index. For that, we performed Genome-Wide Association Study (GWAS) on 1326 landrace progenies developed by the CIMMYT Genetic Resources Program, originating from 20 landraces populations collected in arid regions. Phenotypic data were obtained from early testcross trials conducted in three sites and two contrasting irrigation environments, full irrigation (well-watered) and reduced irrigation (drought). The populations were genotyped using the DArTSeq^® platform, and a final set of 5,695 SNPs markers was used. The genotypic values were estimated using spatial adjustment in a two-stage analysis. First, we performed the individual analysis for each site/irrigation treatment combination. The best linear unbiased estimates (BLUEs) were used to calculate the Harmonic Mean of Relative Performance (HMRP) as a drought tolerance index for each testcross. The second stage was a joint analysis, which was performed using the HMRP to obtain the best linear unbiased predictions (BLUPs) of the index for each genotype. Then, GWAS was performed to determine the marker-index associations and the marker-Grain Yield (GY) associations for the two irrigation treatments. We detected two significant markers associated with the drought-tolerance index, four associated with GY in drought condition, and other four associated with GY in irrigated conditions each. Although each of these markers explained less than 0.1% of the phenotypic variation for the index and GY, we found two genes likely related to the plant response to drought stress. For these markers, alleles from landraces provide a slightly higher yield under drought conditions. Our results indicate that the positive diversity delivered by landraces are still present on the backcrosses and this is a potential breeding strategy for improving maize for drought tolerance and for trait introgression bringing new superior allelic diversity from landraces to breeding populations.

Genome evolution in Oryza allopolyploids of various ages: Insights into the process of diploidization

The prevalence and recurrence of whole-genome duplication in plants and its major role in evolution have been well recognized. Despite great efforts, many aspects of genome evolution, particularly the temporal progression of genomic responses to allopolyploidy and the underlying mechanisms, remain poorly understood. The rice genus Oryza consists of both recently formed and older allopolyploid species, representing an attractive system for studying the genome evolution after allopolyploidy. In this study, through screening BAC libraries and sequencing and annotating the targeted BAC clones, we generated orthologous genomic sequences surrounding the DEP1 locus, a major grain yield QTL in cultivated rice, from four Oryza polyploids of various ages and their likely diploid genome donors or close relatives. Based on sequenced DEP1 region and published data from three other genomic regions, we investigated the temporal evolutionary dynamics of four polyploid genomes at both genetic and expression levels. In the recently formed BBCC polyploid, Oryza minuta, genome dominance was not observed and its short-term responses to allopolyploidy are mainly manifested as a high proportion of homoeologous gene pairs showing unequal expression. This could partly be explained by parental legacy, rewiring of divergent regulatory networks and epigenetic modulation. Moreover, we detected an ongoing diploidization process in this genus, and suggest that the expression divergence driven by changes of selective constraint probably plays a big role in the long-term diploidization. These findings add novel insights into our understanding of genome evolution after allopolyploidy, and could facilitate crop improvements through hybridization and polyploidization.

In situ dissecting the evolution of gene duplication with different histone modification patterns based on high-throughput data analysis in Arabidopsis thaliana

Background

Genetic regulation is known to contribute to the divergent expression of duplicate genes; however, little is known about how epigenetic modifications regulate the expression of duplicate genes in plants.

Methods

The histone modification (HM) profile patterns of different modes of gene duplication, including the whole genome duplication, proximal duplication, tandem duplication and transposed duplication were characterized based on ChIP-chip or ChIP-seq datasets. In this study, 10 distinct HM marks including H2Bub, H3K4me1, H3K4me2, H3K4me3, H3K9ac, H3K9me2, H3K27me1, H3K27me3, H3K36me3 and H3K14ac were analyzed. Moreover, the features of gene duplication with different HM patterns were characterized based on 88 RNA-seq datasets of Arabidopsis thaliana.

Results

This study showed that duplicate genes in Arabidopsis have a more similar HM pattern than single-copy genes in both their promoters and protein-coding regions. The evolution of HM marks is found to be coupled with coding sequence divergence and expression divergence after gene duplication. We found that functionally selective constraints may impose on epigenetic evolution after gene duplication. Furthermore, duplicate genes with distinct functions have more divergence in histone modification compared with the ones with the same function, while higher expression divergence is found with mutations of chromatin modifiers. This study shows the role of epigenetic marks in regulating gene expression and functional divergence after gene duplication in plants based on sequencing data.

Genome-wide characterization of GRAS family genes in Medicago truncatula reveals their evolutionary dynamics and functional diversification

The GRAS gene family is a large plant-specific family of transcription factors that are involved in diverse processes during plant development. Medicago truncatula is an ideal model plant for genetic research in legumes, and specifically for studying nodulation, which is crucial for nitrogen fixation. In this study, 59 MtGRAS genes were identified and classified into eight distinct subgroups based on phylogenetic relationships. Motifs located in the C-termini were conserved across the subgroups, while motifs in the N-termini were subfamily specific. Gene duplication was the main evolutionary force for MtGRAS expansion, especially proliferation of the LISCL subgroup. Seventeen duplicated genes showed strong effects of purifying selection and diverse expression patterns, highlighting their functional importance and diversification after duplication. Thirty MtGRAS genes, including NSP1 and NSP2, were preferentially expressed in nodules, indicating possible roles in the process of nodulation. A transcriptome study, combined with gene expression analysis under different stress conditions, suggested potential functions of MtGRAS genes in various biological pathways and stress responses. Taken together, these comprehensive analyses provide basic information for understanding the potential functions of GRAS genes, and will facilitate further discovery of MtGRAS gene functions.

Impact of transposable elements on polyploid plant genomes

BACKGROUND

The growing wealth of knowledge on whole-plant genome sequences is highlighting the key role of transposable elements (TEs) in plant evolution, as a driver of drastic changes in genome size and as a source of an important number of new coding and regulatory sequences. Together with polyploidization events, TEs should thus be considered the major players in evolution of plants.

SCOPE

This review outlines the major mechanisms by which TEs impact plant genome evolution and how polyploidy events can affect these impacts, and vice versa. These include direct effects on genes, by providing them with new coding or regulatory sequences, an effect on the epigenetic status of the chromatin close to genes, and more subtle effects by imposing diverse evolutionary constraints to different chromosomal regions. These effects are particularly relevant after polyploidization events. Polyploidization often induces bursts of transposition probably due to a relaxation in their epigenetic control, and, in the short term, this can increase the rate of gene mutations and changes in gene regulation due to the insertion of TEs next to or into genes. Over longer times, TE bursts may induce global changes in genome structure due to inter-element recombination including losses of large genome regions and chromosomal rearrangements that reduce the genome size and the chromosome number as part of a process called diploidization.

CONCLUSIONS

TEs play an essential role in genome and gene evolution, in particular after polyploidization events. Polyploidization can induce TE activity that may explain part of the new phenotypes observed. TEs may also play a role in the diploidization that follows polyploidization events. However, the extent to which TEs contribute to diploidization and fractionation bias remains unclear. Investigating the multiple factors controlling TE dynamics and the nature of ancient and recent polyploid genomes may shed light on these processes.

Identification of Major Quantitative Trait Loci for Seed Oil Content in Soybeans by Combining Linkage and Genome-Wide Association Mapping

Soybean oil is the most widely produced vegetable oil in the world and its content in soybean seed is an important quality trait in breeding programs. More than 100 quantitative trait loci (QTLs) for soybean oil content have been identified. However, most of them are genotype specific and/or environment sensitive. Here, we used both a linkage and association mapping methodology to dissect the genetic basis of seed oil content of Chinese soybean cultivars in various environments in the Jiang-Huai River Valley. One recombinant inbred line (RIL) population (NJMN-RIL), with 104 lines developed from a cross between M8108 and NN1138-2, was planted in five environments to investigate phenotypic data, and a new genetic map with 2,062 specific-locus amplified fragment markers was constructed to map oil content QTLs. A derived F2 population between MN-5 (a line of NJMN-RIL) and NN1138-2 was also developed to confirm one major QTL. A soybean breeding germplasm population (279 lines) was established to perform a genome-wide association study (GWAS) using 59,845 high-quality single nucleotide polymorphism markers. In the NJMN-RIL population, 8 QTLs were found that explained a range of phenotypic variance from 6.3 to 26.3% in certain planting environments. Among them, qOil-5-1, qOil-10-1, and qOil-14-1 were detected in different environments, and qOil-5-1 was further confirmed using the secondary F2 population. Three loci located on chromosomes 5 and 20 were detected in a 2-year long GWAS, and one locus that overlapped with qOil-5-1 was found repeatedly and treated as the same locus. qOil-5-1 was further localized to a linkage disequilibrium block region of approximately 440 kb. These results will not only increase our understanding of the genetic control of seed oil content in soybean, but will also be helpful in marker-assisted selection for breeding high seed oil content soybean and gene cloning to elucidate the mechanisms of seed oil content.

Genomic Identification and Comparative Expansion Analysis of the Non-Specific Lipid Transfer Protein Gene Family in Gossypium

Plant non-specific lipid transfer proteins (nsLTPs) are involved in many biological processes. In this study, 51, 47 and 91 nsLTPs were identified in Gossypium arboreum, G. raimondii and their descendant allotetraploid G. hirsutum, respectively. All the nsLTPs were phylogenetically divided into 8 distinct subfamilies. Besides, the recent duplication, which is considered cotton-specific whole genome duplication, may have led to nsLTP expansion in Gossypium. Both tandem and segmental duplication contributed to nsLTP expansion in G. arboreum and G. hirsutum, while tandem duplication was the dominant pattern in G. raimondii. Additionally, the interspecific orthologous gene pairs in Gossypium were identified. Some GaLTPs and GrLTPs lost their orthologs in the At and Dt subgenomes, respectively, of G. hirsutum. The distribution of these GrLTPs and GaLTPs within each subfamily was complementary, suggesting that the loss and retention of nsLTPs in G. hirsutum might not be random. Moreover, the nsLTPs in the At and Dt subgenomes might have evolved symmetrically. Furthermore, both intraspecific and interspecific orthologous genes showed considerable expression variation, suggesting that their functions were strongly differentiated. Our results lay an important foundation for expansion and evolutionary analysis of the nsLTP family in Gossypium, and advance nsLTP studies in other plants, especially polyploid plants.

Multiple origins of BBCC allopolyploid species in the rice genus (Oryza)

In the rice genus (Oryza), about one half of the species are allopolyploids. These species are not only important resources for rice breeding but also provide a unique opportunity for studying evolution of polyploid species. In the present study, we sequenced four biparentally inherited nuclear loci and three maternally inherited chloroplast fragments from all diploid and tetraploid species with the B- and C-genome types in this genus. We detected at least three independent origins of three BC-genome tetraploid species. Specifically, the diploid O. punctata (B-genome) and O. officinalis (C-genome) were the parental progenitors of O. minuta and O. malampuzhaensis with O. punctata being the maternal donors, whereas the diploid O. punctata and O. eichingeri (C-genome) were the progenitors of tetraploid O. punctata with O. punctata being the paternal donor. Our relaxed clock analyses suggest that all the BBCC species originated within the last one million years, which is coincident with the severe climate oscillations occurred during the last ice age, implying the potential impact of climate change on their formations and dispersals. In addition, our results support previous taxonomic arguments that the tetraploid O. punctata might be better treated as a separate species (O. schweinfurthiana).

Evolution of an epigenetic gene ensemble within the genus Anopheles

Epigenetic control of gene expression has important implications for the regulation of developmental processes, for mediating homeostasis and responses to the external environment, and for transgenerational inheritance of gene expression patterns. Genes that mediate epigenetic control have been well-characterized in Drosophila melanogaster, and we have identified and analyzed an orthologous gene ensemble in Anopheles gambiae that comprises 169 orthologs related to a 215-member epigenetic gene ensemble in D. melanogaster. We find that this ensemble is highly conserved among anopheline mosquitoes, as we identify only seven gene family expansion/contraction events within the ensemble among 12 mosquito species we have studied within the genus Anopheles. Comparative analyses of the epigenetic gene expression across the genera Drosophila and Anopheles reveal distinct tissue-associated expression patterns in the two genera, but similar temporal expression patterns. The A. gambiae complex and D. melanogaster subgroup epigenetic gene ensembles exhibit similar evolutionary rates, as assessed by their respective dN/dS values. These differences in tissue-associated expression patterns, in contrast to similarities in evolutionary rates and temporal expression patterns, may imply that some members of the epigenetic gene ensemble have been redeployed within one or both genera, in comparison to the most recent common ancestor of these two clades. Members of this epigenetic gene ensemble may constitute another set of potential targets for vector control and enable further reductions in the burden of human malaria, by analogy to recent success in development of small molecule antagonists for mammalian epigenetic machinery.