Distinct structural variants and repeat landscape shape the genomes of the ancient grapes Aglianico and Falanghina

Mounting evidence recognizes structural variations (SVs) and repetitive DNA sequences as crucial players in shaping the existing grape phenotypic diversity at intra- and inter-species levels. To deepen our understanding on the abundance, diversity, and distribution of SVs and repetitive DNAs, including transposable elements (TEs) and tandemly repeated satellite DNA (satDNAs), we re-sequenced the genomes of the ancient grapes Aglianico and Falanghina. The analysis of large copy number variants (CNVs) detected candidate polymorphic genes that are involved in the enological features of these varieties. In a comparative analysis of Aglianico and Falanghina sequences with 21 publicly available genomes of cultivated grapes, we provided a genome-wide annotation of grape TEs at the lineage level. We disclosed that at least two main clusters of grape cultivars could be identified based on the TEs content. Multiple TEs families appeared either significantly enriched or depleted. In addition, in silico and cytological analyses provided evidence for a diverse chromosomal distribution of several satellite repeats between Aglianico, Falanghina, and other grapes. Overall, our data further improved our understanding of the intricate grape diversity held by two Italian traditional varieties, unveiling a pool of unique candidate genes never so far exploited in breeding for improved fruit quality.

On the Way to the Technological Development of Newly Selected Non-Saccharomyces Yeasts Selected as Innovative Biocontrol Agents in Table Grapes

Post-harvest decay of fresh table grapes causes considerable annual production losses. The main fungal agents of decay both in pre- and post-harvest are B. cinerea, Penicillium spp., Aspergillus spp., Alternaria spp., and Cladosporium spp. To date, the use of agrochemicals and SO2 are the main methods to control grape molds in pre- and postharvest, respectively. Significant improvements, however, have already been made in to apply innovative and more environmentally sustainable control strategies, such as Biological Control Agents (BCAs), which can reduce disease severity in both pre- and post-harvest. In this study, 31 new non-Saccharomyces yeast strains, isolated from berries of native Apulian table grape genotypes, were tested for their in vivo effectiveness against grey mold of table grapes, resulting in two St. bacillaris ('N22_I1' and 'S13_I3'), one S. diversa ('N22_I3'), one A. pullulans ('OLB_9.1_VL') and one H. uvarum ('OLB_9.1_BR') yeast strains that were marked as efficient and good BCAs. Their mechanisms of action were characterized through in vitro assays, and additional characteristics were evaluated to assess the economic feasibility and viability for future technological employment. Their effectiveness was tested by reducing the working concentration, their antagonistic effect on a wide range of fungal pathogens, their ability to survive in formulations with long shelf life, and their safety to human health.

Transcriptomic and methylation analysis of susceptible and tolerant grapevine genotypes following Plasmopara viticola infection

Downy mildew, caused by the biotrophic oomycete Plasmopara viticola, is one of the most economically significant grapevine diseases worldwide. Current strategies to cope with this threat rely on the massive use of chemical compounds during each cultivation season. The economic costs and negative environmental impact associated with these applications increased the urge to search for sustainable strategies of disease control. Improved knowledge of plant mechanisms to counteract pathogen infection may allow the development of alternative strategies for plant protection. Epigenetic regulation, in particular DNA methylation, is emerging as a key factor in the context of plant-pathogen interactions associated with the expression modulation of defence genes. To improve our understanding of the genetic and epigenetic mechanisms underpinning grapevine response to P. viticola, we studied the modulation of both 5-mC methylation and gene expression at 6 and 24 h post-infection (hpi). Leaves of two table grape genotypes (Vitis vinifera), selected by breeding activities for their contrasting level of susceptibility to the pathogen, were analysed. Following pathogen infection, we found variations in the 5-mC methylation level and the gene expression profile. The results indicate a genotype-specific response to pathogen infection. The tolerant genotype (N23/018) at 6 hpi exhibits a lower methylation level compared to the susceptible one (N20/020), and it shows an early modulation (at 6 hpi) of defence and epigenetic-related genes during P. viticola infection. These data suggest that the timing of response is an important mechanism to efficiently counteract the pathogen attack.

VviAGL11 self-regulates and targets hormone- and secondary metabolism-related genes during seed development

VviAGL11, the Arabidopsis SEEDSTICK homolog, has been proposed to have a causative role in grapevine stenospermocarpy. An association between a mutation in the coding sequence (CDS) and the seedless phenotype was reported, however, no working mechanisms have been demonstrated yet. We performed a deep investigation of the full VviAGL11 gene sequence in a collection of grapevine varieties belonging to several seedlessness classes that revealed three different promoter-CDS combinations. By investigating the expression of the three VviAGL11 alleles, and by evaluating their ability to activate the promoter region, we observed that VviAGL11 self-activates in a specific promoter-CDS combination manner. Furthermore, by transcriptomic analyses on ovule and developing seeds in seeded and seedless varieties and co-expression approaches, candidate VviAGL11 targets were identified and further validated through luciferase assay and in situ hybridization. We demonstrated that VviAGL11 Wild Type CDS activates Methyl jasmonate esterase and Indole-3-acetate beta-glucosyltransferase, both involved in hormone signaling and Isoflavone reductase, involved in secondary metabolism. The dominant-negative effect of the mutated CDS was also functionally ectopically validated in target induction. VviAGL11 was shown to co-localize with its targets in the outer seed coat integument, supporting its direct involvement in seed development, possibly by orchestrating the crosstalk among MeJA, auxin, and isoflavonoids synthesis. In conclusion, the VviAGL11 expression level depends on the promoter-CDS allelic combination, and this will likely affect its ability to activate important triggers of the seed coat development. The dominant-negative effect of the mutated VviAGL11 CDS on the target genes activation was molecularly validated. A new regulatory mechanism correlating VviAGL11 haplotype assortment and seedlessness class in grapevine is proposed.

