Natural variation of the Y chromosome suppresses sex ratio distortion and modulates testis-specific gene expression in Drosophila simulans

X-linked sex-ratio distorters that disrupt spermatogenesis can cause a deficiency in functional Y-bearing sperm and a female-biased sex ratio. Y-linked modifiers that restore a normal sex ratio might be abundant and favored when a X-linked distorter is present. Here we investigated natural variation of Y-linked suppressors of sex-ratio in the Winters systems and the ability of these chromosomes to modulate gene expression in Drosophila simulans. Seventy-eight Y chromosomes of worldwide origin were assayed for their resistance to the X-linked sex-ratio distorter gene Dox. Y chromosome diversity caused males to sire ∼63% to ∼98% female progeny. Genome-wide gene expression analysis revealed hundreds of genes differentially expressed between isogenic males with sensitive (high sex ratio) and resistant (low sex ratio) Y chromosomes from the same population. Although the expression of about 75% of all testis-specific genes remained unchanged across Y chromosomes, a subset of post-meiotic genes was upregulated by resistant Y chromosomes. Conversely, a set of accessory gland-specific genes and mitochondrial genes were downregulated in males with resistant Y chromosomes. The D. simulans Y chromosome also modulated gene expression in XXY females in which the Y-linked protein-coding genes are not transcribed. The data suggest that the Y chromosome might exert its regulatory functions through epigenetic mechanisms that do not require the expression of protein-coding genes. The gene network that modulates sex ratio distortion by the Y chromosome is poorly understood, other than that it might include interactions with mitochondria and enriched for genes expressed in post-meiotic stages of spermatogenesis.     

X linkage of AP3A, a homolog of the Y-linked MADS-box gene AP3Y in Silene latifolia and S. dioica

Background

The duplication of autosomal genes onto the Y chromosome may be an important element in the evolution of sexual dimorphism.A previous cytological study reported on a putative example of such a duplication event in a dioecious tribe of Silene (Caryophyllaceae): it was inferred that the Y-linked MADS-box gene AP3Y originated from a duplication of the reportedly autosomal orthologAP3A. However, a recent study, also using cytological methods, indicated that AP3A is X-linked in Silenelatifolia.

Methodology/Principal Findings

In this study, we hybridized S. latifolia and S. dioicato investigate whether the pattern of X linkage is consistent among distinct populations, occurs in both species, and is robust to genetic methods. We found inheritance patterns indicative of X linkage of AP3A in widely distributed populations of both species.

Conclusions/Significance

X linkage ofAP3A and Y linkage of AP3Yin both species indicates that the genes'' ancestral progenitor resided on the autosomes that gave rise to the sex chromosomesand that neither gene has moved between chromosomes since species divergence.Consequently, our results do not support the contention that inter-chromosomal gene transfer occurred in the evolution of SlAP3Y from SlAP3A.     

AZFc somatic microdeletions and copy number polymorphism of the DAZ genes in human males exposed to natural background radiation

Ionizing radiations are known to induce tumors, chromosomal lesions and minisatellite length variations, yet no correlation has been demonstrated between radiation exposure and indels or copy number polymorphism (CNP) of the genes. We studied the impact of natural background radiation (NBR) on the human Y chromosome owing to its haploid status and clonal inheritance. We analyzed the AZFc region using the DNA from blood and semen of 100 males living near the coastal peninsula in Kerala (India), exposed to NBR along with other 50 normal fertile males. STS mapping of AZFc region showed random microdeletions without conclusive gr/gr or b1/b3 phenotypes. Using a highly specific novel Taqman assay based on sY587 sequence, we detected four copies of the DAZ genes in normal males and 4–16 in those exposed to NBR. Amongst NBR exposed males with multiples copies of the DAZ genes, 75% showed varying FISH signals for DAZ genes with cosmid 18E8 whereas 30% showed mosaicism in terms of presence/absence of the signals in 6–8% cells and unexpected number of signals in 9–12% interphase nuclei. Startlingly, all germline samples studied were found to be free from AZFc microdeletions and CNP of the DAZ genes. Since the DAZ genes are heavily implicated with the germ cell development, the cells with DAZ deletion/duplication are unlikely to survive. Alternatively, an innate mechanism may be operative to protect the germline from the effects of NBR.     

Genetic Variants in Nuclear-Encoded Mitochondrial Genes Influence AIDS Progression

Background

