Sequencing and characterization of telomere and subtelomere regions on rice chromosomes 1S, 2S, 2L, 6L, 7S, 7L and 8S

Telomere and subtelomere high polymorphism might contribute to the specificity of homologous recognition and pairing during meiosis in barley in the context of breeding

Barley (Hordeum vulgare) is one of the most popular cereal crops globally. Although it is a diploid species, (2n = 2x = 14) the study of its genome organization is necessary in the framework of plant breeding since barley is often used in crosses with other cereals like wheat to provide them with advantageous characters. We already have an extensive knowledge on different stages of the meiosis, the cell division to generate the gametes in species with sexual reproduction, such as the formation of the synaptonemal complex, recombination, and chromosome segregation. But meiosis really starts with the identification of homologous chromosomes and pairing initiation, and it is still unclear how chromosomes exactly choose a partner to appropriately pair for additional recombination and segregation. In this work we present an exhaustive molecular analysis of both telomeres and subtelomeres of barley chromosome arms 2H-L, 3H-L and 5H-L. As expected, the analysis of multiple features, including transposable elements, repeats, GC content, predicted CpG islands, recombination hotspots, G4 quadruplexes, genes and targeted sequence motifs for key DNA-binding proteins, revealed a high degree of variability both in telomeres and subtelomeres. The molecular basis for the specificity of homologous recognition and pairing occurring in the early chromosomal interactions at the start of meiosis in barley may be provided by these polymorphisms. A more relevant role of telomeres and most distal part of subtelomeres is suggested.

A complete telomere-to-telomere assembly of the maize genome

A complete telomere-to-telomere (T2T) finished genome has been the long pursuit of genomic research. Through generating deep coverage ultralong Oxford Nanopore Technology (ONT) and PacBio HiFi reads, we report here a complete genome assembly of maize with each chromosome entirely traversed in a single contig. The 2,178.6 Mb T2T Mo17 genome with a base accuracy of over 99.99% unveiled the structural features of all repetitive regions of the genome. There were several super-long simple-sequence-repeat arrays having consecutive thymine-adenine-guanine (TAG) tri-nucleotide repeats up to 235 kb. The assembly of the entire nucleolar organizer region of the 26.8 Mb array with 2,974 45S rDNA copies revealed the enormously complex patterns of rDNA duplications and transposon insertions. Additionally, complete assemblies of all ten centromeres enabled us to precisely dissect the repeat compositions of both CentC-rich and CentC-poor centromeres. The complete Mo17 genome represents a major step forward in understanding the complexity of the highly recalcitrant repetitive regions of higher plant genomes.

The RNA binding protein OsLa influences grain and anther development in rice

La proteins are found widely in eukaryotes and play a variety of vital roles. AtLa1 has been identified as an La protein that is necessary for embryogenesis in Arabidopsis; however, the existence and biological functions of La proteins in rice (Oryza sativa L.) remain unclear. In this study, we identified and characterized two La proteins in rice that are homologous to AtLa1 and named them OsLa1 and OsLa2. Both the OsLa1 and OsLa2 genes encode RNA-binding proteins with an La domain and two RNA-binding domains. Mutant OsLa1 reduced grain length and pollen fertility, whereas OsLa1 overexpression caused the opposite phenotypes. Further experiments indicated that OsLa1 modulates grain size by influencing cell expansion. Interestingly, mutant OsLa2 resulted in thin grains with decreased weight and a low seed-setting rate. We also found that OsLa1 interacted with OsLa2 and that both OsLa1 and OsLa2 interacted with OseIF6.1, a eukaryotic translation initiation factor involved in ribosome biogenesis. In addition, OsLa1 was able to bind to OseIF6.1 mRNA to modulate its expression. Complete OseIF6.1 knockout caused lethality and OseIF6.1/oseif6.1 heterozygous plants displayed low fertility and low seed setting. Together, our results enrich our knowledge of the role of La proteins in rice growth and development, as well as the relationship between La and eIF6 in rice.

Rice OsBRCA2 Is Required for DNA Double-Strand Break Repair in Meiotic Cells

The mammalian BREAST CANCER 2 (BRCA2) gene is a tumor suppressor that plays a crucial role in DNA repair and homologous recombination (HR). Here, we report the identification and characterization of OsBRCA2, the rice orthologue of human BRCA2. Osbrca2 mutant plants exhibit normal vegetative growth but experience complete male and female sterility as a consequence of severe meiotic defects. Pairing, synapsis and recombination are impaired in osbrca2 male meiocytes, leading to chromosome entanglements and fragmentation. In the absence of OsBRCA2, localization to the meiotic chromosome axes of the strand-invasion proteins OsRAD51 and OsDMC1 is severely reduced and in vitro OsBRCA2 directly interacts with OsRAD51 and OsDMC1. These results indicate that OsBRCA2 is essential for facilitating the loading of OsRAD51 and OsDMC1 onto resected ends of programmed double-strand breaks (DSB) during meiosis to promote single-end invasions of homologous chromosomes and accurate recombination. In addition, treatment of osbrca2-1 seedlings with mitomycin C (MMC) led to hypersensitivity. As MMC is a genotoxic agent that creates DNA lesions in the somatic cells that can only be repaired by HR, these results suggest that OsBRCA2 has a conserved role in DSB repair and HR in rice.