Novel and emerging biotechnological crop protection approaches

Traditional breeding or genetically modified organisms (GMOs) have for a long time been the sole approaches to effectively cope with biotic and abiotic stresses and implement the quality traits of crops. However, emerging diseases as well as unpredictable climate changes affecting agriculture over the entire globe force scientists to find alternative solutions required to quickly overcome seasonal crises. In this review, we first focus on cisgenesis and genome editing as challenging biotechnological approaches for breeding crops more tolerant to biotic and abiotic stresses. In addition, we take into consideration a toolbox of new techniques based on applications of RNA interference and epigenome modifications, which can be adopted for improving plant resilience. Recent advances in these biotechnological applications are mainly reported for non-model plants and woody crops in particular. Indeed, the characterization of RNAi machinery in plants is fundamental to transform available information into biologically or biotechnologically applicable knowledge. Finally, here we discuss how these innovative and environmentally friendly techniques combined with traditional breeding can sustain a modern agriculture and be of potential contribution to climate change mitigation.

Correction to: SNP genotyping elucidates the genetic diversity of Magna Graecia grapevine germplasm and its historical origin and dissemination

An amendment to this paper has been published and can be accessed via the original article.

Microarray data and pathway analyses of peripheral blood mononuclear cells from healthy subjects after a three weeks grape-rich diet

Using Human Gene Expression Microarrays (Agilent) technologies, we investigated changes of the level of gene expression in peripheral blood mononuclear cells of healthy subjects after 21 days of fresh table grape-rich diet and after an additional 28-day washout. Several hundreds of genes were differentially expressed after grape intake or after washout. The functional analysis of these genes detected significant changes in key processes such as inflammation and immunity, thrombosis, DNA and protein repair, autophagy and mitochondrial biogenesis. Moreover, fresh grape intake was found to influence the expression of many long non-coding RNA genes. The data can be valuable for researchers interested in nutrigenetics and nutrigenomics studies and are related to the research article “Gene expression signature induced by grape intake in healthy subjects reveals wide-spread beneficial effects on PBMCs” [1].

FT-NIR Analysis of Intact Table Grape Berries to Understand Consumer Preference Driving Factors

Fourier-transform near infrared spectroscopy (FT-NIR) is a technique used in the compositional and sensory analysis of foodstuffs. In this work, we have measured the main maturity parameters for grape (sugars and acids) using hundreds of intact berry samples to build models for the prediction of these parameters from berries of two very different varieties: “Victoria” and “Autumn Royal”. Together with the chemical composition in terms of sugar and acidic content, we have carried out a sensory analysis on single berries. Employing the models built for sugars and acids it was possible to learn the sweetness and acidity of each berry before the destructive sensory analysis. The direct correlation of sensory data with FT-NIR spectra is difficult; therefore, spectral data were exported from the spectrometer built-in software and analyzed with R software using a statistical analysis technique (Spearman correlation) which allowed the correlation of berry appreciation data with specific wavelengths that were then related to sugar and acidic content. In this article, we show how it is possible to carry out the analysis of single berries to obtain data on chemical composition parameters and consumer appreciation with a fast, simple, and non-destructive technique with a clear advantage for producers and consumers.

Transcriptomic and genomic structural variation analyses on grape cultivars reveal new insights into the genotype-dependent responses to water stress

Grapevine (Vitis vinifera L.) is importantly cultivated worldwide for table grape and wine production. Its cultivation requires irrigation supply, especially in arid and semiarid areas. Water deficiency can affect berry and wine quality mostly depending on the extent of plant perceived stress, which is a cultivar-specific trait. We tested the physiological and molecular responses to water deficiency of two table grape cultivars, Italia and Autumn royal, and we highlighted their different adaptation. Microarray analyses revealed that Autumn royal reacts involving only 29 genes, related to plant stress response and ABA/hormone signal transduction, to modulate the response to water deficit. Instead, cultivar Italia orchestrates a very broad response (we found 1037 differentially expressed genes) that modifies the cell wall organization, carbohydrate metabolism, response to reactive oxygen species, hormones and osmotic stress. For the first time, we integrated transcriptomic data with cultivar-specific genomics and found that ABA-perception and -signalling are key factors mediating the varietal-specific behaviour of the early response to drought. We were thus able to isolate candidate genes for the genotype-dependent response to drought. These insights will allow the identification of reliable plant stress indicators and the definition of sustainable cultivar-specific protocols for water management.

SNP genotyping elucidates the genetic diversity of Magna Graecia grapevine germplasm and its historical origin and dissemination

BACKGROUND

Magna Graecia is the ancient name for the modern geopolitical region of South Italy extensively populated by Greek colonizers, shown by archeological and historical evidence to be the oldest wine growing region of Italy, crucial for the spread of specialized viticulture around Mediterranean shores. Here, the genetic diversity of Magna Graecia grape germplasm was assessed and its role in grapevine propagation around the Mediterranean basin was underlined.

RESULTS

A large collection of grapevines from Magna Graecia was compared with germplasm from Georgia to the Iberian Peninsula using the 18 K SNP array. A high level of genetic diversity of the analyzed germplasm was determined; clustering, structure analysis and DAPC (Discriminant Analysis of Principal Components) highlighted the genetic relationships among genotypes from South Italy and the Eastern Mediterranean (Greece). Gene flow from east (Georgia) to west (Iberian Peninsula) was identified throughout the large number of detected admixed samples. Pedigree analysis showed a complex and well-structured network of first degree relationships, where the cultivars from Magna Graecia were mainly involved.

CONCLUSIONS

This study provided evidence that Magna Graecia germplasm was shaped by historical events that occurred in the area due to the robust link between South Italian and Greek genotypes, as well as, by the availability of different thermal resources for cultivars growing in such different winegrowing areas. The uniqueness of this ampelographic platform was mainly an outcome of complex natural or human-driven crosses involving elite cultivars.

