Genetic gain from selection for rooting ability and early growth in vegetatively propagated clones of loblolly pine

A successful clonal forestry program for loblolly pine based on rooted cutting technology needs to consider selection for both rooting ability and subsequent field growth. Rooting ability and second-year height were assessed in more than 2,000 clones from 70 full-sib families of loblolly pine. The bivariate analysis of rooting ability from five rooting trials and field growth from six field trials allowed for estimation of the genetic covariance between rooting ability and second-year height for parental effects, full-sib family effects, and the total genetic value of clones within full-sib family. There was a positive genetic relationship between rooting ability and second-year height at all three genetic levels. The genetic correlation at the parental level between rooting ability and second-year height was 0.32. At the full-sib family level, the genetic correlation between traits was 0.39. The correlation of total genetic values of clones for rooting ability and second-year height was 0.29. The genetic gain in rooting ability and second-year height was estimated for a number of deployment options based on various selection scenarios using the best linear unbiased prediction (BLUP) values from the bivariate analysis. The deployment strategies compared were (1) half-sib family deployment, (2) full-sib family deployment, and (3) clonal deployment. Moderate to high family and clonal mean heritabilities, moderate to high type B genetic correlations, and substantial among-family and among-clone genetic variation indicate the potential for increasing rooting efficiency and improving growth.     

Inheritance of foliar stable carbon isotope discrimination and third-year height in Pinus taeda clones on contrasting sites in Florida and Georgia

Quantifying foliar stable carbon isotope discrimination (Δ) is a powerful approach for understanding genetic variation in gas exchange traits in large populations. The genetic architecture of Δ and third-year height is described for more than 1,000 clones of Pinus taeda tested on two contrasting sites. for Δ was 0.14 (±0.03), 0.20 (±0.07), and 0.09 (±0.04) at Florida, Georgia, and across sites, respectively. for stable carbon isotope discrimination ranged from 0.25 (±0.03) at the Florida site to 0.33 (±0.03) at the Georgia site, while the across-site estimate of was 0.19 (±0.02). For third-year height, ranged from 0.13 (±0.05) at the Georgia site to 0.20 (±0.06) at the Florida site with an across-site estimate of 0.09 (±0.05). Broad-sense heritability estimates for third-year height were 0.23 (±0.03), 0.28 (±0.03), and 0.13 (±0.02) at the Florida site, Georgia site, and across sites, respectively. Type B total genetic correlation for Δ was 0.70 ± 0.06, indicating that clonal rankings were relatively stable across sites, while for third-year height, rankings of clones were more unstable across the two trials . Third-year height and Δ were negatively correlated at the parental , full-sib family , and clonal levels, suggesting that genetic variation for Δ in P. taeda may be a result of differences in photosynthetic capacity. We conclude that Δ may be a useful selection trait to improve water-use efficiency and for guiding deployment decisions in P. taeda.     

Genotype by environment interaction for growth of Eucalyptus globulus in Australia

This study assessed the genotype by environment (G × E) interaction for diameter growth in 15 Eucalyptus globulus progeny trials in Australia. Single-site analyses revealed significant subrace and family-within-subrace variance in all trials. Across-site subrace () and family () correlations were estimated by linear mixed model analyses of pairs of trials. Using a factor analytic structure for subrace and family random terms in a multi-environment mixed model analysis, best linear unbiased predictions of subrace effects were obtained for each trial. These were then averaged for each of four states (Victoria, Tasmania, South Australia and Western Australia) and across all sites. Statistically significant G × E interaction was detected, and weighted means across states for and were 0.73 and 0.76, respectively. Nevertheless, the three subraces from the Otway Ranges were both fast growing and relatively stable in their ranks over all sites. We evaluated the sensitivity of subraces to changing environmental conditions, on the basis of random coefficient models regressing subrace performance on selected trial climatic variables. The results suggested differential susceptibility of subraces to water, light and (to a less extent) temperature stresses during summer. Moreover, using multivariate techniques to visualize and interpret the across-site correlation structure for subrace effects, we could identify site clusters of reduced G × E interaction related to soil water availability and evaporative demand during summer. A revised site-type classification using these factors should allow a better capture of genetic gains from breeding and deployment.     

Gains from the clonal and the clonal seed-orchard options compared for tree breeding programs

