Giusi Zaina

Conserved noncoding genomic sequences associated with a flowering-time quantitative trait locus in maize

Flowering time is a fundamental trait of maize adaptation to different agricultural environments. Although a large body of information is available on the map position of quantitative trait loci for flowering time, little is known about the molecular basis of quantitative trait loci. Through positional cloning and association mapping, we resolved the major flowering-time quantitative trait locus, Vegetative to generative transition 1 (Vgt1), to an ≈2-kb noncoding region positioned 70 kb upstream of an Ap2-like transcription factor that we have shown to be involved in flowering-time control. Vgt1 functions as a cis-acting regulatory element as indicated by the correlation of the Vgt1 alleles with the transcript expression levels of the downstream gene. Additionally, within Vgt1, we identified evolutionarily conserved noncoding sequences across the maize–sorghum–rice lineages. Our results support the notion that changes in distant cis-acting regulatory regions are a key component of plant genetic adaptation throughout breeding and evolution.

Disentangling the Roles of History and Local Selection in Shaping Clinal Variation of Allele Frequencies and Gene Expression in Norway Spruce (Picea abies)

Understanding the genetic basis of local adaptation is challenging due to the subtle balance among conflicting evolutionary forces that are involved in its establishment and maintenance. One system with which to tease apart these difficulties is clines in adaptive characters. Here we analyzed genetic and phenotypic variation in bud set, a highly heritable and adaptive trait, among 18 populations of Norway spruce (Picea abies), arrayed along a latitudinal gradient ranging from 47°N to 68°N. We confirmed that variation in bud set is strongly clinal, using a subset of five populations. Genotypes for 137 single-nucleotide polymorphisms (SNPs) chosen from 18 candidate genes putatively affecting bud set and 308 control SNPs chosen from 264 random genes were analyzed for patterns of genetic structure and correlation to environment. Population genetic structure was low (FST = 0.05), but latitudinal patterns were apparent among Scandinavian populations. Hence, part of the observed clinal variation should be attributable to population demography. Conditional on patterns of genetic structure, there was enrichment of SNPs within candidate genes for correlations with latitude. Twenty-nine SNPs were also outliers with respect to FST. The enrichment for clinal variation at SNPs within candidate genes (i.e., SNPs in PaGI, PaPhyP, PaPhyN, PaPRR7, and PaFTL2) indicated that local selection in the 18 populations, and/or selection in the ancestral populations from which they were recently derived, shaped the observed cline. Validation of these genes using expression studies also revealed that PaFTL2 expression is significantly associated with latitude, thereby confirming the central role played by this gene in the control of phenology in plants.

Bud set in poplar – genetic dissection of a complex trait in natural and hybrid populations

• The seasonal timing of growth events is crucial to tree distribution and conservation. The seasonal growth cycle is strongly adapted to the local climate that is changing because of global warming. We studied bud set as one cornerstone of the seasonal growth cycle in an integrative approach. • Bud set was dissected at the phenotypic level into several components, and phenotypic components with most genetic variation were identified. While phenotypic variation resided in the timing of growth cessation, and even so more in the duration from growth cessation to bud set, the timing of growth cessation had a stronger genetic component in both natural and hybrid populations. • Quantitative trait loci (QTL) were identified for the most discriminative phenotypic bud-set components across four poplar pedigrees. The QTL from different pedigrees were recurrently detected in six regions of the poplar genome. • These regions of 1.83-4.25 Mbp in size, containing between 202 and 394 genes, form the basis for further molecular-genetic dissection of bud set.

Nucleotide diversity and linkage disequilibrium in Populus nigra cinnamyl alcohol dehydrogenase (CAD4) gene

Cinnamyl alcohol dehydrogenase (CAD) is involved in the biosynthesis of lignin, a component of plant cell wall which negatively impacts paper pulp processing and biomass fermentation to ethanol. Transgenic poplars with depressed CAD activity show structural alterations of lignin. Natural CAD mutants have been identified in several plants; however, no natural CAD mutants have been identified in poplar. We surveyed the natural genetic variation in CAD4, a gene coding for CAD, in 360 poplar trees from Western Europe. We measured linkage disequilibrium (LD) between single-nucleotide polymorphisms (SNPs), performed neutrality tests and estimated diversity indexes, and investigated their dependence from sample size. We identified 45 SNPs, six of which caused an amino acid substitution. Our results suggest a short span of LD in Populus nigra CAD4 gene. We identified carriers of different nonsynonymous SNPs in CAD4; those subjects are candidate to be used in classical breeding programs to obtain carriers of different combinations of functional polymorphisms. We showed that use of small sample size might lead to biased estimates of LD, neutrality tests, and diversity indexes.