The human mitochondrial genome includes only 13 coding genes while nuclear-encoded genes account for 99% of proteins responsible for mitochondrial morphology, redox regulation, and energetics. Mitochondrial pathogenesis occurs in HIV patients and genetically, mitochondrial DNA haplogroups with presumed functional differences have been associated with differential AIDS progression.

Methodology/Principal Findings

Here we explore whether single nucleotide polymorphisms (SNPs) within 904 of the estimated 1,500 genes that specify nuclear-encoded mitochondrial proteins (NEMPs) influence AIDS progression among HIV-1 infected patients. We examined NEMPs for association with the rate of AIDS progression using genotypes generated by an Affymetrix 6.0 genotyping array of 1,455 European American patients from five US AIDS cohorts. Successfully genotyped SNPs gave 50% or better haplotype coverage for 679 of known NEMP genes. With a Bonferroni adjustment for the number of genes and tests examined, multiple SNPs within two NEMP genes showed significant association with AIDS progression: acyl-CoA synthetase medium-chain family member 4 (ACSM4) on chromosome 12 and peroxisomal D3,D2-enoyl-CoA isomerase (PECI) on chromosome 6.

Conclusions

Our previous studies on mitochondrial DNA showed that European haplogroups with presumed functional differences were associated with AIDS progression and HAART mediated adverse events. The modest influences of nuclear-encoded mitochondrial genes found in the current study add support to the idea that mitochondrial function plays a role in AIDS pathogenesis.     

Active Trachoma and Ocular Chlamydia trachomatis Infection in Two Gambian Regions: On Course for Elimination by 2020?

Background

Trachoma has been endemic in The Gambia for decades. National trachoma control activities have been in place since the mid-1980''s, but with no mass antibiotic treatment campaign. We aimed to assess the prevalence of active trachoma and of actual ocular Chlamydia trachomatis infection as measured by polymerase chain reaction (PCR) in the two Gambian regions that had had the highest prevalence of trachoma in the last national survey in 1996 prior to planned national mass antibiotic treatment distribution in 2006.

Methodology/Principal Findings

Two stage random sampling survey in 61 randomly selected Enumeration Areas (EAs) in North Bank Region (NBR) and Lower River Region (LRR). Fifty randomly selected children aged under 10 years were examined per EA for clinical signs of trachoma. In LRR, swabs were taken to test for ocular C. trachomatis infection. Unadjusted prevalences of active trachoma were calculated, as would be done in a trachoma control programme. The prevalence of trachomatous inflammation, follicular (TF) in the 2777 children aged 1–9 years was 12.3% (95% CI 8.8%–17.0%) in LRR and 10.0% (95% CI 7.7%–13.0%) in NBR, with significant variation within divisions (p<0.01), and a design effect of 3.474. Infection with C. trachomatis was found in only 0.3% (3/940) of children in LRR.

Conclusions/Significance

This study shows a large discrepancy between the prevalence of trachoma clinical signs and ocular C. trachomatis infection in two Gambian regions. Assessment of trachoma based on clinical signs alone may lead to unnecessary treatment, since the prevalence of active trachoma remains high but C. trachomatis infection has all but disappeared. Assuming that repeated infection is required for progression to blinding sequelae, blinding trachoma is on course for elimination by 2020 in The Gambia.     

Interchromosomal duplications on the Bactrocera oleae Y chromosome imply a distinct evolutionary origin of the sex chromosomes compared to Drosophila

Background

Diptera have an extraordinary variety of sex determination mechanisms, and Drosophila melanogaster is the paradigm for this group. However, the Drosophila sex determination pathway is only partially conserved and the family Tephritidae affords an interesting example. The tephritid Y chromosome is postulated to be necessary to determine male development. Characterization of Y sequences, apart from elucidating the nature of the male determining factor, is also important to understand the evolutionary history of sex chromosomes within the Tephritidae. We studied the Y sequences from the olive fly, Bactrocera oleae. Its Y chromosome is minute and highly heterochromatic, and displays high heteromorphism with the X chromosome.

Methodology/Principal Findings

A combined Representational Difference Analysis (RDA) and fluorescence in-situ hybridization (FISH) approach was used to investigate the Y chromosome to derive information on its sequence content. The Y chromosome is strewn with repetitive DNA sequences, the majority of which are also interdispersed in the pericentromeric regions of the autosomes. The Y chromosome appears to have accumulated small and large repetitive interchromosomal duplications. The large interchromosomal duplications harbour an importin-4-like gene fragment. Apart from these importin-4-like sequences, the other Y repetitive sequences are not shared with the X chromosome, suggesting molecular differentiation of these two chromosomes. Moreover, as the identified Y sequences were not detected on the Y chromosomes of closely related tephritids, we can infer divergence in the repetitive nature of their sequence contents.

