Genome-wide quantification of homeolog expression ratio revealed nonstochastic gene regulation in synthetic allopolyploid Arabidopsis

Genome duplication with hybridization, or allopolyploidization, occurs commonly in plants, and is considered to be a strong force for generating new species. However, genome-wide quantification of homeolog expression ratios was technically hindered because of the high homology between homeologous gene pairs. To quantify the homeolog expression ratio using RNA-seq obtained from polyploids, a new method named HomeoRoq was developed, in which the genomic origin of sequencing reads was estimated using mismatches between the read and each parental genome. To verify this method, we first assembled the two diploid parental genomes of Arabidopsis halleri subsp. gemmifera and Arabidopsis lyrata subsp. petraea (Arabidopsis petraea subsp. umbrosa), then generated a synthetic allotetraploid, mimicking the natural allopolyploid Arabidopsis kamchatica. The quantified ratios corresponded well to those obtained by Pyrosequencing. We found that the ratios of homeologs before and after cold stress treatment were highly correlated (r = 0.870). This highlights the presence of nonstochastic polyploid gene regulation despite previous research identifying stochastic variation in expression. Moreover, our new statistical test incorporating overdispersion identified 226 homeologs (1.11% of 20 369 expressed homeologs) with significant ratio changes, many of which were related to stress responses. HomeoRoq would contribute to the study of the genes responsible for polyploid-specific environmental responses.     

Genomic compatibility between two phyllotine rodent species evaluated through their hybrids

In order to investigate the genomic compatibility between allopatric rodent species, Phyllotis darwini and Phyllotis magister, we have studied several cytogenetic and reproductive features of their laboratory hybrids. Of thirty-one pairings between species, only five were successful, producing eleven newborns. Like parents, hybrids had 38 metacentric chromosomes, except for the subtelocentric Y chromosome inherited from P. magister. There was almost total C and G band correspondence between homeologous autosomes. However, parental sex chromosomes had different morphology, C and G bands. Ag-NOR bands appeared as small telomeric Ag+ regions, distributed in four chromosomal pairs of darwini, three of magister and four homeologous chromosomes of the hybrids. The three forms had similar indexes of NOR activity per cell, in spite of the variability in NOR expression which was always detected. Usually, only one member of parental homologous chromosomes showed AgNOR+; nevertheless, both homeologous chromosomes were active in many hybrid cells. The frequencies of cells that expressed their ribosomal genes in the two homologous or homeologous NOR chromosomes were similar in parental and hybrid cells. These results strongly suggest that ribosomal genes of both parental genomes would function codominantly in the hybrids. The gonad histological and morphometric analyses showed that hybrids conformed to Haldane's rule, since females were fertile and males were infertile. Our results indicate that P. darwini and P. magister genomes can function in relative harmony and compatibility when they are placed together in their laboratory generated hybrids, suggesting that these species have few genetic differences, probably because they have recently diverged.     

Comprehensive analyses of hox gene expression in Xenopus laevis embryos and adult tissues

From whole genome sequencing of an allotetraploid frog, Xenopus laevis, two homeologous sets (L and S) of four Hox clusters A through D (HoxA.L/S, HoxB.L/S, HoxC.L/S, and HoxD.L/S) and 13 paralogous groups (PGs) with 76 genes were identified, allowing us to carry out the first comprehensive analyses of hox gene expression in vertebrates. Expression of all hox genes during development and in adult tissues was analyzed by RNA‐sequencing. The expression levels of most hox genes were similar between homeologs, but in some pairs, large differences were observed and several of these were confirmed by RT‐PCR and whole mount in situ hybridization experiments. These results indicate that subfunctionalization of hox genes may have occurred since allotetraploidization. Furthermore, comprehensive analysis of hox gene expression during early development did not agree with the hypothesis of temporal collinearity especially in genes belonging to PG2 to PG10 .     

The alpha-tubulin gene family in wheat (Triticum aestivum L.) and differential gene expression during cold acclimation.

The alpha-tubulins and beta-tubulins are the major constituents of microtubules, which have been recognized as important structural elements in cell growth and morphogenesis, and, recently, for their role in regulation and signal transduction. We have identified 15 full-length cDNAs for the members of the alpha-tubulin gene family in hexaploid bread wheat (Triticum aestivum L.). The genes were clustered into 5 homeologous groups of 3 genes. Representatives of the 5 homeologous groups were mapped to different chromosome arms, and the genome of origin was determined for each gene. Changes in mRNA levels were observed for the paralogous members of the gene family during cold acclimation. Three members of the family had initial decreases in mRNA levels in response to cold treatment, which were followed by increases, each with a different pattern of reinduction. One gene-family member showed increased mRNA for up to 14 d during cold acclimation and had decreased levels after 36 d of cold treatment; a fifth paralogous member of the gene family had slowly declining mRNA levels up to 36 d. Subtle differences in the level of gene expression among homeologs and large differences among paralogs were detected by comparing the relative abundance of wheat alpha-tubulin expressed sequence tags (ESTs) in public databases.     

Recombinant genome of cereals: The pattern of formation and the role in evolution of polyploid species

The pivotal-differential model of evolution of polyploid species of cereals has been experimentally reproduced, and the pattern of the formation of a recombinant genome has been analyzed. It has been found that mutual substitution of chromosomes of the original genomes is subjected to selection pressure and, hence, is nonrandom. The selection occurs at the level of homeologs, whose selective advantages are determined by interactions between the genotype and the environment. If a homeolog has distinct selective advantages, the chromosomal composition of the corresponding homeologous group is completed rapidly, which leads to the formation of intergenomic recombination at the level of whole chromosomes. If homeologs have the same competitiveness, the composition of the group is stabilized more slowly. Domination of the genetic systems of the basic genome ensures a high rate of pairing of homeologous chromosomes of the recombinant genome during meiosis, which leads to recombinations at the level of chromosomal segments. It has been demonstrated that different combinations of chromosomes from original genomes are selected at different conditions of plant growth.     