Inter-varietal structural variation in grapevine genomes

Grapevine (Vitis vinifera L.) is one of the world's most important crop plants, which is of large economic value for fruit and wine production. There is much interest in identifying genomic variations and their functional effects on inter-varietal, phenotypic differences. Using an approach developed for the analysis of human and mammalian genomes, which combines high-throughput sequencing, array comparative genomic hybridization, fluorescent in situ hybridization and quantitative PCR, we created an inter-varietal atlas of structural variations and single nucleotide variants (SNVs) for the grapevine genome analyzing four economically and genetically relevant table grapevine varieties. We found 4.8 million SNVs and detected 8% of the grapevine genome to be affected by genomic variations. We identified more than 700 copy number variation (CNV) regions and more than 2000 genes subjected to CNV as potential candidates for phenotypic differences between varieties.

The molecular characterization by SSRs reveals a new South Italian kinship and the origin of the cultivar Uva di Troia

Vitis vinifera L. varieties were spread through cuttings following historic migrations of people, trades, or after biological crises due to pests outbreaks. Some today’s varieties could be more than a 1000 years old and, although over the centuries these varieties generated most of the remaining cultivars, their origin could be impossible to track back. The Italian grapevine biodiversity is one of most important, most likely due to its strategic position in the middle of the Mediterranean sea. Unravelling of its structure is challenging because of its complexity and the lack of historical documentation. In this paper molecular data are compared with historical documentations. Simple Sequence Repeats fingerprinting are molecular markers best suited to investigate genetic relationships and identify pedigrees. South-Italian germplasm was studied with 54 nuclear microsatellites. A family was identified, consisting of two parents and three siblings and further genetically characterized with six nuclear and five chloroplast microsatellites and described with ampelographic and phylometric analysis. Although these latter were not informative for the kinship identification. The common Bombino bianco was the female parent and the previously unknown Uva rosa antica was the male parent. Bombino nero, Impigno and the popular Uva di Troia, all typical of the south-east Italy, were the offspring. Further research showed that the Uva rosa antica was a synonym of Quagliano and Bouteillan noir, both minor varieties. Quagliano was considered to be autochthonous of some alpine valleys in the north-west of Italy and Bouteillan noir is a neglected variety of Vancluse in France. This finding uncovers the intricate nature of Italian grape cultivars, considered peculiar of an area, but possibly being the remains of ancient latin founding varieties. Consequently, intriguing new hypotheses are discussed and some conclusions are drawn, based on the peculiar geographical origin of the parents, on the distribution of the offspring, on the chance of a single, and perhaps intentional, crossing event.

Morphological variability in leaves and molecular characterization of novel table grape candidate cultivars (Vitis vinifera L.)

The present work report the characterization of twenty-one table grapes candidate cultivars plus five registered ones included as reference, by means of 47 ampelographic traits, 23 ampelometric measurements and six microsatellite loci. The final goal of the research was to analyse the possibility of reducing the number of morphological and molecular tools required for a precise and effective description of a grape genotype or cultivar. This would be of great help for future biodiversity description on a larger sample of more than 300 table grapes accessions today grown at the 'Consiglio per la Ricerca e la sperimentazione in Agricoltura (C.R.A.)-Unità di ricerca per l'uva da tavola e la vitivinicoltura in ambiente mediterraneo (Bari-Italy)'. OIV ampelographic traits showed a clear distinction among all twenty-six genotypes analysed, suggesting the relevant morphological variability investigated. Principal component analysis based on ampelometric traits revealed main veins ON(3), ON(4) and O(3)N(4); ratios between main veins; angles between main veins and of petiolar sinus, to be the most effective records in differentiating cultivars, for a total variation of 69.9 % described by the first three components. Molecular analysis based on six microsatellite loci was performed on all genotypes, providing a detailed molecular profile and a dendrogram of genetic similarity, in which all genotypes were clearly distinguishable. Finally, with the goal of using the minimum possible number of markers to differentiate genotypes, microsatellites VVMD5 and VVMD27 were selected to be sufficient to distinguish among all the candidate cultivars included in the analysis, representing a possible 'step by step' approach when a molecular characterization has to be undertaken on a large number of genotypes, by first testing few markers and increasing their number only if necessary.

Evolutionary history of linked D4Z4 and Beta satellite clusters at the FSHD locus (4q35)

We performed a detailed genomic investigation of the chimpanzee locus syntenic to human chromosome 4q35.2, associated to the facioscapulohumeral dystrophy. Two contigs of approximately 150 kb and 200 kb were derived from PTR chromosomes 4q35 and 3p12, respectively: both regions showed a very similar sequence organization, including D4Z4 and Beta satellite linked clusters. Starting from these findings, we derived a hypothetical evolutionary history of human 4q35, 10q26 and 3p12 chromosome regions focusing on the D4Z4–Beta satellite linked organization. The D4Z4 unit showed an open reading frame (DUX4) at both PTR 4q35 and 3p12 regions; furthermore some subregions of the Beta satellite unit showed a high degree of conservation between chimpanzee and humans. In conclusion, this paper provides evidence that at the 4q subtelomere the linkage between D4Z4 and Beta satellite arrays is a feature that appeared late during evolution and is conserved between chimpanzee and humans.

Copy number variation of individual cattle genomes using next-generation sequencing