Conclusions/Significance

The identification of Y-linked sequences may tell us much about the repetitive nature, the origin and the evolution of Y chromosomes. We hypothesize how these repetitive sequences accumulated and were maintained on the Y chromosome during its evolutionary history. Our data reinforce the idea that the sex chromosomes of the Tephritidae may have distinct evolutionary origins with respect to those of the Drosophilidae and other Dipteran families.     

Frequent MAGE Mutations in Human Melanoma

Background

Cancer/testis (CT) genes are expressed only in the germ line and certain tumors and are most frequently located on the X-chromosome (the CT-X genes). Amongst the best studied CT-X genes are those encoding several MAGE protein families. The function of MAGE proteins is not well understood, but several have been shown to potentially influence the tumorigenic phenotype.

Methodology/Principal Findings

We undertook a mutational analysis of coding regions of four CT-X MAGE genes, MAGEA1, MAGEA4, MAGEC1, MAGEC2 and the ubiquitously expressed MAGEE1 in human melanoma samples. We first examined cell lines established from tumors and matching blood samples from 27 melanoma patients. We found that melanoma cell lines from 37% of patients contained at least one mutated MAGE gene. The frequency of mutations in the coding regions of individual MAGE genes varied from 3.7% for MAGEA1 and MAGEA4 to 14.8% for MAGEC2. We also examined 111 fresh melanoma samples collected from 86 patients. In this case, samples from 32% of the patients exhibited mutations in one or more MAGE genes with the frequency of mutations in individual MAGE genes ranging from 6% in MAGEA1 to 16% in MAGEC1.

Significance

These results demonstrate for the first time that the MAGE gene family is frequently mutated in melanoma.     

Sex-Related Differences in Gene Expression Following Coxiella burnetii Infection in Mice: Potential Role of Circadian Rhythm

Background

Q fever, a zoonosis due to Coxiella burnetii infection, exhibits sexual dimorphism; men are affected more frequently and severely than women for a given exposure. Here we explore whether the severity of C. burnetii infection in mice is related to differences in male and female gene expression profiles.

Methodology/Principal Findings

Mice were infected with C. burnetii for 24 hours, and gene expression was measured in liver cells using microarrays. Multiclass analysis identified 2,777 probes for which expression was specifically modulated by C. burnetti infection. Only 14% of the modulated genes were sex-independent, and the remaining 86% were differentially expressed in males and females. Castration of males and females showed that sex hormones were responsible for more than 60% of the observed gene modulation, and this reduction was most pronounced in males. Using functional annotation of modulated genes, we identified four clusters enriched in males that were related to cell-cell adhesion, signal transduction, defensins and cytokine/Jak-Stat pathways. Up-regulation of the IL-10 and Stat-3 genes may account for the high susceptibility of men with Q fever to C. burnetii infection and autoantibody production. Two clusters were identified in females, including the circadian rhythm pathway, which consists of positive (Clock, Arntl) and negative (Per) limbs of a feedback loop. We found that Clock and Arntl were down-modulated whereas Per was up-regulated; these changes may be associated with efficient bacterial elimination in females but not in males, in which an exacerbated host response would be prominent.

Conclusion

This large-scale study revealed for the first time that circadian rhythm plays a major role in the anti-infectious response of mice, and it provides a new basis for elucidating the role of sexual dimorphism in human infections.     

Identifying new sex-linked genes through BAC sequencing in the dioecious plant Silene latifolia

Background

Silene latifolia represents one of the best-studied plant sex chromosome systems. A new approach using RNA-seq data has recently identified hundreds of new sex-linked genes in this species. However, this approach is expected to miss genes that are either not expressed or are expressed at low levels in the tissue(s) used for RNA-seq. Therefore other independent approaches are needed to discover such sex-linked genes.

Results