Summary Gains expected from clonal propagation of selections for plantation from a breeding population were compared with those expected from seed propagation via clonal seed-orchards of selections from the same breeding population. Assumptions were made about numbers of clones selected, size of the breeding population, relative sizes of additive and dominance genetic variance components and time required for various operations. Even when dominance variance is zero, considerable extra gain is obtained by the clonal option over the seed-orchard option; mostly due to the shorter time between selection in the breeding population and field planting. When dominance variance equals additive variance, the advantage of the clonal option due to time saved is approximately equal to the advantage due to genetics (i.e. use of more of the additive variance, use of non-additive variance and greater precision of selection). This means that there is a substantial gain to be made simply by getting superior genotypes into plantations more quickly via the clonal option. The gains obtainable through the use of clonal forestry may also be obtained through seed orchards, but some decades later. In no case was the seed-orchard option superior to the clonal option in terms of the gains obtained. No clonal propagation program can advance without a strong sexually-based breeding program to supply it with improved genotypes. The opportunity for improvement comes from genetic recombination.     

Genetic variation in nine Shorea species (Dipterocarpaceae) in Indonesia revealed by AFLPs

Shorea is the largest and most important genus of the Dipterocarpaceae. The genetic diversity and structure of nine Shorea species from two different locations, namely Nanjak Makmur in Sumatra and Sumalindo in Borneo, were evaluated using amplified fragment length polymorphism (AFLP) markers. A total of 274 trees were investigated at 85 polymorphic AFLP loci. Levels of genetic diversity of these species ranged from  = 0.100 for S. acuminata to  = 0.165 for S. blumutensis. The population of rare species S. blumutensis possessed the highest genetic diversity suggesting that geographically restricted species can have levels of genetic variation comparable to closely related widespread common congeners. Analyses of molecular variance revealed that the genetic variation was mainly found among species in both locations (57.7% in Sumatra; 56.3% in Borneo). The unweighted pairgroup method using arithmetic averages dendrogram of all samples revealed an almost complete separation of species. Thus, AFLP markers proved appropriate for phylogenetic studies of Shorea species. Specific markers have been detected showing high-frequency differences among species and between regions within species. Sequence information of these markers can be used to develop specific polymerase chain reaction markers for wood identification. The possibility of interspecific hybridization was discussed.     

Patterns of genetic parameters for height in field genetic tests of Picea abies and Pinus sylvestris in Sweden

We reviewed the genetic parameter estimates carried out from 1992 to 2006 for height increment in genetic tests of Norway spruce and Scots pine, to describe patterns of genetic variation, heritability, and genetic correlations. The material included seedling and clonal tests in Sweden, aged between 5 and 20 years. Multiple regression was used to explore relationships between parameter values and test environments. Results showed moderate narrow-sense heritabilities ([^(h)]² {\hat{h}^2} : mean =0.29 in Norway spruce; mean =0.23 in Scots pine) that decreased with test site latitude for both species. In Norway spruce, [^(h)]² {\hat{h}^2} increased with better growth and decreased with tree age, while for Scots pine, [^(h)]² {\hat{h}^2} increased with tree age and southward transfer. The additive genetic coefficient of variation (; mean 15%), in Norway spruce, decreased with growth as well as site latitude. in Scots pine (mean =8.5%) increased with southward transfer and more southerly test latitude. Additive and genotypic within-site genetic age-age correlations in Norway spruce were high, with mean r _A and r _G of 0.92 and 0.85, respectively. Corresponding across-sites estimates were on average lower. Genetic parameters were better expressed on favorable sites, at younger ages in Norway spruce and at older ages in Scots pine. The results imply that gain calculations should be based on different parameters in the two species. For maximizing genetic gain in the Swedish breeding program, testing times could be shorter for Norway spruce than for Scots pine. The investigation showed a large variation in parameter estimates from different field experiments, highlighting the importance of testing over multiple sites.     

Genotypic variation in wood properties and growth traits of triploid hybrid clones of Populus tomentosa at three clonal trials

To better understand the genetic control of growth traits (tree height, dbh, and stem volume) and wood properties (basic wood density and fiber length) in triploid hybrid clones of Populus tomentosa, genetic relationships among selected wood properties with growth traits were examined in 5-year-old clonal field trials located in Yanzhou, Gaotang, and Xiangfen, northern China. In total, 180 trees from 10 clones were sampled from the three sites. The site had a moderate effect on basic wood density (BWD), stem wood dry weight (DWT), and tree growth and had a highly significant effect on fiber length (FL) (P?<?0.001). Clonal effects were also significant (P?<?0.05) for all studied traits (except for diameter at breast height (DBH) and stem volume (SV)). Clone × site interaction was significant for all the studied traits except for FL. The estimated repeatability of clonal means for FL (0.91) was higher than for BWD (0.71), DWT (0.62), tree height (0.62), DBH (0.61), and SV (0.55). Intersite genetic correlation estimates indicated that wood properties were more stable than growth traits. Phenotypic correlation estimates between SV and BWD ranged from ?0.29 to ?0.10.     

