Body condition scoring facilitates healthcare monitoring in Hermann's Tortoises (Testudo hermanni ssp.).

Clinical assessment of body condition is crucial in captive and free-ranging reptiles, since a large percentage of diseases result from inadequate nutrition. However, preventive health care is restricted by the lack of a practical method for the assessment in tortoises. Pre-existing evaluation systems based on weight and shell measurements are laborious and ignore the clinical presentation of the animal. The present study aimed to facilitate the assessment by establishing a body condition score. A total of 373 Hermann's Tortoises (Testudo hermanni) (n = 281 tortoises kept as pets in Germany and n = 92 tortoises originating from a free-ranging population (68) or a rearing station (24) in France) were examined and data (weight (g), carapace length (cm), width (cm), height (cm)) were recorded in a standard protocol between October 2020 and October 2021. A modified version of a body condition score for Mojave Desert Tortoises (Gopherus agassizii) (1 = cachectic, 3 = ideal, 5 = obese) was utilized and tested against pre-existing shell measurement systems (Jackson's ratio, body condition index, volume condition index, circumferential product). German captive tortoises were significantly heavier and larger than French specimens. In the Spearman's correlation matrix, the body condition score showed a statistically significant correlation with all measurement methods in the total population of captive tortoises (Testudo hermanni boettgeri), with a medium correlation strength, and a lack of correlation in free-ranging tortoises (Testudo hermanni hermanni). However, individual animal data suggested misleading results of mathematical equations in terms of body condition. Clinical evaluation of tortoises, including a body condition score, should be considered essential to provide good healthcare and should be an integral part of general examination.

Reduce, Replace, Refine: Determining A Posteriori Reference Intervals for Biochemistry in Hermann's Tortoise (Testudo hermanni).

Biochemical and hematological analyses are important for the assessment of animal health. However, for most wild species their use is hindered by the scarcity of reliable reference intervals. Indeed, collecting body fluids (e.g., blood, urine) in free-ranging animals is often technically challenging. Further, sampling many individuals would be essential to consider major sources of variations, such as species, populations, sex, age, and seasons. One alternative, according to the reduction, refinement, and replacement framework, is to establish reference intervals a posteriori using literature survey and unpublished data. We produced reference intervals for free-ranging Hermann's tortoises (Testudo hermanni), using analyses performed on blood samples collected in previous studies and conservation programs conducted in the field between 2010 and 2016 in southern France (n=195 individuals). Thirteen parameters were analyzed: packed-cell volume, blood concentrations of corticosterone, testosterone, glycemia, cholesterol, triglyceride, urea, uric acid, calcium, sodium, potassium, asparagine aminotransferase (AST), and alanine aminotransferases (ALT). Reference intervals for subgroups defined by sex and season were relevant for corticosterone, triglyceride, and calcium (sex) and cholesterol (season). Comparing our results with those obtained in captive individuals in Germany, except for urea and AST levels the intervals from both free-ranging versus captive tortoises were similar, suggesting that reference intervals established from captive individuals may be suitable for free-ranging populations in this species.

Proteomic insight into the venom composition of the largest European rear-fanged snake, Malpolon monspessulanus monspessulanus.

Snake envenomations constitute a worldwide neglected tropical disease, with the vast majority of lethal bites inflicted by front-fanged snakes from the viperid and elapid groups. Rear-fanged snakes (colubrids) were often considered harmless and as a result, are much less studied, but several documented deaths have suggested potent venom in this group as well. The largest European snake (Malpolon monspessulanus monspessulanus), known as the "Montpellier snake", is such a rear-fanged snake that belongs to the Lamprophiidae family. Its venom remains largely unknown but cases of envenomation with neurological symptoms have been reported. Here, we provide the first insights into the composition of its venom using mass spectrometry methods. First, liquid chromatography coupled mass spectrometry analysis of the manually collected venom samples reveals a complex profile, with the majority of masses encompassing the range 500-3000 Da, 4000-8000 Da, and 10 000-30 000 Da. Next, shotgun proteomics allowed the identification of a total of 42 different known families of proteins, including snake venom metalloproteinases, peptidase M1, and cysteine-rich secretory proteins, as the most prominent. Interestingly, three-finger toxins were not detected, suggesting that neurotoxicity may occur via other, yet to be determined, toxin types. Overall, our results provide the basis for a better understanding of the effects of a peculiar snake venom on human symptomatology, but also on the main prey consumed by this species.

