Duplicated genes evolve independently after polyploid formation in cotton

Allopolyploid subgenome identification and implications for evolutionary analysis

Whole-genome duplications (WGDs) are widespread genomic events in eukaryotes that are hypothesized to contribute to the evolutionary success of many lineages, including flowering plants, Saccharomyces yeast, and vertebrates. WGDs generally can be classified into autopolyploids (ploidy increase descended from one species) or allopolyploids (ploidy increase descended from multiple species). Assignment of allopolyploid progenitor species (called subgenomes in the polyploid) is important to understanding the biology and evolution of polyploids, including the asymmetric subgenome evolution following hybridization (biased fractionation). Here, I review the different methodologies used to identify the ancestors of allopolyploid subgenomes, discuss the advantages and disadvantages of these methods, and outline the implications of how these methods affect the subsequent evolutionary analysis of these genomes.

Phylogenetic and functional analysis of tiller angle control homeologs in allotetraploid cotton

Introduction

Plants can adapt their growth to optimize light capture in competitive environments, with branch angle being a crucial factor influencing plant phenotype and physiology. Decreased branch angles in cereal crops have been shown to enhance productivity in high-density plantings. The Tiller Angle Control (TAC1) gene, known for regulating tiller inclination in rice and corn, has been found to control branch angle in eudicots. Manipulating TAC1 in field crops like cotton offers the potential for improving crop productivity.

Methods

Using a homolog-based methodology, we examined the distribution of TAC1-related genes in cotton compared to other angiosperms. Furthermore, tissue-specific qPCR analysis unveiled distinct expression patterns of TAC1 genes in various cotton tissues. To silence highly expressed specific TAC1 homeologs in the stem, we applied CRISPR-Cas9 gene editing and Agrobacterium-mediated transformation, followed by genotyping and subsequent phenotypic validation of the mutants.

Results

Gene duplication events of TAC1 specific to the Gossypium lineage were identified, with 3 copies in diploid progenitors and 6 copies in allotetraploid cottons. Sequence analysis of the TAC1 homeologs in Gossypium hirsutum revealed divergence from other angiosperms with 1-2 copies, suggesting possible neo- or sub-functionalization for the duplicated copies. These TAC1 homeologs exhibited distinct gene expression patterns in various tissues over developmental time, with elevated expression of A11G109300 and D11G112200, specifically in flowers and stems, respectively. CRISPR-mediated loss of these TAC1 homeologous genes resulted in a reduction in branch angle and altered petiole angles, and a 5 to 10-fold reduction in TAC1 expression in the mutants, confirming their role in controlling branch and petiole angles. This research provides a promising strategy for genetically engineering branch and petiole angles in commercial cotton varieties, potentially leading to increased productivity.

Genome-Wide Identification and Expression Analysis of RLCK-VII Subfamily Genes Reveal Their Roles in Stress Responses of Upland Cotton

Receptor-like cytoplasmic kinase VII (RLCK-VII) subfamily members are vital players in plant innate immunity and are also involved in plant development and abiotic stress tolerance. As a widely cultivated textile crop, upland cotton (Gossypium hirsutum) attaches great importance to the cotton industry worldwide. To obtain details of the composition, phylogeny, and putative function of RLCK-VII genes in upland cotton, genome-wide identification, evolutionary event analysis, and expression pattern examination of RLCK-VII subfamily genes in G. hirsutum were performed. There are 129 RLCK-VII members in upland cotton (GhRLCKs) and they were divided into nine groups based on their phylogenetic relationships. The gene structure and sequence features are relatively conserved within each group, which were divided based on their phylogenetic relationships, and consistent with those in Arabidopsis. The phylogenetic analysis results showed that RLCK-VII subfamily genes evolved in plants before the speciation of Arabidopsis and cotton, and segmental duplication was the major factor that caused the expansion of GhRLCKs. The diverse expression patterns of GhRLCKs in response to abiotic stresses (temperature, salt, and drought) and V. dahliae infection were observed. The candidates that may be involved in cotton's response to these stresses are highlighted. GhRLCK7 (GhRLCK7A and D), which is notably induced by V. dahliae infection, was demonstrated to positively regulate cotton defense against V. dahliae by the loss-of-function assay in cotton. These findings shed light on the details of the RLCK-VII subfamily in cotton and provide a scaffold for the further function elucidation and application of GhRLCKs for the germplasm innovation of cotton.

Recent progression and future perspectives in cotton genomic breeding

Upland cotton is an important global cash crop for its long seed fibers and high edible oil and protein content. Progress in cotton genomics promotes the advancement of cotton genetics, evolutionary studies, functional genetics, and breeding, and has ushered cotton research and breeding into a new era. Here, we summarize high-impact genomics studies for cotton from the last 10 years. The diploid Gossypium arboreum and allotetraploid Gossypium hirsutum are the main focus of most genetic and genomic studies. We next review recent progress in cotton molecular biology and genetics, which builds on cotton genome sequencing efforts, population studies, and functional genomics, to provide insights into the mechanisms shaping abiotic and biotic stress tolerance, plant architecture, seed oil content, and fiber development. We also suggest the application of novel technologies and strategies to facilitate genome-based crop breeding. Explosive growth in the amount of novel genomic data, identified genes, gene modules, and pathways is now enabling researchers to utilize multidisciplinary genomics-enabled breeding strategies to cultivate "super cotton", synergistically improving multiple traits. These strategies must rise to meet urgent demands for a sustainable cotton industry.

Identification and Characterization of Malate Dehydrogenases in Tomato (Solanum lycopersicum L.)

Malate dehydrogenase, which facilitates the reversible conversion of malate to oxaloacetate, is essential for energy balance, plant growth, and cold and salt tolerance. However, the genome-wide study of the MDH family has not yet been carried out in tomato (Solanum lycopersicum L.). In this study, 12 MDH genes were identified from the S. lycopersicum genome and renamed according to their chromosomal location. The tomato MDH genes were split into five groups based on phylogenetic analysis and the genes that clustered together showed similar lengths, and structures, and conserved motifs in the encoded proteins. From the 12 tomato MDH genes on the chromosomes, three pairs of segmental duplication events involving four genes were found. Each pair of genes had a Ka/Ks ratio < 1, indicating that the MDH gene family of tomato was purified during evolution. Gene expression analysis exhibited that tomato MDHs were differentially expressed in different tissues, at various stages of fruit development, and differentially regulated in response to abiotic stresses. Molecular docking of four highly expressed MDHs revealed their substrate and co-factor specificity in the reversible conversion process of malate to oxaloacetate. Further, co-localization of tomato MDH genes with quantitative trait loci (QTL) of salt stress-related phenotypes revealed their broader functions in salt stress tolerance. This study lays the foundation for functional analysis of MDH genes and genetic improvement in tomato.

Global phylogeography of a pantropical mangrove genus Rhizophora

Rhizophora is a key genus for revealing the formation process of the pantropical distribution of mangroves. In this study, in order to fully understand the historical scenario of Rhizophora that achieved pantropical distribution, we conducted phylogeographic analyses based on nucleotide sequences of chloroplast and nuclear DNA as well as microsatellites for samples collected worldwide. Phylogenetic trees suggested the monophyly of each AEP and IWP lineages respectively except for R. samoensis and R. × selala. The divergence time between the two lineages was 10.6 million years ago on a dated phylogeny, and biogeographic stochastic mapping analyses supported these lineages separated following a vicariant event. These data suggested that the closure of the Tethys Seaway and the reduction in mangrove distribution followed by Mid-Miocene cooling were key factors that caused the linage diversification. Phylogeographic analyses also suggested the formation of the distinctive genetic structure at the AEP region across the American continents around Pliocene. Furthermore, long-distance trans-pacific dispersal occurred from the Pacific coast of American continents to the South Pacific and formed F1 hybrid, resulting in gene exchange between the IWP and AEP lineages after 11 million years of isolation. Considering the phylogeny and phylogeography with divergence time, a comprehensive picture of the historical scenario behind the pantropical distribution of Rhizophora is updated.

Genome-wide identification and characterization of olfactory receptor genes in common carp (Cyprinus carpio)