Sequence analysis of wheat subtelomeres reveals a high polymorphism among homoeologous chromosomes

Bread wheat, Triticum aestivum L., is one of the most important crops in the world. Understanding its genome organization (allohexaploid; AABBDD; 2n = 6x = 42) is essential for geneticists and plant breeders. Particularly, the knowledge of how homologous chromosomes (equivalent chromosomes from the same genome) specifically recognize each other to pair at the beginning of meiosis, the cellular process to generate gametes in sexually reproducing organisms, is fundamental for plant breeding and has a big influence on the fertility of wheat plants. Initial homologous chromosome interactions contribute to specific recognition and pairing between homologues at the onset of meiosis. Understanding the molecular basis of these critical processes can help to develop genetic tools in a breeding context to promote interspecific chromosome associations in hybrids or interspecific genetic crosses to facilitate the transfer of desirable agronomic traits from related species into a crop like wheat. The terminal regions of chromosomes, which include telomeres and subtelomeres, participate in chromosome recognition and pairing. We present a detailed molecular analysis of subtelomeres of wheat chromosome arms 1AS, 4AS, 7AS, 7BS and 7DS. Results showed a high polymorphism in the subtelomeric region among homoeologues (equivalent chromosomes from related genomes) for all the features analyzed, including genes, transposable elements, repeats, GC content, predicted CpG islands, recombination hotspots and targeted sequence motifs for relevant DNA-binding proteins. These polymorphisms might be the molecular basis for the specificity of homologous recognition and pairing in initial chromosome interactions at the beginning of meiosis in wheat.

Bioinformatic and Molecular Analysis of Satellite Repeat Diversity in Vaccinium Genomes

Bioinformatic and molecular characterization of satellite repeats was performed to understand the impact of their diversification on Vaccinium genome evolution. Satellite repeat diversity was evaluated in four cultivated and wild species, including the diploid species Vaccinium myrtillus and Vaccinium uliginosum, as well as the tetraploid species Vaccinium corymbosum and Vaccinium arctostaphylos. We comparatively characterized six satellite repeat families using in total 76 clones with 180 monomers. We observed that the monomer units of VaccSat1, VaccSat2, VaccSat5, and VaccSat6 showed a higher order repeat (HOR) structure, likely originating from the organization of two adjacent subunits with differing similarity, length and size. Moreover, VaccSat1, VaccSat3, VaccSat6, and VaccSat7 were found to have sequence similarity to parts of transposable elements. We detected satellite-typical tandem organization for VaccSat1 and VaccSat2 in long arrays, while VaccSat5 and VaccSat6 distributed in multiple sites over all chromosomes of tetraploid V. corymbosum, presumably in long arrays. In contrast, very short arrays of VaccSat3 and VaccSat7 are dispersedly distributed over all chromosomes in the same species, likely as internal parts of transposable elements. We provide a comprehensive overview on satellite species specificity in Vaccinium, which are potentially useful as molecular markers to address the taxonomic complexity of the genus, and provide information for genome studies of this genus.

Cloning and preliminary verification of telomere-associated sequences in upland cotton

Telomeres are structures enriched in repetitive sequences at the end of chromosomes. In this study, using the telomere primer AA(CCCTAAA)₃CCC for the single primer PCR, two DNA sequences were obtained from Gossypium hirsutum (Linnaeus, 1753) accession (acc.) TM-1. Sequence analysis showed that the two obtained sequences were all rich in A/T base, which was consistent with the characteristic of the telomere-associated sequence (TAS). They were designated as GhTAS1 and GhTAS2 respectively. GhTAS1 is 489 bp long, with 57.6% of A/T, and GhTAS2 is 539 bp long, with 63.9% of A/T. Fluorescence in situ hybridization results showed that both of the cloned TASs were located at the ends of the partial chromosomes of G. hirsutum, with the strong signals, which further confirmed that GhTAS1 and GhTAS2 were telomere-associated sequences including highly tandemly repetitive sequences. Results of blast against the assembled genome of G. hirsutum showed that GhTAS sequences may be missed on some assembled chromosomes. The results provide important evidence for the evaluation of the integrity of assembled chromosome end sequences, and will also contribute to the further perfection of the draft genomes of cotton.

Genome-wide analysis of MIKC-type MADS-box genes in wheat: pervasive duplications, functional conservation and putative neofunctionalization

Wheat (Triticum aestivum) is one of the most important crops worldwide. Given a growing global population coupled with increasingly challenging cultivation conditions, facilitating wheat breeding by fine-tuning important traits is of great importance. MADS-box genes are prime candidates for this, as they are involved in virtually all aspects of plant development. Here, we present a detailed overview of phylogeny and expression of 201 wheat MIKC-type MADS-box genes. Homoeolog retention is significantly above the average genome-wide retention rate for wheat genes, indicating that many MIKC-type homoeologs are functionally important and not redundant. Gene expression is generally in agreement with the expected subfamily-specific expression pattern, indicating broad conservation of function of MIKC-type genes during wheat evolution. We also found extensive expansion of some MIKC-type subfamilies, especially those potentially involved in adaptation to different environmental conditions like flowering time genes. Duplications are especially prominent in distal telomeric regions. A number of MIKC-type genes show novel expression patterns and respond, for example, to biotic stress, pointing towards neofunctionalization. We speculate that conserved, duplicated and neofunctionalized MIKC-type genes may have played an important role in the adaptation of wheat to a diversity of conditions, hence contributing to the importance of wheat as a global staple food.