The large repertoire of conifer NLR resistance genes includes drought responsive and highly diversified RNLs.

The NLRs or NBS-LRRs (nucleotide-binding, leucine-rich-repeat) form the largest resistance gene family in plants, with lineage-specific contingents of TNL, CNL and RNL subfamilies and a central role in resilience to stress. The origin, evolution and distribution of NLR sequences has been unclear owing in part to the variable size and diversity of the RNL subfamily and a lack of data in Gymnosperms. We developed, searched and annotated transcriptomes assemblies of seven conifers and identified a resource of 3816 expressed NLR sequences. Our analyses encompassed sequences data spanning the major groups of land plants and determinations of NLR transcripts levels in response to drought in white spruce. We showed that conifers have among the most diverse and numerous RNLs in tested land plants. We report an evolutionary swap in the formation of RNLs, which emerged from the fusion of an RPW8 domain to a NB-ARC domain of CNL. We uncovered a quantitative relationship between RNLs and TNLs across all land plants investigated, with an average ratio of 1:10. The conifer RNL repertoire harbours four distinct groups, with two that differ from Angiosperms, one of which contained several upregulated sequences in response to drought while the majority of responsive NLRs are downregulated.

Gene copy number variations in adaptive evolution: The genomic distribution of gene copy number variations revealed by genetic mapping and their adaptive role in an undomesticated species, white spruce (Picea glauca).

Gene copy number variation (CNV) has been associated with phenotypic variability in animals and plants, but a genomewide understanding of their impacts on phenotypes is largely restricted to human and agricultural systems. As such, CNVs have rarely been considered in investigations of the genomic architecture of adaptation in wild species. Here, we report on the genetic mapping of gene CNVs in white spruce, which lacks a contiguous assembly of its large genome (~20 Gb), and their relationships with adaptive phenotypic variation. We detected 3,911 gene CNVs including de novo structural variations using comparative genome hybridization on arrays (aCGH) in a large progeny set. We inferred the heterozygosity at CNV loci within parents by comparing haploid and diploid tissues and genetically mapped 82 gene CNVs. Our analysis showed that CNVs were distributed over 10 linkage groups and identified four CNV hotspots that we predict to occur in other species of the Pinaceae. Significant relationships were found between 29 of the gene CNVs and adaptive traits based on regression analyses with timings of bud set and bud flush, and height growth, suggesting a role for CNVs in climate adaptation. The importance of CNVs in adaptive evolution of white spruce was also indicated by functional gene annotations and the clustering of 31% of the mapped adaptive gene CNVs in CNV hotspots. Taken together, these results illustrate the feasibility of studying CNVs in undomesticated species and represent a major step towards a better understanding of the roles of CNVs in adaptive evolution.

CNVs into the wild: screening the genomes of conifer trees (Picea spp.) reveals fewer gene copy number variations in hybrids and links to adaptation.

BACKGROUND: Copy number variations (CNVs) have been linked to different phenotypes in human, including many diseases. A genome-scale understanding of CNVs is available in a few plants but none are wild species, leaving a knowledge gap regarding their genome biology and evolutionary role. We developed a reliable CNV detection method for species lacking contiguous reference genome. We selected multiple probes within 14,078 gene sequences and developed comparative genome hybridization on arrays. Gene CNVs were assessed in three full-sib families from species with 20 Gb genomes, i.e., white and black spruce, and interior spruce - a natural hybrid. RESULTS: We discovered hundreds of gene CNVs in each species, 3612 in total, which were enriched in functions related to stress and defense responses and narrow expression profiles, indicating a potential role in adaptation. The number of shared CNVs was in accordance with the degree of relatedness between individuals and species. The genetically mapped subset of these genes showed a wide distribution across the genome, implying numerous structural variations. The hybrid family presented significantly fewer CNVs, suggesting that the admixture of two species within one genome reduces the occurrence of CNVs. CONCLUSIONS: The approach we developed is of particular interest in non-model species lacking a reference genome. Our findings point to a role for CNVs in adaptation. Their reduced abundance in the hybrid may limit genetic variability and evolvability of hybrids. We propose that CNVs make a qualitatively distinct contribution to adaptation which could be important for short term change.

Sex-specific ecophysiological responses to environmental fluctuations of free-ranging Hermann's tortoises: implication for conservation.