Here we used 10 well-characterized S. latifolia sex-linked genes and their homologs in Silene vulgaris, a species without sex chromosomes, to screen BAC libraries of both species. We isolated and sequenced 4 Mb of BAC clones of S. latifolia X and Y and S. vulgaris genomic regions, which yielded 59 new sex-linked genes (with S. vulgaris homologs for some of them). We assembled sequences that we believe represent the tip of the Xq arm. These sequences are clearly not pseudoautosomal, so we infer that the S. latifolia X has a single pseudoautosomal region (PAR) on the Xp arm. The estimated mean gene density in X BACs is 2.2 times lower than that in S. vulgaris BACs, agreeing with the genome size difference between these species. Gene density was estimated to be extremely low in the Y BAC clones. We compared our BAC-located genes with the sex-linked genes identified in previous RNA-seq studies, and found that about half of them (those with low expression in flower buds) were not identified as sex-linked in previous RNA-seq studies. We compiled a set of ~70 validated X/Y genes and X-hemizygous genes (without Y copies) from the literature, and used these genes to show that X-hemizygous genes have a higher probability of being undetected by the RNA-seq approach, compared with X/Y genes; we used this to estimate that about 30 % of our BAC-located genes must be X-hemizygous. The estimate is similar when we use BAC-located genes that have S. vulgaris homologs, which excludes genes that were gained by the X chromosome.

Conclusions

Our BAC sequencing identified 59 new sex-linked genes, and our analysis of these BAC-located genes, in combination with RNA-seq data suggests that gene losses from the S. latifolia Y chromosome could be as high as 30 %, higher than previous estimates of 10-20 %.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1698-7) contains supplementary material, which is available to authorized users.     

Association between Common Germline Genetic Variation in 94 Candidate Genes or Regions and Risks of Invasive Epithelial Ovarian Cancer

Background

Recent studies have identified several single nucleotide polymorphisms (SNPs) in the population that are associated with variations in the risks of many different diseases including cancers such as breast, prostate and colorectal. For ovarian cancer, the known highly penetrant susceptibility genes (BRCA1 and BRCA2) are probably responsible for only 40% of the excess familial ovarian cancer risks, suggesting that other susceptibility genes of lower penetrance exist.

Methods

We have taken a candidate approach to identifying moderate risk susceptibility alleles for ovarian cancer. To date, we have genotyped 340 SNPs from 94 candidate genes or regions, in up to 1,491 invasive epithelial ovarian cancer cases and 3,145 unaffected controls from three different population based studies from the UK, Denmark and USA.

Results

After adjusting for population stratification by genomic control, 18 SNPs (5.3%) were significant at the 5% level, and 5 SNPs (1.5%) were significant at the 1% level. The most significant association was for the SNP rs2107425, located on chromosome 11p15.5, which has previously been identified as a susceptibility allele for breast cancer from a genome wide association study (P-trend = 0.0012). When SNPs/genes were stratified into 7 different pathways or groups of validation SNPs, the breast cancer associated SNPs were the only group of SNPs that were significantly associated with ovarian cancer risk (P-heterogeneity = 0.0003; P-trend = 0.0028; adjusted (for population stratification) P-trend = 0.006). We did not find statistically significant associations when the combined data for all SNPs were analysed using an admixture maximum likelihood (AML) experiment-wise test for association (P-heterogeneity = 0.051; P-trend = 0.068).

Conclusion

These data suggest that a proportion of the SNPs we evaluated were associated with ovarian cancer risk, but that the effect sizes were too small to detect associations with individual SNPs.     

Trichomonas Transmembrane Cyclases Result from Massive Gene Duplication and Concomitant Development of Pseudogenes

Background

Trichomonas vaginalis has an unusually large genome (∼160 Mb) encoding ∼60,000 proteins. With the goal of beginning to understand why some Trichomonas genes are present in so many copies, we characterized here a family of ∼123 Trichomonas genes that encode transmembrane adenylyl cyclases (TMACs).

Methodology/Principal Findings

The large family of TMACs genes is the result of recent duplications of a small set of ancestral genes that appear to be unique to trichomonads. Duplicated TMAC genes are not closely associated with repetitive elements, and duplications of flanking sequences are rare. However, there is evidence for TMAC gene replacements by homologous recombination. A high percentage of TMAC genes (∼46%) are pseudogenes, as they contain stop codons and/or frame shifts, or the genes are truncated. Numerous stop codons present in the genome project G3 strain are not present in orthologous genes of two other Trichomonas strains (S1 and B7RC2). Each TMAC is composed of a series of N-terminal transmembrane helices and a single C-terminal cyclase domain that has adenylyl cyclase activity. Multiple TMAC genes are transcribed by Trichomonas cloned by limiting dilution.

Conclusions/Significance

