Homeologous Pairs Research Articles

Subfunctionalization of Parental Polyamine Oxidase (PAO) Genes in the Allopolyploid Tobacco Nicotiana tabacum (L.).

Polyamines play an important role in developmental and environmental stress responses in plants. Polyamine oxidases (PAOs) are flavin-adenine-dinucleotide-dependent enzymes associated with polyamine catabolism. In this study, 14 genes were identified in the tobacco genome that code for PAO proteins being named based on their sequence homology with Arabidopsis PAOs (AtPAO1-5): NtPAO1A-B; NtPAO2A-C, NtPAO4A-D, and NtPAO5A-E. Sequence analysis confirmed that the PAO gene family of the allopolyploid hybrid Nicotiana tabacum is not an exact combination of the PAO genes of the maternal Nicotiana sylvestris and paternal Nicotiana tomentosiformis ones. The loss of the N. sylvestris homeolog of NtPAO5E and the gain of an extra NtPAO2 copy, likely of Nicotiana othophora origin, was revealed. The latter adds to the few pieces of evidence suggesting that the paternal parent of N. tabacum was an introgressed hybrid of N. tomentosiformis and N. othophora. Gene expression analysis indicated that all 14 PAO genes kept their expression following the formation of the hybrid species. The homeologous gene pairs showed similar or opposite regulation depending on the investigated organ, applied stress, or hormone treatment. The data indicate that the expression pattern of the homeologous genes is diversifying in a process of subfunctionalization.

A near-complete genome assembly of the allotetrapolyploid Cenchrus fungigraminus (JUJUNCAO) provides insights into its evolution and C4 photosynthesis

JUJUNCAO (Cenchrus fungigraminus; 2n = 4x = 28) is a Cenchrus grass with the highest biomass production among cultivated plants, and it can be used for mushroom cultivation, animal feed, and biofuel production. Here, we report a nearly complete genome assembly of JUJUNCAO and reveal that JUJUNCAO is an allopolyploid that originated ∼2.7 million years ago (mya). Its genome consists of two subgenomes, and subgenome A shares high collinear synteny with pearl millet. We also investigated the genome evolution of JUJUNCAO and suggest that the ancestral karyotype of Cenchrus split into the A and B ancestral karyotypes of JUJUNCAO. Comparative transcriptome and DNA methylome analyses revealed functional divergence of homeologous gene pairs between the two subgenomes, which was a further indication of asymmetric DNA methylation. The three types of centromeric repeat in the JUJUNCAO genome (CEN137, CEN148, and CEN156) may have evolved independently within each subgenome, with some introgressions of CEN156 from the B to the A subgenome. We investigated the photosynthetic characteristics of JUJUNCAO, revealing its typical C4 Kranz anatomy and high photosynthetic efficiency. NADP-ME and PEPCK appear to cooperate in the major C4 decarboxylation reaction of JUJUNCAO, which is different from other C4 photosynthetic subtypes and may contribute to its high photosynthetic efficiency and biomass yield. Taken together, our results provide insights into the highly efficient photosynthetic mechanism of JUJUNCAO and provide a valuable reference genome for future genetic and evolutionary studies, as well as genetic improvement of Cenchrus grasses.

Conservation and divergence of expression of GA2-oxidase homeologs in apple (Malus x domestica Borkh.).

In domesticated apple (Malus x domestica Borkh.) and other woody perennials, floral initiation can be repressed by gibberellins (GAs). The associated mechanism is a major unanswered question in plant physiology, and understanding organismal aspects of GA signaling in apple has important commercial applications. In plants, the major mechanism for elimination of GAs and resetting of GA signaling is through catabolism by GA2-oxidases (GA2ox). We found that the GA2ox gene family in apple comprises 16 genes representing eight, clearly defined homeologous pairs, which were named as MdGA2ox1A/1B to MdGA2ox8A/8B. Expression of the genes was analyzed in the various structures of the spur, where flowers are initiated, as well as in various structures of seedlings over one diurnal cycle and in response to water-deficit and salt stress. Among the results, we found that MdGA2ox2A/2B dominated expression in the shoot apex and were strongly upregulated in the apex after treatment with exogenous GA3, suggesting potential involvement in repression of flowering. Several MdGA2ox genes also showed preferential expression in the leaf petiole, fruit pedicel, and the seed coat of developing seeds, potentially representing mechanisms to limit diffusion of GAs across these structures. In all contexts studied, we documented both concerted and distinct expression of individual homeologs. This work introduces an accessible woody plant model for studies of GA signaling, GA2ox gene regulation, and conservation/divergence of expression of homeologous genes, and should find application in development of new cultivars of apple and other tree fruits.