The environment contains a large extent of chemical information, which could be detected as olfactory sense. Olfactory in vertebrates plays important roles on many aspects during life time, including localizing prey or food, avoiding predators, mating behavior and social communication. Considering the essential role of olfactory receptors in the specific recognition of diverse stimuli, understanding the evolutionary dynamics of olfactory receptors in teleost means a lot, especially in the allotetraploid common carp, who has undergone the fourth whole-genome duplication event. Here, we identified the whole set of olfactory receptor genes in representative teleosts and found a significant contraction in common carp when compared with other teleosts. Odorant receptor genes (OR) occupy the most among four groups of olfactory receptors, including 33 functional genes and 16 pseudogenes. Furthermore, 6 trace amine-associated receptor (TAAR) genes (including 1 pseudogene), 7 odorant-related-A receptor genes, and 10 olfactory C family receptor genes (including 3 pseudogenes) were identified in common carp. Phylogenetic and motif analysis were performed to illustrate the phylogenetic relationship and structural conservation of teleost olfactory receptors. Selection pressure analysis suggested that olfactory receptor groups in common carp were all under relaxed purifying-selection. Additionally, gene expression divergences for olfactory receptor genes were investigated during embryonic development stages of common carp. We aim to determine the abundance of common carp olfactory receptor genes, explore the evolutionary fate and expression dynamics, and provide some genomic clues for the evolution of polyploid olfactory after whole-genome duplication and for future studies of teleost olfactory.

GhTULP34, a member of tubby-like proteins, interacts with GhSKP1A to negatively regulate plant osmotic stress

Tubby-like protein genes (TULPs), present in the form of large multigene families, play important roles in environmental stress. However, little is known regarding the TULP family genes in cotton. In this study, we systematically identified and analyzed the membership, characterization, and evolutionary relationship of TULPs in four species of cotton. Transcriptome analysis indicated that GhTULPs participate in environmental stress and cotton tissue development. At the same time, we also predicted and analyzed the potential molecular regulatory mechanisms and functions of TULPs. GhTULP34, as a candidate gene, significantly reduced the germination rate of transgenic Arabidopsis plants under salt stress, and inhibited root development and stomatal closure under mannitol stress. The yeast two-hybrid and luciferase (LUC) systems showed that GhTULP34 can interact with GhSKP1A, a subunit of the SCF-type (Skp1-Cullin-1-F-box) complex. This study will provide a basis and reference for future research on their roles in stress tolerance.

Genome-wide identification and evolution of Dof transcription factor family in cultivated and ancestral cotton species

The DNA-binding with one finger (Dof) proteins are transcription factors involved in many biological processes in plants. To predict the evolutionary pattern, a genome-wide in-silico analysis of Dof TFs family in diploid (Gossypium arboreum and Gossypium raimondii) and allotetraploid (Gossypium hirsutum and Gossypium barbadense) cotton species were carried out. In G. arboreum, we have identified 58 non-redundant genes encoding Dof proteins renamed as GaDof (G. arboreum Dof), 55 Dof genes were identified in G. raimondii (GrDof), 89 were predicted ffrom G. hirsutum (GhDof) and the highest, 110 Dof genes were identified in G. barbadense (GbDof). The phylogenetic analysis, physical location, gene structure, conserved domain analyses were also investigated for G. arboreum, G. raimondii, and G. hirsutum. The gene expression pattern in G. hirsutum, at different growth stages, revealing the probable involvement of some GhDof genes in growth and development. These genes may improve seed germination and growth in cotton.

The effects of hybridization and genome doubling in plant evolution via allopolyploidy

Polyploidy is a pervasive and recurring phenomenon across the tree of life, which occurred at variable time scales, ecological amplitudes and cell types, and is especially prominent in the evolutionary histories of plants. Importantly, many of the world's most important crops and noxious invasive weeds are recent polyploids. Polyploidy includes two major types, autopolyploidy, referring to doubling of a single species genome, and allopolyploidy referring to doubling of two or more merged genomes via biological hybridization of distinct but related species. The prevalence of both types of polyploidy implies that both genome doubling alone and doubling coupled with hybridization confer selective advantages over their diploid progenitors under specific circumstances. In cases of allopolyploidy, the two events, genome doubling and hybridization, have both advantages and disadvantages. Accumulated studies have established that, in allopolyploidy, some advantage(s) of doubling may compensate for the disadvantage(s) of hybridity and vice versa, although further study is required to validate generality of this trend. Some studies have also revealed a variety of non-Mendelian genetic and genomic consequences induced by doubling and hybridization separately or concertedly in nascent allopolyploidy; however, the significance of which to the immediate establishment and longer-term evolutionary success of allopolyploid species remain to be empirically demonstrated and ecologically investigated. This review aims to summarize recent advances in our understanding of the roles of hybridization and genome doubling, in separation and combination, in the evolution of allopolyploid genomes, as well as fruitful future research directions that are emerging from these studies.

Genome-Wide Identification, Sequence Variation, and Expression of the Glycerol-3-Phosphate Acyltransferase (GPAT) Gene Family in Gossypium

Cotton is an economically important crop grown for natural fiber and seed oil production. Cottonseed oil ranks third after soybean oil and colza oil in terms of edible oilseed tonnage worldwide. Glycerol-3-phosphate acyltransferase (GPAT) genes encode enzymes involved in triacylglycerol biosynthesis in plants. In the present study, 85 predicted GPAT genes were identified from the published genome data in Gossypium. Among them, 14, 16, 28, and 27 GPAT homologs were identified in G. raimondii, G. arboreum, G. hirsutum, and G. barbadense, respectively. Phylogenetic analysis revealed that a total of 108 GPAT genes from cotton, Arabidopsis and cacao could be classified into three groups. Furthermore, through comparison, the gene structure analyses indicated that GPAT genes from the same group were highly conserved between Arabidopsis and cotton. Segmental duplication could be the major driver for GPAT gene family expansion in the four cotton species above. Expression patterns of GhGPAT genes were diverse in different tissues. Most GhGPAT genes were induced or suppressed after salt or cold stress in Upland cotton. Eight GhGPAT genes were co-localized with oil and protein quantitative trait locus (QTL) regions. Thirty-two single nucleotide polymorphisms (SNPs) were detected from 12 GhGPAT genes, sixteen of which in nine GhGPAT genes were classified as synonymous, and sixteen SNPs in ten GhGPAT genes non-synonymous. Two SNP markers of the GhGPAT16 and GhGPAT26 genes were significantly correlated with cotton oil content in one of the three field tests. This study shed lights on the molecular evolutionary properties of GPAT genes in cotton, and provided reference for improvement of cotton response to abiotic stress and the genetic improvement of cotton oil content.

Insights into the Evolution of the New World Diploid Cottons (Gossypium, Subgenus Houzingenia) Based on Genome Sequencing

We employed phylogenomic methods to study molecular evolutionary processes and phylogeny in the geographically widely dispersed New World diploid cottons (Gossypium, subg. Houzingenia). Whole genome resequencing data (average of 33× genomic coverage) were generated to reassess the phylogenetic history of the subgenus and provide a temporal framework for its diversification. Phylogenetic analyses indicate that the subgenus likely originated following transoceanic dispersal from Africa about 6.6 Ma, but that nearly all of the biodiversity evolved following rapid diversification in the mid-Pleistocene (0.5-2.0 Ma), with multiple long-distance dispersals required to account for range expansion to Arizona, the Galapagos Islands, and Peru. Comparative analyses of cpDNAversus nuclear data indicate that this history was accompanied by several clear cases of interspecific introgression. Repetitive DNAs contribute roughly half of the total 880 Mb genome, but most transposable element families are relatively old and stable among species. In the genic fraction, pairwise synonymous mutation rates average 1% per Myr, with nonsynonymous changes being about seven times less frequent. Over 1.1 million indels were detected and phylogenetically polarized, revealing a 2-fold bias toward deletions over small insertions. We suggest that this genome down-sizing bias counteracts genome size growth by TE amplification and insertions, and helps explain the relatively small genomes that are restricted to this subgenus. Compared with the rate of nucleotide substitution, the rate of indel occurrence is much lower averaging about 17 nucleotide substitutions per indel event.

Origin and parental genome characterization of the allotetraploid Stylosanthes scabra Vogel (Papilionoideae, Leguminosae), an important legume pasture crop

BACKGROUNDS AND AIMS

The genus Stylosanthes includes nitrogen-fixing and drought-tolerant species of considerable economic importance for perennial pasture, green manure and land recovery. Stylosanthes scabra is adapted to variable soil conditions, being cultivated to improve pastures and soils worldwide. Previous studies have proposed S. scabra as an allotetraploid species (2n = 40) with a putative diploid A genome progenitor S. hamata or S. seabrana (2n = 20) and the B genome progenitor S. viscosa (2n = 20). We aimed to provide conclusive evidence for the origin of S. scabra.