N6-Methyladenine DNA Methylation in Japonica and Indica Rice Genomes and Its Association with Gene Expression, Plant Development, and Stress Responses

N⁶-Methyladenine (6mA) DNA methylation has recently been implicated as a potential new epigenetic marker in eukaryotes, including the dicot model Arabidopsis thaliana. However, the conservation and divergence of 6mA distribution patterns and functions in plants remain elusive. Here we report high-quality 6mA methylomes at single-nucleotide resolution in rice based on substantially improved genome sequences of two rice cultivars, Nipponbare (Nip; Japonica) and 93-11 (Indica). Analysis of 6mA genomic distribution and its association with transcription suggest that 6mA distribution and function is rather conserved between rice and Arabidopsis. We found that 6mA levels are positively correlated with the expression of key stress-related genes, which may be responsible for the difference in stress tolerance between Nip and 93-11. Moreover, we showed that mutations in DDM1 cause defects in plant growth and decreased 6mA level. Our results reveal that 6mA is a conserved DNA modification that is positively associated with gene expression and contributes to key agronomic traits in plants.

The Nipponbare genome and the next-generation of rice genomics research in Japan

The map-based genome sequence of the japonica rice cultivar Nipponbare remains to date as the only monocot genome that has been sequenced to a high-quality level. It has become the reference sequence for understanding the diversity among thousands of rice cultivars and its wild relatives as well as the major cereal crops that comprised the food source for the entire human race. This review focuses on the accomplishments in rice genomics in Japan encompassing the last 10 years which have led into deeper understanding of the genome, characterization of many agronomic traits, comprehensive analysis of the transcriptome, and the map-based cloning of many genes associated with agronomic traits.

Telomere-binding protein regulates the chromosome ends through the interaction with histone deacetylases in Arabidopsis thaliana

Telomeres are nucleoprotein complexes at the end of eukaryotic chromosomes. Many telomere-binding proteins bind to telomeric repeat sequences and further generate T-loops in animals. However, it is not clear if they regulate telomere organization using epigenetic mechanisms and how the epigenetic molecules are involved in regulating the telomeres. Here, we show direct interactions between the telomere-binding protein, AtTRB2 and histone deacetylases, HDT4 and HDA6, in vitro and in vivo AtTRB2 mediates the associations of HDT4 and HDA6 with telomeric repeats. Telomere elongation is found in AtTRB2, HDT4 and HDA6 mutants over generations, but also in met1 and cmt3 DNA methyltransferases mutants. We also characterized HDT4 as an Arabidopsis H3K27 histone deacetylase. HDT4 binds to acetylated peptides at residue K27 of histone H3 in vitro, and deacetylates this residue in vivo Our results suggest that AtTRB2 also has a role in the regulation of telomeric chromatin as a possible scaffold protein for recruiting the epigenetic regulators in Arabidopsis, in addition to its telomere binding and length regulation activity. Our data provide evidences that epigenetic molecules associate with telomeres by direct physical interaction with telomere-binding proteins and further regulate homeostasis of telomeres in Arabidopsis thaliana.

Interstitial telomeric sites and Robertsonian translocations in species of Ipheion and Nothoscordum (Amaryllidaceae)

The genera Nothoscordum and Ipheion (Allioideae, Amaryllidaceae) are cytologically characterized by a dysploid series with variable numbers of metacentric and acrocentric chromosomes typical of karyotypes rearranged by Robertsonian translocations (RT). Since they have large chromosomes, low diploid numbers, and possess two telomeric motifs [the vertebrate-type (TTAGGG) n and the Arabidopsis-type (TTTAGGG) n ] they are suitable for investigating the occurrence and possible role of interstitial telomeric sites (ITS) associated with RT. We analyzed the distributions of telomeric sites in 12 species of Nothoscordum and Ipheion and found that both telomeric probes colocalized in all chromosome termini. Cloning and sequencing PCR products obtained using both telomeric primers simultaneously revealed long stretches of (TTAGGG) n and (TTTAGGG) n sequences together with degenerated telomeric sequences. Most acrocentric chromosomes have a 45S rDNA site at the terminal region of the short arms adjacent to the most distal telomeric sites. Telomeric signals were found at all chromosome ends, but ITS were also detected in a few proximal and subterminal regions in some Nothoscordum species. Although RT are common in this group of plants, our findings suggest that proximal positioning of telomeric motifs are not necessarily related to that kind of rearrangement. Rather, transposition of telomeric sequences followed by amplification, could better explain the presence of (TTAGGG) n and (TTTAGGG) n repeats at those sites. Furthermore, a few small interstitial sites found in some Nothoscordum species indicate that dispersion of these sequences was not restricted to the proximal region.

Extending the model of Arabidopsis telomere length and composition across Brassicaceae