Dioecy and chromosomal sex determination are maintained through allopolyploid speciation in the plant genus Mercurialis.

Polyploidization may precipitate dramatic changes to the genome, including chromosome rearrangements, gene loss, and changes in gene expression. In dioecious plants, the sex-determining mechanism may also be disrupted by polyploidization, with the potential evolution of hermaphroditism. However, while dioecy appears to have persisted through a ploidy transition in some species, it is unknown whether the newly formed polyploid maintained its sex-determining system uninterrupted, or whether dioecy re-evolved after a period of hermaphroditism. Here, we develop a bioinformatic pipeline using RNA-sequencing data from natural populations to demonstrate that the allopolyploid plant Mercurialis canariensis directly inherited its sex-determining region from one of its diploid progenitor species, M. annua, and likely remained dioecious through the transition. The sex-determining region of M. canariensis is smaller than that of its diploid progenitor, suggesting that the non-recombining region of M. annua expanded subsequent to the polyploid origin of M. canariensis. Homeologous pairs show partial sexual subfunctionalization. We discuss the possibility that gene duplicates created by polyploidization might contribute to resolving sexual antagonism.

A GWAS identified a major QTL for resistance to Fusarium wilt (Fusarium oxysporum f. sp. vasinfectum) race 4 in a MAGIC population of Upland cotton and a meta-analysis of QTLs for Fusarium wilt resistance.

A major QTL conferring resistance to Fusarium wilt race 4 in a narrow region of chromosome D02 was identified in a MAGIC population of 550 RILs of Upland cotton. Numerous studies have been conducted to investigate the genetic basis of Fusarium wilt (FW, caused by Fusarium oxysporum f. sp. vasinfectum, FOV) resistance using bi-parental and association mapping populations in cotton. In this study, a multi-parent advanced generation inter-cross (MAGIC) population of 550 recombinant inbred lines (RILs), together with their 11 Upland cotton (Gossypium hirsutum) parents, was used to identify QTLs for FOV race 4 (FOV4) resistance. Among the parents, Acala Ultima, M-240 RNR, and Stoneville 474 were the most resistant, while Deltapine Acala 90, Coker 315, and Stoneville 825 were the most susceptible. Twenty-two MAGIC lines were consistently resistant to FOV4. Through a genome-wide association study (GWAS) based on 473,516 polymorphic SNPs, a major FOV4 resistance QTL within a narrow region on chromosomes D02 was detected, allowing identification of 14 candidate genes. Additionally, a meta-analysis of 133 published FW resistance QTLs showed a D subgenome and individual chromosome bias and no correlation between homeologous chromosome pairs. This study represents the first GWAS study using a largest genetic population and the most comprehensive meta-analysis for FW resistance in cotton. The results illustrated that 550 lines were not enough for high resolution mapping to pinpoint a candidate gene, and experimental errors in phenotyping cotton for FW resistance further compromised the accuracy and precision in QTL localization and identification of candidate genes. This study identified FOV4-resistant parents and MAGIC lines, and the first major QTL for FOV4 resistancein Upland cotton, providing useful information for developing FOV4-resistant cultivars and further genomic studies towards identification of causal genes for FOV4 resistance in cotton.

Comparative Cytogenetics Analysis Among Peckoltia Species (Siluriformes, Loricariidae): Insights on Karyotype Evolution and Biogeography in the Amazon Region.