Linkage arrangement of Na,K-ATPase genes in the tetraploid-derived genome of the rainbow trout (Oncorhynchus mykiss)

As part of our efforts to characterize Na,K-ATPase isoforms in salmonid fish, we investigated the linkage arrangement of genes coding for the alpha and beta-subunits of the enzyme complex in the tetraploid-derived genome of the rainbow trout (Oncorhynchus mykiss). Genetic markers were developed from four of five previously characterized alpha-subunit isoforms (alpha1b, alpha1c, alpha2 and alpha3) and four expressed sequence tags derived from yet undescribed beta-subunit isoforms (beta1a, beta1b, beta3a and beta3b). Sex-specific linkage analysis of polymorphic loci in a reference meiotic panel revealed that Na,K-ATPase genes are generally dispersed throughout the rainbow trout genome. A notable exception was the colocalization of two alpha-subunit genes and one beta-subunit gene on linkage group RT-12, which may thus share a conserved orthologous segment with linkage group 1 in zebrafish (Danio rerio). Consistent with previously reported homeologous relationships among the chromosomes of the rainbow trout, primers designed from the alpha3-isoform detected a pair of duplicated genes on linkage groups RT-27 and RT-31. Similarly, the evolutionary conservation of homeologous regions on linkage groups RT-12 and RT-16 was further supported by the map localization of gene duplicates for the beta1b isoform. The detection of homeologs within each gene family also raises the possibility that novel isoforms may be discovered as functional duplicates.     

Identification and Molecular Characterization of Homeologous Δ9-Stearoyl Acyl Carrier Protein Desaturase 3 Genes from the Allotetraploid Peanut (Arachis hypogaea)

The stearoyl–acyl carrier protein (ACP) desaturase (SAD) is a nuclear-encoded, plastid-localized soluble desaturase that catalyzes the conversion of stearoyl-ACP to oleoyl-ACP and plays a key role in the determination of the properties of the majority of cellular glycerolipids. Sad genes from a variety of plant species have been cloned and characterized. However, in peanut (Arachis hypogaea), an important edible and oilseed crop, these genes have not yet been characterized. By searching peanut expressed sequence tag (EST) and parallel sequencing (454) libraries, we have identified three members of the ahSad gene family. Among them, only one gene, ahSad3, was exclusively expressed during seed development and in a manner fully corresponding to oil accumulation. Both ahSad3 homeologous genes (ahSad3A and ahSad3B) were recovered from the allotetraploid peanut, and their mRNA expression levels were characterized. The open reading frames for ahSad3A and ahSad3B are 98% identical and consist of 1,158 bp, encoding a 386-full-amino-acid protein, with one intron in the coding sequence. Comparisons of the sequences of these two homeologous genes revealed seven single-nucleotide polymorphisms and one triplet insertion in the coding region. Southern blot analysis indicated that there are only two copies of the ahSad3 gene in the peanut genome. Homeolog-specific gene expression analysis showed that both ahSad3 homeologs are expressed in developing seeds, but gene expression is significantly biased toward the B genome. Our results point to ahSad3 as a possible target gene for manipulation of fatty acid saturation in A. hypogaea.     

Recombinant genome of cereals: the pattern of formation and the role in evolution of polyploid species

The pivotal-differential model of evolution of polyploid species of cereals has been experimentally reproduced, and the pattern of the formation of a recombinant genome has been analyzed. It has been found that mutual substitution of chromosomes of the original genomes is subjected to selection pressure and, hence, is nonrandom. The selection occurs at the level of homeologs, whose selective advantages are determined by interactions between the genotype and the environment. If a homeolog has distinct selective advantages, the chromosomal composition of the corresponding homeologous group is completed rapidly, which leads to the formation of intergenomic recombination at the level of whole chromosomes. If homeologs have the same competitiveness, the composition of the group is stabilized more slowly. Domination of the genetic systems of the basic genome ensures a high rate of pairing of homeologous chromosomes of the recombinant genome during meiosis, which leads to recombinations at the level of chromosomal segments. It has been demonstrated that different combinations of chromosomes from original genomes are selected at different conditions of plant growth.     

Impairing both HMA4 homeologs is required for cadmium reduction in tobacco

In tobacco, the heavy metal P1B‐ATPases HMA4.1 and HMA4.2 function in root‐to‐shoot zinc and cadmium transport. We present greenhouse and field data that dissect the possibilities to impact the two homeologous genes in order to define the best strategy for leaf cadmium reduction. In a first step, both genes were silenced using an RNAi approach leading to >90% reduction of leaf cadmium content. To modulate HMA4 function more precisely, mutant HMA4.1 and HMA4.2 alleles of a Targeting Induced Local Lesions IN Genomes (TILLING) population were combined. As observed with RNAi plants, knockout of both homeologs decreased cadmium root‐to‐shoot transfer by >90%. Analysis of plants with segregating null and wild‐type alleles of both homeologs showed that one functional HMA4 allele is sufficient to maintain wild‐type cadmium levels. Plant development was affected in HMA4 RNAi and double knockout plants that included retarded growth, necrotic lesions, altered leaf morphology and increased water content. The combination of complete functional loss (nonsense mutation) in one homeologous HMA4 gene and the functional reduction in the other HMA4 gene (missense mutation) is proposed as strategy to limit cadmium leaf accumulation without developmental effects.