Copy number variations (CNVs) affect a wide range of phenotypic traits; however, CNVs in or near segmental duplication regions are often intractable. Using a read depth approach based on next-generation sequencing, we examined genome-wide copy number differences among five taurine (three Angus, one Holstein, and one Hereford) and one indicine (Nelore) cattle. Within mapped chromosomal sequence, we identified 1265 CNV regions comprising ~55.6-Mbp sequence--476 of which (~38%) have not previously been reported. We validated this sequence-based CNV call set with array comparative genomic hybridization (aCGH), quantitative PCR (qPCR), and fluorescent in situ hybridization (FISH), achieving a validation rate of 82% and a false positive rate of 8%. We further estimated absolute copy numbers for genomic segments and annotated genes in each individual. Surveys of the top 25 most variable genes revealed that the Nelore individual had the lowest copy numbers in 13 cases (~52%, χ(2) test; P-value <0.05). In contrast, genes related to pathogen- and parasite-resistance, such as CATHL4 and ULBP17, were highly duplicated in the Nelore individual relative to the taurine cattle, while genes involved in lipid transport and metabolism, including APOL3 and FABP2, were highly duplicated in the beef breeds. These CNV regions also harbor genes like BPIFA2A (BSP30A) and WC1, suggesting that some CNVs may be associated with breed-specific differences in adaptation, health, and production traits. By providing the first individualized cattle CNV and segmental duplication maps and genome-wide gene copy number estimates, we enable future CNV studies into highly duplicated regions in the cattle genome.

Real-time PCR for the detection of precise transgene copy number in durum wheat

Recent results obtained in various crops indicate that real-time PCR could be a powerful tool for the detection and characterization of transgene locus structures. The determination of transgenic locus number through real-time PCR overcomes the problems linked to phenotypic segregation analysis (i.e. lack of detectable expression even when the transgenes are present) and can analyse hundreds of samples in a day, making it an efficient method for estimating gene copy number. Despite these advantages, many authors speak of "estimating" copy number by real-time PCR, and this is because the detection of a precise number of transgene depends on how well real-time PCR performs.This study was conducted to determine transgene copy number in transgenic wheat lines and to investigate potential variability in sensitivity and resolution of real-time chemistry by TaqMan probes. We have applied real-time PCR to a set of four transgenic durum wheat lines previously obtained. A total of 24 experiments (three experiments for two genes in each transgenic line) were conducted and standard curves were obtained from serial dilutions of the plasmids containing the genes of interest. The correlation coefficients ranged from 0.95 to 0.97. By using TaqMan quantitative real-time PCR we were able to detect 1 to 41 copies of transgenes per haploid genome in the DNA of homozygous T4 transformants. Although a slight variability was observed among PCR experiments, in our study we found real-time PCR to be a fast, sensitive and reliable method for the detection of transgene copy number in durum wheat, and a useful adjunct to Southern blot and FISH analyses to detect the presence of transgenic DNA in plant material.

Genome-wide characterization of centromeric satellites from multiple mammalian genomes

Despite its importance in cell biology and evolution, the centromere has remained the final frontier in genome assembly and annotation due to its complex repeat structure. However, isolation and characterization of the centromeric repeats from newly sequenced species are necessary for a complete understanding of genome evolution and function. In recent years, various genomes have been sequenced, but the characterization of the corresponding centromeric DNA has lagged behind. Here, we present a computational method (RepeatNet) to systematically identify higher-order repeat structures from unassembled whole-genome shotgun sequence and test whether these sequence elements correspond to functional centromeric sequences. We analyzed genome datasets from six species of mammals representing the diversity of the mammalian lineage, namely, horse, dog, elephant, armadillo, opossum, and platypus. We define candidate monomer satellite repeats and demonstrate centromeric localization for five of the six genomes. Our analysis revealed the greatest diversity of centromeric sequences in horse and dog in contrast to elephant and armadillo, which showed high-centromeric sequence homogeneity. We could not isolate centromeric sequences within the platypus genome, suggesting that centromeres in platypus are not enriched in satellite DNA. Our method can be applied to the characterization of thousands of other vertebrate genomes anticipated for sequencing in the near future, providing an important tool for annotation of centromeres.

Analysis of copy number variations among diverse cattle breeds

Genomic structural variation is an important and abundant source of genetic and phenotypic variation. Here, we describe the first systematic and genome-wide analysis of copy number variations (CNVs) in modern domesticated cattle using array comparative genomic hybridization (array CGH), quantitative PCR (qPCR), and fluorescent in situ hybridization (FISH). The array CGH panel included 90 animals from 11 Bos taurus, three Bos indicus, and three composite breeds for beef, dairy, or dual purpose. We identified over 200 candidate CNV regions (CNVRs) in total and 177 within known chromosomes, which harbor or are adjacent to gains or losses. These 177 high-confidence CNVRs cover 28.1 megabases or approximately 1.07% of the genome. Over 50% of the CNVRs (89/177) were found in multiple animals or breeds and analysis revealed breed-specific frequency differences and reflected aspects of the known ancestry of these cattle breeds. Selected CNVs were further validated by independent methods using qPCR and FISH. Approximately 67% of the CNVRs (119/177) completely or partially span cattle genes and 61% of the CNVRs (108/177) directly overlap with segmental duplications. The CNVRs span about 400 annotated cattle genes that are significantly enriched for specific biological functions, such as immunity, lactation, reproduction, and rumination. Multiple gene families, including ULBP, have gone through ruminant lineage-specific gene amplification. We detected and confirmed marked differences in their CNV frequencies across diverse breeds, indicating that some cattle CNVs are likely to arise independently in breeds and contribute to breed differences. Our results provide a valuable resource beyond microsatellites and single nucleotide polymorphisms to explore the full dimension of genetic variability for future cattle genomic research.

New insights into centromere organization and evolution from the white-cheeked gibbon and marmoset