Peckoltia is widely distributed genus in the Amazon and Orinoco basins and the Guiana Shield, containing 18 valid species, and distinct morphotypes still needing description in the scientific literature due to its great taxonomic complexity. This study performed a comparative chromosomal analysis of two undescribed Peckoltia species (Peckoltia sp. 3 Jarumã and Peckoltia sp. 4 Caripetuba) from the Brazilian Amazon using conventional chromosome bands methods and in situ localization of the repetitive DNA (5S and 18S rRNA and U1 snRNA genes and telomeric sequences). Both species presented 2n = 52 but differed in their karyotype formula, probably due to inversions or translocations. The nucleolus organizer regions (NORs) showed distal location on a probably homeologous submetacentric pair in both species, besides an extra signal in a subtelocentric chromosome in Peckoltia sp. 4 Caripetuba. Heterochromatin occurred in large blocks, with different distributions in the species. The mapping of the 18S and 5S rDNA, and U1 snDNA showed differences in locations and number of sites. No interstitial telomeric sites were detected using the (TTAGGG)n probes. Despite 2n conservationism in Peckoltia species, the results showed variation in karyotype formulas, chromosomal bands, and locations of repetitive sites, demonstrating great chromosomal diversity. A proposal for Peckoltia karyotype evolution was inferred in this study based on the diversity of location and number of chromosomal markers analyzed. A comparative analysis with other Peckoltia karyotypes described in the literature, their biogeography patterns, and molecular phylogeny led to the hypothesis that the derived karyotype was raised in the left bank of the Amazon River.

Introgression and genetic mapping of leaf rust and stripe rust resistance in Aegilops triuncialis

The growing and cultivating resistant wheat crop varieties is important to meet the demands of the growing population and minimizing the yield losses due to foliar diseases. More important is the identification of novel resistance sources and transfer of resistance in ready to use form. In the current study, leaf rust (LR) and stripe rust (YR) resistant tetraploid nonprogenitors of wheat Aegilops triuncialis (UtUtCtCt) acc pau 3462 was crossed and backcrossed susceptible cultivar WL711(NN) by inducing homeologous pairing using CS ph1. Recurrent parent type plants were selected in subsequent generation with resistance to LR and YR and BC2F7 introgression line (2n=42) named ILtri have been developed. To understand the nature and inheritance of LR and YR resistance genes and to map their genomic location, F2 and F2:3 mapping populations were developed by crossing ILtri with WL711(NN). In F2 and F2:3, the seedlings and adult plants segregated into 3R:1S and 1HR:2Seg:1HS ratios, respectively for both LR and YR, indicating inheritance of single dominant all stage resistance gene working against both the rusts. These genes were temporary designated as Lrtri and Yrtri and were inherited independently.Molecular mapping of 614 SSR markers mapped the Lrtri at a distance of 11.2 cM from SSR marker Xwmc606.

Asymmetric epigenome maps of subgenomes reveal imbalanced transcription and distinct evolutionary trends in Brassica napus

The complexity of the epigenome landscape and transcriptional regulation is significantly increased during plant polyploidization, which drives genome evolution and contributes to the increased adaptability to diverse environments. However, a comprehensive epigenome map of Brassica napus is still unavailable. In this study, we performed integrative analysis of five histone modifications, RNA polymerase II occupancy, DNA methylation, and transcriptomes in two B. napus lines (2063A and B409), and established global maps of regulatory elements, chromatin states, and their dynamics for the whole genome (including the An and Cn subgenomes) in four tissue types (young leaf, flower bud, silique, and root) of these two lines. Approximately 65.8% of the genome was annotated with different epigenomic signals. Compared with the Cn subgenome, the An subgenome possesses a higher level of active epigenetic marks and lower level of repressive epigenetic marks. Genes from subgenome-unique regions contribute to the major differences between the An and Cn subgenomes. Asymmetric histone modifications between homeologous gene pairs reflect their biased expression patterns. We identified a novel bivalent chromatin state (with H3K4me1 and H3K27me3) in B. napus that is associated with tissue-specific gene expression. Furthermore, we observed that different types of duplicated genes have discrepant patterns of histone modification and DNA methylation levels. Collectively, our findings provide a valuable epigenetic resource for allopolyploid plants.

Evolution of Homeologous Gene Expression in Polyploid Wheat.