The potential of aspen clonal forestry in alberta: breeding regions and estimates of genetic gain from selection

Background

Aspen naturally grows in large, single-species, even-aged stands that regenerate clonally after fire disturbance. This offers an opportunity for an intensive clonal forestry system that closely emulates the natural life history of the species. In this paper, we assess the potential of genetic tree improvement and clonal deployment to enhance the productivity of aspen forests in Alberta. We further investigate geographic patterns of genetic variation in aspen and infer forest management strategies under uncertain future climates.

Methodology/Principal Findings

Genetic variation among 242 clones from Alberta was evaluated in 13 common garden trials after 5–8 growing seasons in the field. Broad-sense heritabilities for height and diameter at breast height (DBH) ranged from 0.36 to 0.64, allowing 5–15% genetic gains in height and 9–34% genetic gains in DBH. Geographic partitioning of genetic variance revealed predominant latitudinal genetic differentiation. We further observed that northward movement of clones almost always resulted in increased growth relative to local planting material, while southward movement had a strong opposite effect.

Conclusion/Significance

Aspen forests are an important natural resource in western Canada that is used for pulp and oriented strandboard production, accounting for ∼40% of the total forest harvest. Moderate to high broad-sense heritabilities in growth traits suggest good potential for a genetic tree improvement program with aspen. Significant productivity gains appear possible through clonal selection from existing trials. We propose two breeding regions for Alberta, and suggest that well-tested southern clones may be used in the northern breeding region, accounting for a general warming trend observed over the last several decades in Alberta.     

Genetic transformation of <Emphasis Type="Italic">Gentiana macrophylla</Emphasis> with <Emphasis Type="Italic">Agrobacterium rhizogenes</Emphasis>: growth and production of secoiridoid glucoside gentiopicroside in transformed hairy root cultures

Hairy root cultures of Gentiana macrophylla were established by infecting the different explants four Agrobacterium rhizogenes strains namely A₄GUS, R1000, LBA 9402 and ATCC11325, and hairy root lines were established with A. rhizogenes strain R1000 in 1/2 MS + B₅ medium. Initially, 42 independent hairy root clones were maintained and seven clones belongs to different category were evaluated for growth, morphology, integration and expression of Ri T-DNA genes, and alkaloid contents in dry root samples. On the basis of total root elongation, lateral root density and biomass accumulation on solid media, hairy root clones were separated into three categories. PCR and Southern hybridization analysis revealed both left and right T-DNA integration in the root clones and RT-PCR analysis confirmed the expression of hairy root inducible gene. GUS assay was also performed to confirm the integration of left T-DNA. The accumulation of considerable amounts of the root-specific secoiridoid glucosides gentiopicroside was observed in GM1 ( and ) and the GM2 ( and DNA) type clones in considerably higher amount whether as two but callus-type clones (GM3) accumulated much less or only very negligible amounts of gentiopicroside. Out of four media composition the 1/2 MS + B₅ vitamin media was found most suitable. We found that initial establishment of root cultures largely depends on root:media ratio. Maximum growth rate was recorded in 1:50 root:media ratio. The maximum biomass in terms of fresh weight (33-fold) was achieved in 1/2 MS + B₅ media composition after 35 days in comparison to sixfold increase in control. The biomass increase was most abundant maximum from 15 to 30 days. Influence of A. rhizogenes strains and Ri plasmid of hairy root induction, the possible role of the T_L-DNA and T_R-DNA genes on growth pattern of hairy root, initial root inoculum:media ratio and effect of media composition is discussed.     

Effects of common soil anions and pH on the uptake and accumulation of perchlorate in lettuce