Physiological parameters provide indicators to evaluate how organisms respond to conservation actions. For example, individuals translocated during reinforcement programmes may not adapt to their novel host environment and may exhibit elevated chronic levels of stress hormones and/or decreasing body condition. Conversely, successful conservation actions should be associated with a lack of detrimental physiological perturbation. However, physiological references fluctuate over time and are influenced by various factors (e.g. sex, age, reproductive status). It is therefore necessary to determine the range of natural variations of the selected physiological metrics to establish useful baselines. This study focuses on endangered free-ranging Hermann's tortoises (Testudo hermanni hermanni), where conservation actions have been preconized to prevent extinction of French mainland populations. The influence of sex and of environmental factors (site, year and season) on eight physiological parameters (e.g. body condition, corticosterone concentrations) was assessed in 82 individuals from two populations living in different habitats. Daily displacements were monitored by radio-tracking. Most parameters varied between years and seasons and exhibited contrasting sex patterns but with no or limited effect of site. By combining behavioural and physiological traits, this study provides sex-specific seasonal baselines that can be used to monitor the health status of Hermann's tortoises facing environmental threats (e.g. habitat changes) or during conservation actions (e.g. translocation). These results might also assist in selection of the appropriate season for translocation.

Modular organization of the white spruce (Picea glauca) transcriptome reveals functional organization and evolutionary signatures.

Transcript profiling has shown the molecular bases of several biological processes in plants but few studies have developed an understanding of overall transcriptome variation. We investigated transcriptome structure in white spruce (Picea glauca), aiming to delineate its modular organization and associated functional and evolutionary attributes. Microarray analyses were used to: identify and functionally characterize groups of co-expressed genes; investigate expressional and functional diversity of vascular tissue preferential genes which were conserved among Picea species, and identify expression networks underlying wood formation. We classified 22 857 genes as variable (79%; 22 coexpression groups) or invariant (21%) by profiling across several vegetative tissues. Modular organization and complex transcriptome restructuring among vascular tissue preferential genes was revealed by their assignment to coexpression groups with partially overlapping profiles and partially distinct functions. Integrated analyses of tissue-based and temporally variable profiles identified secondary xylem gene networks, showed their remodelling over a growing season and identified PgNAC-7 (no apical meristerm (NAM), Arabidopsis transcription activation factor (ATAF) and cup-shaped cotyledon (CUC) transcription factor 007 in Picea glauca) as a major hub gene specific to earlywood formation. Reference profiling identified comprehensive, statistically robust coexpressed groups, revealing that modular organization underpins the evolutionary conservation of the transcriptome structure.

Opposite action of R2R3-MYBs from different subgroups on key genes of the shikimate and monolignol pathways in spruce.

Redundancy and competition between R2R3-MYB activators and repressors on common target genes has been proposed as a fine-tuning mechanism for the regulation of plant secondary metabolism. This hypothesis was tested in white spruce [Picea glauca (Moench) Voss] by investigating the effects of R2R3-MYBs from different subgroups on common targets from distinct metabolic pathways. Comparative analysis of transcript profiling data in spruces overexpressing R2R3-MYBs from loblolly pine (Pinus taeda L.), PtMYB1, PtMYB8, and PtMYB14, defined a set of common genes that display opposite regulation effects. The relationship between the closest MYB homologues and 33 putative target genes was explored by quantitative PCR expression profiling in wild-type P. glauca plants during the diurnal cycle. Significant Spearman's correlation estimates were consistent with the proposed opposite effect of different R2R3-MYBs on several putative target genes in a time-related and tissue-preferential manner. Expression of sequences coding for 4CL, DHS2, COMT1, SHM4, and a lipase thio/esterase positively correlated with that of PgMYB1 and PgMYB8, but negatively with that of PgMYB14 and PgMYB15. Complementary electrophoretic mobility shift assay (EMSA) and transactivation assay provided experimental evidence that these different R2R3-MYBs are able to bind similar AC cis-elements in the promoter region of Pg4CL and PgDHS2 genes but have opposite effects on their expression. Competitive binding EMSA experiments showed that PgMYB8 competes more strongly than PgMYB15 for the AC-I MYB binding site in the Pg4CL promoter. Together, the results bring a new perspective to the action of R2R3-MYB proteins in the regulation of distinct but interconnecting metabolism pathways.

Potential link between biotic defense activation and recalcitrance to induction of somatic embryogenesis in shoot primordia from adult trees of white spruce (Picea glauca).