We conclude that one reason for the unusually large genome of Trichomonas is the presence of unstable families of genes such as those encoding TMACs that are undergoing massive gene duplication and concomitant development of pseudogenes.     

Diversity of 23S rRNA Genes within Individual Prokaryotic Genomes

Background

The concept of ribosomal constraints on rRNA genes is deduced primarily based on the comparison of consensus rRNA sequences between closely related species, but recent advances in whole-genome sequencing allow evaluation of this concept within organisms with multiple rRNA operons.

Methodology/Principal Findings

Using the 23S rRNA gene as an example, we analyzed the diversity among individual rRNA genes within a genome. Of 184 prokaryotic species containing multiple 23S rRNA genes, diversity was observed in 113 (61.4%) genomes (mean 0.40%, range 0.01%–4.04%). Significant (1.17%–4.04%) intragenomic variation was found in 8 species. In 5 of the 8 species, the diversity in the primary structure had only minimal effect on the secondary structure (stem versus loop transition). In the remaining 3 species, the diversity significantly altered local secondary structure, but the alteration appears minimized through complex rearrangement. Intervening sequences (IVS), ranging between 9 and 1471 nt in size, were found in 7 species. IVS in Deinococcus radiodurans and Nostoc sp. encode transposases. T. tengcongensis was the only species in which intragenomic diversity >3% was observed among 4 paralogous 23S rRNA genes.

Conclusions/Significance

These findings indicate tight ribosomal constraints on individual 23S rRNA genes within a genome. Although classification using primary 23S rRNA sequences could be erroneous, significant diversity among paralogous 23S rRNA genes was observed only once in the 184 species analyzed, indicating little overall impact on the mainstream of 23S rRNA gene-based prokaryotic taxonomy.     

Genome Sequence of Cronobacter sakazakii BAA-894 and Comparative Genomic Hybridization Analysis with Other Cronobacter Species

Background

The genus Cronobacter (formerly called Enterobacter sakazakii) is composed of five species; C. sakazakii, C. malonaticus, C. turicensis, C. muytjensii, and C. dublinensis. The genus includes opportunistic human pathogens, and the first three species have been associated with neonatal infections. The most severe diseases are caused in neonates and include fatal necrotizing enterocolitis and meningitis. The genetic basis of the diversity within the genus is unknown, and few virulence traits have been identified.

Methodology/Principal Findings

We report here the first sequence of a member of this genus, C. sakazakii strain BAA-894. The genome of Cronobacter sakazakii strain BAA-894 comprises a 4.4 Mb chromosome (57% GC content) and two plasmids; 31 kb (51% GC) and 131 kb (56% GC). The genome was used to construct a 387,000 probe oligonucleotide tiling DNA microarray covering the whole genome. Comparative genomic hybridization (CGH) was undertaken on five other C. sakazakii strains, and representatives of the four other Cronobacter species. Among 4,382 annotated genes inspected in this study, about 55% of genes were common to all C. sakazakii strains and 43% were common to all Cronobacter strains, with 10–17% absence of genes.

Conclusions/Significance

CGH highlighted 15 clusters of genes in C. sakazakii BAA-894 that were divergent or absent in more than half of the tested strains; six of these are of probable prophage origin. Putative virulence factors were identified in these prophage and in other variable regions. A number of genes unique to Cronobacter species associated with neonatal infections (C. sakazakii, C. malonaticus and C. turicensis) were identified. These included a copper and silver resistance system known to be linked to invasion of the blood-brain barrier by neonatal meningitic strains of Escherichia coli. In addition, genes encoding for multidrug efflux pumps and adhesins were identified that were unique to C. sakazakii strains from outbreaks in neonatal intensive care units.     

Correlating Global Gene Regulation to Angiogenesis in the Developing Chick Extra-Embryonic Vascular System

Background

Formation of blood vessels requires the concerted regulation of an unknown number of genes in a spatial-, time- and dosage-dependent manner. Determining genes, which drive vascular maturation is crucial for the identification of new therapeutic targets against pathological angiogenesis.

Methology/Principal Findings