The evolutionary history of alpha-satellite DNA, the major component of primate centromeres, is hardly defined because of the difficulty in its sequence assembly and its rapid evolution when compared with most genomic sequences. By using several approaches, we have cloned, sequenced, and characterized alpha-satellite sequences from two species representing critical nodes in the primate phylogeny: the white-cheeked gibbon, a lesser ape, and marmoset, a New World monkey. Sequence analyses demonstrate that white-cheeked gibbon and marmoset alpha-satellite sequences are formed by units of approximately 171 and approximately 342 bp, respectively, and they both lack the high-order structure found in humans and great apes. Fluorescent in situ hybridization characterization shows a broad dispersal of alpha-satellite in the white-cheeked gibbon genome including centromeric, telomeric, and chromosomal interstitial localizations. On the other hand, centromeres in marmoset appear organized in highly divergent dimers roughly of 342 bp that show a similarity between monomers much lower than previously reported dimers, thus representing an ancient dimeric structure. All these data shed light on the evolution of the centromeric sequences in Primates. Our results suggest radical differences in the structure, organization, and evolution of alpha-satellite DNA among different primate species, supporting the notion that 1) all the centromeric sequence in Primates evolved by genomic amplification, unequal crossover, and sequence homogenization using a 171 bp monomer as the basic seeding unit and 2) centromeric function is linked to relatively short repeated elements, more than higher-order structure. Moreover, our data indicate that complex higher-order repeat structures are a peculiarity of the hominid lineage, showing the more complex organization in humans.

Evolutionary toggling of the MAPT 17q21.31 inversion region

Using comparative sequencing approaches, we investigated the evolutionary history of the European-enriched 17q21.31 MAPT inversion polymorphism. We present a detailed, BAC-based sequence assembly of the inverted human H2 haplotype and compare it to the sequence structure and genetic variation of the corresponding 1.5-Mb region for the noninverted H1 human haplotype and that of chimpanzee and orangutan. We found that inversion of the MAPT region is similarly polymorphic in other great ape species, and we present evidence that the inversions occurred independently in chimpanzees and humans. In humans, the inversion breakpoints correspond to core duplications with the LRRC37 gene family. Our analysis favors the H2 configuration and sequence haplotype as the likely great ape and human ancestral state, with inversion recurrences during primate evolution. We show that the H2 architecture has evolved more extensive sequence homology, perhaps explaining its tendency to undergo microdeletion associated with mental retardation in European populations.

Hominoid chromosomal rearrangements on 17q map to complex regions of segmental duplication

BACKGROUND

Chromosomal rearrangements, such as translocations and inversions, are recurrent phenomena during evolution, and both of them are involved in reproductive isolation and speciation. To better understand the molecular basis of chromosome rearrangements and their part in karyotype evolution, we have investigated the history of human chromosome 17 by comparative fluorescence in situ hybridization (FISH) and sequence analysis.

RESULTS

Human bacterial artificial chromosome/p1 artificial chromosome probes spanning the length of chromosome 17 were used in FISH experiments on great apes, Old World monkeys and New World monkeys to study the evolutionary history of this chromosome. We observed that the macaque marker order represents the ancestral organization. Human, chimpanzee and gorilla homologous chromosomes differ by a paracentric inversion that occurred specifically in the Homo sapiens/Pan troglodytes/Gorilla gorilla ancestor. Detailed analyses of the paracentric inversion revealed that the breakpoints mapped to two regions syntenic to human 17q12/21 and 17q23, both rich in segmental duplications.

CONCLUSION

Sequence analyses of the human and macaque organization suggest that the duplication events occurred in the catarrhine ancestor with the duplication blocks continuing to duplicate or undergo gene conversion during evolution of the hominoid lineage. We propose that the presence of these duplicons has mediated the inversion in the H. sapiens/P. troglodytes/G. gorilla ancestor. Recently, the same duplication blocks have been shown to be polymorphic in the human population and to be involved in triggering microdeletion and duplication in human. These results further support a model where genomic architecture has a direct role in both rearrangement involved in karyotype evolution and genomic instability in human.

Characterization and evolution of the novel gene family FAM90A in primates originated by multiple duplication and rearrangement events

Genomic plasticity of human chromosome 8p23.1 region is highly influenced by two groups of complex segmental duplications (SDs), termed REPD and REPP, that mediate different kinds of rearrangements. Part of the difficulty to explain the wide range of phenotypes associated with 8p23.1 rearrangements is that REPP and REPD are not yet well characterized, probably due to their polymorphic status. Here, we describe a novel primate-specific gene family, named FAM90A (family with sequence similarity 90), found within these SDs. According to the current human reference sequence assembly, the FAM90A family includes 24 members along 8p23.1 region plus a single member on chromosome 12p13.31, showing copy number variation (CNV) between individuals. These genes can be classified into subfamilies I and II, which differ in their upstream and 5'-untranslated region sequences, but both share the same open reading frame and are ubiquitously expressed. Sequence analysis and comparative fluorescence in situ hybridization studies showed that FAM90A subfamily II suffered a big expansion in the hominoid lineage, whereas subfamily I members were likely generated sometime around the divergence of orangutan and African great apes by a fusion process. In addition, the analysis of the Ka/Ks ratios provides evidence of functional constraint of some FAM90A genes in all species. The characterization of the FAM90A gene family contributes to a better understanding of the structural polymorphism of the human 8p23.1 region and constitutes a good example of how SDs, CNVs and rearrangements within themselves can promote the formation of new gene sequences with potential functional consequences.

Evolutionary history of chromosome 11 featuring four distinct centromere repositioning events in Catarrhini

Panels of BAC clones used in FISH experiments allow a detailed definition of chromosomal marker arrangement and orientation during evolution. This approach has disclosed the centromere repositioning phenomenon, consisting in the activation of a novel, fully functional centromere in an ectopic location, concomitant with the inactivation of the old centromere. In this study, appropriate panels of BAC clones were used to track the chromosome 11 evolutionary history in primates and nonprimate boreoeutherian mammals. Chromosome 11 synteny was found to be highly conserved in both primate and boreoeutherian mammalian ancestors. Amazingly, we detected four centromere repositioning events in primates (in Old World monkeys, in gibbons, in orangutans, and in the Homo-Pan-Gorilla (H-P-G) clade ancestor), and one in Equidae. Both H-P-G and Lar gibbon novel centromeres were flanked by large duplicons with high sequence similarity. Outgroup species analysis revealed that this duplicon was absent in phylogenetically more distant primates. The chromosome 11 ancestral centromere was probably located near the HSA11q telomere. The domain of this inactivated centromere, in humans, is almost devoid of segmental duplications. An inversion occurred in chromosome 11 in the common ancestor of H-P-G. A large duplicon, again absent in outgroup species, was found located adjacent to the inversion breakpoints. In Hominoidea, almost all the five largest duplicons of this chromosome appeared involved in significant evolutionary architectural changes.