BACKGROUND: Among the many commercial opportunities afforded by somatic embryogenesis (SE), it is the ability to clonally propagate individual plants with rare or elite traits that has some of the most significant implications. This is particularly true for many long-lived species, such as conifers, but whose long generation times pose substantive challenges, including increased recalcitrance for SE as plants age. Identification of a clonal line of somatic embryo-derived trees whose shoot primordia have remained responsive to SE induction for over a decade, provided a unique opportunity to examine the molecular aspects underpinning SE within shoot tissues of adult white spruce trees. RESULTS: Microarray analysis was used to conduct transcriptome-wide expression profiling of shoot explants taken from this responsive genotype following one week of SE induction, which when compared with that of a nonresponsive genotype, led to the identification of four of the most differentially expressed genes within each genotype. Using absolute qPCR to expand the analysis to three weeks of induction revealed that differential expression of all eight candidate genes was maintained to the end of the induction treatment, albeit to differing degrees. Most striking was that both the magnitude and duration of candidate gene expression within the nonresponsive genotype was indicative of an intense physiological response. Examining their putative identities further revealed that all four encoded for proteins with similarity to angiosperm proteins known to play prominent roles in biotic defense, and that their high-level induction over an extended period is consistent with activation of a biotic defense response. In contrast, the more temperate response within the responsive genotype, including induction of a conifer-specific dehydrin, is more consistent with elicitation of an adaptive stress response. CONCLUSIONS: While additional evidence is required to definitively establish an association between SE responsiveness and a specific physiological response, these results suggest that biotic defense activation may be antagonistic, likely related to the massive transcriptional and metabolic reprogramming that it elicits. A major issue for future work will be to determine how and if suppressing biotic defense activation could be used to promote a physiological state more conducive to SE induction.

Transcriptome profiling in conifers and the PiceaGenExpress database show patterns of diversification within gene families and interspecific conservation in vascular gene expression.

BACKGROUND: Conifers have very large genomes (13 to 30 Gigabases) that are mostly uncharacterized although extensive cDNA resources have recently become available. This report presents a global overview of transcriptome variation in a conifer tree and documents conservation and diversity of gene expression patterns among major vegetative tissues. RESULTS: An oligonucleotide microarray was developed from Picea glauca and P. sitchensis cDNA datasets. It represents 23,853 unique genes and was shown to be suitable for transcriptome profiling in several species. A comparison of secondary xylem and phelloderm tissues showed that preferential expression in these vascular tissues was highly conserved among Picea spp. RNA-Sequencing strongly confirmed tissue preferential expression and provided a robust validation of the microarray design. A small database of transcription profiles called PiceaGenExpress was developed from over 150 hybridizations spanning eight major tissue types. In total, transcripts were detected for 92% of the genes on the microarray, in at least one tissue. Non-annotated genes were predominantly expressed at low levels in fewer tissues than genes of known or predicted function. Diversity of expression within gene families may be rapidly assessed from PiceaGenExpress. In conifer trees, dehydrins and late embryogenesis abundant (LEA) osmotic regulation proteins occur in large gene families compared to angiosperms. Strong contrasts and low diversity was observed in the dehydrin family, while diverse patterns suggested a greater degree of diversification among LEAs. CONCLUSION: Together, the oligonucleotide microarray and the PiceaGenExpress database represent the first resource of this kind for gymnosperm plants. The spruce transcriptome analysis reported here is expected to accelerate genetic studies in the large and important group comprised of conifer trees.

Association genetics of wood physical traits in the conifer white spruce and relationships with gene expression.

Marker-assisted selection holds promise for highly influencing tree breeding, especially for wood traits, by considerably reducing breeding cycles and increasing selection accuracy. In this study, we used a candidate gene approach to test for associations between 944 single-nucleotide polymorphism markers from 549 candidate genes and 25 wood quality traits in white spruce. A mixed-linear model approach, including a weak but nonsignificant population structure, was implemented for each marker-trait combination. Relatedness among individuals was controlled using a kinship matrix estimated either from the known half-sib structure or from the markers. Both additive and dominance effect models were tested. Between 8 and 21 single-nucleotide polymorphisms (SNPs) were found to be significantly associated (P ≤ 0.01) with each of earlywood, latewood, or total wood traits. After controlling for multiple testing (Q ≤ 0.10), 13 SNPs were still significant across as many genes belonging to different families, each accounting for between 3 and 5% of the phenotypic variance in 10 wood characters. Transcript accumulation was determined for genes containing SNPs associated with these traits. Significantly different transcript levels (P ≤ 0.05) were found among the SNP genotypes of a 1-aminocyclopropane-1-carboxylate oxidase, a ß-tonoplast intrinsic protein, and a long-chain acyl-CoA synthetase 9. These results should contribute toward the development of efficient marker-assisted selection in an economically important tree species.

