Sex-specific polymorphism of MET1 and ARR17 genes in Populus × sibirica

Sex-determining region complements traditionally used in phylogenetic studies nuclear and chloroplast sequences in investigation of Aigeiros Duby and Tacamahaca Spach poplars (genus Populus L., Salicaceae)

Members of the genus Populus L. play an important role in the formation of forests in the northern hemisphere and are used in urban landscaping and timber production. Populus species of closely related sections show extensive hybridization. Therefore, the systematics of the genus is rather complicated, especially for poplars of hybrid origin. We aimed to assess the efficiency of application of the sex-determining region (SDR) in addition to the nuclear and chloroplast genome loci traditionally used in phylogenetic studies of poplars to investigate relationships in sections Aigeiros Duby and Tacamahaca Spach. Targeted deep sequencing of NTS 5S rDNA, ITS, DSH 2, DSH 5, DSH 8, DSH 12, DSH 29, 6, 15, 16, X18, trnG-psbK-psbI, rps2-rpoC2, rpoC2-rpoC1, as well as SDR and ARR17 gene was performed for 379 poplars. The SDR and ARR17 gene together with traditionally used multicopy and single-copy loci of nuclear and chloroplast DNA allowed us to obtain a clustering that is most consistent with poplar systematics based on morphological data and to shed light on several controversial hypotheses about the origin of the studied taxa (for example, the inexpediency of separating P. koreana, P. maximowiczii, and P. suaveolens into different species). We present a scheme of relationships between species and hybrids of sections Aigeiros and Tacamahaca based on molecular genetic, morphological, and geographical data. The geographical proximity of species and, therefore, the possibility of hybridization between them appear to be more important than the affiliation of species to the same section. We speculate that sections Aigeiros and Tacamahaca are distinguished primarily on an ecological principle (plain and mountain poplars) rather than on a genetic basis. Joint analysis of sequencing data for the SDR and chloroplast genome loci allowed us to determine the ancestors of P. × petrovskoe - P. laurifolia (female tree) × P. × canadensis (male tree), and P. × rasumovskoe - P. nigra (female tree) × P. suaveolens (male tree). Thus, the efficiency of using the SDR for the study of poplars of sections Aigeiros and Tacamahaca and the prospects of its use for the investigation of species of the genus Populus were shown.

Cloning and Functional Analysis of the VfRR17 Gene from tung tree (Vernicia fordii)

Tung tree (Vernicia fordii) is one of the four major woody oilseed species in China. However, its fruit yield is severely hampered by the low number of female flowers and the imbalanced male-to-female flower ratio, which is a problem that restricts the development of the oilseed industry. Previous research has demonstrated that the exogenous application of cytokinins can significantly augment the number of female flowers, although the underlying regulatory mechanism remains elusive. To elucidate the involvement of VfRR17, a member of the A-type ARRs family, in the exogenous cytokinin regulation of flower sexual differentiation in tung tree, this study conducted a comprehensive bioinformatic analysis of the physicochemical properties, structural characteristics, and evolutionary relationships of the protein encoded by VfRR17. Additionally, gene function analysis was performed using subcellular localization, qRT-PCR, and genetic transformation techniques. The findings revealed that the VfRR17 gene's coding region spanned 471 bp, encoding an unstable protein of 156 amino acids with a relative molecular mass of 17.4 kDa. Phylogenetic analysis demonstrated a higher similarity between VfRR17 of the tung tree and the RR17 gene of Jatropha curcas, Ricinus communis, Hevea brasiliensis, and other species within the Euphorbiaceae family compared to other species, with the greatest similarity of 86% observed with the RR17 gene of Jatropha curcas. The qRT-PCR analysis indicated that VfRR17 exhibited high expression levels during the early stage of tung tree inflorescence buds following 6-BA treatment, peaking at 24 h and displaying a 3.47-fold increase compared to that at 0 h. In female and male flowers of the tung tree, the expression in female flowers during the 1 DBF period was significantly higher than in male flowers, exhibiting a difference of approximately 47.91-fold. Furthermore, notable differential expression was observed in the root, leaf, and petiole segments of the tung tree under low-temperature stress at the 12-h time point. In transgenic Arabidopsis, the VfRR17 lines and wild-type lines exhibited significantly different flowering times under an exogenous 6-BA treatment at a concentration of 2 mg/L, with the VfRR17 lines experiencing an 11-day delay compared to the wild-type lines. Additionally, the number of fruit pods in VfRR17 transgenic Arabidopsis lines was significantly reduced by 28 compared to the wild-type lines at a 6-BA concentration of 3 mg/L. These findings suggest that VfRR17 likely plays a critical role in regulating flower development in response to exogenous 6-BA, providing valuable insights into the mechanisms underlying exogenous 6-BA-mediated regulation of female flower development in the tung tree.