METHODS

We performed fluorescence in situ hybridization (FISH) and genomic in situ hybridization (GISH) experiments and Illumina paired-end sequencing of S. scabra, S. hamata and S. viscosa genomic DNA, to assemble and compare complete ribosomal DNA (rDNA) units and chloroplast genomes. Plastome- and genome-wide single nucleotide variation detection was also performed.

KEY RESULTS

GISH and phylogenetic analyses of plastid DNA and rDNA sequences support that S. scabra is an allotetraploid formed from 0.63 to 0.52 million years ago (Mya), from progenitors with a similar genome structure to the maternal donor S. hamata and the paternal donor S. viscosa. FISH revealed a non-additive number of 35S rDNA sites in S. scabra compared with its progenitors, indicating the loss of one locus from A genome origin. In S. scabra, most 5S rDNA units were similar to S. viscosa, while one 5S rDNA site of reduced size most probably came from an A genome species as revealed by GISH and in silico analysis.

CONCLUSIONS

Our approach combined whole-plastome and rDNA assembly with additional cytogenetic analysis to shed light successfully on the allotetraploid origin of S. scabra. We propose a Middle Pleistocene origin for S. scabra involving species with maternal A and paternal B genomes. Our data also suggest that variation found in rDNA units in S. scabra and its progenitors reveals differences that can be explained by homogenization, deletion and amplification processes that have occurred since its origin.

Evolutionary dynamics of 3D genome architecture following polyploidization in cotton

The formation of polyploids significantly increases the complexity of transcriptional regulation, which is expected to be reflected in sophisticated higher-order chromatin structures. However, knowledge of three-dimensional (3D) genome structure and its dynamics during polyploidization remains poor. Here, we characterize 3D genome architectures for diploid and tetraploid cotton, and find the existence of A/B compartments and topologically associated domains (TADs). By comparing each subgenome in tetraploids with its extant diploid progenitor, we find that genome allopolyploidization has contributed to the switching of A/B compartments and the reorganization of TADs in both subgenomes. We also show that the formation of TAD boundaries during polyploidization preferentially occurs in open chromatin, coinciding with the deposition of active chromatin modification. Furthermore, analysis of inter-subgenomic chromatin interactions has revealed the spatial proximity of homoeologous genes, possibly associated with their coordinated expression. This study advances our understanding of chromatin organization in plants and sheds new light on the relationship between 3D genome evolution and transcriptional regulation.

Comparative genomic and transcriptomic analyses of Family-1 UDP glycosyltransferase in three Brassica species and Arabidopsis indicates stress-responsive regulation

In plants, UGTs (UDP-glycosyltransferases) glycosylate various phytohormones and metabolites in response to biotic and abiotic stresses. Little is known about stress-responsive glycosyltransferases in plants. Therefore, it is important to understand the genomic and transcriptomic portfolio of plants with regard to biotic and abiotic stresses. Here, we identified 140, 154, and 251 putative UGTs in Brassica rapa, Brassica oleracea, and Brassica napus, respectively, and clustered them into 14 major phylogenetic groups (A–N). Fourteen major KEGG pathways and 24 biological processes were associated with the UGTs, highlighting them as unique modulators against environmental stimuli. Putative UGTs from B. rapa and B. oleracea showed a negative selection pressure and biased gene fractionation pattern during their evolution. Polyploidization increased the intron proportion and number of UGT-containing introns among Brassica. The putative UGTs were preferentially expressed in developing tissues and at the senescence stage. Differential expression of up- and down-regulated UGTs in response to phytohormone treatments, pathogen responsiveness and abiotic stresses, inferred from microarray and RNA-Seq data in Arabidopsis and Brassica broaden the glycosylation impact at the molecular level. This study identifies unique candidate UGTs for the manipulation of biotic and abiotic stress pathways in Brassica and Arabidopsis.

Identification and Characterization of miRNA Transcriptome in Asiatic Cotton (Gossypium arboreum) Using High Throughput Sequencing

MicroRNAs (miRNAs) are small 20-24nt molecules that have been well studied over the past decade due to their important regulatory roles in different cellular processes. The mature sequences are more conserved across vast phylogenetic scales than their precursors and some are conserved within entire kingdoms, hence, their loci and function can be predicted by homology searches. Different studies have been performed to elucidate miRNAs using de novo prediction methods but due to complex regulatory mechanisms or false positive in silico predictions, not all of them express in reality and sometimes computationally predicted mature transcripts differ from the actual expressed ones. With the availability of a complete genome sequence of Gossypium arboreum, it is important to annotate the genome for both coding and non-coding regions using high confidence transcript evidence, for this cotton species that is highly resistant to various biotic and abiotic stresses. Here we have analyzed the small RNA transcriptome of G. arboreum leaves and provided genome annotation of miRNAs with evidence from miRNA/miRNA∗ transcripts. A total of 446 miRNAs clustered into 224 miRNA families were found, among which 48 families are conserved in other plants and 176 are novel. Four short RNA libraries were used to shortlist best predictions based on high reads per million. The size, origin, copy numbers and transcript depth of all miRNAs along with their isoforms and targets has been reported. The highest gene copy number was observed for gar-miR7504 followed by gar-miR166, gar-miR8771, gar-miR156, and gar-miR7484. Altogether, 1274 target genes were found in G. arboreum that are enriched for 216 KEGG pathways. The resultant genomic annotations are provided in UCSC, BED format.

Genome-wide identification and characterization of phospholipase C gene family in cotton (Gossypium spp.)

Phospholipase C (PLC) are important regulatory enzymes involved in several lipid and Ca²⁺-dependent signaling pathways. Previous studies have elucidated the versatile roles of PLC genes in growth, development and stress responses of many plants, however, the systematic analyses of PLC genes in the important fiber-producing plant, cotton, are still deficient. In this study, through genome-wide survey, we identified twelve phosphatidylinositol-specific PLC (PI-PLC) and nine non-specific PLC (NPC) genes in the allotetraploid upland cotton Gossypium hirsutum and nine PI-PLC and six NPC genes in two diploid cotton G. arboretum and G.raimondii, respectively. The PI-PLC and NPC genes of G. hirsutum showed close phylogenetic relationship with their homologous genes in the diploid cottons and Arabidopsis. Segmental and tandem duplication contributed greatly to the formation of the gene family. Expression profiling indicated that few of the PLC genes are constitutely expressed, whereas most of the PLC genes are preferentially expressed in specific tissues and abiotic stress conditions. Promoter analyses further implied that the expression of these PLC genes might be regulated by MYB transcription factors and different phytohormones. These results not only suggest an important role of phospholipase C members in cotton plant development and abiotic stress response but also provide good candidate targets for future molecular breeding of superior cotton cultivars.

Genome-Wide Identification and Expression Analysis of the Biotin Carboxyl Carrier Subunits of Heteromeric Acetyl-CoA Carboxylase in Gossypium

Acetyl-CoA carboxylase is an important enzyme, which catalyzes acetyl-CoA's carboxylation to produce malonyl-CoA and to serve as a committed step for de novo fatty acid biosynthesis in plastids. In this study, 24 putative cotton BCCP genes were identified based on the lately published genome data in Gossypium. Among them, 4, 4, 8, and 8 BCCP homologs were identified in Gossypium raimondii, G. arboreum, G. hirsutum, and G. barbadense, respectively. These genes were divided into two classes based on a phylogenetic analysis. In each class, these homologs were relatively conserved in gene structure and motifs. The chromosomal distribution pattern revealed that all the BCCP genes were distributed equally on corresponding chromosomes or scaffold in the four cotton species. Segmental duplication was a predominant duplication event in both of G. hirsutum and G. barbadense. The analysis of the expression profile showed that 8 GhBCCP genes expressed in all the tested tissues with changed expression levels, and GhBCCP genes belonging to class II were predominantly expressed in developing ovules. Meanwhile, the expression analysis for the 16 cotton BCCP genes from G. raimondii, G. arboreum and G. hirsutum showed that they were induced or suppressed by cold or salt stress, and their expression patterns varied among different tissues. These findings will help to determine the functional and evolutionary characteristics of the BCCP genes in Gossypium species.

Polyploidy and the relationship between leaf structure and function: implications for correlated evolution of anatomy, morphology, and physiology in Brassica

BACKGROUND

