Genome-wide analysis of genomic imprinting in the endosperm and allelic variation in flax

Genomic imprinting is an epigenetic phenomenon that causes biased expression of maternally and paternally inherited alleles. In flowering plants, genomic imprinting predominantly occurs in the triploid endosperm and plays a vital role in seed development. In this study, we identified 248 candidate imprinted genes including 114 maternally expressed imprinted genes (MEGs) and 134 paternally expressed imprinted genes (PEGs) in flax (Linum usitatissimum L.) endosperm using deep RNA sequencing. These imprinted genes were neither clustered in specific chromosomal regions nor well conserved among flax and other plant species. MEGs tended to be expressed specifically in the endosperm, whereas the expression of PEGs was not tissue-specific. Imprinted single nucleotide polymorphisms differentiated 200 flax cultivars into the oil flax, oil-fiber dual purpose flax and fiber flax subgroups, suggesting that genomic imprinting contributed to intraspecific variation in flax. The nucleotide diversity of imprinted genes in the oil flax subgroup was significantly higher than that in the fiber flax subgroup, indicating that some imprinted genes underwent positive selection during flax domestication from oil flax to fiber flax. Moreover, imprinted genes that underwent positive selection were related to flax functions. Thirteen imprinted genes related to flax seed size and weight were identified using a candidate gene-based association study. Therefore, our study provides information for further exploration of the function and genomic variation of imprinted genes in the flax population.     

Identification of a Pair of Phospholipid:Diacylglycerol Acyltransferases from Developing Flax (Linum usitatissimum L.) Seed Catalyzing the Selective Production of Trilinolenin

The oil from flax (Linum usitatissimum L.) has high amounts of α-linolenic acid (ALA; 18:3^cisΔ9,12,15) and is one of the richest sources of omega-3 polyunsaturated fatty acids (ω-3-PUFAs). To produce ∼57% ALA in triacylglycerol (TAG), it is likely that flax contains enzymes that can efficiently transfer ALA to TAG. To test this hypothesis, we conducted a systematic characterization of TAG-synthesizing enzymes from flax. We identified several genes encoding acyl-CoA:diacylglycerol acyltransferases (DGATs) and phospholipid:diacylglycerol acyltransferases (PDATs) from the flax genome database. Due to recent genome duplication, duplicated gene pairs have been identified for all genes except DGAT2-2. Analysis of gene expression indicated that two DGAT1, two DGAT2, and four PDAT genes were preferentially expressed in flax embryos. Yeast functional analysis showed that DGAT1, DGAT2, and two PDAT enzymes restored TAG synthesis when produced recombinantly in yeast H1246 strain. The activity of particular PDAT enzymes (LuPDAT1 and LuPDAT2) was stimulated by the presence of ALA. Further seed-specific expression of flax genes in Arabidopsis thaliana indicated that DGAT1, PDAT1, and PDAT2 had significant effects on seed oil phenotype. Overall, this study indicated the existence of unique PDAT enzymes from flax that are able to preferentially catalyze the synthesis of TAG containing ALA acyl moieties. The identified LuPDATs may have practical applications for increasing the accumulation of ALA and other polyunsaturated fatty acids in oilseeds for food and industrial applications.     

Developing Fiber Specific Promoter-Reporter Transgenic Lines to Study the Effect of Abiotic Stresses on Fiber Development in Cotton

Cotton is one of the most important cash crops in US agricultural industry. Environmental stresses, such as drought, high temperature and combination of both, not only reduce the overall growth of cotton plants, but also greatly decrease cotton lint yield and fiber quality. The impact of environmental stresses on fiber development is poorly understood due to technical difficulties associated with the study of developing fiber tissues and lack of genetic materials to study fiber development. To address this important question and provide the need for scientific community, we have generated transgenic cotton lines harboring cotton fiber specific promoter (CFSP)-reporter constructs from six cotton fiber specific genes (Expansin, E6, Rac13, CelA1, LTP, and Fb late), representing genes that are expressed at different stages of fiber development. Individual CFSP::GUS or CFSP::GFP construct was introduced into Coker 312 via Agrobacterium mediated transformation. Transgenic cotton lines were evaluated phenotypically and screened for the presence of selectable marker, reporter gene expression, and insertion numbers. Quantitative analysis showed that the patterns of GUS reporter gene activity during fiber development in transgenic cotton lines were similar to those of the native genes. Greenhouse drought and heat stress study showed a correlation between the decrease in promoter activities and decrease in fiber length, increase in micronaire and changes in other fiber quality traits in transgenic lines grown under stressed condition. These newly developed materials provide new molecular tools for studying the effects of abiotic stresses on fiber development and may be used in study of cotton fiber development genes and eventually in the genetic manipulation of fiber quality.     