Development and Complex Application of Methods for the Identification of Mutations in the FAD3A and FAD3B Genes Resulting in the Reduced Content of Linolenic Acid in Flax Oil

Flax is grown worldwide for seed and fiber production. Linseed varieties differ in their oil composition and are used in pharmaceutical, food, feed, and industrial production. The field of application primarily depends on the content of linolenic (LIN) and linoleic (LIO) fatty acids. Inactivating mutations in the FAD3A and FAD3B genes lead to a decrease in the LIN content and an increase in the LIO content. For the identification of the three most common low-LIN mutations in flax varieties (G-to-A in exon 1 of FAD3A substituting tryptophan with a stop codon, C-to-T in exon 5 of FAD3A leading to arginine to a stop codon substitution, and C-to-T in exon 2 of FAD3B resulting in histidine to tyrosine substitution), three approaches were proposed: (1) targeted deep sequencing, (2) high resolution melting (HRM) analysis, (3) cleaved amplified polymorphic sequences (CAPS) markers. They were tested on more than a thousand flax samples of various types and showed promising results. The proposed approaches can be used in marker-assisted selection to choose parent pairs for crosses, separate heterogeneous varieties into biotypes, and select genotypes with desired homozygous alleles of the FAD3A and FAD3B genes at the early stages of breeding for the effective development of varieties with a particular LIN and LIO content, as well as in basic studies of the molecular mechanisms of fatty acid synthesis in flax seeds to select genotypes adequate to the tasks.

Genome Assembly and Sex-Determining Region of Male and Female Populus × sibirica

The genus Populus is presented by dioecious species, and it became a promising object to study the genetics of sex in plants. In this work, genomes of male and female Populus × sibirica individuals were sequenced for the first time. To achieve high-quality genome assemblies, we used Oxford Nanopore Technologies and Illumina platforms. A protocol for the isolation of long and pure DNA from young poplar leaves was developed, which enabled us to obtain 31 Gb (N50 = 21 kb) for the male poplar and 23 Gb (N50 = 24 kb) for the female one using the MinION sequencer. Genome assembly was performed with different tools, and Canu provided the most complete and accurate assemblies with a length of 818 Mb (N50 = 1.5 Mb) for the male poplar and 816 Mb (N50 = 0.5 Mb) for the female one. After polishing with Racon and Medaka (Nanopore reads) and then with POLCA (Illumina reads), assembly completeness was 98.45% (87.48% duplicated) for the male and 98.20% (76.77% duplicated) for the female according to BUSCO (benchmarking universal single-copy orthologs). A high proportion of duplicated BUSCO and the increased genome size (about 300 Mb above the expected) pointed at the separation of haplotypes in a large part of male and female genomes of P. × sibirica. Due to this, we were able to identify two haplotypes of the sex-determining region (SDR) in both assemblies; and one of these four SDR haplotypes, in the male genome, contained partial repeats of the ARR17 gene (Y haplotype), while the rest three did not (X haplotypes). The analysis of the male P. × sibirica SDR suggested that the Y haplotype originated from P. nigra, while the X haplotype is close to P. trichocarpa and P. balsamifera species. Moreover, we revealed a Populus-specific repeat that could be involved in translocation of the ARR17 gene or its part to the SDR of P. × sibirica and other Populus species. The obtained results expand our knowledge on SDR features in the genus Populus and poplar phylogeny.

Genome and Transcriptome Sequencing of Populus × sibirica Identified Sex-Associated Allele-Specific Expression of the CLC Gene

Transcriptome sequencing of leaves, catkin axes, and flowers from male and female trees of Populus × sibirica and genome sequencing of the same plants were performed for the first time. The availability of both genome and transcriptome sequencing data enabled the identification of allele-specific expression. Such an analysis was performed for genes from the sex-determining region (SDR). P. × sibirica is an intersectional hybrid between species from sections Aigeiros (Populus nigra) and Tacamahaca (Populus laurifolia, Populus suaveolens, or Populus × moskoviensis); therefore, a significant number of heterozygous polymorphisms were identified in the SDR that allowed us to distinguish between alleles. In the SDR, both allelic variants of the TCP (T-complex protein 1 subunit gamma), CLC (Chloride channel protein CLC-c), and MET1 (DNA-methyltransferase 1) genes were expressed in females, while in males, two allelic variants were expressed for TCP and MET1 but only one allelic variant prevailed for CLC. Targeted sequencing of TCP, CLC, and MET1 regions on a representative set of trees confirmed the sex-associated allele-specific expression of the CLC gene in generative and vegetative tissues of P. × sibirica. Our study brings new knowledge on sex-associated differences in Populus species.