Polyploidy is well studied from a genetic and genomic perspective, but the morphological, anatomical, and physiological consequences of polyploidy remain relatively uncharacterized. Whether these potential changes bear on functional integration or are idiosyncratic remains an open question. Repeated allotetraploid events and multiple genomic combinations as well as overlapping targets of artificial selection make the Brassica triangle an excellent system for exploring variation in the connection between plant structure (anatomy and morphology) and function (physiology). We examine phenotypic integration among structural aspects of leaves including external morphology and internal anatomy with leaf-level physiology among several species of Brassica. We compare diploid and allotetraploid species to ascertain patterns of phenotypic correlations among structural and functional traits and test the hypothesis that allotetraploidy results in trait disintegration allowing for transgressive phenotypes and additional evolutionary and crop improvement potential.

RESULTS

Among six Brassica species, we found significant effects of species and ploidy level for morphological, anatomical and physiological traits. We identified three suites of intercorrelated traits in both diploid parents and allotetraploids: Morphological traits (such as leaf area and perimeter) anatomic traits (including ab- and ad- axial epidermis) and aspects of physiology. In general, there were more correlations between structural and functional traits for allotetraploid hybrids than diploid parents. Parents and hybrids did not have any significant structure-function correlations in common. Of particular note, there were no significant correlations between morphological structure and physiological function in the diploid parents. Increased phenotypic integration in the allotetraploid hybrids may be due, in part, to increased trait ranges or simply different structure-function relationships.

CONCLUSIONS

Genomic and chromosomal instability in early generation allotetraploids may allow Brassica species to explore new trait space and potentially reach higher adaptive peaks than their progenitor species could, despite temporary fitness costs associated with unstable genomes. The trait correlations that disappear after hybridization as well as the novel trait correlations observed in allotetraploid hybrids may represent relatively evolutionarily labile associations and therefore could be ideal targets for artificial selection and crop improvement.

Comparative Genome-Wide Analysis of the Malate Dehydrogenase Gene Families in Cotton

Malate dehydrogenases (MDHs) play crucial roles in the physiological processes of plant growth and development. In this study, 13 and 25 MDH genes were identified from Gossypium raimondii and Gossypium hirsutum, respectively. Using these and 13 previously reported Gossypium arboretum MDH genes, a comparative molecular analysis between identified MDH genes from G. raimondii, G. hirsutum, and G. arboretum was performed. Based on multiple sequence alignments, cotton MDHs were divided into five subgroups: mitochondrial MDH, peroxisomal MDH, plastidial MDH, chloroplastic MDH and cytoplasmic MDH. Almost all of the MDHs within the same subgroup shared similar gene structure, amino acid sequence, and conserved motifs in their functional domains. An analysis of chromosomal localization suggested that segmental duplication played a major role in the expansion of cotton MDH gene families. Additionally, a selective pressure analysis indicated that purifying selection acted as a vital force in the evolution of MDH gene families in cotton. Meanwhile, an expression analysis showed the distinct expression profiles of GhMDHs in different vegetative tissues and at different fiber developmental stages, suggesting the functional diversification of these genes in cotton growth and fiber development. Finally, a promoter analysis indicated redundant but typical cis-regulatory elements for the potential functions and stress activity of many MDH genes. This study provides fundamental information for a better understanding of cotton MDH gene families and aids in functional analyses of the MDH genes in cotton fiber development.

Genome-wide characterization and comparative analysis of the MLO gene family in cotton

In plants, MLO (Mildew Locus O) gene encodes a plant-specific seven transmembrane (TM) domain protein involved in several cellular processes, including susceptibility to powdery mildew (PM). In this study, a genome-wide characterization of the MLO gene family in G. raimondii L., G. arboreum L. and G. hirsutum L. was performed. In total, 22, 17 and 38 homologous sequences were identified for each species, respectively. Gene organization, including chromosomal location, gene clustering and gene duplication, was investigated. Homologues related to PM susceptibility in upland cotton were inferred by phylogenetic relationships with functionally characterized MLO proteins. To conduct a comparative analysis between MLO candidate genes from G. raimondii L., G. arboreum L. and G. hirsutum L., orthologous relationships and conserved synteny blocks were constructed. The transcriptional variation of 38 GhMLO genes in response to exogenous application of salt, mannitol (Man), abscisic acid (ABA), ethylene (ETH), jasmonic acid (JA) and salicylic acid (SA) was monitored. Further studies should be conducted to elucidate the functions of MLO genes in PM susceptibility and phytohormone signalling pathways.

Molecular evolution and phylogenetic analysis of eight COL superfamily genes in group I related to photoperiodic regulation of flowering time in wild and domesticated cotton (Gossypium) species

Flowering time is an important ecological trait that determines the transition from vegetative to reproductive growth. Flowering time in cotton is controlled by short-day photoperiods, with strict photoperiod sensitivity. As the CO-FT (CONSTANS-FLOWER LOCUS T) module regulates photoperiodic flowering in several plants, we selected eight CONSTANS genes (COL) in group I to detect their expression patterns in long-day and short-day conditions. Further, we individually cloned and sequenced their homologs from 25 different cotton accessions and one outgroup. Finally, we studied their structures, phylogenetic relationship, and molecular evolution in both coding region and three characteristic domains. All the eight COLs in group I show diurnal expression. In the orthologous and homeologous loci, each gene structure in different cotton species is highly conserved, while length variation has occurred due to insertions/deletions in intron and/or exon regions. Six genes, COL2 to COL5, COL7 and COL8, exhibit higher nucleotide diversity in the D-subgenome than in the A-subgenome. The Ks values of 98.37% in all allotetraploid cotton species examined were higher in the A-D and At-Dt comparison than in the A-At and D-Dt comparisons, and the Pearson’s correlation coefficient (r) of Ks between A vs. D and At vs. Dt also showed positive, high correlations, with a correlation coefficient of at least 0.797. The nucleotide polymorphism in wild species is significantly higher compared to G. hirsutum and G. barbadense, indicating a genetic bottleneck associated with the domesticated cotton species. Three characteristic domains in eight COLs exhibit different evolutionary rates, with the CCT domain highly conserved, while the B-box and Var domain much more variable in allotetraploid species. Taken together, COL1, COL2 and COL8 endured greater selective pressures during the domestication process. The study improves our understanding of the domestication-related genes/traits during cotton evolutionary process.

Construction of cytogenetic map of Gossypium herbaceum chromosome 1 and its integration with genetic maps

Background

Cytogenetic map can provide not only information of the genome structure, but also can build a solid foundation for genetic research. With the developments of molecular and cytogenetic studies in cotton (Gossypium), the construction of cytogenetic map is becoming more and more imperative.

Results

A cytogenetic map of chromosome 1 (A₁01) of Gossypium herbaceum (A₁) which includes 10 bacterial artificial chromosome (BAC) clones was constructed by using fluorescent in situ hybridization (FISH). Meanwhile, comparison and analysis were made for the cytogenetic map of chromosome 1 (A₁01) of G. herbaceum with four genetic linkage maps of chromosome 1 (A_h01) of G. hirsutum ((AD)₁) and one genetic linkage map of chromosome 1 of (A₁01) G. arboreum (A₂). The 10 BAC clones were also used to be localized on G. raimondii (D₅) chromosome 1 (D₅01), and 2 of them showed clear unique hybridized signals. Furthermore, these 2 BAC clones were also shown localized on chromosome 1 of both A sub-genome and D sub-genome of G. hirsutum.

Conclusion

The comparison of the cytogenetic map with genetic linkage maps showed that most of the identified marker-tagged BAC clones appearing same orders in different maps except three markers showing different positions, which might indicate chromosomal segmental rearrangements. The positions of the 2 BAC clones which were localized on A_h01 and D_h01 chromosomes were almost the same as that on A₁01 and D₅01 chromosomes. The corresponding anchored SSR markers of these 2 BAC clones were firstly found to be localized on chromosome D₅01 (D_h01) as they were not seen mapped like this in any genetic map reported.

Mapping genomic loci for cotton plant architecture, yield components, and fiber properties in an interspecific (Gossypium hirsutum L. × G. barbadense L.) RIL population