We accessed global gene regulation throughout maturation of the chick chorio-allantoic membrane (CAM), a highly vascularized tissue, using pan genomic microarrays. Seven percent of analyzed genes showed a significant change in expression (>2-fold, FDR<5%) with a peak occurring from E7 to E10, when key morphogenetic and angiogenic genes such as BMP4, SMO, HOXA3, EPAS1 and FGFR2 were upregulated, reflecting the state of an activated endothelium. At later stages, a general decrease in gene expression occurs, including genes encoding mitotic factors or angiogenic mediators such as CYR61, EPAS1, MDK and MYC. We identified putative human orthologs for 77% of significantly regulated genes and determined endothelial cell enrichment for 20% of the orthologs in silico. Vascular expression of several genes including ENC1, FSTL1, JAM2, LDB2, LIMS1, PARVB, PDE3A, PRCP, PTRF and ST6GAL1 was demonstrated by in situ hybridization. Up to 9% of the CAM genes were also overexpressed in human organs with related functions, such as placenta and lung or the thyroid. 21–66% of CAM genes enriched in endothelial cells were deregulated in several human cancer types (P<.0001). Interfering with PARVB (encoding parvin, beta) function profoundly changed human endothelial cell shape, motility and tubulogenesis, suggesting an important role of this gene in the angiogenic process.

Conclusions/Significance

Our study underlines the complexity of gene regulation in a highly vascularized organ during development. We identified a restricted number of novel genes enriched in the endothelium of different species and tissues, which may play crucial roles in normal and pathological angiogenesis.     

Differences in the number of intrinsically disordered regions between yeast duplicated proteins, and their relationship with functional divergence

Background

Intrinsically disordered regions are enriched in short interaction motifs that play a critical role in many protein-protein interactions. Since new short interaction motifs may easily evolve, they have the potential to rapidly change protein interactions and cellular signaling. In this work we examined the dynamics of gain and loss of intrinsically disordered regions in duplicated proteins to inspect if changes after genome duplication can create functional divergence. For this purpose we used Saccharomyces cerevisiae and the outgroup species Lachancea kluyveri.

Principal Findings

We find that genes duplicated as part of a genome duplication (ohnologs) are significantly more intrinsically disordered than singletons (p<2.2^e-16, Wilcoxon), reflecting a preference for retaining intrinsically disordered proteins in duplicate. In addition, there have been marked changes in the extent of intrinsic disorder following duplication. A large number of duplicated genes have more intrinsic disorder than their L. kluyveri ortholog (29% for duplicates versus 25% for singletons) and an even greater number have less intrinsic disorder than the L. kluyveri ortholog (37% for duplicates versus 25% for singletons). Finally, we show that the number of physical interactions is significantly greater in the more intrinsically disordered ohnolog of a pair (p = 0.003, Wilcoxon).

Conclusion

This work shows that intrinsic disorder gain and loss in a protein is a mechanism by which a genome can also diverge and innovate. The higher number of interactors for proteins that have gained intrinsic disorder compared with their duplicates may reflect the acquisition of new interaction partners or new functional roles.     

Sequence composition and gene content of the short arm of rye (Secale cereale) chromosome 1

Background

The purpose of the study is to elucidate the sequence composition of the short arm of rye chromosome 1 (Secale cereale) with special focus on its gene content, because this portion of the rye genome is an integrated part of several hundreds of bread wheat varieties worldwide.

Methodology/Principal Findings

Multiple Displacement Amplification of 1RS DNA, obtained from flow sorted 1RS chromosomes, using 1RS ditelosomic wheat-rye addition line, and subsequent Roche 454FLX sequencing of this DNA yielded 195,313,589 bp sequence information. This quantity of sequence information resulted in 0.43× sequence coverage of the 1RS chromosome arm, permitting the identification of genes with estimated probability of 95%. A detailed analysis revealed that more than 5% of the 1RS sequence consisted of gene space, identifying at least 3,121 gene loci representing 1,882 different gene functions. Repetitive elements comprised about 72% of the 1RS sequence, Gypsy/Sabrina (13.3%) being the most abundant. More than four thousand simple sequence repeat (SSR) sites mostly located in gene related sequence reads were identified for possible marker development. The existence of chloroplast insertions in 1RS has been verified by identifying chimeric chloroplast-genomic sequence reads. Synteny analysis of 1RS to the full genomes of Oryza sativa and Brachypodium distachyon revealed that about half of the genes of 1RS correspond to the distal end of the short arm of rice chromosome 5 and the proximal region of the long arm of Brachypodium distachyon chromosome 2. Comparison of the gene content of 1RS to 1HS barley chromosome arm revealed high conservation of genes related to chromosome 5 of rice.

Conclusions

The present study revealed the gene content and potential gene functions on this chromosome arm and demonstrated numerous sequence elements like SSRs and gene-related sequences, which can be utilised for future research as well as in breeding of wheat and rye.