Integrative genomics reveals paths to sex dimorphism in Salix purpurea L

Sex dimorphism and gene expression were studied in developing catkins in 159 F2 individuals from the bioenergy crop Salix purpurea, and potential mechanisms and pathways for regulating sex development were explored. Differential expression, eQTL, bisulfite sequencing, and network analysis were used to characterize sex dimorphism, detect candidate master regulator genes, and identify pathways through which the sex determination region (SDR) may mediate sex dimorphism. Eleven genes are presented as candidates for master regulators of sex, supported by gene expression and network analyses. These include genes putatively involved in hormone signaling, epigenetic modification, and regulation of transcription. eQTL analysis revealed a suite of transcription factors and genes involved in secondary metabolism and floral development that were predicted to be under direct control of the sex determination region. Furthermore, data from bisulfite sequencing and small RNA sequencing revealed strong differences in expression between males and females that would implicate both of these processes in sex dimorphism pathways. These data indicate that the mechanism of sex determination in Salix purpurea is likely different from that observed in the related genus Populus. This further demonstrates the dynamic nature of SDRs in plants, which involves a multitude of mechanisms of sex determination and a high rate of turnover.

A General Model to Explain Repeated Turnovers of Sex Determination in the Salicaceae

Dioecy, the presence of separate sexes on distinct individuals, has evolved repeatedly in multiple plant lineages. However, the specific mechanisms by which sex systems evolve and their commonalities among plant species remain poorly understood. With both XY and ZW sex systems, the family Salicaceae provides a system to uncover the evolutionary forces driving sex chromosome turnovers. In this study, we performed a genome-wide association study to characterize sex determination in two Populus species, P. euphratica and P. alba. Our results reveal an XY system of sex determination on chromosome 14 of P. euphratica, and a ZW system on chromosome 19 of P. alba. We further assembled the corresponding sex-determination regions, and found that their sex chromosome turnovers may be driven by the repeated translocations of a Helitron-like transposon. During the translocation, this factor may have captured partial or intact sequences that are orthologous to a type-A cytokinin response regulator gene. Based on results from this and other recently published studies, we hypothesize that this gene may act as a master regulator of sex determination for the entire family. We propose a general model to explain how the XY and ZW sex systems in this family can be determined by the same RR gene. Our study provides new insights into the diversification of incipient sex chromosomes in flowering plants by showing how transposition and rearrangement of a single gene can control sex in both XY and ZW systems.

Genetic diversity of SAD and FAD genes responsible for the fatty acid composition in flax cultivars and lines

BACKGROUND

Flax (Linum usitatissimum L.) is grown for fiber and seed in many countries. Flax cultivars differ in the oil composition and, depending on the ratio of fatty acids, are used in pharmaceutical, food, or paint industries. It is known that genes of SAD (stearoyl-ACP desaturase) and FAD (fatty acid desaturase) families play a key role in the synthesis of fatty acids, and some alleles of these genes are associated with a certain composition of flax oil. However, data on genetic polymorphism of these genes are still insufficient.

RESULTS

On the basis of the collection of the Institute for Flax (Torzhok, Russia), we formed a representative set of 84 cultivars and lines reflecting the diversity of fatty acid composition of flax oil. An approach for the determination of full-length sequences of SAD1, SAD2, FAD2A, FAD2B, FAD3A, and FAD3B genes using the Illumina platform was developed and deep sequencing of the 6 genes in 84 flax samples was performed on MiSeq. The obtained high coverage (about 400x on average) enabled accurate assessment of polymorphisms in SAD1, SAD2, FAD2A, FAD2B, FAD3A, and FAD3B genes and evaluation of cultivar/line heterogeneity. The highest level of genetic diversity was observed for FAD3A and FAD3B genes - 91 and 62 polymorphisms respectively. Correlation analysis revealed associations between particular variants in SAD and FAD genes and predominantly those fatty acids whose conversion they catalyze: SAD - stearic and oleic acids, FAD2 - oleic and linoleic acids, FAD3 - linoleic and linolenic acids. All except one low-linolenic flax cultivars/lines contained both the substitution of tryptophan to stop codon in the FAD3A gene and histidine to tyrosine substitution in the FAD3B gene, while samples with only one of these polymorphisms had medium content of linolenic acid and cultivars/lines without them were high-linolenic.

CONCLUSIONS