A quantitative trait locus (QTL) mapping was conducted to better understand the genetic control of plant architecture (PA), yield components (YC), and fiber properties (FP) in the two cultivated tetraploid species of cotton (Gossypium hirsutum L. and G. barbadense L.). One hundred and fifty-nine genomic regions were identified on a saturated genetic map of more than 2,500 SSR and SNP markers, constructed with an interspecific recombinant inbred line (RIL) population derived from the genetic standards of the respective cotton species (G. hirsutum acc. TM-1 × G. barbadense acc. 3-79). Using the single nonparametric and MQM QTL model mapping procedures, we detected 428 putative loci in the 159 genomic regions that confer 24 cotton traits in three diverse production environments [College Station F&B Road (FB), TX; Brazos Bottom (BB), TX; and Shafter (SH), CA]. These putative QTL loci included 25 loci for PA, 60 for YC, and 343 for FP, of which 3, 12, and 60, respectively, were strongly associated with the traits (LOD score ≥ 3.0). Approximately 17.7 % of the PA putative QTL, 32.9 % of the YC QTL, and 48.3 % of the FP QTL had trait associations under multiple environments. The At subgenome (chromosomes 1-13) contributed 72.7 % of loci for PA, 46.2 % for YC, and 50.4 % for FP while the Dt subgenome (chromosomes 14-26) contributed 27.3 % of loci for PA, 53.8 % for YC, and 49.6 % for FP. The data obtained from this study augment prior evidence of QTL clusters or gene islands for specific traits or biological functions existing in several non-homoeologous cotton chromosomes. DNA markers identified in the 159 genomic regions will facilitate further dissection of genetic factors underlying these important traits and marker-assisted selection in cotton.

Consequences of Whole-Genome Triplication as Revealed by Comparative Genomic Analyses of the Wild Radish Raphanus raphanistrum and Three Other Brassicaceae Species

Polyploidization events are frequent among flowering plants, and the duplicate genes produced via such events contribute significantly to plant evolution. We sequenced the genome of wild radish (Raphanus raphanistrum), a Brassicaceae species that experienced a whole-genome triplication event prior to diverging from Brassica rapa. Despite substantial gene gains in these two species compared with Arabidopsis thaliana and Arabidopsis lyrata, ∼70% of the orthologous groups experienced gene losses in R. raphanistrum and B. rapa, with most of the losses occurring prior to their divergence. The retained duplicates show substantial divergence in sequence and expression. Based on comparison of A. thaliana and R. raphanistrum ortholog floral expression levels, retained radish duplicates diverged primarily via maintenance of ancestral expression level in one copy and reduction of expression level in others. In addition, retained duplicates differed significantly from genes that reverted to singleton state in function, sequence composition, expression patterns, network connectivity, and rates of evolution. Using these properties, we established a statistical learning model for predicting whether a duplicate would be retained postpolyploidization. Overall, our study provides new insights into the processes of plant duplicate loss, retention, and functional divergence and highlights the need for further understanding factors controlling duplicate gene fate.

Upland cotton gene GhFPF1 confers promotion of flowering time and shade-avoidance responses in Arabidopsis thaliana

Extensive studies on floral transition in model species have revealed a network of regulatory interactions between proteins that transduce and integrate developmental and environmental signals to promote or inhibit the transition to flowering. Previous studies indicated FLOWERING PROMOTING FACTOR 1 (FPF1) gene was involved in the promotion of flowering, but the molecular mechanism was still unclear. Here, FPF1 homologous sequences were screened from diploid Gossypium raimondii L. (D-genome, n = 13) and Gossypium arboreum L. genome (A-genome, n = 13) databases. Orthologous genes from the two species were compared, suggesting that distinctions at nucleic acid and amino acid levels were not equivalent because of codon degeneracy. Six FPF1 homologous genes were identified from the cultivated allotetraploid Gossypium hirsutum L. (AD-genome, n = 26). Analysis of relative transcripts of the six genes in different tissues revealed that this gene family displayed strong tissue-specific expression. GhFPF1, encoding a 12.0-kDa protein (Accession No: KC832319) exerted more transcripts in floral apices of short-season cotton, hinting that it could be involved in floral regulation. Significantly activated APETALA 1 and suppressed FLOWERING LOCUS C expression were induced by over-expression of GhFPF1 in the Arabidopsis Columbia-0 ecotype. In addition, transgenic Arabidopsis displayed a constitutive shade-avoiding phenotype that is characterized by long hypocotyls and petioles, reduced chlorophyll content, and early flowering. We propose that GhFPF1 may be involved in flowering time control and shade-avoidance responses.

Phylogeny of a genomically diverse group of elymus (poaceae) allopolyploids reveals multiple levels of reticulation

The grass tribe Triticeae (=Hordeeae) comprises only about 300 species, but it is well known for the economically important crop plants wheat, barley, and rye. The group is also recognized as a fascinating example of evolutionary complexity, with a history shaped by numerous events of auto- and allopolyploidy and apparent introgression involving diploids and polyploids. The genus Elymus comprises a heterogeneous collection of allopolyploid genome combinations, all of which include at least one set of homoeologs, designated St, derived from Pseudoroegneria. The current analysis includes a geographically and genomically diverse collection of 21 tetraploid Elymus species, and a single hexaploid species. Diploid and polyploid relationships were estimated using four molecular data sets, including one that combines two regions of the chloroplast genome, and three from unlinked nuclear genes: phosphoenolpyruvate carboxylase, β-amylase, and granule-bound starch synthase I. Four gene trees were generated using maximum likelihood, and the phylogenetic placement of the polyploid sequences reveals extensive reticulation beyond allopolyploidy alone. The trees were interpreted with reference to numerous phenomena known to complicate allopolyploid phylogenies, and introgression was identified as a major factor in their history. The work illustrates the interpretation of complicated phylogenetic results through the sequential consideration of numerous possible explanations, and the results highlight the value of careful inspection of multiple independent molecular phylogenetic estimates, with particular focus on the differences among them.

Molecular characterisation of two homoeologous elicitor-responsive lipin genes in cotton

The identification and molecular characterisation of two lipin-like gene copies (GhLIPN) in cotton, Gossypium hirsutum, an allotetraploid derived from two progenitor diploid Gossypium species, is described. Sequence analyses of the GhLIPN copies, designated GhLIPN-1 and -2, revealed that they contain 11 exons, separated by ten introns. They each have a 2,643 bp open reading frame that encodes 880 aa proteins, and share a 97.7 and 95.5 % sequence similarity at the translated nucleotide and amino acid level, respectively. The GhLIPN genes have a distinct domain architecture consisting of an archetypical N-terminal lipin domain, followed by a haloacid dehalogenase (HAD) domain towards the C-terminus. A Southern blot did not distinguish between the two gene copies, which suggests that they may be homoeologs rather than paralogs. GhLIPN-2 is more similar to a homoeologous sequence from G. raimondii, representing the ancestral D-genome, compared to GhLIPN-1 that matches G. herbaceum and that represents the A-genome. Our data indicate that GhLIPN-1 and GhLIPN-2 are homoeologs that derive from the A- and the D-diploid genomes, respectively. The promoter sequences of GhLIPN-1 and -2 differ by 56 %, as a result of multiple indels. In silico analysis of the promoter regions revealed that both genes contain numerous putative defence-related and elicitor-responsive cis-elements that support a role for GhLIPN in defence responses. Relative quantification real-time PCR confirmed the up-regulation in response to a cell-wall-derived V. dahliae elicitor, which supported the association of GhLIPN with defence signalling. The results add a new dimension to the proposed roles of lipins in plants by suggesting that lipins may have a role in defence signalling.

C, Percy R, Stewart JM. Genetic diversity and population structure of cotton (Gossypiumspp.) of the New World assessed by SSR markers

A global analysis of cotton (Gossypium spp.) genetic diversity is the first step to understanding its geographical distribution, dissemination, genetic relatedness, and population structure. To assess the genetic diversity and population structure in Gossypium species, 111 cotton accessions representing five allotetraploids (AD1–AD5genomes), 23 Asiatic diploids of the Old World (A1and A2genomes), and 82 diploids of the New World subgenus Houzingenia (D1–D11genomes) species were assessed using simple sequence repeats (SSR) markers with wide genome coverage. The mean genetic distance (GD) between the two most important New World tetraploid cottons (Upland (Gossypium hirsutum L.) and Pima (Gossypium barbadense L.)) was 0.39. Among the three shrub type sections (Houzingenia, Integrifolia, and Caducibracteolata) and three arborescent sections (Erioxylum, Selera, and Austroamericana), the GD ranged between 0.19 and 0.41. Phylogenetic analyses clustered all species into distinct phylogenetic groups, which were consistent with genomic origin, evolutionary history, and geographic distribution or ecotypes of these accessions, suggesting the existence of clear structured strata. With all of the genomes, the highest statistical analysis of Structure test through measurements of ad hoc (ΔK) occurred at K = 2, with group Q1 with the Old World diploid A genomes and with group Q2 with all the New World diploids of the D genome. AD genome accessions shared nearly equal alleles from both Q1 and Q2 groups. With all of the diploids of the New World D genomes, the highest value of ΔK occurred at K = 5. These results are consistent with the fundamental knowledge of tetraploid AD-genome formation and the rapid radiation of the American diploid cotton linage that took place somewhere in southwestern Mexico, followed by a differentiation–speciation during angiosperm evolution. In addition, SSR markers provide an alternative solution for distinguishing phylogenetic relationships between accessions of different ecotypes and for elucidating population structure of cottons of the New World.