Evolutionary Formation of New Centromeres in Macaque

A systematic fluorescence in situ hybridization comparison of macaque and human synteny organization disclosed five additional macaque evolutionary new centromeres (ENCs) for a total of nine ENCs. To understand the dynamics of ENC formation and progression, we compared the ENC of macaque chromosome 4 with the human orthologous region, at 6q24.3, that conserves the ancestral genomic organization. A 250-kilobase segment was extensively duplicated around the macaque centromere. These duplications were strictly intrachromosomal. Our results suggest that novel centromeres may trigger only local duplication activity and that the absence of genes in the seeding region may have been important in ENC maintenance and progression.

Evolutionary genomic remodelling of the human 4q subtelomere (4q35.2)

BACKGROUND

In order to obtain insights into the functionality of the human 4q35.2 domain harbouring the facioscapulohumeral muscular dystrophy (FSHD) locus, we investigated in African apes genomic and chromatin organisations, and the nuclear topology of orthologous regions.

RESULTS

A basic block consisting of short D4Z4 arrays (10-15 repeats), 4q35.2 specific sequences, and approximately 35 kb of interspersed repeats from different LINE subfamilies was repeated at least twice in the gorilla 4qter. This genomic organisation has undergone evolutionary remodelling, leading to the single representation of both the D4Z4 array and LINE block in chimpanzee, and the loss of the LINE block in humans. The genomic remodelling has had an impact on 4qter chromatin organisation, but not its interphase nuclear topology. In comparison with humans, African apes show very low or undetectable levels of FRG1 and FRG2 histone 4 acetylation and gene transcription, although histone deacetylase inhibition restores gene transcription to levels comparable with those of human cells, thus indicating that the 4qter region is capable of acquiring a more open chromatin structure. Conversely, as in humans, the 4qter region in African apes has a very peripheral nuclear localisation.

CONCLUSION

The 4q subtelomere has undergone substantial genomic changes during evolution that have had an impact on chromatin condensation and the region's transcriptional regulation. Consequently, the 4qter genes in African apes and humans seem to be subjected to a different strategy of regulation in which LINE and D4Z4 sequences may play a pivotal role. However, the effect of peripheral nuclear anchoring of 4qter on these regulation mechanisms is still unclear. The observed differences in the regulation of 4qter gene expression between African apes and humans suggest that the human 4q35.2 locus has acquired a novel functional relevance.

Independent centromere formation in a capricious, gene-free domain of chromosome 13q21 in Old World monkeys and pigs

BACKGROUND

Evolutionary centromere repositioning and human analphoid neocentromeres occurring in clinical cases are, very likely, two stages of the same phenomenon whose properties still remain substantially obscure. Chromosome 13 is the chromosome with the highest number of neocentromeres. We reconstructed the mammalian evolutionary history of this chromosome and characterized two human neocentromeres at 13q21, in search of information that could improve our understanding of the relationship between evolutionarily new centromeres, inactivated centromeres, and clinical neocentromeres.

RESULTS

Chromosome 13 evolution was studied, using FISH experiments, across several diverse superordinal phylogenetic clades spanning >100 million years of evolution. The analysis revealed exceptional conservation among primates (hominoids, Old World monkeys, and New World monkeys), Carnivora (cat), Perissodactyla (horse), and Cetartiodactyla (pig). In contrast, the centromeres in both Old World monkeys and pig have apparently repositioned independently to a central location (13q21). We compared these results to the positions of two human 13q21 neocentromeres using chromatin immunoprecipitation and genomic microarrays.

CONCLUSION

We show that a gene-desert region at 13q21 of approximately 3.9 Mb in size possesses an inherent potential to form evolutionarily new centromeres over, at least, approximately 95 million years of mammalian evolution. The striking absence of genes may represent an important property, making the region tolerant to the extensive pericentromeric reshuffling during subsequent evolution. Comparison of the pericentromeric organization of chromosome 13 in four Old World monkey species revealed many differences in sequence organization. The region contains clusters of duplicons showing peculiar features.

The boundary of macaque rDNA is constituted by low-copy sequences conserved during evolution

In Macaca mulatta, the single rDNA array is flanked by a patchwork of sequences including subregions of human Yp11.2, 4q35.2, and 10p15.3. This composite DNA region is characterized by unique or low-copy sequences, resembling a potentially transcribed region. The analysis of Cercopithecus aethiops, Presbytis cristata, and Hylobates lar suggests that this complex sequence organization could be shared by Old World monkey and lesser ape species. After the lesser apes/great apes divergence, the unique or nonduplicated DNA region underwent amplification and spreading, preferentially marking the p arm of acrocentric chromosomes bearing the rDNA. The molecular analysis of human acrocentric chromosomes revealed some extent of remodeling of the rDNA boundary: near the human NOR, a large 4q35.2 duplication partially resembles that found in MMU; conversely, infrequently represented Yp11.2 sequences totally differed from those of the macaque, and 10p15.3 sequences were lacking. Thus, although evolutionary events modified the sequence organization of the MMU rDNA boundary, its overall sequence feature and the preferential location in vicinity to the NOR have been conserved.

Evolutionary movement of centromeres in horse, donkey, and zebra