Genetic polymorphism of SAD and FAD genes was evaluated in the collection of flax cultivars and lines with diverse oil composition, and associations between particular polymorphisms and the ratio of fatty acids were revealed. The achieved results are the basis for the development of marker-assisted selection and DNA-based certification of flax cultivars.

Sequencing and Analysis of the Sex Determination Region of Populus trichocarpa

The ages and sizes of a sex-determination region (SDR) are difficult to determine in non-model species. Due to the lack of recombination and enrichment of repetitive elements in SDRs, the quality of assembly with short sequencing reads is universally low. Unique features present in the SDRs help provide clues about how SDRs are established and how they evolve in the absence of recombination. Several Populus species have been reported with a male heterogametic configuration of sex (XX/XY system) mapped on chromosome 19, but the exact location of the SDR has been inconsistent among species, and thus far, none of these SDRs has been fully assembled in a genomic context. Here we identify the Y-SDR from a Y-linked contig directly from a long-read PacBio assembly of a Populus trichocarpa male individual. We also identified homologous gene sequences in the SDR of P. trichocarpa and the SDR of the W chromosome in Salix purpurea. We show that inverted repeats (IRs) found in the Y-SDR and the W-SDR are lineage-specific. We hypothesize that, although the two IRs are derived from the same orthologous gene within each species, they likely have independent evolutionary histories. Furthermore, the truncated inverted repeats in P. trichocarpa may code for small RNAs that target the homologous gene for RNA-directed DNA methylation. These findings support the hypothesis that diverse sex-determining systems may be achieved through similar evolutionary pathways, thereby providing a possible mechanism to explain the lability of sex-determination systems in plants in general.

Data on genetic polymorphism of flax (Linum usitatissimum L.) pathogenic fungi of Fusarium, Colletotrichum, Aureobasidium, Septoria, and Melampsora genera

Being a valuable agricultural plant, flax (Linum usitatissimum L.) is used for oil and fiber production. However, the cultivation of this agriculture faces an urgent problem of flax susceptibility to fungal diseases. The most destructive ones are caused by the representatives of Fusarium, Colletotrichum, Aureobasidium, Septoria, and Melampsora genera, reducing flax yields significantly. To combat such pathogens effectively, it is of high importance to assess their genetic diversity that can be used to develop molecular markers to distinguish fungal genera and species. Morphological analysis traditionally carried out for fungal identification requires a given amount of time and tends to be difficult. In the present work, we determined the DNA sequences that are frequently used for phylogenetic studies in fungi - internal transcribed spacer (ITS) and beta-tubulin (tub2), translation elongation factor 1-alpha (tef1), RNA polymerase II largest subunit (RPB1), RNA polymerase II second largest subunit (RPB2), and minichromosome maintenance protein (MCM7) genes - for 203 flax fungal pathogens of Fusarium oxysporum, F. avenaceum, F. solani, F. sporotrichiella, F. moniliforme, F. culmorum, F. semitectum, F. gibbosum, Colletotrichum lini, Aureobasidium pullulans, Septoria linicola, and Melampsora lini species. The sequencing was performed using the Illumina MiSeq platform with a 300+300 bp kit, and on average, about 2350 reads per sample were obtained that allows accurate identification of the genetic polymorphism. Raw data are stored at the Sequence Read Archive under the accession number PRJNA596387. The obtained data can be used for fungal phylogenetic studies and the development of a PCR-based test system for flax pathogen identification.

A willow sex chromosome reveals convergent evolution of complex palindromic repeats

Background

Sex chromosomes have arisen independently in a wide variety of species, yet they share common characteristics, including the presence of suppressed recombination surrounding sex determination loci. Mammalian sex chromosomes contain multiple palindromic repeats across the non-recombining region that show sequence conservation through gene conversion and contain genes that are crucial for sexual reproduction. In plants, it is not clear if palindromic repeats play a role in maintaining sequence conservation in the absence of homologous recombination.

Results

Here we present the first evidence of large palindromic structures in a plant sex chromosome, based on a highly contiguous assembly of the W chromosome of the dioecious shrub Salix purpurea. The W chromosome has an expanded number of genes due to transpositions from autosomes. It also contains two consecutive palindromes that span a region of 200 kb, with conspicuous 20-kb stretches of highly conserved sequences among the four arms that show evidence of gene conversion. Four genes in the palindrome are homologous to genes in the sex determination regions of the closely related genus Populus, which is located on a different chromosome. These genes show distinct, floral-biased expression patterns compared to paralogous copies on autosomes.

Conclusion

The presence of palindromes in sex chromosomes of mammals and plants highlights the intrinsic importance of these features in adaptive evolution in the absence of recombination. Convergent evolution is driving both the independent establishment of sex chromosomes as well as their fine-scale sequence structure.