Analysis of complete nucleotide sequences of 12 Gossypium chloroplast genomes: origin and evolution of allotetraploids

BACKGROUND

Cotton (Gossypium spp.) is a model system for the analysis of polyploidization. Although ascertaining the donor species of allotetraploid cotton has been intensively studied, sequence comparison of Gossypium chloroplast genomes is still of interest to understand the mechanisms underlining the evolution of Gossypium allotetraploids, while it is generally accepted that the parents were A- and D-genome containing species. Here we performed a comparative analysis of 13 Gossypium chloroplast genomes, twelve of which are presented here for the first time.

METHODOLOGY/PRINCIPAL FINDINGS

The size of 12 chloroplast genomes under study varied from 159,959 bp to 160,433 bp. The chromosomes were highly similar having >98% sequence identity. They encoded the same set of 112 unique genes which occurred in a uniform order with only slightly different boundary junctions. Divergence due to indels as well as substitutions was examined separately for genome, coding and noncoding sequences. The genome divergence was estimated as 0.374% to 0.583% between allotetraploid species and A-genome, and 0.159% to 0.454% within allotetraploids. Forty protein-coding genes were completely identical at the protein level, and 20 intergenic sequences were completely conserved. The 9 allotetraploids shared 5 insertions and 9 deletions in whole genome, and 7-bp substitutions in protein-coding genes. The phylogenetic tree confirmed a close relationship between allotetraploids and the ancestor of A-genome, and the allotetraploids were divided into four separate groups. Progenitor allotetraploid cotton originated 0.43-0.68 million years ago (MYA).

CONCLUSION

Despite high degree of conservation between the Gossypium chloroplast genomes, sequence variations among species could still be detected. Gossypium chloroplast genomes preferred for 5-bp indels and 1-3-bp indels are mainly attributed to the SSR polymorphisms. This study supports that the common ancestor of diploid A-genome species in Gossypium is the maternal source of extant allotetraploid species and allotetraploids have a monophyletic origin. G. hirsutum AD1 lineages have experienced more sequence variations than other allotetraploids in intergenic regions. The available complete nucleotide sequences of 12 Gossypium chloroplast genomes should facilitate studies to uncover the molecular mechanisms of compartmental co-evolution and speciation of Gossypium allotetraploids.

Molecular evolution and phylogenetic analysis of genes related to cotton fibers development from wild and domesticated cotton species in Gossypium

The domestication of both diploid and tetraploid cotton species was carried out for fiber utilization. To understand the origin and domestication of fibers, 18 genes related to fiber development were individually cloned and sequenced from 22 different cotton species. Their structures, phylogenetic relationship and molecular evolution were further studied. In the orthologous and homeologous loci of the 18 genes, the sequence and structure of 72.22% were conserved and 27.78% were diverse. Tree topologies constructed based on the combined sequences showed that all 13 D-genome species were congruent with Fryxell's subsection taxonomy, the A- and D-subgenomes independently evolved in the allopolyploid after polyploid formation, and Gossypium raimondii had the closest relationship with all allotetraploids of D-subgenomes. The molecular evolutionary rates revealed approximately equivalent rates among different D-genome species, and purifying selection acted on all genes in the wild D-genome species. Among orthologs and homeologs, the D-subgenomes had higher evolutionary rates than the A-subgenomes in tetraploid cotton species, and the cultivars had higher evolutionary rates than either the semi-domesticated or wild species. Our study revealed that human domestication altered the molecular evolutionary pattern of genes related to fiber development, and Gossypium hirsutum endured greater selective pressures than Gossypium barbadense during the domestication process.

Cloning and characterization of homeologous cellulose synthase catalytic subunit 2 genes from allotetraploid cotton (Gossypium hirsutum L.)

Cellulose synthase catalytic subunits (CesAs) are the catalytic sites within a multisubunit complex for cellulose biosynthesis in plants. CesAs have been extensively studied in diploid plants, but are not well characterized in polyploid plants. Gossypium hirsutum is an allotetraploid cotton species producing over 90% of the world's cotton fibers. Although G. hirsutum CesAs (GhCesAs) are responsible for cellulose production in cotton fiber, very limited numbers of GhCesA genes have been identified. Here, we report isolating and characterizing a pair of homeologous CesA2 genes and their full-length cDNAs from allotetraploid cotton. The GhCesA2-A(T) gene from the A-subgenome and GhCesA2-D(T) gene from the D-subgenome were screened from a G. hirsutum BAC library. These genes shared 92% sequence similarity throughout the entire sequence. The coding sequences were nearly identical, and the deduced amino acid sequences from GhCesA2-A(T) (1,039 amino acids) and GhCesA2-D(T) (1,040 amino acids) were identical except four amino acids, whereas the noncoding sequences showed divergence. Sequence analyses showed that all exons of GhCesA2-A(T) contained consensus splice donor dinucleotides, but one exon in GhCesA2-D(T) contained nonconsensus splice donor dinucleotides. Although the nonconsensus splice donor dinucleotides were previously suggested to be involved in alternative splice or pseudogenization, our results showed that a majority of GhCesA2-A(T) and GhCesA2-D(T) transcripts consisted of functional and full-length transcripts with little evidence for alternative mRNA isoforms in developing cotton fibers. Expression analyses showed that GhCesA2-A(T) and GhCesA2-D(T) shared common temporal and spatial expression patterns, and they were highly and preferentially expressed during the cellulose biosynthesis stage in developing cotton fibers. The observations of higher expression levels of both GhCesA2-A(T) and GhCesA2-D(T) in developing fibers of one near-isogenic line (NIL) with higher fiber bundle strength over the other NIL with lower fiber bundle strength suggested that the differential expression of genes associated with secondary cell wall cellulose biosynthesis in developing fiber might affect cotton fiber properties.

Molecular evolution of the clustered MIC-3 multigene family of Gossypium species

The Gossypium MIC-3 (Meloidogyne Induced Cotton-3) gene family is of great interest for molecular evolutionary studies because of its uniqueness to Gossypium species, multi-gene content, clustered localization, and root-knot nematode resistance-associated features. Molecular evolution of the MIC-3 gene family was studied in 15 tetraploid and diploid Gossypium genotypes that collectively represent seven phylogenetically distinct genomes. Synonymous (d(S)) and non-synonymous (d(N)) nucleotide substitution rates suggest that the second of the two exons of the MIC-3 genes has been under strong positive selection pressure, while the first exon has been under strong purifying selection to preserve function. Based on nucleotide substitution rates, we conclude that MIC-3 genes are evolving by a birth-and-death process and that a 'gene amplification' mechanism has helped to retain all duplicate copies, which best fits with the "bait and switch" model of R-gene evolution. The data indicate MIC-3 gene duplication events occurred at various rates, once per 1 million years (MY) in the allotetraploids, once per ~2 MY in the A/F genome clade, and once per ~8 MY in the D-genome clade. Variations in the MIC-3 gene family seem to reflect evolutionary selection for increased functional stability, while also expanding the capacity to develop novel "switch" pockets for responding to diverse pests and pathogens. Such evolutionary roles are congruent with the hypothesis that members of this unique resistance gene family provide fitness advantages in Gossypium.

Local synteny and codon usage contribute to asymmetric sequence divergence of Saccharomyces cerevisiae gene duplicates

Background