Centromere repositioning (CR) is a recently discovered biological phenomenon consisting of the emergence of a new centromere along a chromosome and the inactivation of the old one. After a CR, the primary constriction and the centromeric function are localized in a new position while the order of physical markers on the chromosome remains unchanged. These events profoundly affect chromosomal architecture. Since horses, asses, and zebras, whose evolutionary divergence is relatively recent, show remarkable morphological similarity and capacity to interbreed despite their chromosomes differing considerably, we investigated the role of CR in the karyotype evolution of the genus Equus. Using appropriate panels of BAC clones in FISH experiments, we compared the centromere position and marker order arrangement among orthologous chromosomes of Burchelli's zebra (Equus burchelli), donkey (Equus asinus), and horse (Equus caballus). Surprisingly, at least eight CRs took place during the evolution of this genus. Even more surprisingly, five cases of CR have occurred in the donkey after its divergence from zebra, that is, in a very short evolutionary time (approximately 1 million years). These findings suggest that in some species the CR phenomenon could have played an important role in karyotype shaping, with potential consequences on population dynamics and speciation.

A preliminary comparative analysis of primate segmental duplications shows elevated substitution rates and a great-ape expansion of intrachromosomal duplications

Compared with other sequenced animal genomes, human segmental duplications appear larger, more interspersed, and disproportionately represented as high-sequence identity alignments. Global sequence divergence estimates of human duplications have suggested an expansion relatively recently during hominoid evolution. Based on primate comparative sequence analysis of 37 unique duplication-transition regions, we establish a molecular clock for their divergence that shows a significant increase in their effective substitution rate when compared with unique genomic sequence. Fluorescent in situ hybridization (FISH) analyses from 1053 random nonhuman primate BACs indicate that great-ape species have been enriched for interspersed segmental duplications compared with representative Old World and New World monkeys. These findings support computational analyses that show a 12-fold excess of recent (>98%) intrachromosomal duplications when compared with duplications between nonhomologous chromosomes. These architectural shifts in genomic structure and elevated substitution rates have important implications for the emergence of new genes, gene-expression differences, and structural variation among humans and great apes.

Evolutionary history of chromosome 20

The evolutionary history of human chromosome 20 in primates was investigated using a panel of human BAC/PAC probes spaced along the chromosome. Oligonucleotide primers derived from the sequence of each human clone were used to screen horse, cat, pig, and black lemur BAC libraries to assemble, for each species, a panel of probes mapping to chromosomal loci orthologous to the loci encompassed by the human BACs. This approach facilitated marker-order comparison aimed at defining marker arrangement in primate ancestor. To this goal, we also took advantage of the mouse and rat draft sequences. The almost perfect colinearity of chromosome 20 sequence in humans and mouse could be interpreted as evidence that their form was ancestral to primates. Contrary to this view, we found that horse, macaque, and two New World monkeys share the same marker-order arrangement from which the human and mouse forms can be derived, assuming similar but distinct inversions that fully account for the small difference in marker arrangement between humans and mouse. The evolutionary history of this chromosome unveiled also two centromere repositioning events in New World monkey species.

Recurrent sites for new centromere seeding

Using comparative FISH and genomics, we have studied and compared the evolution of chromosome 3 in primates and two human neocentromere cases on the long arm of this chromosome. Our results show that one of the human neocentromere cases maps to the same 3q26 chromosomal region where a new centromere emerged in a common ancestor of the Old World monkeys approximately 25-40 million years ago. Similarly, the locus in which a new centromere was seeded in the great apes' ancestor was orthologous to the site in which a new centromere emerged in the New World monkeys' ancestor. These data suggest the recurrent use of longstanding latent centromeres and that there is an inherent potential of these regions to form centromeres. The second human neocentromere case (3q24) revealed unprecedented features. The neocentromere emergence was not accompanied by any chromosomal rearrangement that usually triggers these events. Instead, it involved the functional inactivation of the normal centromere, and was present in an otherwise phenotypically normal individual who transmitted this unusual chromosome to the next generation. We propose that the formation of neocentromeres in humans and the emergence of new centromeres during the course of evolution share a common mechanism.

Evolution of beta satellite DNA sequences: evidence for duplication-mediated repeat amplification and spreading

In this article, we report studies on the evolutionary history of beta satellite repeats (BSR) in primates. In the orangutan genome, the bulk of BSR sequences was found organized as very short stretches of approximately 100 to 170 bp, embedded in a 60-kb to 80-kb duplicated DNA segment. The estimated copy number of the duplicon that carries BSR sequences ranges from 70 to 100 per orangutan haploid genome. In both macaque and gibbon, the duplicon mapped to a single chromosomal region at the boundary of the rDNA on the marker chromosome (chromosome 13 and 12, respectively). However, only in the gibbon, the duplicon comprised 100 bp of beta satellite. Thus, the ancestral copy of the duplicon appeared in Old World monkeys ( approximately 25 to approximately 35 MYA), whereas the prototype of beta satellite repeats took place in a gibbon ancestor, after apes/Old World monkeys divergence ( approximately 25 MYA). Subsequently, a burst in spreading of the duplicon that carries the beta satellite was observed in the orangutan, after lesser apes divergence from the great apes-humans lineage ( approximately 18 MYA). The analysis of the orangutan genome also indicated the existence of two variants of the duplication that differ for the length (100 or 170 bp) of beta satellite repeats. The latter organization was probably generated by nonhomologous recombination between two 100-bp repeated regions, and it likely led to the duplication of the single Sau3A site present in the 100-bp variant, which generated the prototype of Sau3A 68-bp beta satellite tandem organization. The two variants of the duplication, although with a different ratios, characterize the hominoid genomes from the orangutan to humans, preferentially involving acrocentric chromosomes. At variance to alpha satellite, which appeared before the divergence of New World and Old World monkeys, the beta satellite evolutionary history began in apes ancestor, where we have first documented a low-copy, nonduplicated BSR sequence. The first step of BSR amplification and spreading occurred, most likely, because the BSR was part of a large duplicon, which underwent a burst dispersal in great apes' ancestor after the lesser apes' branching. Then, after orangutan divergence, BSR acquired the clustered structural organization typical of satellite DNA.