Polyploidization has played a prominent role in the evolutionary history of plants. Two recent and sequential allopolyploidization events have resulted in the formation of wheat species with different ploidies, and which provide a model to study the effects of polyploidization on the evolution of gene expression. In this study, we identified differentially expressed genes (DEGs) between four BBAA tetraploid wheats of three different ploidy backgrounds. DEGs were found to be unevenly distributed among functional categories and duplication modes. We observed more DEGs in the extracted tetraploid wheat (ETW) than in natural tetraploid wheats (TD and TTR13) as compared to a synthetic tetraploid (AT2). Furthermore, DEGs showed higher Ka/Ks ratios than those that did not show expression changes (non-DEGs) between genotypes, indicating DEGs and non-DEGs experienced different selection pressures. For A-B homeolog pairs with DEGs, most of them had only one differentially expressed copy, however, when both copies of a homeolog pair were DEGs, the A and B copies were more likely to be regulated to the same direction. Our results suggest that both cis- and inter-subgenome trans-regulatory changes are important drivers in the evolution of homeologous gene expression in polyploid wheat, with ploidy playing a significant role in the process.

Targeted inactivation of the AGO1 homeologues of Nicotiana benthamiana reveals their distinct roles in development and antiviral defence.

The Solanaceae family includes numerous highly valuable crops. Understanding the viral diseases that affect them is of great importance. Nicotiana benthamiana has contributed greatly to unravelling antiviral RNA interference, and can also be regarded as an adequate model for studying viral diseases of solanaceous crops. This species, however, as with many of its relatives, possesses an allopolyploid genome, in which homeologous gene pairs frequently occur. AGO1 is a pivotal component of most plant RNA silencing pathways. The Nicotiana benthamiana genome encodes two highly similar AGO1 homeologues: AGO1A and AGO1B. To understand their roles in planta, their genes were selectively inactivated. Given the inherent limitations of RNA interference-based techniques, we used genome editing to achieve this goal. We found that AGO1A was not required for normal development, while AGO1B was indispensable for that. By contrast, the two homeologues both contributed to antiviral defence. Additionally, we observed that AGO1B utilised miR168 poorly, which may help to retain a significant level of antiviral RNA interference during viral infection. Our results have important implications for the better understanding of viral diseases of economically important solanaceous crops.

Comparative Cytogenetics of Four Sea Turtle Species (Cheloniidae): G-Banding Pattern and in situ Localization of Repetitive DNA Units

Sea turtles are considered flagship species for marine biodiversity conservation and are considered to be at varying risk of extinction globally. Cases of hybridism have been reported in sea turtles, but chromosomal analyses are limited to classical karyotype descriptions and a few molecular cytogenetic studies. In order to compare karyotypes and understand evolutive mechanisms related to chromosome dif­ferentiation in this group, Chelonia mydas, Caretta caretta, Eretmochelys imbricata, and Lepidochelys olivacea were cytogenetically characterized in the present study. When the obtained cytogenetic data were compared with the putative ancestral Cryptodira karyotype, the studied species showed the same diploid number (2n) of 56 chromosomes, with some variations in chromosomal morphology (karyotypic formula) and minor changes in longitudinal band locations. In situ localization using a 18S ribosomal DNA probe indicated a homeologous microchromosome pair bearing a 45S ribosomal DNA locus and size heteromorphism in all 4 species. Interstitial telomeric sites were identified in a microchromosome pair in C. mydas and C. caretta. The data showed that interspecific variations occurred in chromosomal sets among the Cheloniidae species, in addition to other Cryptodira karyotypes. These variations generated lineage-specific karyotypic diversification in sea turtles, which will have considerable implications for hybrid recognition and for the study, the biology, ecology, and evolutionary history of regional and global populations. Furthermore, we demonstrated that some chromosome rearrangements occurred in sea turtle species, which is in conflict with the hypothesis of conserved karyotypes in this group.

A Recently Formed Triploid Cardamine insueta Inherits Leaf Vivipary and Submergence Tolerance Traits of Parents.