Duplicated genes frequently experience asymmetric rates of sequence evolution. Relaxed selective constraints and positive selection have both been invoked to explain the observation that one paralog within a gene-duplicate pair exhibits an accelerated rate of sequence evolution. In the majority of studies where asymmetric divergence has been established, there is no indication as to which gene copy, ancestral or derived, is evolving more rapidly. In this study we investigated the effect of local synteny (gene-neighborhood conservation) and codon usage on the sequence evolution of gene duplicates in the S. cerevisiae genome. We further distinguish the gene duplicates into those that originated from a whole-genome duplication (WGD) event (ohnologs) versus small-scale duplications (SSD) to determine if there exist any differences in their patterns of sequence evolution.

Results

For SSD pairs, the derived copy evolves faster than the ancestral copy. However, there is no relationship between rate asymmetry and synteny conservation (ancestral-like versus derived-like) in ohnologs. mRNA abundance and optimal codon usage as measured by the CAI is lower in the derived SSD copies relative to ancestral paralogs. Moreover, in the case of ohnologs, the faster-evolving copy has lower CAI and lowered expression.

Conclusions

Together, these results suggest that relaxation of selection for codon usage and gene expression contribute to rate asymmetry in the evolution of duplicated genes and that in SSD pairs, the relaxation of selection stems from the loss of ancestral regulatory information in the derived copy.

Structure, expression differentiation and evolution of duplicated fiber developmental genes in Gossypium barbadense and G. hirsutum

BACKGROUND

Both Gossypium hirsutum and G. barbadense probably originated from a common ancestor, but they have very different agronomic and fiber quality characters. Here we selected 17 fiber development-related genes to study their structures, tree topologies, chromosomal location and expression patterns to better understand the interspecific divergence of fiber development genes in the two cultivated tetraploid species.

RESULTS

The sequence and structure of 70.59% genes were conserved with the same exon length and numbers in different species, while 29.41% genes showed diversity. There were 15 genes showing independent evolution between the A- and D-subgenomes after polyploid formation, while two evolved via different degrees of colonization. Chromosomal location showed that 22 duplicate genes were located in which at least one fiber quality QTL was detected. The molecular evolutionary rates suggested that the D-subgenome of the allotetraploid underwent rapid evolutionary differentiation, and selection had acted at the tetraploid level. Expression profiles at fiber initiation and early elongation showed that the transcripts levels of most genes were higher in Hai7124 than in TM-1. During the primary-secondary transition period, expression of most genes peaked earlier in TM-1 than in Hai7124. Homeolog expression profile showed that A-subgenome, or the combination of A- and D-subgenomes, played critical roles in fiber quality divergence of G. hirsutum and G. barbadense. However, the expression of D-subgenome alone also played an important role.

CONCLUSION

Integrating analysis of the structure and expression to fiber development genes, suggests selective breeding for certain desirable fiber qualities played an important role in divergence of G. hirsutum and G. barbadense.

To be or not to be the odd one out--allele-specific transcription in pentaploid dogroses (Rosa L. sect. Caninae (DC.) Ser)

BACKGROUND

Multiple hybridization events gave rise to pentaploid dogroses which can reproduce sexually despite their uneven ploidy level by the unique canina meiosis. Two homologous chromosome sets are involved in bivalent formation and are transmitted by the haploid pollen grains and the tetraploid egg cells. In addition the egg cells contain three sets of univalent chromosomes which are excluded from recombination. In this study we investigated whether differential behavior of chromosomes as bivalents or univalents is reflected by sequence divergence or transcription intensity between homeologous alleles of two single copy genes (LEAFY, cGAPDH) and one ribosomal DNA locus (nrITS).

RESULTS

We detected a maximum number of four different alleles of all investigated loci in pentaploid dogroses and identified the respective allele with two copies, which is presumably located on bivalent forming chromosomes. For the alleles of the ribosomal DNA locus and cGAPDH only slight, if any, differential transcription was determined, whereas the LEAFY alleles with one copy were found to be significantly stronger expressed than the LEAFY allele with two copies. Moreover, we found for the three marker genes that all alleles have been under similar regimes of purifying selection.

CONCLUSIONS

Analyses of both molecular sequence evolution and expression patterns did not support the hypothesis that unique alleles probably located on non-recombining chromosomes are less functional than duplicate alleles presumably located on recombining chromosomes.

The Rad50 genes of diploid and polyploid wheat species. Analysis of homologue and homoeologue expression and interactions with Mre11

The MRN complex plays a central role in the DNA repair pathways of eukaryotic cells and takes part in many other processes, including cell cycle checkpoint signalling, meiosis, DNA replication and telomere maintenance. This complex is formed by the interaction of the products of the Mre11, Rad50 and Nbs1 genes. This paper reports the molecular characterization, expression and interactions of the Rad50 gene in several wheat species with different levels of ploidy. The homoeologous Rad50 wheat genes were found to show a high level of conservation. Most of the RAD50 domains and motifs previously described in other species were also present in wheat RAD50; these proteins are therefore likely to have similar functions. Interactions between the RAD50 wheat proteins and their MRE11 counterparts in the MRN complex were observed. The level of expression of Rad50 in each of the species examined was determined and compared with those previously reported for the Mre11 genes. In some cases similar levels of expression were seen, as expected. The expression of the RAD50 homoeologous genes was assessed in two polyploid wheat species using quantitative PCR. In both cases, an overexpression of the Rad50B gene was detected. Although the results indicate the maintenance of function of these species' three homoeologous Rad50 genes, the biased expression of Rad50B might indicate ongoing silencing of one or both other homoeologues in polyploid wheat. To assess the consequences of such silencing on the formation of the MRN complex, the interactions between individual homoeologues of Rad50 and their genomic counterpart Mre11 genes were examined. The results indicate the inexistence of genomic specificity in the interactions between these genes. This would guarantee the formation of an MRN complex in wheat.

Spatio-temporal patterns of genome evolution in allotetraploid species of the genus Oryza

Despite knowledge that polyploidy is widespread and a major evolutionary force in flowering plant diversification, detailed comparative molecular studies on polyploidy have been confined to only a few species and families. The genus Oryza is composed of 23 species that are classified into ten distinct 'genome types' (six diploid and four polyploid), and is emerging as a powerful new model system to study polyploidy. Here we report the identification, sequence and comprehensive comparative annotation of eight homoeologous genomes from a single orthologous region (Adh1-Adh2) from four allopolyploid species representing each of the known Oryza genome types (BC, CD, HJ and KL). Detailed comparative phylogenomic analyses of these regions within and across species and ploidy levels provided several insights into the spatio-temporal dynamics of genome organization and evolution of this region in 'natural' polyploids of Oryza. The major findings of this study are that: (i) homoeologous genomic regions within the same nucleus experience both independent and parallel evolution, (ii) differential lineage-specific selection pressures do not occur between polyploids and their diploid progenitors, (iii) there have been no dramatic structural changes relative to the diploid ancestors, (iv) a variation in the molecular evolutionary rate exists between the two genomes in the BC complex species even though the BC and CD polyploid species appear to have arisen <2 million years ago, and (v) there are no clear distinctions in the patterns of genome evolution in the diploid versus polyploid species.

Identification of the family of aquaporin genes and their expression in upland cotton (Gossypium hirsutum L.)

BACKGROUND

Cotton (Gossypium spp.) is produced in over 30 countries and represents the most important natural fiber in the world. One of the primary factors affecting both the quantity and quality of cotton production is water. A major facilitator of water movement through cell membranes of cotton and other plants are the aquaporin proteins. Aquaporin proteins are present as diverse forms in plants, where they function as transport systems for water and other small molecules. The plant aquaporins belong to the large major intrinsic protein (MIP) family. In higher plants, they consist of five subfamilies including plasma membrane intrinsic proteins (PIP), tonoplast intrinsic proteins (TIP), NOD26-like intrinsic proteins (NIP), small basic intrinsic proteins (SIP), and the recently discovered X intrinsic proteins (XIP). Although a great deal is known about aquaporins in plants, very little is known in cotton.

RESULTS

From a molecular cloning effort, together with a bioinformatic homology search, 71 upland cotton (G. hirsutum) aquaporin genes were identified. The cotton aquaporins consist of 28 PIP and 23 TIP members with high sequence similarity. We also identified 12 NIP and 7 SIP members that showed more divergence. In addition, one XIP member was identified that formed a distinct 5th subfamily. To explore the physiological roles of these aquaporin genes in cotton, expression analyses were performed for a select set of aquaporin genes from each subfamily using semi-quantitative reverse transcription (RT)-PCR. Our results suggest that many cotton aquaporin genes have high sequence similarity and diverse roles as evidenced by analysis of sequences and their expression.