Localization of beta-defensin genes in non human primates

Defensins are a family of host defence peptides that play an important role in the innate immunity of mammalian and avian species. In humans, four beta-defensins have been isolated so far, corresponding to the products of the genes DEFB1 (h-BD1, GenBank accession number NM_005218); DEFB4 (h-Bd2, NM_004942.2), DEFB103 (h-BD3, NM_018661); and DEFB104 (hBD4, NM_080389) mapping on chromosome 8p23.22. We have localized beta-defensin genes on metaphasic chromosomes of great apes and several non-human primate species to determine their physical mapping. Using fluorescent in situ hybridization and BAC probes containing the four beta-defensin genes, we have mapped the homologous regions to the beta-defensin genes on chromosome 8p23-p.22 in non-human primates, while no signals were detected on prosimians chromosomes.

Refinement of a chimpanzee pericentric inversion breakpoint to a segmental duplication cluster

BACKGROUND

Pericentric inversions are the most common euchromatic chromosomal differences among humans and the great apes. The human and chimpanzee karyotype differs by nine such events, in addition to several constitutive heterochromatic increases and one chromosomal fusion event. Reproductive isolation and subsequent speciation are thought to be the potential result of pericentric inversions, as reproductive boundaries form as a result of hybrid sterility.

RESULTS

Here we employed a comparative fluorescence in situ hybridization approach, using probes selected from a combination of physical mapping, genomic sequence, and segmental duplication analyses to narrow the breakpoint interval of a pericentric inversion in chimpanzee involving the orthologous human 15q11-q13 region. We have refined the inversion breakpoint of this chimpanzee-specific rearrangement to a 600 kilobase (kb) interval of the human genome consisting of entirely duplicated material. Detailed analysis of the underlying sequence indicated that this region comprises multiple segmental duplications, including a previously characterized duplication of the alpha7 neuronal nicotinic acetylcholine receptor subunit gene (CHRNA7) in 15q13.3 and several Golgin-linked-to-PML, or LCR15, duplications.

CONCLUSIONS

We conclude that, on the basis of experimental data excluding the CHRNA7 duplicon as the site of inversion, and sequence analysis of regional duplications, the most likely rearrangement site is within a GLP/LCR15 duplicon. This study further exemplifies the genomic plasticity due to the presence of segmental duplications and highlights their importance for a complete understanding of genome evolution.

Analysis of chromosome conservation in Lemur catta studied by chromosome paints and BAC/PAC probes

A panel of human chromosome painting probes and bacterial and P1 artificial chromosome (BAC/PAC) clones were used in fluorescence in situ hybridization (FISH) experiments to investigate the chromosome conservation of the ring-tailed lemur (Lemur catta, LCA) with respect to human. Whole chromosome paints specific for human chromosomes 7, 9, 11, 13, 14, 17, 18, 20, 21, and X were found to identify a single chromosome or an uninterrupted chromosomal region in LCA. A large set of partial chromosome paints and BAC/PAC probes were then used to refine the characterization of the rearrangements differentiating the two karyotypes. The results were also used to reconstruct the ancestral Lemuridae karyotype. Lemur catta, indeed, can be used as an outgroup, allowing symplesiomorphic (ancestral) rearrangements to be distinguished from apomorphic (derived) rearrangements in lemurs. Some LCA chromosomes are difficult to distinguish morphologically. The 'anchorage' of most LCA chromosomes to specific probes will contribute to the standardization of the karyotype of this species.

Localization of a new highly repeated DNA sequence of Lemur cafta (Lemuridae, Strepsirhini)

We have isolated and cloned an 800-bp highly repeated DNA (HRDNA) sequence from Lemur catta (LCA) and described its localization on LCA chromosomes. Lemur catta HRDNA sequences were localized by performing FISH experiments on standard and elongated metaphasic chromosomes using an LCA HRDNA probe (LCASAT). A complex hybridization pattern was detected. A strong pericentromeric hybridization signal was observed on most LCA chromosomes. Chromosomes 7 and 13 were lit in pericentromeric regions, as well as in the interspersed heterochromatin. Chromosomes 1, 3, 4, 17, 19, X, and microchromosomes (20, 25, 26, and 27) showed no signals in the pericentromeric region, but chromosomes 3 and 4 showed a positive hybridization in heterochromatic regions. The 800-bp L catta HRDNA was species specific. We performed FISH experiments with the LCASAT probe on Eulemur macaco macaco (EMA) and Eulemur fulvus fulvus (EFU) metaphases and no positive signal of hybridization was detected. These findings were also confirmed by Southern blot analysis and PCR.

Characterization of a highly repeated DNA sequence family in five species of the genus Eulemur

The karyotypes of Eulemur species exhibit a high degree of variation, as a consequence of the Robertsonian fusion and/or centromere fission. Centromeric and pericentromeric heterochromatin of eulemurs is constituted by highly repeated DNA sequences (including some telomeric TTAGGG repeats) which have so far been investigated and used for the study of the systematic relationships of the different species of the genus Eulemur. In our study, we have cloned a set of repetitive pericentromeric sequences of five Eulemur species: E. fulvus fulvus (EFU), E. mongoz (EMO), E. macaco (EMA), E. rubriventer (ERU), and E. coronatus (ECO). We have characterized these clones by sequence comparison and by comparative fluorescence in situ hybridization analysis in EMA and EFU. Our results showed a high degree of sequence similarity among Eulemur species, indicating a strong conservation, within the five species, of these pericentromeric highly repeated DNA sequences.