Contemporary speciation provides a unique opportunity to directly observe the traits and environmental responses of a new species. Cardamine insueta is an allotriploid species that appeared within the past 150 years in a Swiss village, Urnerboden. In contrast to its two progenitor species, Cardamine amara and Cardamine rivularis that live in wet and open habitats, respectively, C. insueta is found in-between their habitats with temporal water level fluctuation. This triploid species propagates clonally and serves as a triploid bridge to form higher ploidy species. Although niche separation is observed in field studies, the mechanisms underlying the environmental robustness of C. insueta are not clear. To characterize responses to a fluctuating environment, we performed a time-course analysis of homeolog gene expression in C. insueta in response to submergence treatment. For this purpose, the two parental (C. amara and C. rivularis) genome sequences were assembled with a reference-guided approach, and homeolog-specific gene expression was quantified using HomeoRoq software. We found that C. insueta and C. rivularis initiated vegetative propagation by forming ectopic meristems on leaves, while C. amara did not. We examined homeolog-specific gene expression of three species at nine time points during the treatment. The genome-wide expression ratio of homeolog pairs was 2:1 over the time-course, consistent with the ploidy number. By searching the genes with high coefficient of variation of expression over time-course transcriptome data, we found many known key transcriptional factors related to meristem development and formation upregulated in both C. rivularis and rivularis-homeolog of C. insueta, but not in C. amara. Moreover, some amara-homeologs of these genes were also upregulated in the triploid, suggesting trans-regulation. In turn, Gene Ontology analysis suggested that the expression pattern of submergence tolerant genes in the triploid was inherited from C. amara. These results suggest that the triploid C. insueta combined advantageous patterns of parental transcriptomes to contribute to its establishment in a new niche along a water-usage gradient.

Sequence-based genetic mapping of Cynodon dactylon Pers. reveals new insights into genome evolution in Poaceae

Bermudagrass (Cynodon dactylon Pers.) is an important warm-season perennial used extensively for turf, forage, soil conservation and remediation worldwide. However, limited genomic information has hindered the application of molecular tools towards understanding genome evolution and in breeding new cultivars. We genotype a first-generation selfed population derived from the tetraploid (4x = 36) ‘A12359’ using genotyping-by-sequencing. A high-density genetic map of 18 linkage groups (LGs) is constructed with 3,544 markers. Comparative genomic analyses reveal that each of nine homeologous LG pairs of C. dactylon corresponds to one of the first nine chromosomes of Oropetium thomaeum. Two nested paleo-ancestor chromosome fusions (ρ6-ρ9-ρ6, ρ2-ρ10-ρ2) may have resulted in a 12-to-10 chromosome reduction. A segmental dissemination of the paleo-chromosome ρ12 (ρ1-ρ12-ρ1, ρ6-ρ12-ρ6) leads to the 10-to-9 chromosome reduction in C. dactylon genome. The genetic map will assist in an ongoing whole genome sequence assembly and facilitate marker-assisted selection (MAS) in developing new cultivars.

First report for availability of HRT-like genes in Eragrostis tef and in silico analysis for elucidating their potential functions

Tef [Eragrostis tef (Zucc.) Troetter] is called ‘super grain’ because of its healthy and nutritious seeds along with compatibility to suboptimal environment. Plant Transcription factors (TFs) are key regulators of gene expression which affect transcription through interaction with proteins and DNA binding. Hordeum Repressor Transcription (HRT) are one of TF gene families with too scanty of information on their abundance, distribution, motifs, phylogenetic relationship among plants, promoter and association network analysis on cereals including tef. CoGeBlast, MEME, MEGA X and PLANTCARE were employed accordingly. Two HRT-like TFs (ET) were identified residing on chromosome 5A and 5B, the proteins exhibit CH residue at the C- terminal end and GIY-YIG domain in the N-terminal. The phylogenetic tree obtained depicted orthology in HRT genes in S. bicolor and Z. mays; O. sativa and O. glaberrima, S. italica and S. virdis; and H. vulgare and T. aestivum and homeologous pairs in E. tef. From the promoter analysis of EtHRTs, CRE elements like ABRE, MYBs, G-box, GC-motif, TC-rich repeats, O2-site were identified known to involve in the control of hormones, sugars and abiotic stresses (ABA, GA, zein and α-amylase) confirming the involvement/role of HRT-like TFs in growth, development and stress response of tef. The study contributes for the understanding of structure and function of HRT-like TFs among plants.

Introgression of powdery mildew resistance from Aegilops triuncialis into wheat through induced homeologous pairing