A mechanistic understanding of perchlorate () entry into plants is important for establishing the human health risk associated with consumption of contaminated produce and for assessing the effectiveness of phytoremediation. To determine whether common soil anions affect uptake and accumulation in higher plants, a series of competition experiments using lettuce (Lactuca sativa L.) were conducted between (50 nM) and (4–12 mM), (1–10 mM), or Cl⁻ (5–15 mM) in hydroponic solution. The effects of (0–5 mM) and pH (5.5–7.5) on uptake were also examined. Increasing in solution significantly reduced the amount of taken up by green leaf, butter head, and crisphead lettuces. Sulfate and Cl⁻ had no significant effects on uptake in lettuce over the concentrations tested. Increasing pH significantly reduced the amount of taken up by crisphead and green leaf lettuces, whereas increasing significantly reduced uptake in butter head lettuce. The inhibition by across all lettuce genotypes suggests that may share an ion carrier with , and the decrease in uptake with increasing pH or provides macroscopic evidence for cotransport across the plasma membrane.     

Influence of inbreeding depression on a lake population of Nymphoides peltata after restoration from the soil seed bank

The negative effects of inbreeding depression on fragmented small populations are likely to be expressed more strongly after restoration efforts if regeneration processes have been highly restricted in degraded habitats. We examined the potential influences of inbreeding depression on a population of Nymphoides peltata (Menyanthaceae) restored from the remnant soil seed bank. A hand-pollination experiment demonstrated self-compatibility of a single remaining homostyle genet and significant inbreeding depression in selfed progeny, especially in parameters related to seedling growth (–0.6 for biomass, and –0.4 for relative growth rate). Our genetic analysis indicated that the presumed number of parents contributing to the current soil seed bank was only 2–8 genets and that a single sib-family dominated at each of three sampling sites. The results also showed that the selfed progeny of the homostyle genet were overwhelmingly dominant at two sites (86.8 and 94.7%). As a result, the growth performance of the seed bank seedlings was significantly reduced, to a level as low as that of the selfed progeny. Active restoration efforts to minimize the negative effects of the genetic bottleneck and continuous monitoring based on genetic and demographic study are recommended.     

Use of isotopic signatures to assess the food web in a tropical shallow marine ecosystem of Southeastern Brazil