Phenotypic plasticity, QTL mapping and genomic characterization of bud set in black poplar

BackgroundThe genetic control of important adaptive traits, such as bud set, is still poorly understood in most forest trees species. Poplar is an ideal model tree to study bud set because of its indeterminate shoot growth. Thus, a full-sib family derived from an intraspecific cross of P. nigra with 162 clonally replicated progeny was used to assess the phenotypic plasticity and genetic variation of bud set in two sites of contrasting environmental conditions.ResultsSix crucial phenological stages of bud set were scored. Night length appeared to be the most important signal triggering the onset of growth cessation. Nevertheless, the effect of other environmental factors, such as temperature, increased during the process. Moreover, a considerable role of genotype × environment (G × E) interaction was found in all phenological stages with the lowest temperature appearing to influence the sensitivity of the most plastic genotypes.Descriptors of growth cessation and bud onset explained the largest part of phenotypic variation of the entire process. Quantitative trait loci (QTL) for these traits were detected. For the four selected traits (the onset of growth cessation (date2.5), the transition from shoot to bud (date1.5), the duration of bud formation (subproc1) and bud maturation (subproc2)) eight and sixteen QTL were mapped on the maternal and paternal map, respectively. The identified QTL, each one characterized by small or modest effect, highlighted the complex nature of traits involved in bud set process. Comparison between map location of QTL and P. trichocarpa genome sequence allowed the identification of 13 gene models, 67 bud set-related expressional and six functional candidate genes (CGs). These CGs are functionally related to relevant biological processes, environmental sensing, signaling, and cell growth and development. Some strong QTL had no obvious CGs, and hold great promise to identify unknown genes that affect bud set.ConclusionsThis study provides a better understanding of the physiological and genetic dissection of bud set in poplar. The putative QTL identified will be tested for associations in P. nigra natural populations. The identified QTL and CGs will also serve as useful targets for poplar breeding.

Characterization of the Poplar Pan-Genome by Genome-Wide Identification of Structural Variation

Many recent studies have emphasized the important role of structural variation (SV) in determining human genetic and phenotypic variation. In plants, studies aimed at elucidating the extent of SV are still in their infancy. Evidence has indicated a high presence and an active role of SV in driving plant genome evolution in different plant species. With the aim of characterizing the size and the composition of the poplar pan-genome, we performed a genome-wide analysis of structural variation in three intercrossable poplar species: Populus nigra, Populus deltoides, and Populus trichocarpa. We detected a total of 7,889 deletions and 10,586 insertions relative to the P. trichocarpa reference genome, covering respectively 33.2 Mb and 62.9 Mb of genomic sequence, and 3,230 genes affected by copy number variation (CNV). The majority of the detected variants are inter-specific in agreement with a recent origin following separation of species. Insertions and deletions (INDELs) were preferentially located in low-gene density regions of the poplar genome and were, for the majority, associated with the activity of transposable elements. Genes affected by SV showed lower-than-average expression levels and higher levels of dN/dS, suggesting that they are subject to relaxed selective pressure or correspond to pseudogenes. Functional annotation of genes affected by INDELs showed over-representation of categories associated with transposable elements activity, while genes affected by genic CNVs showed enrichment in categories related to resistance to stress and pathogens. This study provides a genome-wide catalogue of SV and the first insight on functional and structural properties of the poplar pan-genome.

A Picea abies linkage map based on SNP markers identifies QTLs for four aspects of resistance to Heterobasidion parviporum infection.

A consensus linkage map of Picea abies, an economically important conifer, was constructed based on the segregation of 686 SNP markers in a F1 progeny population consisting of 247 individuals. The total length of 1889.2 cM covered 96.5% of the estimated genome length and comprised 12 large linkage groups, corresponding to the number of haploid P. abies chromosomes. The sizes of the groups (from 5.9 to 9.9% of the total map length) correlated well with previous estimates of chromosome sizes (from 5.8 to 10.8% of total genome size). Any locus in the genome has a 97% probability to be within 10 cM from a mapped marker, which makes the map suited for QTL mapping. Infecting the progeny trees with the root rot pathogen Heterobasidion parviporum allowed for mapping of four different resistance traits: lesion length at the inoculation site, fungal spread within the sapwood, exclusion of the pathogen from the host after initial infection, and ability to prevent the infection from establishing at all. These four traits were associated with two, four, four and three QTL regions respectively of which none overlapped between the traits. Each QTL explained between 4.6 and 10.1% of the respective traits phenotypic variation. Although the QTL regions contain many more genes than the ones represented by the SNP markers, at least four markers within the confidence intervals originated from genes with known function in conifer defence; a leucoanthocyanidine reductase, which has previously been shown to upregulate during H. parviporum infection, and three intermediates of the lignification process; a hydroxycinnamoyl CoA shikimate/quinate hydroxycinnamoyltransferase, a 4-coumarate CoA ligase, and a R2R3-MYB transcription factor.

