Dynamic gene copy number variation in collinear regions of grass genomes

A salient feature of genomes of higher organisms is the birth and death of gene copies. An example is the alpha prolamin genes, which encode seed storage proteins in grasses (Poaceae) and represent a medium-size gene family. To better understand the mechanism, extent, and pace of gene amplification, we compared prolamin gene copies in the genomes of two different tribes in the Panicoideae, the Paniceae and the Andropogoneae. We identified alpha prolamin (setarin) gene copies in the diploid foxtail millet (Paniceae) genome (490 Mb) and compared them with orthologous regions in diploid sorghum (730 Mb) and ancient allotetraploid maize (2,300 Mb) (Andropogoneae). Because sequenced genomes of other subfamilies of Poaceae like rice (389 Mb) (Ehrhartoideae) and Brachypodium (272 Mb) (Pooideae) do not have alpha prolamin genes, their collinear regions can serve as "empty" reference sites. A pattern emerged, where genes were copied and inserted into other chromosomal locations followed by additional tandem duplications (clusters). We observed both recent (species-specific) insertion events and older ones that are shared by these tribes. Many older copies were deleted by unequal crossing over of flanking sequences or damaged by truncations. However, some remain intact with active and inactive alleles. These results indicate that genomes reflect only a snapshot of the gene content of a species and are far less static than conventional genetics has suggested. Nucleotide substitution rates for active alpha prolamins genes were twice as high as for low copy number beta, gamma, and delta prolamin genes, suggesting that gene amplification accelerates the pace of divergence.     

Matrix attachment regions and structural colinearity in the genomes of two grass species.

In order to gain insights into the relationship between spatial organization of the genome and genome function we have initiated studies of the co-linear Sh2/A1- homologous regions of rice (30 kb) and sorghum (50 kb). We have identified the locations of matrix attachment regions (MARs) in these homologous chromosome segments, which could serve as anchors for individual structural units or loops. Despite the fact that the nucleotide sequences serving as MARs were not detectably conserved, the general organizational patterns of MARs relative to the neighboring genes were preserved. All identified genes were placed in individual loops that were of comparable size for homologous genes. Hence, gene composition, gene orientation, gene order and the placement of genes into structural units has been evolutionarily conserved in this region. Our analysis demonstrated that the occurrence of various 'MAR motifs' is not indicative of MAR location. However, most of the MARs discovered in the two genomic regions were found to co-localize with miniature inverted repeat transposable elements (MITEs), suggesting that MITEs preferentially insert near MARs and/or that they can serve as MARs.     

Colinearity and gene density in grass genomes

Grasses are the single most important plant family in agriculture. In the past years, comparative genetic mapping has revealed conserved gene order (colinearity) among many grass species. Recently, the first studies at gene level have demonstrated that microcolinearity of genes is less conserved: small scale rearrangements and deletions complicate the microcolinearity between closely related species, such as sorghum and maize, but also between rice and other crop plants. In spite of these problems, rice remains the model plant for grasses as there is limited useful colinearity between Arabidopsis and grasses. However, studies in rice have to be complemented by more intensive genetic work on grass species with large genomes (maize, Triticeae). Gene-rich chromosomal regions in species with large genomes, such as wheat, have a high gene density and are ideal targets for partial genome sequencing.     

Structural organization of the barley D-hordein locus in comparison with its orthologous regions of wheat genomes.

D hordein, a prolamin storage protein of barley endosperms, is highly homologous to the high molecular weight (HWM) glutenin subunits, which are the major determinants of bread-making quality in wheat flour. In hexaploid wheat (AABBDD), each genome contains two paralogous copies of HMW-glutenin genes that encode the x- and y-type HMW-glutenin subunits. Previously, we reported the sequence analysis of a 102-kb genomic region that contains the HMW-glutenin locus of the D genome from Aegilops tauschii, the donor of the D genome of hexaploid wheat. Here, we present the sequence analysis of a 120-kb D-hordein region of the barley genome, a more distantly related member of the Triticeae grass tribe. Comparative sequence analysis revealed that gene content and order are generally conserved. Genes included in both of these orthologous regions are arranged in the following order: a Xa21-like receptor kinase, an endosperm globulin, an HMW prolamin, and a serine (threonine) protein kinase. However, in the wheat D genome, a region containing both the globulin and HMW-glutenin gene was duplicated, indicating that this duplication event occurred after the separation of the wheat and barley genomes. The intergenic regions are divergent with regard to the sequence and structural organization. It was found that different types of retroelements are responsible for the intergenic structure divergence in the wheat and barley genomes. In the barley region, we identified 16 long terminal repeat (LTR) retrotransposons in three distinct nested clusters. These retroelements account for 63% of the contig sequence. In addition, barley D hordein was compared with wheat HMW glutenins in terms of cysteine residue conservation and repeat domain organization.     