Telomeres are repetitive TG-rich DNA elements essential for maintaining the stability of genomes and replicative capacity of cells in almost all eukaryotes. Most of what is known about telomeres in plants comes from the angiosperm Arabidopsis thaliana, which has become an important comparative model for telomere biology. Arabidopsis tolerates numerous insults to its genome, many of which are catastrophic or lethal in other eukaryotic systems such as yeast and vertebrates. Despite the importance of Arabidopsis in establishing a model for the structure and regulation of plant telomeres, only a handful of studies have used this information to assay components of telomeres from across land plants, or even among the closest relatives of Arabidopsis in the plant family Brassicaceae. Here, we determined how well Arabidopsis represents Brassicaceae by comparing multiple aspects of telomere biology in species that represent major clades in the family tree. Specifically, we determined the telomeric repeat sequence, measured bulk telomere length, and analyzed variation in telomere length on syntenic chromosome arms. In addition, we used a phylogenetic approach to infer the evolutionary history of putative telomere-binding proteins, CTC1, STN1, TEN1 (CST), telomere repeat-binding factor like (TRFL), and single Myb histone (SMH). Our analyses revealed conservation of the telomeric DNA repeat sequence, but considerable variation in telomere length among the sampled species, even in comparisons of syntenic chromosome arms. We also found that the single-stranded and double-stranded telomeric DNA-binding complexes CST and TRFL, respectively, differ in their pattern of gene duplication and loss. The TRFL and SMH gene families have undergone numerous duplication events, and these duplicate copies are often retained in the genome. In contrast, CST components occur as single-copy genes in all sampled genomes, even in species that experienced recent whole genome duplication events. Taken together, our results place the Arabidopsis model in the context of other species in Brassicaceae, making the family the best characterized plant group in regard to telomere architecture.

The Subtelomeric khipu Satellite Repeat from Phaseolus vulgaris: Lessons Learned from the Genome Analysis of the Andean Genotype G19833

Subtelomeric regions in eukaryotic organisms are known for harboring species-specific tandemly repeated satellite sequences. However, studies on the molecular organization and evolution of subtelomeric repeats are scarce, especially in plants. Khipu is a satellite DNA of 528-bp repeat unit, specific of the Phaseolus genus, with a subtelomeric distribution in common bean, P. vulgaris. To investigate the genomic organization and the evolution of khipu, we performed genome-wide analysis on the complete genome sequence of the common bean genotype G19833. We identified 2,460 khipu units located at most distal ends of the sequenced regions. Khipu units are arranged in discrete blocks of 2-55 copies and are heterogeneously distributed among the different chromosome ends of G19833 (from 0 to 555 khipus units per chromosome arm). Phylogenetically related khipu units are spread between numerous chromosome ends, suggesting frequent exchanges between non-homologous subtelomeres. However, most subclades contain numerous khipu units from only one or few chromosome ends indicating that local duplication is also driving khipu expansion. Unexpectedly, we also identified 81 khipu units located at centromeres. All the centromeric khipu units belong to a single divergent clade also comprised of a few units from several subtelomeres, suggesting that a few sequence exchanges between centromeres and subtelomeres took place in the common bean genome. The divergence and low copy number of these centromeric units from the subtelomeric units could explain why they were not detected by FISH (Fluorescence in situ Hybridization) although it can not be excluded that these centromeric units may have resulted from errors in the pseudomolecule assembly. Altogether our data highlight extensive sequence exchanges in subtelomeres between non-homologous chromosomes in common bean and confirm that subtelomeres represent one of the most dynamic and rapidly evolving regions in eukaryotic genomes.

Systematic application of DNA fiber-FISH technique in cotton

Fluorescence in situ hybridization on extended DNA (fiber-FISH) is a powerful tool in high-resolution physical mapping. To introduce this technique into cotton, we developed the technique and tested it by deliberately mapping of telomere and 5S rDNA. Results showed that telomere-length ranged from 0.80 kb to 37.86 kb in three species, G. hirsutum, G. herbaceum and G. arboreum. However, most of the telomeres (>91.0%) were below 10 kb. The length of 5S rDNA was revealed as 964 kb in G. herbaceum whereas, in G. arboreum, it was approximately three times longer (3.1 Mb). A fiber-FISH based immunofluorescence method was also described to assay the DNA methylation. Using this technique, we revealed that both telomere and 5S rDNA were methylated at different levels. In addition, we developed a BAC molecule-based fiber-FISH technique. Using this technique, we can precisely map BAC clones on each other and evaluated the size and location of overlapped regions. The development and application of fiber-FISH technique will facilitate high-resolution physical mapping and further directed sequencing projects for cotton.

Construction of rice mini-chromosomes by telomere-mediated chromosomal truncation

Telomere truncation has been shown to be an efficient technology for the creation of mini-chromosomes that can be used as artificial chromosome platforms for genetic engineering. Artificial chromosome-based genetic engineering is considered to be superior to the existing techniques of randomized gene integration by Agrobacterium or biolistic-mediated genetic transformation. It organizes multiple transgenes as a unique genetic linkage block for subsequent manipulations in breeding. Telomere truncation technology relies on three components: the telomere sequence that mediates chromosomal truncation, a selection marker that allows the selection of transgenic events, and a site-specific recombination system that can be used to accept future genes into the mini-chromosome by gene targeting. These elements are usually pre-assembled before transformation, a process that is both time and labor consuming. We found in this research that the three elements could be mixed to transform plant cells in a biolistic transformation, and produced efficient chromosomal truncations and mini-chromosomes in rice. This system will allow rapid construction of mini-chromosomes with a flexible selection of resistant markers, site-specific recombination systems and other desirable elements. In addition, a rice telotrisomic line was used as the starting material for chromosomal truncations. Mini-chromosomes from the truncations of both the telocentric chromosome and other chromosomes were recovered. The mini-chromosomes remained stable during 2 years of subculture. The construction of mini-chromosomes in rice, an economically important crop, will provide a platform for future artificial chromosome-based genetic engineering of rice for stacking multiple genes.