Comprehensive Analysis of the COBRA-Like (COBL) Gene Family in Gossypium Identifies Two COBLs Potentially Associated with Fiber Quality

COBRA-Like (COBL) genes, which encode a plant-specific glycosylphosphatidylinositol (GPI) anchored protein, have been proven to be key regulators in the orientation of cell expansion and cellulose crystallinity status. Genome-wide analysis has been performed in A. thaliana, O. sativa, Z. mays and S. lycopersicum, but little in Gossypium. Here we identified 19, 18 and 33 candidate COBL genes from three sequenced cotton species, diploid cotton G. raimondii, G. arboreum and tetraploid cotton G. hirsutum acc. TM-1, respectively. These COBL members were anchored onto 10 chromosomes in G. raimondii and could be divided into two subgroups. Expression patterns of COBL genes showed highly developmental and spatial regulation in G. hirsutum acc. TM-1. Of them, GhCOBL9 and GhCOBL13 were preferentially expressed at the secondary cell wall stage of fiber development and had significantly co-upregulated expression with cellulose synthase genes GhCESA4, GhCESA7 and GhCESA8. Besides, GhCOBL9 D_t and GhCOBL13 D_t were co-localized with previously reported cotton fiber quality quantitative trait loci (QTLs) and the favorable allele types of GhCOBL9 D_t had significantly positive correlations with fiber quality traits, indicating that these two genes might play an important role in fiber development.     

Genome-wide identification,characterization and expression analysis of the amino acid permease gene family in tea plants (Camellia sinensis)

Amino acid permeases (AAPs) are involved in transporting a broad spectrum of amino acids and regulating physiological processes in plants. In this study, 19 AAP genes were identified from the tea plants genome database and named CsAAP1–19. Based on phylogenetic analysis, the CsAAP genes were classified into three groups, having significantly different structures and conserved motifs. In addition, an expression analysis revealed that most of CsAAP genes were specifically expressed in different tissues, especially CsAAP19 was expressed only in root. These genes also were significantly expressed in the Baiye 1 and Huangjinya cultivars. Nitrogen treatments indicated that the CsAAPs were obviously expressed in root. CsAAP2, −6, −12, −13 and − 16 were significantly expressed at 6 d after the glutamate treatment, while the expression trend at 24 h after contained the ammonium. These results improve our understanding of the CsAAP genes and their functions in nitrogen utilization in tea plants.     

Isolation and characterization of genes encoding lipid transfer proteins in Linum usitatissimum

Very little is known about lipid transfer proteins from flax (Linum usitatissimum L.). In the present work, three genes encoding a lipid transfer protein (LTP) were isolated from flax, two of which encoded Type-1 and one Type-2 LTPs with molecular masses of about 9 and 7 kDa, respectively. The analysis of deduced amino acid sequence reveals that only Type 2 of the L. usitatissimum leaf specific LTP (LuLTP_Ls) had an N terminal signal peptide consisting of 23 amino acids. The phylogenetic analyses of LuLTP_Ls suggest their closest relatedness with respective proteins from Dimocarpus longan and Vitis vinifera. The gene expression analysis shows that LTP Type 1 genes, which include LuLTP_Ls1 and LuLTP_Ls3, were progressively expressed during leaf development, whereas LuLTP_Ls4 (Type 2) was expressed only at initial and terminal senescence stages of cotyledons. The results suggest that both types of LuLTP_Ls were differentially yet significantly expressed in cotyledons implicating their function in transport and scavenging lipidic skeletons for the benefit of other developing parts of the plant.     