Powdery mildew is a serious fungal disease of wheat caused by Blumeria graminis f. sp. tritici. Chromosome 5U of Aegilops triuncialis carrying powdery mildew resistance was induced to recombine with chromosome 5A of wheat using ph1b deletion line. The resulting recombinant progeny lines were tested for resistance and screened with molecular markers. Among 367 progeny lines 8 were immune to powdery mildew, 53 were moderately resistant, 231 moderately susceptible and 75 highly susceptible. Screening of 5U/5A recombinant lines with 5A specific SSR markers indicated that a major powdery mildew resistance gene from A. triuncialis was introgressed on the long arm of wheat chromosome 5A, along with two other loci, one located in the terminal end of the same arm and the other on the short arm of 5A. Meiotic metaphase-I pairing indicated that chromosome pairing in the resistant recombinant lines was regular with 21 bivalents and high seed set. The DNA markers associated with the introgressions can be used in marker assisted selection for improved powdery mildew resistance in wheat.

Transfer of grain softness from 5U-5A wheat-Aegilops triuncialis substitution line to bread wheat through induced homeologous pairing

Bread wheat sustains genes for grain softness on “Ha” locus of short arm of 5D chromosome. Pina-D1 and Pinb-D1 alleles of “Ha” locus together are responsible for soft grain texture of wheat. The wheat cultivar WL711 had wild type Pinb-D1a and mutated Pina allele, imparting hard kernel texture. A 5U-5A substitution line of Aegilops triuncialis had extra soft grains. Efforts were made to transfer puroindolines from Aegilops triuncialis 5U as 5U-5A substitution line BTC17 to its homeologous chromosome 5A of bread wheat. A total of 367 5U-5A substitution recombinants developed through induced homeologous pairing by ph1b deletion of Ph1 locus were screened for the presence of Pina and Pinb. Only 23 lines were found to have both Pina and Pinb. Grain hardness index of recombinants with functional puroindolines was reduced from 6 to 67%. Deletion mapping of wheat 5A chromosome specific SSR markers revealed the maintenance of synteny in 5U and 5A homeologous chromosomes and puroindolines were found to be translocated at the telomeric end of 5A short arm. Translocation of puroindolines from 5U was confirmed by sequencing of Pina and Pinb genes. The Pina nucleotide sequence was found 100% similar to Aegilops markgrafii (CC). Puroindoline proteins extracted from the selected softer lines exhibited significant antimicrobial activity against gram positive and gram negative bacteria. Chromosome 5U puroindolines transferred to 5A may be pyramided with 5D of soft wheat to develop extra soft varieties for superior biscuit making quality.

Tissue-specific gene expression and protein abundance patterns are associated with fractionation bias in maize

BackgroundMaize experienced a whole-genome duplication event approximately 5 to 12 million years ago. Because this event occurred after speciation from sorghum, the pre-duplication subgenomes can be partially reconstructed by mapping syntenic regions to the sorghum chromosomes. During evolution, maize has had uneven gene loss between each ancient subgenome. Fractionation and divergence between these genomes continue today, constantly changing genetic make-up and phenotypes and influencing agronomic traits.ResultsHere we regenerate the subgenome reconstructions for the most recent maize reference genome assembly. Based on both expression and abundance data for homeologous gene pairs across multiple tissues, we observed functional divergence of genes across subgenomes. Although the genes in the larger maize subgenome are often expressing more highly than their homeologs in the smaller subgenome, we observed cases where homeolog expression dominance switches in different tissues. We demonstrate for the first time that protein abundances are higher in the larger subgenome, but they also show tissue-specific dominance, a pattern similar to RNA expression dominance. We also find that pollen expression is uniquely decoupled from protein abundance.ConclusionOur study shows that the larger subgenome has a greater range of functional assignments and that there is a relative lack of overlap between the subgenomes in terms of gene functions than would be suggested by similar patterns of gene expression and protein abundance. Our study also revealed that some reactions are catalyzed uniquely by the larger and smaller subgenomes. The tissue-specific, nonequivalent expression-level dominance pattern observed here implies a change in regulatory control which favors differentiated selective pressure on the retained duplicates leading to eventual change in gene functions.