Preparations of meiotic pachytene chromosomes and extended DNA fibers from cotton suitable for fluorescence in situ hybridization

Fluorescence in situ hybridization (FISH) has become one of the most important techniques applied in plant molecular cytogenetics. However, the application of this technique in cotton has lagged behind because of difficulties in chromosome preparation. The focus of this article was FISH performed not only on cotton pachytene chromosomes, but also on cotton extended DNA fibers. The cotton pollen mother cells (PMCs) instead of buds or anthers were directly digested in enzyme to completely breakdown the cell wall. Before the routine acetic acid treatment, PMCs were incubated in acetic acid and enzyme mixture to remove the cytoplasm and clear the background. The method of ice-cold Carnoy's solution spreading chromosome was adopted instead of nitrogen removed method to avoid chromosomes losing and fully stretch chromosome. With the above-improved steps, the high-quality well-differentiated pachytene chromosomes with clear background were obtained. FISH results demonstrated that a mature protocol of cotton pachytene chromosomes preparation was presented. Intact and no debris cotton nuclei were obtained by chopping from etiolation cotyledons instead of the conventional liquid nitrogen grinding method. After incubating the nuclei with nucleus lysis buffer on slide, the parallel and clear background DNA fibers were acquired along the slide. This method overcomes the twist, accumulation and fracture of DNA fibers compared with other methods. The entire process of DNA fibers preparation requires only 30 min, in contrast, it takes 3 h with routine nitrogen grinding method. The poisonous mercaptoethanol in nucleus lysis buffer is replaced by nonpoisonous dithiothreitol. PVP40 in nucleus isolation buffer is used to prevent oxidation. The probability of success in isolating nuclei for DNA fiber preparation is almost 100% tested with this method in cotton. So a rapid, safe, and efficient method for the preparation of cotton extended DNA fibers suitable for FISH was established.

Microevolution of Cryptococcus neoformans driven by massive tandem gene amplification

The subtelomeric regions of organisms ranging from protists to fungi undergo a much higher rate of rearrangement than is observed in the rest of the genome. While characterizing these ~40-kb regions of the human fungal pathogen Cryptococcus neoformans, we have identified a recent gene amplification event near the right telomere of chromosome 3 that involves a gene encoding an arsenite efflux transporter (ARR3). The 3,177-bp amplicon exists in a tandem array of 2-15 copies and is present exclusively in strains with the C. neoformans var. grubii subclade VNI A5 MLST profile. Strains bearing the amplification display dramatically enhanced resistance to arsenite that correlates with the copy number of the repeat; the origin of increased resistance was verified as transport-related by functional complementation of an arsenite transporter mutant of Saccharomyces cerevisiae. Subsequent experimental evolution in the presence of increasing concentrations of arsenite yielded highly resistant strains with the ARR3 amplicon further amplified to over 50 copies, accounting for up to ~1% of the whole genome and making the copy number of this repeat as high as that seen for the ribosomal DNA. The example described here therefore represents a rare evolutionary intermediate-an array that is currently in a state of dynamic flux, in dramatic contrast to relatively common, static relics of past tandem duplications that are unable to further amplify due to nucleotide divergence. Beyond identifying and engineering fungal isolates that are highly resistant to arsenite and describing the first reported instance of microevolution via massive gene amplification in C. neoformans, these results suggest that adaptation through gene amplification may be an important mechanism that C. neoformans employs in response to environmental stresses, perhaps including those encountered during infection. More importantly, the ARR3 array will serve as an ideal model for further molecular genetic analyses of how tandem gene duplications arise and expand.

Organization and evolution of subtelomeric satellite repeats in the potato genome

Subtelomeric domains immediately adjacent to telomeres represent one of the most dynamic and rapidly evolving regions in eukaryotic genomes. A common feature associated with subtelomeric regions in different eukaryotes is the presence of long arrays of tandemly repeated satellite sequences. However, studies on molecular organization and evolution of subtelomeric repeats are rare. We isolated two subtelomeric repeats, CL14 and CL34, from potato (Solanum tuberosum). The CL14 and CL34 repeats are organized as independent long arrays, up to 1-3 Mb, of 182 bp and 339 bp monomers, respectively. The CL14 and CL34 repeat arrays are directly connected with the telomeric repeats at some chromosomal ends. The CL14 repeat was detected at the subtelomeric regions among highly diverged Solanum species, including tomato (Solanum lycopersicum). In contrast, CL34 was only found in potato and its closely related species. Interestingly, the CL34 repeat array was always proximal to the telomeres when both CL14 and CL34 were found at the same chromosomal end. In addition, the CL34 repeat family showed more sequence variability among monomers compared with the CL14 repeat family. We conclude that the CL34 repeat family emerged recently from the subtelomeric regions of potato chromosomes and is rapidly evolving. These results provide further evidence that subtelomeric domains are among the most dynamic regions in eukaryotic genomes.

Analysis of BAC clones containing homologous sequences on the end of the Xq arm and on chromosome 7 in the dioecious plant Silene latifolia