Endonuclease genes up-regulated in tissues undergoing programmed cell death are expressed during male gametogenesis in barley

In the process of programmed cell death (PCD), a key role has been attributed to endonucleases capable to cleave nuclear DNA at internucleosomal sites. In barley (Hordeum vulgare L.), two such nucleases (Bnuc1 and BEN1) were individually identified in unrelated tissues. In the present work, we demonstrate that their genes are also expressed in immature anthers at different stages of pollen development. Further experiments carried out on RNA extracted from immature barley anthers led to discover a novel endonuclease gene, namely Bnuc2 (AJ311603 in the EMBL/GenBank/DDBJ databases), eventually found up-regulated at the tetrad stage. The protein encoded was found to conserve large sequence portions of Bnuc1 and BEN1 endonucleases, including the domain regions involved in secretion and DNA/RNA binding. A survey conducted on barley EST libraries showed that Bnuc2 and BEN1 mRNAs are jointly present also in the transcriptome of 20 DAP spike and that other endonuclease ESTs are co-expressed with Bnuc1 or BEN1 in tissues where PCD has been recorded. Therefore, it can be concluded that during the PCD process, a set of S1-type endonucleases is synthesised regardless of the tissue considered.

New resources for genetic studies in Populus nigra: genome-wide SNP discovery and development of a 12k Infinium array.

Whole genome resequencing of 51 Populus nigra (L.) individuals from across Western Europe was performed using Illumina platforms. A total number of 1 878 727 SNPs distributed along the P. nigra reference sequence were identified. The SNP calling accuracy was validated with Sanger sequencing. SNPs were selected within 14 previously identified QTL regions, 2916 expressional candidate genes related to rust resistance, wood properties, water‐use efficiency and bud phenology and 1732 genes randomly spread across the genome. Over 10 000 SNPs were selected for the construction of a 12k Infinium Bead‐Chip array dedicated to association mapping. The SNP genotyping assay was performed with 888 P. nigra individuals. The genotyping success rate was 91%. Our high success rate was due to the discovery panel design and the stringent parameters applied for SNP calling and selection. In the same set of P. nigra genotypes, linkage disequilibrium throughout the genome decayed on average within 5–7 kb to half of its maximum value. As an application test, ADMIXTURE analysis was performed with a selection of 600 SNPs spread throughout the genome and 706 individuals collected along 12 river basins. The admixture pattern was consistent with genetic diversity revealed by neutral markers and the geographical distribution of the populations. These newly developed SNP resources and genotyping array provide a valuable tool for population genetic studies and identification of QTLs through natural‐population based genetic association studies in P. nigra.

Identifying Genetic Signatures of Natural Selection Using Pooled Population Sequencing in Picea abies

The joint inference of selection and past demography remain a costly and demanding task. We used next generation sequencing of two pools of 48 Norway spruce mother trees, one corresponding to the Fennoscandian domain, and the other to the Alpine domain, to assess nucleotide polymorphism at 88 nuclear genes. These genes are candidate genes for phenological traits, and most belong to the photoperiod pathway. Estimates of population genetic summary statistics from the pooled data are similar to previous estimates, suggesting that pooled sequencing is reliable. The nonsynonymous SNPs tended to have both lower frequency differences and lower FST values between the two domains than silent ones. These results suggest the presence of purifying selection. The divergence between the two domains based on synonymous changes was around 5 million yr, a time similar to a recent phylogenetic estimate of 6 million yr, but much larger than earlier estimates based on isozymes. Two approaches, one of them novel and that considers both FST and difference in allele frequencies between the two domains, were used to identify SNPs potentially under diversifying selection. SNPs from around 20 genes were detected, including genes previously identified as main target for selection, such as PaPRR3 and PaGI.