The Chromosome-Level Genome Assembly of European Grayling Reveals Aspects of a Unique Genome Evolution Process Within Salmonids.

Salmonids represent an intriguing taxonomical group for investigating genome evolution in vertebrates due to their relatively recent last common whole genome duplication event, which occurred between 80 and 100 million years ago. Here, we report on the chromosome-level genome assembly of European grayling (Thymallus thymallus), which represents one of the earliest diverged salmonid subfamilies. To achieve this, we first generated relatively long genomic scaffolds by using a previously published draft genome assembly along with long-read sequencing data and a linkage map. We then merged those scaffolds by applying synteny evidence from the Atlantic salmon (Salmo salar) genome. Comparisons of the European grayling genome assembly to the genomes of Atlantic salmon and Northern pike (Esox lucius), the latter used as a nonduplicated outgroup, detailed aspects of the characteristic chromosome evolution process that has taken place in European grayling. While Atlantic salmon and other salmonid genomes are portrayed by the typical occurrence of numerous chromosomal fusions, European grayling chromosomes were confirmed to be fusion-free and were characterized by a relatively large proportion of paracentric and pericentric inversions. We further reported on transposable elements specific to either the European grayling or Atlantic salmon genome, on the male-specific sdY gene in the European grayling chromosome 11A, and on regions under residual tetrasomy in the homeologous European grayling chromosome pairs 9A-9B and 25A-25B. The same chromosome pairs have been observed under residual tetrasomy in Atlantic salmon and in other salmonids, suggesting that this feature has been conserved since the subfamily split.

Paracentric Inversions Differentiate the Conservative Karyotypes in Two Centropomus Species (Teleostei: Centropomidae)

Centropomus is the sole genus of the Centropomidae family (Teleostei), comprising 12 species widely distributed throughout the Western Atlantic and Eastern Pacific, with 6 of them occurring in the Western Atlantic in extensive sympatry. Their life history and phylogenetic relationships are well characterized; however, aspects of chromosomal evolution are still unknown. Here, cytogenetic analyses of 2 Centropomus species of great economic value (C. undecimalis and C. mexicanus) were performed using conventional (Giemsa, Ag-NOR, and fluorochrome staining, C- and replication banding) and molecular (chromosomal mapping of 18S and 5S rDNA, H2A-H2B and H3 hisDNA, and (TTAGGG)_n repeats) approaches. The karyotypes of both species were composed of 48 solely acrocentric chromosomes (2n = 48; FN = 48), but the single ribosomal site was located in varying positions in the long arms of the second largest chromosome pair. Replication bands were generally similar, although conspicuous differences were observed in some chromosome regions. In both species, the histone H3 genes were located on 3 apparently homeologous chromosome pairs, but the exact position of these clusters differed slightly. Interspecific hisDNA and rDNA site displacements can indicate the occurrence of multiple paracentric inversions during the evolutionary diversification of the Centropomus genomes. Although the karyotypes remained similar in both species, our data demonstrate an unsuspected microstructural reorganization between them, driven most likely by a series of paracentric inversions.

A likelihood ratio test for changes in homeolog expression bias

BackgroundGene duplications are a major source of raw material for evolution and a likely contributor to the diversity of life on earth. Duplicate genes (i.e., homeologs, in the case of a whole genome duplication) may retain their ancestral function, sub- or neofunctionalize, or be lost entirely. A primary way that duplicate genes evolve new functions is by altering their expression patterns. Comparing the expression patterns of duplicate genes gives clues as to whether any of these evolutionary processes have occurred.ResultsWe develop a likelihood ratio test for the analysis of the expression ratios of duplicate genes across two conditions (e.g., tissues). We demonstrate an application of this test by comparing homeolog expression patterns of 1448 homeologous gene pairs using RNA-seq data generated from leaves and petals of an allotetraploid monkeyflower (Mimulus luteus). We assess the sensitivity of this test to different levels of homeolog expression bias and compare the method to several alternatives.ConclusionsThe likelihood ratio test derived here is a direct, transparent, and easily implemented method for detecting changes in homeolog expression bias that outperforms alternative approaches. While our method was derived with homeolog analysis in mind, this method can be used to analyze changes in the ratio of expression levels between any two genes in any two conditions.