Silene latifolia is a model dioecious plant with morphologically distinguishable XY sex chromosomes. The end of the Xq arm is quite different from that of the Yp arm, although both are located at opposite ends of their respective chromosomes relative to a pseudo-autosomal region. The Xq arm does not seem to originate from the same autosome as the Yp arm. Bacterial artificial chromosome clone #15B12 has an insert containing a 130-kb stretch in which a 313-bp satellite DNA is repeated 420 times. PCR with a single primer revealed that this 130-kb stretch consists of three reversals of the orientation of the satellite DNA. A non-long terminal repeat retroelement and two sequences that share homology with an Oryza sativa RING zinc finger and a putative Arabidopsis thaliana protein, respectively, were found in the sequences that flank the satellite DNA. Fluorescence in situ hybridization carried out using this low-copy region of #15B12 as a probe confirmed that these sequences originated from the X chromosome and that homologous sequences exist at the end of chromosome 7. The region distal to DD44X on the Xq arm is postulated to have recombined with a region containing satellite DNA on chromosome 7 during the process of sex chromosome evolution.

Systematic overrepresentation of DNA termini and underrepresentation of subterminal regions among sequencing templates prepared from hydrodynamically sheared linear DNA molecules

Background

Analysis of fungal genome sequence assemblies reveals that telomeres are poorly represented even though telomeric reads tend to be superabundant. We surmised that the problem might lie in the DNA shearing conditions used to create clone libraries for genome sequencing.

Results

A shotgun strategy was used to sequence and assemble circular and linear cosmid DNAs sheared using conditions typical for a genome project. The DNA sheared in circular form assembled into a single sequence contig. However, the linearized cosmid produced an incomplete assembly because the two DNA termini, though greatly overrepresented in the clone library used for sequencing, were separated from neighboring sequences by gaps of ~1.4 and 1.8 kb. These gap sizes were reduced, but not eliminated, by shearing the linear cosmid into smaller fragments. Mapping of shearing breakpoints revealed a paucity of breaks in the subterminal regions of the linearized cosmid and also near chromosome ends of the fungus Neurospora crassa.

Conclusion

Together, our data indicate that the ends of linear DNA molecules are recalcitrant to hydrodynamic shearing. We propose that this causes DNA termini to be overrepresented in the resulting fragment population but ultimately prevents their incorporation into sequence assemblies.

Recent advances in rice genome and chromosome structure research by fluorescence in situ hybridization (FISH)

Fluorescence in situ hybridization (FISH) is an effective method for the physical mapping of genes and repetitive DNA sequences on chromosomes. Physical mapping of unique nucleotide sequences on specific rice chromosome regions was performed using a combination of chromosome identification and highly sensitive FISH. Increases in the detection sensitivity of smaller DNA sequences and improvements in spatial resolution have ushered in a new phase in FISH technology. Thus, it is now possible to perform in situ hybridization on somatic chromosomes, pachytene chromosomes, and even on extended DNA fibers (EDFs). Pachytene-FISH allows the integration of genetic linkage maps and quantitative chromosome maps. Visualization methods using FISH can reveal the spatial organization of the centromere, heterochromatin/euchromatin, and the terminal structures of rice chromosomes. Furthermore, EDF-FISH and the DNA combing technique can resolve a spatial distance of 1 kb between adjacent DNA sequences, and the detection of even a 300-bp target is now feasible. The copy numbers of various repetitive sequences and the sizes of various DNA molecules were quantitatively measured using the molecular combing technique. This review describes the significance of these advances in molecular cytology in rice and discusses future applications in plant studies using visualization techniques.

Isolation and sequence analysis of the wheat B genome subtelomeric DNA

Background

Telomeric and subtelomeric regions are essential for genome stability and regular chromosome replication. In this work, we have characterized the wheat BAC (bacterial artificial chromosome) clones containing Spelt1 and Spelt52 sequences, which belong to the subtelomeric repeats of the B/G genomes of wheats and Aegilops species from the section Sitopsis.

Results

The BAC library from Triticum aestivum cv. Renan was screened using Spelt1 and Spelt52 as probes. Nine positive clones were isolated; of them, clone 2050O8 was localized mainly to the distal parts of wheat chromosomes by in situ hybridization. The distribution of the other clones indicated the presence of different types of repetitive sequences in BACs. Use of different approaches allowed us to prove that seven of the nine isolated clones belonged to the subtelomeric chromosomal regions. Clone 2050O8 was sequenced and its sequence of 119 737 bp was annotated. It is composed of 33% transposable elements (TEs), 8.2% Spelt52 (namely, the subfamily Spelt52.2) and five non-TE-related genes. DNA transposons are predominant, making up 24.6% of the entire BAC clone, whereas retroelements account for 8.4% of the clone length. The full-length CACTA transposon Caspar covers 11 666 bp, encoding a transposase and CTG-2 proteins, and this transposon accounts for 40% of the DNA transposons. The in situ hybridization data for 2050O8 derived subclones in combination with the BLAST search against wheat mapped ESTs (expressed sequence tags) suggest that clone 2050O8 is located in the terminal bin 4BL-10 (0.95-1.0). Additionally, four of the predicted 2050O8 genes showed significant homology to four putative orthologous rice genes in the distal part of rice chromosome 3S and confirm the synteny to wheat 4BL.