Subgroup 4 R2R3-MYBs in conifer trees: gene family expansion and contribution to the isoprenoid- and flavonoid-oriented responses.

Transcription factors play a fundamental role in plants by orchestrating temporal and spatial gene expression in response to environmental stimuli. Several R2R3-MYB genes of the Arabidopsis subgroup 4 (Sg4) share a C-terminal EAR motif signature recently linked to stress response in angiosperm plants. It is reported here that nearly all Sg4 MYB genes in the conifer trees Picea glauca (white spruce) and Pinus taeda (loblolly pine) form a monophyletic clade (Sg4C) that expanded following the split of gymnosperm and angiosperm lineages. Deeper sequencing in P. glauca identified 10 distinct Sg4C sequences, indicating over-representation of Sg4 sequences compared with angiosperms such as Arabidopsis, Oryza, Vitis, and Populus. The Sg4C MYBs share the EAR motif core. Many of them had stress-responsive transcript profiles after wounding, jasmonic acid (JA) treatment, or exposure to cold in P. glauca and P. taeda, with MYB14 transcripts accumulating most strongly and rapidly. Functional characterization was initiated by expressing the P. taeda MYB14 (PtMYB14) gene in transgenic P. glauca plantlets with a tissue-preferential promoter (cinnamyl alcohol dehydrogenase) and a ubiquitous gene promoter (ubiquitin). Histological, metabolite, and transcript (microarray and targeted quantitative real-time PCR) analyses of PtMYB14 transgenics, coupled with mechanical wounding and JA application experiments on wild-type plantlets, allowed identification of PtMYB14 as a putative regulator of an isoprenoid-oriented response that leads to the accumulation of sesquiterpene in conifers. Data further suggested that PtMYB14 may contribute to a broad defence response implicating flavonoids. This study also addresses the potential involvement of closely related Sg4C sequences in stress responses and plant evolution.

High nitrogen fertilization and stem leaning have overlapping effects on wood formation in poplar but invoke largely distinct molecular pathways.

Previous studies indicated that high nitrogen fertilization may impact secondary xylem development and alter fibre anatomy and composition. The resulting wood shares some resemblance with tension wood, which has much thicker cell walls than normal wood due to the deposition of an additional layer known as the G-layer. This report compares the short-term effects of high nitrogen fertilization and tree leaning to induce tension wood, either alone or in combination, upon wood formation in young trees of Populus trichocarpa (Torr. & Gray) × P. deltoides Bartr. ex Marsh. Fibre anatomy, chemical composition and transcript profiles were examined in newly formed secondary xylem. Each of the treatments resulted in thicker cell walls relative to the controls. High nitrogen and tree leaning had overlapping effects on chemical composition based on Fourier transform infrared analysis, specifically indicating that secondary cell wall composition was shifted in favour of cellulose and hemicelluloses relative to lignin content. In contrast, the high-nitrogen trees had shorter fibres, whilst the leaning trees had longer fibres that the controls. Microarray transcript profiling carried out after 28 days of treatment identified 180 transcripts that accumulated differentially in one or more treatments. Only 10% of differentially expressed transcripts were affected in all treatments relative to the controls. Several of the affected transcripts were related to carbohydrate metabolism, secondary cell wall formation, nitrogen metabolism and osmotic stress. RT-qPCR analyses at 1, 7 and 28 days showed that several transcripts followed very different accumulation profiles in terms of rate and level of accumulation, depending on the treatment. Our findings suggest that high nitrogen fertilization and tension wood induction elicit largely distinct and molecular pathways with partial overlap. When combined, the two types of environmental cue yielded additive effects.

Identification of conserved core xylem gene sets: conifer cDNA microarray development, transcript profiling and computational analyses.