Genetic variation of six desaturase genes in flax and their impact on fatty acid composition

Flax (Linum usitatissimum L.) is one of the richest plant sources of omega-3 fatty acids praised for their health benefits. In this study, the extent of the genetic variability of genes encoding stearoyl-ACP desaturase (SAD), and fatty acid desaturase 2 (FAD2) and 3 (FAD3) was determined by sequencing the six paralogous genes from 120 flax accessions representing a broad range of germplasm including some EMS mutant lines. A total of 6 alleles for sad1 and sad2, 21 for fad2a, 5 for fad2b, 15 for fad3a and 18 for fad3b were identified. Deduced amino acid sequences of the alleles predicted 4, 2, 3, 4, 6 and 7 isoforms, respectively. Allele frequencies varied greatly across genes. Fad3a, with 110 SNPs and 19 indels, and fad3b, with 50 SNPs and 5 indels, showed the highest levels of genetic variations. While most of the SNPs and all the indels were silent mutations, both genes carried nonsense SNP mutations resulting in premature stop codons, a feature not observed in sad and fad2 genes. Some alleles and isoforms discovered in induced mutant lines were absent in the natural germplasm. Correlation of these genotypic data with fatty acid composition data of 120 flax accessions phenotyped in six field experiments revealed statistically significant effects of some of the SAD and FAD isoforms on fatty acid composition, oil content and iodine value. The novel allelic variants and isoforms identified for the six desaturases will be a resource for the development of oilseed flax with unique and useful fatty acid profiles.     

Development of specific and degenerate primers for CesA genes encoding cellulose synthase in flax (Linum usitatissimum L.)

Representatives of the CesA multigene family that control the synthesis of the catalytic subunits of the cellulose synthase complex were described for a number of higher plants. It has been established that the HVR2 region of these genes is class-specific and determines the involvement of the gene product in the synthesis of either the primary or secondary cell wall. The purpose of the current research was to develop degenerate and specific primers for parts of the CesA genes to allow the construction of molecular markers for the class-specific HVR2 region. Two pairs of specific primers for the CesA-1 and CesA-6 genes as well as a pair of degenerate primers for the HVR2 region of all flax CesA genes were developed, based on analysis of the CesA ESTs as well as the full-length cDNA sequences of the CesA genes in Arabidopsis, poplar, maize, and cotton that are available in the GenBank. Fragments of the expected size were amplified using flax cDNA as a template (201 bp for CesA-1, 300 bp for CesA-6, and 600 bp for HVR2). The markers developed in this research can be used for CesA gene cloning and sequencing, analysis of gene copy numbers as well as characterization of tissue- and development specific gene expression.     

Imprinting analysis of the mouse chromosome 7C region in DNMT1-null embryos

The mouse chromosome 7C, orthologous to the human 15q11–q13 has an imprinted domain, where most of the genes are expressed only from the paternal allele. The imprinted domain contains paternally expressed genes, Snurf/Snrpn, Ndn, Magel2, Mkrn3, and Frat3, C/D-box small nucleolar RNAs (snoRNAs), and the maternally expressed gene, Ube3a. Imprinted expression in this large (approximately 3–4 Mb) domain is coordinated by a bipartite cis-acting imprinting center (IC), located upstream of the Snurf/Snrpn gene. The molecular mechanism how IC regulates gene expression of the whole domain remains partially understood. Here we analyzed the relationship between imprinted gene expression and DNA methylation in the mouse chromosome 7C using DNA methyltransferase 1 (DNMT1)-null mutant embryos carrying Dnmt1^ps alleles, which show global loss of DNA methylation and embryonic lethality. In the DNMT1-null embryos at embryonic day 9.5, the paternally expressed genes were biallelically expressed. Bisulfite DNA methylation analysis revealed loss of methylation on the maternal allele in the promoter regions of the genes. These results demonstrate that DNMT1 is necessary for monoallelic expression of the imprinted genes in the chromosome 7C domain, suggesting that DNA methylation in the secondary differentially methylated regions (DMRs), which are acquired during development serves primarily to control the imprinted expression from the maternal allele in the mouse chromosome 7C.