Conclusion

Satellite DNA sequences from the subtelomeric regions of diploid wheat progenitor can be used for selecting the BAC clones from the corresponding regions of hexaploid wheat chromosomes. It has been demonstrated for the first time that Spelt52 sequences were involved in the evolution of terminal regions of common wheat chromosomes. Our research provides new insights into the microcollinearity in the terminal regions of wheat chromosomes 4BL and rice chromosome 3S.

Characterization of a tandemly repeated subtelomeric sequence with inverted telomere repeats in maize

In this study, two complementary telomere primers were applied to a single-primer PCR. A clear amplification band was obtained with one primer, while a smear pattern was seen with the other primer. Sequence analysis of the isolated clones from this specific amplification band revealed that a 412 bp clone designated as MTAS1 shared high homology with a reported subtelomeric sequence (382 bp) from maize (Zea mays L.), which indicated that this clone was possibly present at subtelomeric regions. The clone MTAS1 displayed a novel structural feature flanked by the forward and inverted telomere repeats. Southern hybridization revealed a ladder of hybridization bands, suggesting that MTAS1 was a tandemly repeated sequence. Fluorescence in situ hybridization results showed that the strong MTAS1 signals were present at the ends of short arms of several long chromosomes, confirming that MTAS1 was a subtelomeric sequence and the high brightness of signals further indicated this cloned sequence was a highly and tandemly repetitive sequence in maize. Fluorescence in situ hybridization with telomeric DNA and MTAS1 as probes on metaphase chromosomes and extended genomic DNA fibers showed that hybridization signals of this clone located adjacent to or overlapped with signals of telomere tandem repeats distributed heterogeneously in subtelomeric regions of several chromosomes and even exhibited differences in two subtelomeres of a single chromosome.

Analysis of Carica papaya Telomeres and Telomere-Associated Proteins: Insights into the Evolution of Telomere Maintenance in Brassicales

Telomeres are terminal regions of linear eukaryotic chromosomes that are critical for genome stability and continued cell proliferation. The draft assembly of the papaya genome provides an opportunity to analyze and compare the evolution of telomeric DNA sequence composition and telomere maintenance machinery in this and other organisms of the Brassicales Order, which includes Arabidopsis. Here we investigate telomere size and sequence variation at papaya chromosome ends. As with most other plant species, papaya telomeres consist of TTTAGGG repeats. However, in contrast to members of the closely related Brassicaceae family, telomeres in papaya are ~10-fold longer. Sequence analysis reveals that many centromereproximal telomere repeats in papaya harbor nucleotide substitutions and insertions of Gs and Ts. In contrast, we found very few N-to-C substitutions, and even fewer instances of nucleotide deletion, suggesting that a six-nucleotide telomere repeat is not well tolerated. The papaya genome encodes single-copy sequence homologues of several genes involved in telomere maintenance and chromosome end protection, including the Telomerase Reverse Transcriptase (TERT) and Protection Of Telomeres (POT1). Notably, unlike Arabidopsis, which encodes six Telomere Repeat binding Factor-like (TRFL) proteins that bind double-stranded telomere DNA, papaya appears to encode only two such proteins. Thus, the more streamlined genome of papaya will provide an excellent resource for comparative and functional analysis of telomeres in plants.

The subtelomere of Oryza sativa chromosome 3 short arm as a hot bed of new gene origination in rice

Despite general observations of non-random genomic distribution of new genes, it is unclear whether or not new genes preferentially occur in certain genomic regions driven by related molecular mechanisms. Using 1.5 Mb of genomic sequences from short arms of chromosome 3 of Oryza glaberrima and O. punctata, we conducted a comparative genomic analysis with the reference O. sativa ssp. japonica genome. We identified a 60-kb segment located in the middle of the subtelomeric region of chromosome 3, which is unique to the species O. sativa. The region contained gene duplicates that occurred in Asian cultivated rice species that diverged from the ancestor of Asian and African cultivated rice one million years ago (MYA). For the 12 genes and one complete retrotransposon identified in this segment in O. sativa ssp. japonica, we searched for their parental genes. The high similarity between duplicated paralogs further supports the recent origination of these genes. We found that this segment was recently generated through multiple independent gene recombination and transposon insertion events. Among the 12 genes, we found that five had chimeric gene structures derived from multiple parental genes. Nine out of the 12 new genes seem to be functional, as suggested by Ka/Ks analysis and the presence of cDNA and/or MPSS data. Furthermore, for the eight transcribed genes, at least two genes could be classified as defense or stress response-related genes. Given these findings, and the fact that subtelomeres are associated with high rates of recombination and transcription, it is likely that subtelomeres may facilitate gene recombination and transposon insertions and serve as hot spots for new gene origination in rice genomes.

Cytological characterization of the tandem repetitive sequences and their methylation status in the Antirrhinum majus genome