PCR-based identification of wheat genomes

We describe a PCR system that distinguishes the A, B and D genomes in wheat DNA extracts. PCRs were directed at the ‘non-transcribed spacer’ regions of the rDNA loci. The spacers within the D genome locus have a 71-bp insertion that is absent from the corresponding A and B loci PCR product sizes therefore enable D- and D+ genomes to be distinguished. The A and B genomes can be differentiated by PCR with an internal primer which does not anneal to A genome sequences. This work is relevant to the ancient ecology of wheat, as it is often difficult to determine ploidy level from morphological examination of archaeobotanical remains.     

Evolution of the defensin-like gene family in grass genomes

Plant defensins are small, diverse, cysteine-rich peptides, belonging to a group of pathogenesis-related defense mechanism proteins, which can provide a barrier against a broad range of pathogens. In this study, 51 defensin-like (DEFL) genes in Gramineae, including brachypodium, rice, maize and sorghum were identified based on bioinformatics methods. Using the synteny analysis method, we found that 21 DEFL genes formed 30 pairs of duplicated blocks that have undergone large-scale duplication events, mostly occurring between species. In particular, some chromosomal regions are highly conserved in the four grasses. Using mean K_s values, we estimated the approximate time of divergence for each pair of duplicated regions and found that these regions generally diverged more than 40 million years ago (Mya). Selection pressure analysis showed that the DEFL gene family is subjected to purifying selection. However, sliding window analysis detected partial regions of duplicated genes under positive selection. The evolutionary patterns within DEFL gene families among grasses can be used to explore the subsequent functional divergence of duplicated genes and to further analyse the antimicrobial effects of defensins during plant development.     

Origin of the duplicated regions in the yeast genomes

The genome of Saccharomyces cerevisiae contains several duplicated regions. The recent sequencing results of several yeast species suggest that the duplicated regions found in the modern Saccharomyces species are probably the result of a single gross duplication, as well as a series of sporadic independent short-segment duplications. The gross duplication might coincide with the origin of the ability to grow under anaerobic conditions.     

Phylogeny of the a genomes of wild and cultivated wheat species

Diploid species of the genus Triticum L. are its most ancient representatives and have the A genome, which was more recently inherited by all polyploid species. Studies of the phylogenetic relationships among diploid and polyploid wheat species help to identify the donors of elementary genomes and to examine the species specificity of genomes. In this study, molecular analysis of the variable sequences of three nuclear genes (Acc-1, Pgk-1, and Vrn-1) was performed for wild and cultivated wheat species, including both diploids and polyploids. Based on the sequence variations found in the genes, clear differences were observed among elementary genomes, but almost no polymorphism was detected within each genome in polyploids. At the same time, the regions of the three genes proved to be rather heterogeneous in the diploid species Triticum boeoticum Boiss., T. urartu Thum. ex Gandil., and T. monococcum L., thus representing mixed populations. A genome variant identical to the A genome of polyploid species was observed only in T. urartu. Species-specific molecular markers discriminating the diploid species were not found. Analysis of the inheritance of morphological characters also failed to identify a species-specific character for the three diploid wheat species apart from the hairy leaf blade type, described previously.     

Relative suitability of crested wheatgrass and other perennial grass hosts for the Russian wheat aphid (Homoptera: Aphididae)

The Russian wheat aphid, Diuraphis noxia (Mordvilko) (Homoptera: Aphididae), reproduces parthenogenetically in North America and must survive year-round on host plants, including in late summer when small grains are not in cultivation. During this time, cool-season perennial wheatgrasses (Poaceae: Triticeae) contribute substantially to aphid survival, crested wheatgrass (Agropyron spp.) particularly. In greenhouse studies, the number of aphids per plant was measured after four infestation periods on unvernalized and vernalized wheatgrasses. Before placement on these test plant species, aphids were reared either on winter wheat or on the grass host species on which aphid progeny were counted. On vernalized plants, aphids reared on wheat resulted in more aphids per test plant than when the aphids were reared on wheatgrasses, but on unvernalized plants the number of aphids per test plant did not differ significantly regardless of rearing host. Aphids on crested wheatgrass were similar in number to the other grasses when plants were unvernalized. However, when plants were vernalized, crested wheatgrass supported significantly more aphids than some of the other hosts. Aphid numbers increased on all test species as infestation period lengthened, and plant growth was largely unaffected by aphid feeding. These results suggest if sufficient moisture is available during summer when small grains are not in cultivation, all host species observed are capable of sustaining aphids. Crested wheatgrass is an abundant and important host of the Russian wheat aphid in its northern range of the western United States, but other less prevalent wheatgrasses also may contribute to aphid survival during late summer when small grains are not in cultivation.     