CONCLUSION

This study presents a comprehensive identification of 71 cotton aquaporin genes. Phylogenetic analysis of amino acid sequences divided the large and highly similar multi-gene family into the known 5 aquaporin subfamilies. Together with expression and bioinformatic analyses, our results support the idea that the genes identified in this study represent an important genetic resource providing potential targets to modify the water use properties of cotton.

Meta-analysis of cotton fiber quality QTLs across diverse environments in a Gossypium hirsutum x G. barbadense RIL population

BACKGROUND

Cotton fibers (produced by Gossypium species) are the premier natural fibers for textile production. The two tetraploid species, G. barbadense (Gb) and G. hirsutum (Gh), differ significantly in their fiber properties, the former having much longer, finer and stronger fibers that are highly prized. A better understanding of the genetics and underlying biological causes of these differences will aid further improvement of cotton quality through breeding and biotechnology. We evaluated an inter-specific Gh x Gb recombinant inbred line (RIL) population for fiber characteristics in 11 independent experiments under field and glasshouse conditions. Sites were located on 4 continents and 5 countries and some locations were analyzed over multiple years.

RESULTS

The RIL population displayed a large variability for all major fiber traits. QTL analyses were performed on a per-site basis by composite interval mapping. Among the 651 putative QTLs (LOD > 2), 167 had a LOD exceeding permutation based thresholds. Coincidence in QTL location across data sets was assessed for the fiber trait categories strength, elongation, length, length uniformity, fineness/maturity, and color. A meta-analysis of more than a thousand putative QTLs was conducted with MetaQTL software to integrate QTL data from the RIL and 3 backcross populations (from the same parents) and to compare them with the literature. Although the global level of congruence across experiments and populations was generally moderate, the QTL clustering was possible for 30 trait x chromosome combinations (5 traits in 19 different chromosomes) where an effective co-localization of unidirectional (similar sign of additivity) QTLs from at least 5 different data sets was observed. Most consistent meta-clusters were identified for fiber color on chromosomes c6, c8 and c25, fineness on c15, and fiber length on c3.

CONCLUSIONS

Meta-analysis provided a reliable means of integrating phenotypic and genetic mapping data across multiple populations and environments for complex fiber traits. The consistent chromosomal regions contributing to fiber quality traits constitute good candidates for the further dissection of the genetic and genomic factors underlying important fiber characteristics, and for marker-assisted selection.

Duplication, divergence and persistence in the Phytochrome photoreceptor gene family of cottons (Gossypium spp.)

BACKGROUND

Phytochromes are a family of red/far-red photoreceptors that regulate a number of important developmental traits in cotton (Gossypium spp.), including plant architecture, fiber development, and photoperiodic flowering. Little is known about the composition and evolution of the phytochrome gene family in diploid (G. herbaceum, G. raimondii) or allotetraploid (G. hirsutum, G. barbadense) cotton species. The objective of this study was to obtain a preliminary inventory and molecular-evolutionary characterization of the phytochrome gene family in cotton.

RESULTS

We used comparative sequence resources to design low-degeneracy PCR primers that amplify genomic sequence tags (GSTs) for members of the PHYA, PHYB/D, PHYC and PHYE gene sub-families from A- and D-genome diploid and AD-genome allotetraploid Gossypium species. We identified two paralogous PHYA genes (designated PHYA1 and PHYA2) in diploid cottons, the result of a Malvaceae-specific PHYA gene duplication that occurred approximately 14 million years ago (MYA), before the divergence of the A- and D-genome ancestors. We identified a single gene copy of PHYB, PHYC, and PHYE in diploid cottons. The allotetraploid genomes have largely retained the complete gene complements inherited from both of the diploid genome ancestors, with at least four PHYA genes and two genes encoding PHYB, PHYC and PHYE in the AD-genomes. We did not identify a PHYD gene in any cotton genomes examined.

CONCLUSIONS

Detailed sequence analysis suggests that phytochrome genes retained after duplication by segmental duplication and allopolyploidy appear to be evolving independently under a birth-and-death-process with strong purifying selection. Our study provides a preliminary phytochrome gene inventory that is necessary and sufficient for further characterization of the biological functions of each of the cotton phytochrome genes, and for the development of 'candidate gene' markers that are potentially useful for cotton improvement via modern marker-assisted selection strategies.

Nuclear glutamine synthetase evolution in Nicotiana: phylogenetics and the origins of allotetraploid and homoploid (diploid) hybrids

Interspecies relationships in Nicotiana (Solanaceae) are complex because 40 species are diploid (two sets of chromosomes) and 35 species are allotetraploid (four sets of chromosomes, two from each progenitor diploid species). We sequenced a fragment (containing four introns) of the nuclear gene 'chloroplast-expressed glutamine synthetase' (ncpGS) in 65 species of Nicotiana. Here we present the first phylogenetic analysis based on a low-copy nuclear gene for this well studied and important genus. Diploid species have a single-copy of ncpGS, and allotetraploids as expected have two homeologous copies, each derived from their progenitor diploid. Results were particularly useful for determining the paternal lineage of previously enigmatic taxa (for which our previous analyses had revealed only the maternal progenitors). In particular, we were able to shed light on the origins of the two oldest and largest allotetraploid sections, N. sects. Suaveolentes and Repandae. All homeologues have an intact reading frame and apparently similar rates of divergence, suggesting both remain functional. Difficulties in fitting certain diploid species into the sectional classification of Nicotiana on morphological grounds, coupled with discordance between the ncpGS data and previous trees (i.e. plastid, nuclear ribosomal DNA), indicate a number of homoploid (diploid) hybrids in the genus. We have evidence for Nicotiana glutinosa and Nicotiana linearis being of hybrid origin and patterns of intra-allelic recombination also indicate the possibility of reticulate origins for other diploid species.

Homoeologous nonreciprocal recombination in polyploid cotton

Polyploid formation and processes that create partial genomic duplication generate redundant genomic information, whose fate is of particular interest to evolutionary biologists. Different processes can lead to diversification among duplicate genes, which may be counterbalanced by mechanisms that retard divergence, including gene conversion via nonreciprocal homoeologous exchange. Here, we used genomic resources in diploid and allopolyploid cotton (Gossypium) to detect homoeologous single nucleotide polymorphisms provided by expressed sequence tags from G. arboreum (A genome), G. raimondii (D genome) and G. hirsutum (AD genome), allowing us to identify homoeo-single nucleotide polymorphism patterns indicative of potential homoeologous exchanges. We estimated the proportion of contigs in G. hirsutum that have experienced nonreciprocal homoeologous exchanges since the origin of polyploid cotton 1-2 million years ago (Mya) to be between 1.8% and 1.9%. To address the question of when the intergenomic exchange occurred, we assayed six of the genes affected by homoeo-recombination in all five Gossypium allopolyploids using a phylogenetic approach. This analysis revealed that nonreciprocal homoeologous exchanges have occurred throughout polyploid divergence and speciation, as opposed to saltationally with polyploid formation. In addition, some genomic regions show multiple patterns of homoeologous recombination among species.

Structure and size variations between 12A and 12D homoeologous chromosomes based on high-resolution cytogenetic map in allotetraploid cotton

Cotton is a model system for studying polyploidization, genomic organization, and genome-size variation because the allotetraploid was formed 1-2 million years ago, which is old enough for sequence divergence but relatively recent to maintain genome stability. In spite of characterizing random genomic sequences in many polyploidy plants, the cytogenetic and sequence data that decipher homoeologous chromosomes are very limited in allopolyploid species. Here, we reported comprehensive analyses of integrated cytogenetic and linkage maps of homoeologous chromosomes 12A and 12D in allotetraploid cotton using fluorescence in situ hybridization and a large number of bacterial artificial chromosomes that were anchored by simple sequence repeat markers in the corresponding linkage maps. Integration of genetic loci into physical localizations showed considerable variation of genome organization, structure, and size between 12A and 12D homoeologous chromosomes. The distal regions of the chromosomes displayed relatively lower levels of structural and size variation than other regions of the chromosomes. The highest level of variation was found in the pericentric regions in the long arms of the two homoeologous chromosomes. The genome-size difference between A and D sub-genomes in allotetraploid cotton was mainly associated with uneven expansion or contraction between different regions of homoeologous chromosomes. As an attempt for studying on the polyploidy homoeologous chromosomes, these results are of general interest to the understanding and future sequencing of complex genomes in plant species.