One approach for investigating the molecular basis of wood formation is to integrate microarray profiling data sets and sequence analyses, comparing tree species with model plants such as Arabidopsis. Conifers may be included in comparative studies thanks to large-scale expressed sequence tag (EST) analyses, which enable the development of cDNA microarrays with very significant genome coverage. A microarray of 10,400 low-redundancy sequences was designed starting from white spruce (Picea glauca (Moench.) Voss) cDNAs. Computational procedures that were developed to ensure broad transcriptome coverage and efficient PCR amplification were used to select cDNA clones, which were re-sequenced in the microarray manufacture process. White spruce transcript profiling experiments that compared secondary xylem to phloem and needles identified 360 xylem-preferential gene sequences. The functional annotations of all differentially expressed sequences were highly consistent with the results of similar analyses carried out in angiosperm trees and herbaceous plants. Computational analyses comparing the spruce microarray sequences and core xylem gene sets from Arabidopsis identified 31 transcripts that were highly conserved in angiosperms and gymnosperms, in terms of both sequence and xylem expression. Several other spruce sequences have not previously been linked to xylem differentiation (including genes encoding TUBBY-like domain proteins (TLPs) and a gibberellin insensitive (gai) gene sequence) or were shown to encode proteins of unknown function encompassing diverse conserved domains of unknown function.

Involvement of Pinus taeda MYB1 and MYB8 in phenylpropanoid metabolism and secondary cell wall biogenesis: a comparative in planta analysis.

The involvement of two R2R3-MYB genes from Pinus taeda L., PtMYB1 and PtMYB8, in phenylpropanoid metabolism and secondary cell wall biogenesis was investigated in planta. These pine MYBs were constitutively overexpressed (OE) in Picea glauca (Moench) Voss, used as a heterologous conifer expression system. Morphological, histological, chemical (lignin and soluble phenols), and transcriptional analyses, i.e. microarray and reverse transcription quantitative PCR (RT-qPCR) were used for extensive phenotyping of MYB-overexpressing spruce plantlets. Upon germination of somatic embryos, root growth was reduced in both transgenics. Enhanced lignin deposition was also a common feature but ectopic secondary cell wall deposition was more strongly associated with PtMYB8-OE. Microarray and RT-qPCR data showed that overexpression of each MYB led to an overlapping up-regulation of many genes encoding phenylpropanoid enzymes involved in lignin monomer synthesis, while misregulation of several cell wall-related genes and other MYB transcription factors was specifically associated with PtMYB8-OE. Together, the results suggest that MYB1 and MYB8 may be part of a conserved transcriptional network involved in secondary cell wall deposition in conifers.

A 6374 unigene set corresponding to low abundance transcripts expressed following fertilization in Solanum chacoense Bitt, and characterization of 30 receptor-like kinases.

In order to characterize regulatory genes that are expressed in ovule tissues after fertilization we have undertaken an EST sequencing project in Solanum chacoense, a self-incompatible wild potato species. Two cDNA libraries made from ovule tissues covering embryo development from zygote to late torpedo-stage were constructed and plated at high density on nylon membranes. To decrease EST redundancy and enrich for transcripts corresponding to weakly expressed genes a self-probe subtraction method was used to select the colonies harboring the genes to be sequenced. 7741 good sequences were obtained and, from these, 6374 unigenes were isolated. Thus, the self-probe subtraction resulted in a strong enrichment in singletons, a decrease in the number of clones per contigs, and concomitantly, an enrichment in the total number of unigenes obtained (82%). To gain insights into signal transduction events occurring during embryo development all the receptor-like kinases (or protein receptor kinases) were analyzed by quantitative real-time RT-PCR. Interestingly, 28 out of the 30 RLK isolated were predominantly expressed in ovary tissues or young developing fruits, and 23 were transcriptionaly induced following fertilization. Thus, the self-probe subtraction did not select for genes weakly expressed in the target tissue while being highly expressed elsewhere in the plant. Of the receptor-like kinases (RLK) genes isolated, the leucine-rich repeat (LRR) family of RLK was by far the most represented with 25 members covering 11 LRR classes.

Characterization of five RALF-like genes from Solanum chacoense provides support for a developmental role in plants.

Five RALF (rapid alkalinization factor)-like genes, named ScRALF1 to 5, were isolated from fertilized ovule and ovary cDNA libraries of Solanum chacoense. They showed high sequence similarities with the RALF protein sequence from Nicotiana tabacum, and exhibited the characteristic architecture of RALF polypeptides. All ScRALFs were moderately to highly expressed at some stage of fruit maturation. ScRALF1 and ScRALF3 were predominantly expressed in ovaries and larger fruits, while ScRALF2, ScRALF4, and ScRALF5 were also expressed in other tissues, indicating that while some RALFs may be involved in fruit maturation, others could be involved in other developmental processes. Wounding or treatment of plants with growth regulators involved in plant defense responses had no significant impact on the mRNA level of any of these genes. These results suggest and support previous data showing that RALF peptides are more likely to act as a small peptide involved in plant development than in defense responses.