Orthology between genomes of Brachypodium, wheat and rice

Background

Differential expression of perforin (PRF1), a gene with a pivotal role in immune surveillance, can be attributed to differential methylation of CpG sites in its promoter region. A reproducible method for quantitative and CpG site-specific determination of perforin methylation is required for molecular epidemiologic studies of chronic diseases with immune dysfunction.

Findings

We developed a pyrosequencing based method to quantify site-specific methylation levels in 32 out of 34 CpG sites in the PRF1 promoter, and also compared methylation pattern in DNAs extracted from whole blood drawn into PAXgene blood DNA tubes (whole blood DNA) or DNA extracted from peripheral blood mononuclear cells (PBMC DNA) from the same normal subjects. Sodium bisulfite treatment of DNA and touchdown PCR were highly reproducible (coefficient of variation 1.63 to 2.18%) to preserve methylation information. Application of optimized pyrosequencing protocol to whole blood DNA revealed that methylation level varied along the promoter in normal subjects with extremely high methylation (mean 86%; range 82–92%) in the distal enhancer region (CpG sites 1–10), a variable methylation (range 49%–83%) in the methylation sensitive region (CpG sites 11–17), and a progressively declining methylation level (range 12%–80%) in the proximal promoter region (CpG sites 18–32) of PRF1. This pattern of methylation remained the same between whole blood and PBMC DNAs, but the absolute values of methylation in 30 out of 32 CpG sites differed significantly, with higher values for all CpG sites in the whole blood DNA.

Conclusion

This reproducible, site-specific and quantitative method for methylation determination of PRF1 based on pyrosequencing without cloning is well suited for large-scale molecular epidemiologic studies of diseases with immune dysfunction. PBMC DNA may be better suited than whole blood DNA for examining methylation levels in genes associated with immune function.     

OWEN: aligning long collinear regions of genomes

OWEN is an interactive tool for aligning two long DNA sequences that represents similarity between them by a chain of collinear local similarities. OWEN employs several methods for constructing and editing local similarities and for resolving conflicts between them. Alignments of sequences of lengths over 10(6) can often be produced in minutes. OWEN requires memory below 20 L, where L is the sum of lengths of the compared sequences.     

Birth-and-death evolution of protein-coding regions and concerted evolution of non-coding regions in the multi-component genomes of nanoviruses

Genomes of the four plant viruses of the genus Nanovirus consist of multiple circular single-stranded DNA components, each of which encodes a single protein. Protein phylogenies supported the hypothesis that faba bean necrotic yellows virus (FBNYV) and milk vetch disease virus (MDV) are sister taxa; that subterranean clover stunt virus (SCSV) branched next; and that banana bunchy top virus (BBTV) is an outgroup to the three other species. The phylogeny of replication (Rep) proteins indicate that this small viral multi-gene family has evolved by a process of duplication and subsequent loss of Rep-encoding genome components, analogous to the "birth-and-death" process of evolution which has been described in eukaryotic multi-gene families. By contrast, repeated recombinational events between components were found to have homogenized the non-coding portions of several components encoding unrelated components. For example, as result of recent recombination a portion of the non-coding region is virtually identical among SCSV components 1, 3, 4, 5, and 7. Thus, there is a process of concerted evolution of non-coding regions of Nanovirus genome components, which raises the possibility that certain non-coding regions are subject to functional constraint.     

Selective constraint on noncoding regions of hominid genomes

An important challenge for human evolutionary biology is to understand the genetic basis of human–chimpanzee differences. One influential idea holds that such differences depend, to a large extent, on adaptive changes in gene expression. An important step in assessing this hypothesis involves gaining a better understanding of selective constraint on noncoding regions of hominid genomes. In noncoding sequence, functional elements are frequently small and can be separated by large nonfunctional regions. For this reason, constraint in hominid genomes is likely to be patchy. Here we use conservation in more distantly related mammals and amniotes as a way of identifying small sequence windows that are likely to be functional. We find that putatively functional noncoding elements defined in this manner are subject to significant selective constraint in hominids.