A dual isotope approach was used to assess the relative importance of terrestrial vegetation detritus and other primary producers in the trophic web of Flamengo Sound (Ubatuba, SP), SE Brazil, surrounded by the Atlantic Rain Forest. Primary producers showed distinct C signatures and the observed values suggest that little terrestrial or bulk sediment organic matter enter the food web of the sound. Suspended particulate organic matter (POM, supports the bulk of the consumers, with some contribution by macroalgae . Consumers C values ranged from −17.4 to . At least three trophic levels were detectable in the food web. The N value of POM was , while that of sediment and detritus was . The N values of suspension feeding benthic invertebrates were 8.2–, deposit feeders 8.3–, and carnivores 10.7–. Values for fishes were for detritivore, 11.4– for benthic feeders, 12.4– for zooplanktivores, and for piscivores/benthic invertebrate feeders. Squid mean value was . There is a reasonable agreement between feeding habits information from the literature and N values from this study. In the sound, the first and second trophic steps seem to be about 1– higher than those of similar organisms studied in temperate waters and this may reflect an input of allochtonous anthropogenic nitrogen enriched in ¹⁵N from human activities.     

Genetic differentiation of Rosa canina (L.) at regional and continental scales

We investigated the genetic differentiation of Rosa canina at regional and continental scales to delineate provenance regions for seed collection. Using RAPD technique, we examined eight populations in Brandenburg (Germany) and 13 samples from Bavaria and Hungary. AMOVA and $\upphiWe investigated the genetic differentiation of Rosa canina at regional and continental scales to delineate provenance regions for seed collection. Using RAPD technique, we examined eight populations in Brandenburg (Germany) and 13 samples from Bavaria and Hungary. AMOVA and -statistics revealed a significant differentiation ( = 0.23, P < 0.0001) between accessions from Brandenburg and the two other regions. The interpopulation differentiation in Brandenburg was low ( = 0.13). We found the greatest interpopulation distance between Hungary and Brandenburg (Hoher Fl?ming, = 0.47) and the lowest between Dahmeland and Hoher Fl?ming (both Brandenburg, = 0.01). UPGMA dendrogram and PCA showed a clear disjunction between Brandenburg and other regions. Due to the outcrossing mating and seed dispersal system, the genetic variation within R. canina populations is moderately high (87% for Brandenburg, 76% for Bavaria and Hungary). The significant genetic differentiation at the continental scale suggests lower gene flow with increasing geographical distance and underlines the usefulness of regional provenances for planting.     

Genetic structure of an insect-pollinated and bird-dispersed tropical tree in vegetation fragments and corridors: implications for conservation

In the vegetation corridors that connect small remnants of undisturbed primary forest in the Lavras landscape (Brazil), Protium spruceanum is a representative of a mass-flowering insect-pollinated and bird-dispersed tree. Allozyme variation was quantified from five forest remnants (N = 150) from secondary vegetation corridors linking them (N = 80) to generate information for genetic conservation. The species adhered to H-W equilibrium in all fragments in most of the loci. The results indicated high gene diversity in the fragments and corridors positively correlated with the plant density (r = 0.742, R ² = 0.551, d.f. = 4). We did not find evidence of inbreeding within fragments nor overall The genetic differentiation among remnants was low Evidence of recent bottlenecks by anthropogenic disturbance was detected in fragments (P < 0.05, Wilcoxon sign-rank test). The minimal viable population was estimated for conservation in situ, indicating fragments with possibilities of maintaining genetic equilibrium diversity in the short term (except F3) and in the long term (only F5). The ratios was also calculated to contribute to vegetation enrichment, area recovery or creation of new vegetation corridors. We found high levels of gene diversity in the vegetation corridors, genetic identity with the fragments and absence of inbreeding. Thus, our results suggest that landscape management strategies should therefore consider both the creation of new vegetation corridors and the protection of extant ones.     

Modeling the HIV/AIDS Epidemic Among Injecting Drug Users and Sex Workers in Kunming, China

This paper presents a mathematical model of the HIV/AIDS epidemic in Kunming,the provincial capital of Yunnan, China. The population is divided into several groups, with individuals possibly changing group. Two transmission routes of HIV are considered: needle sharing betweeninjecting drug users (IDUs) and commercial sex between female sex worker(FSWs) and clients. The model includes male IDUs who are also clients and female IDUs who are also FSWs. Groups are split in two—risky and safe—according to condom use and needle sharing. A system of partialdifferential equations is derived to describe the spread of the disease. For the simulation, parameters are chosen to fit as much as possibledata publicly available for Kunming. Some mathematical properties of the model—in particular the epidemic threshold R ⁰ which determines the goal of public health interventions—are also presented. Though the model couples two transmission routes of HIV, the approximation , with closed formulas for and , appears to be quite good. The critical levels of condom use and clean needle use necessary to stop both the sexual transmission and the transmission among IDUs can therefore be determined independently.     

Long-term changes in nutrient concentrations of the Changjiang River and principal tributaries

We present long-term nutrient data on the Changjiang River (Yangtze River) at six hydrological stations and eight principal tributaries during the period 1958–1985. Three patterns of temporal changes were observed in nitrate and nitrite : minimal variations in the upper catchment area, rapid increases in the middle watershed towards the end of the 1970s, and a gradual increase in the lower drainage basin. Prior to the 1970s, the level of throughout the Changjiang River system remained fairly constant. In the 1980s, however, this changed, with the lowest values in the upper Changjiang changing rapidly to the highest in the middle reaches and then declining slowly but steadily in the lower courses. Compared to and ammonium and soluble reactive phosphorus (SRP) showed smaller increases or no long-term variations, while dissolved silica (DSi) concentration generally decreased at most stations. These three patterns of and changes in the Changjiang River system were reflective of the difference in chemical fertilizer use and landscape features (e.g., slope, soil type and water body area) of the drainage basins of the primary tributaries. The decreases in DSi were most likely attributed to a reduction in suspended sediment loading due to dam constructions and increasing diatom consumption. The increase in and with a reduction in DSi concentrations in the Changjiang River could have significant effects on the stoichiometric balance of nutrients delivered to the East China Sea and the ecosystem in this dynamic region.     

Regional Assessment of N Saturation using Foliar and Root \varvec {\delta}^{\bf 15}{\bf N}

N saturation induced by atmospheric N deposition can have serious consequences for forest health in many regions. In order to evaluate whether foliar may be a robust, regional-scale measure of the onset of N saturation in forest ecosystems, we assembled a large dataset on atmospheric N deposition, foliar and root and N concentration, soil C:N, mineralization and nitrification. The dataset included sites in northeastern North America, Colorado, Alaska, southern Chile and Europe. Local drivers of N cycling (net nitrification and mineralization, and forest floor and soil C:N) were more closely coupled with foliar than the regional driver of N deposition. Foliar increased non-linearly with nitrification:mineralization ratio and decreased with forest floor C:N. Foliar was more strongly related to nitrification rates than was foliar N concentration, but concentration was more strongly correlated with N deposition. Root was more tightly coupled to forest floor properties than was foliar . We observed a pattern of decreasing foliar values across the following species: American beech>yellow birch>sugar maple. Other factors that affected foliar included species composition and climate. Relationships between foliar and soil variables were stronger when analyzed on a species by species basis than when many species were lumped. European sites showed distinct patterns of lower foliar , due to the importance of ammonium deposition in this region. Our results suggest that examining values of foliage may improve understanding of how forests respond to the cascading effects of N deposition.     

Determination of the chemical composition of ulvan, a cell wall polysaccharide from Ulva spp. (Ulvales, Chlorophyta) by FT-IR and chemometrics

The water-soluble cell wall polysaccharides from green seaweeds of Ulva spp. (Ulvales, Chlorophyta), referred to as ulvan, demonstrate composition- and structure-related functional properties. Mid-infrared spectroscopy combined with chemometric techniques was investigated as a means to rapidly predict the chemical composition of ulvan extracts. A calibration was realized with 41 ulvan extracts from two Ulva species. The variables studied included the constituent sugars (rhamnose, xylose, glucose, galactose, glucuronic acid, iduronic acid), protein, and sulfate contents. The correlation between Fourier transform infrared and chemical data was developed using partial least squares (PLS) regression with full cross-validation (leave one out). The coefficients of determination in cross-validation () and the standard error in cross-validation were determined for each variable. The PLS model validation resulted in a coefficient of determination () and a standard error in prediction. Good predictions were obtained for rhamnose ( = 0.9244), xylose ( = 0.8758), glucuronic acid ( = 0.9415), and sulfate ( = 0.9218), which are the main ulvan constituents. However, minor components such as proteins, glucose, galactose, and iduronic acid were not correctly predicted. This study showed that mid-infrared spectroscopy combined with PLS regression is a reliable and fast method for the quantification of the main chemical constituents of ulvan extracts.     

Regional Assessment of N Saturation using Foliar and Root $\varvec {\delta}^{\bf 15}{\bf N}$\varvec {\delta}^{\bf 15}{\bf N}

N saturation induced by atmospheric N deposition can have serious consequences for forest health in many regions. In order to evaluate whether foliar may be a robust, regional-scale measure of the onset of N saturation in forest ecosystems, we assembled a large dataset on atmospheric N deposition, foliar and root and N concentration, soil C:N, mineralization and nitrification. The dataset included sites in northeastern North America, Colorado, Alaska, southern Chile and Europe. Local drivers of N cycling (net nitrification and mineralization, and forest floor and soil C:N) were more closely coupled with foliar than the regional driver of N deposition. Foliar increased non-linearly with nitrification:mineralization ratio and decreased with forest floor C:N. Foliar was more strongly related to nitrification rates than was foliar N concentration, but concentration was more strongly correlated with N deposition. Root was more tightly coupled to forest floor properties than was foliar . We observed a pattern of decreasing foliar values across the following species: American beech>yellow birch>sugar maple. Other factors that affected foliar included species composition and climate. Relationships between foliar and soil variables were stronger when analyzed on a species by species basis than when many species were lumped. European sites showed distinct patterns of lower foliar , due to the importance of ammonium deposition in this region. Our results suggest that examining values of foliage may improve understanding of how forests respond to the cascading effects of N deposition.     

Topological research on diamagnetic susceptibilities of organic compounds

A novel molecular connectivity index, , based on the adjacency matrix of molecular graphs and novel atomic valence connectivities, , for predicting the molar diamagnetic susceptibilities of organic compounds is proposed. The is defined as: , where and E_i are the atomic valence connectivity and the valence orbital energy of atom i, respectively. A good QSPR model for molar diamagnetic susceptibilities can be constructed from and using multivariate linear regression (MLR). The correlation coefficient r, standard error, and average absolute deviation of the MLR model are 0.9918, 5.56 cgs, and 4.26 cgs, respectively, for the 721 organic compounds tested (training set). Cross-validation using the leave-one-out method demonstrates that the MLR model is highly reliable statistically. Using the MLR model, the average absolute deviations of the predicted values of molar diamagnetic susceptibility of another 360 organic compounds (test set) is 4.34 cgs. The results show that the current method is more effective than literature methods for estimating the molar diamagnetic susceptibility of an organic compound. The MLR method thus provides an acceptable model for the prediction of molar diamagnetic susceptibilities of organic compounds. Figure Plot of calculated vs experimental values of molar diamagnetic susceptibilities using the multivariate linear regression (MLR) model (Eq. 8)