Tandem repetitive sequences are DNA motifs common in the genomes of eukaryotic species and are often embedded in heterochromatic regions. In most eukaryotes, ribosomal genes, as well as centromeres and telomeres or subtelomeres, are associated with abundant tandem arrays of repetitive sequences and typically represent the final barriers to completion of whole-genome sequencing. The nature of these repeats makes it difficult to estimate their actual sizes. In this study, combining the two cytological techniques DNA fiber-FISH and pachytene chromosome FISH allowed us to characterize the tandem repeats distributed genome wide in Antirrhinum majus and identify four types of tandem repeats, 45S rDNA, 5S rDNA, CentA1, and CentA2, representing the major tandem repetitive components, which were estimated to have a total length of 18.50 Mb and account for 3.59% of the A. majus genome. FISH examination revealed that all the tandem repeats correspond to heterochromatic knobs along the pachytene chromosomes. Moreover, the methylation status of the tandem repeats was investigated in both somatic cells and pollen mother cells from anther tissues using an antibody against 5-methylcytosine combined with sequential FISH analyses. Our results showed that these repeats were hypomethylated in anther tissues, especially in the pollen mother cells at pachytene stage.

Characterization of chromosome ends on the basis of the structure of TrsA subtelomeric repeats in rice (Oryza sativa L.)

Subtelomeres contain species-specific repetitive sequences. We characterized rice chromosome ends on the basis of the structure of TrsA, a subtelomeric repetitive sequence of rice. Among the 24 chromosome arms, TrsA was arrayed in tandem on the ends of five: 5L, 6S, 8L, 9L, and 12L. TrsA sequences were arranged in discrete clusters of 3-106 copies in a chromosome-specific manner, instead of being distributed uniformly throughout the subtelomeric regions. The clusters were located at the distal-most end of the sequenced region in 5L, 6S, 8L, and 9L, but in 12L expressed genes were present distal to the clusters. Thus, rice subtelomeres are composed of discrete clusters of a TrsA-rich region and a gene-rich region with high transcriptional activity. Intra-chromosomal duplications have resulted in a striking degree of variation in the number and distribution of TrsAs, suggesting that the areas near the ends of the chromosomes are dynamic and variable.

Gene organization in rice revealed by full-length cDNA mapping and gene expression analysis through microarray

Rice (Oryza sativa L.) is a model organism for the functional genomics of monocotyledonous plants since the genome size is considerably smaller than those of other monocotyledonous plants. Although highly accurate genome sequences of indica and japonica rice are available, additional resources such as full-length complementary DNA (FL-cDNA) sequences are also indispensable for comprehensive analyses of gene structure and function. We cross-referenced 28.5K individual loci in the rice genome defined by mapping of 578K FL-cDNA clones with the 56K loci predicted in the TIGR genome assembly. Based on the annotation status and the presence of corresponding cDNA clones, genes were classified into 23K annotated expressed (AE) genes, 33K annotated non-expressed (ANE) genes, and 5.5K non-annotated expressed (NAE) genes. We developed a 60mer oligo-array for analysis of gene expression from each locus. Analysis of gene structures and expression levels revealed that the general features of gene structure and expression of NAE and ANE genes were considerably different from those of AE genes. The results also suggested that the cloning efficiency of rice FL-cDNA is associated with the transcription activity of the corresponding genetic locus, although other factors may also have an effect. Comparison of the coverage of FL-cDNA among gene families suggested that FL-cDNA from genes encoding rice- or eukaryote-specific domains, and those involved in regulatory functions were difficult to produce in bacterial cells. Collectively, these results indicate that rice genes can be divided into distinct groups based on transcription activity and gene structure, and that the coverage bias of FL-cDNA clones exists due to the incompatibility of certain eukaryotic genes in bacteria.

Recombination: an underappreciated factor in the evolution of plant genomes

Our knowledge of recombination rates and patterns in plants is far from being comprehensive. However, compelling evidence indicates a central role for recombination, through its influences on mutation and selection, in the evolution of plant genomes. Furthermore, recombination seems to be generally higher and more variable in plants than in animals, which could be one of the primary reasons for differences in genome lability between these two kingdoms. Much additional study of recombination in plants is needed to investigate these ideas further.

Sequence comparison of distal and proximal ribosomal DNA arrays in rice (Oryza sativa L.) chromosome 9S and analysis of their flanking regions

Rice (Oryza sativa ssp. japonica cv. Nipponbare) harbors a ribosomal RNA gene (rDNA) cluster in the nucleolar-organizing region at the telomeric end of the short arm of chromosome 9. We isolated and sequenced two genomic clones carrying rice rDNA fragments from this region. The rice rDNA repeat units could be classified into three types based on length, which ranged from 7,928 to 8,934 bp. This variation was due to polymorphism in the number of 254-bp subrepeats in the intergenic spacer (IGS). Polymerase chain reaction (PCR) analysis suggested that the rDNA units in rice vary widely in length and that the copy number of the subrepeats in the IGS ranges from 1 to 12 in the rice genome. PCR and Southern blot analyses showed that most rDNA units have three intact and one truncated copies of the subrepeats in the IGS, and distal (telomere-side) rDNA units have more subrepeats than do proximal (centromere-side) ones. Both genomic clones we studied contained rDNA-flanking DNA sequences of either telomeric repeats (5'-TTTAGGG-3') or a chromosome-specific region, suggesting that they were derived from the distal or proximal end, respectively, of the rDNA cluster. A similarity search indicated that retrotransposons appeared more frequently in a 500-kb portion of the proximal rDNA-flanking region than in other subtelomeric regions or sequenced regions of the genome. This study reveals the repetitive nature of the telomeric end of the short arm of chromosome 9, which consists of telomeric repeats, an rDNA array, and a retrotransposon-rich chromosomal region.