Similar results were found for HT soybeans at the time of their introduction

China will likely be one of the first countries in the world to commercialize GM rice. In the United States, the two most widely visible, potentially commercially viable transgenic rice cultivars are Roundup Ready® rice by Monsanto and LibertyLink® by Bayer CropScience . Both are HT varieties—the former is resistant to Roundup® and the latter to Liberty® , both nonselective herbicides able to control a broad spectrum of weeds . Glyphosate is currently registered for rice in California but not widely utilized while glufosinate is not registered [California Department of Pesticide Regulation ]. As such, it is unlikely that local weeds have developed a natural resistance to these chemicals, unlike, for example, bensulfuron methyl . In 1999, LibertyLink® rice cleared biosafety tests by USDA’s Animal and Plant Health Inspection Service but is not commercially available at this time . The primary direct effects of HT transgenic-rice adoption on the cost structure of California rice growers are reductions in herbicide material and application costs and the likely increased cost of transgenic seed. An HT cultivar differs from conventional seed in that a particular gene has been inserted into the rice plant that renders the species relatively unharmed by a particular active chemical ingredient, thus allowing application of broad-spectrum herbicides directly to the entire planting area . This has the potential to simplify overall weed management strategies and to decrease both the number of active ingredients applied to a particular acreage and the number of applications of any one herbicide, blueberry container thus decreasing weed-management costs.

Reduced chemical use provides the major cost saving for growers. Similarly, herbicide application costs per acre depend on the specific chemical involved and the means of application. Typically, application by ground is 60 to 80 percent more expensive than aerial applications . For this study, other pest-management practices and fertilizer applications are assumed not to change with adoption of HT rice. The cost of transgenic rice seed will be greater than that of conventional seed because companies that sell transgenic varieties typically charge a premium to recoup their research investment costs.8 Based on Roundup Ready® corn and soybeans as a reference point, the technology fee is approximately 30 to 60 percent of conventional seed costs per acre . Seed price premiums are in a similar range for Bt corn varieties . In addition to the technology fee, seed costs for transgenic rice will likely change as a result of the California Rice Certification Act of 2000 signed by Governor Gray Davis in September 2000. With the full support of CRC, the CRCA provides the framework for a voluntary certification program run by the industry, offering assurances of varietal purity, area of origin, and certification of non-GM rice . A second, mandatory provision of the CRCA involves classification of rice varieties that have “characteristics of commercial impact,” defined as “characteristics that may adversely affect the marketability of rice in the event of commingling with other rice and may include, but are not limited to, those characteristics that cannot be visually identified without the aid of specialized equipment or testing, those characteristics that create a significant economic impact in their removal from commingled rice, and those characteristics whose removal from commingled rice is infeasible” . Under this legislation, any person selling seed deemed to have characteristics of commercial impact, which would include anytransgenic cultivars, must pay an assessment “not to exceed five dollars per hundredweight.”

This fee is currently assessed at $0.33 per cwt with specific conditions for planting and handling divided into two tiers .10 In addition, the first handler of rice having these characteristics will pay an assessment of $0.10 per cwt . The $0.33 seed assessment is approximately 2.4 percent of average seed costs while the $0.10 fee represents 1.5 percent of average output price. A portion of these assessments is likely to be passed to the grower, depending on the relative elasticities of supply and demand in the seed and milling markets. In addition to generating cost savings, cultivation of HT rice will affect revenues as well. Net returns will be positively correlated with transgenic yield improvements. HT crops are not engineered to increase yields; rather, they are designed to prevent yield losses arising from pest or weed infestation. As such, potential yield gains depend on the degree of the pest and/or weed problem and the efficacy of the HT treatment relative to the alternatives. Many adopters of transgenic corn, cotton, canola, and soybeans have experienced positive yield effects on the order of 0 to 20 percent . However, under more ideal conditions, a yield drag may occur if the cultivar exhibiting the genetic trait is not the highest-yielding variety or if the gene or gene-insertion process affects potential yields . Field tests of LibertyLink® in California have generally found a yield drag of between 5 and 10 percent relative to traditional medium-grain M-202 varieties . To the extent that a yield drag actually exists in the field, it is expected to quickly dissipate over time as a greater number of varieties with the HT trait become available.Another effect of GM rice cultivation on California growers’ returns is the potential development of price premia for conventional medium-grain rice varieties in world rice markets. Despite the predictions and evidence of producer financial benefits from transgenic crops, there is demand uncertainty in world grain markets, especially in the European Union and Japan . Although challenged by many of the major transgenic-crop producing countries , the EU has prohibited imports of new GM crops.

Many other countries have varying GMcrop threshold labeling regulations, including China, Japan, the Republic of Korea, the Russian Federation, and Thailand . These regulations have the potential to ensure that there is some demand for non-GM grain. Due to segregation requirements and the higher unit cost of production of non-GM crops, this introduces the potential for a price premium for non-GM rice. As a result, nonadopters may indirectly benefit from the introduction of transgenic rice. There is good evidence that foreign regulations have affected export demand for transgenic crops, but there is mixed evidence of price premia for traditional non-GM grains. For example, after the United States started growing GM corn, EU corn imports from the United States dropped from 2.1 million metric tons in 1995 to just under 22,000 metric tons by 2002 [USDA, Foreign Agricultural Service 2003b]. Notably, however, the gap in U.S. corn sales to the EU was filled by Argentina, a transgenic producer that only grows varieties approved by the EU . On the other hand, imports of U.S. corn byproducts to the EU have dropped only slightly since 1995 . The U.S. GM soybean export share in Europe has suffered as well, declining by more than 50 percent since 1997 . Price premia exist for non-U.S. corn in Japan and the Republic of Korea, traditional soybeans in Japan, and non-transgenic corn at elevators in the U.S., typically ranging from 3 to 8 percent . However, there is little evidence for price differentials between the GM and non-GM product in the canola market . The global market for rice differs from the market for soybeans in that the majority of rice sold is for human consumption rather than for animal feed. As a result, the market-acceptance issue is likely to be a key determinant of the success of transgenic rice adoption in California . As can be seen in Table 1, the export market for California rice accounts for approximately one-third to one-half of total annual production with Japan and Turkey as the major destinations. California Japonica rice imported by Japan is channeled through a quota system that was negotiated at the Uruguay Round in 1995. Most of California’s rice exports are purchased by the Japanese government and used for food aid and for other industrial uses, including food and beverage processing . Only a small portion of this imported high-quality rice is released into the domestic Japanese market .

Turkey is reportedly attempting to severely restrict imports of transgenic crops through health regulations, despite importing corn and soybeans from the United States , growing blueberries in containers while Japan requires labeling of 44 crop products that contain more than 5 percent transgenic material as one of the top three ingredients . Currently, several varieties of HT and viral resistant rice have entered the Japanese regulatory system for testing but have not yet been approved for food or feed use . As an illustration of potential market resistance, Monsanto suffered setbacks in Japan in December 2002 when local prefecture authorities withdrew from a collaborative study to develop a transgenicrice cultivar after being presented with a petition from 580,000 Japanese citizens . In 2002, China imposed additional restrictions on transgenic crops, including safety tests and import labeling . However, this action may be nothing more than a trade barrier to reduce soybean imports from the United States. In addition, China is worried that introducing biotech food crops may jeopardize trade with the EU. Nevertheless, China is not taking a back seat in transgenic crop research, as it has a major ongoing research program on biotech rice and other crops and is predicted to be an early adopter . There is also some skepticism in the United States with regard to GM crops. Aventis was sued in 2000 over accidental contamination of taco shells by transgenic corn that was not approved for human consumption, resulting in an expensive food recall. The company subsequently decided to destroy its 2001 LibertyLink® rice crop rather than risk its potential export to hostilenations . Kellogg Company and Coors Brewing Company have publicly stated that they have no plans to use transgenic rice in their products due to fears of consumer rejection, and several consumer and environmental groups favor labeling of foods made from transgenic crops . For most food and beverage products manufactured by these companies, however, rice accounts for a small input cost share, resulting in little financial incentive to support GM crop technology. In May 2004, Monsanto announced that it was pulling out of GM wheat research in North America, partly due to consumer resistance. This has important implications for commercialization of GM rice because both grains are predominantly food crops. Many California rice farmers are concerned over the confusion regarding GM crops and do not want to jeopardize export market sales. This fear has been exacerbated by Measure D on the November 2004 ballot in a major rice-producing county that would have prohibited farmers from growing GM crops. A 2001 survey of California growers performed by the University of California Cooperative Extension showed that, of the respondents, 24 percent planned to use transgenic varieties, 37 percent would not, and the remainder were undecided . Of those growers who answered “no,” 78 percent responded that market concerns were a reason. Nevertheless, if profitability at the farm level increases, it is likely that a subset of California producers will adopt the technology . Presumably, those with the most significant weed problems and hence the highest costs would be the first to adopt.UCCE produces detailed cost and return studies for a wide variety of crops produced in California, including “Rice Only” and “Rice in Rotation.” The studies are specific to the Sacramento Valley region where virtually all California rice is produced. Figures on herbicide applications are based on actual use data as reported by DPR and UC Integrated Pest Management Guidelines . The most recent study completed for rice is by Williams et al. and is used as the basis for this study. As the potential adoption of transgenic rice is unlikely to significantly change farm overhead expenses on average, we focus on returns and operating costs per acre as reported in the sample-costs document. However, given weed-resistance evolution, changing regulations from DPR, and changes in the 2002 Farm Bill, the baseline cost scenario is adjusted here to account for changes in herbicide-use patterns, prices of herbicides and rice, and projected government payments. Using information from the 1999 pesticide use report compiled by DPR, the 2001 sample costs assume applications of bensulfuron and triclopyr, both broadleaf herbicides, on 25 and 30 percent of the acreage, respectively, and applications of the grass herbicides molinate and methyl parathion on 75 and 45 percent, respectively, of the acreage. These figures are updated using data from Rice Pesticide Use and Surface Water Monitoring, a 2002 report by DPR, as interpreted by the authors.

The neural activity can be measured using invasive or noninvasive techniques

In the second experiment, a factorial design was employed with one factor being the dosage of JA, and the second being the dosage of ACC. If JA and ACC did not interact to affect a variable, then only the main effects of JA and ACC were considered. If JA and ACC interacted to affect a variable, then the nature of the interaction was determined, and the data were summarized in two-way tables.The Internet of Things is increasingly used by normal people. There will be 50 billion Internet of Thing devices by the 2030. More and more people have started to adopt and use the Internet of Thing devices in everyday life. This thesis aims to explore and study the possibility of implementing and using electroencephalography as controller in the Internet of Thing environment. Also, this thesis intends to study and integrate the human emotion with the Internet of Things Framework . This chapter introduces what Brain Computer Interface is and discusses the components of the Brain Computer Interface. In addition, this chapter explores some of the techniques used to measure the brain activity. Finally, this chapter discusses this research questions of this thesis.Brain Computer Interface is a communication method that depends on the neural activity generated by the brain regardless of peripheral nerves and muscles. BCI aims to provide a new channel of output for the brain controlled and adapted by the user. There are many Brain Computer Interface applications that can be implemented, such as applications for disabled people to interact with computational devices, blueberry container applications for gamers to play games with their thoughts, social applications to capture feelings and emotions, and application for human brain activities.

It is a medical imaging technique used to capture high quality pictures of the anatomy and physiological processes of the body. It uses powerful magnet, radio waves, and field gradients to generate images of the body. It is non-invasive, painless and does not use radiation . It can provide very detailed high resolution images of a body parts. In particular, it can capture very detailed high resolution images for the brain compared to other imaging techniques such as CT and X-ray because of its ability to differentiate between soft tissues of the body. However, due to the magnet effects, metallic items are not allowed during the scan which because they limit its applications.It is a special MRI technology that measures brain activity by detecting changes associated with blood flow. This technique relies on coupled cerebral blood flow and neuronal activation. The blood flow increases in a region when this region is in use. The idea of this technique lies in the amount of oxygenated and deoxygenated hemoglobin changes in the blood flow during the neural activity. The most common one is Blood Oxygenation Level Dependent fMRI which measures the ratio of Oxy-Hb to Deoxy-Hb in order to measure the oxygen consumption of active neurons. It is also invasive and has an excellent spatial resolution compared to EEG, and records signals from all the brain regions. However, this technique has the same limitations of the MRI technique.It is the technique used to measure the magnetic field over the head generated by the electric current in the brain. The most commonly used technology of MEG currently is SQUIDs. This technique allows capturing MEG of the head efficiently and rapidly. Also, this technique is non-invasive and can be used as complement for other techniques such as EEG and fMRI. Due to the fact that MEG uses magnetic fields,this technique makes less distortion than the electric fields. However, the same restriction applied on fMRI and MRI can be applied to MEG due the to its sensitivity for ferromagnetic.

An electroencephalogram is a method monitoring the electrical activity of the brain using small flat metal discs placed on the scalp. EEG measures voltage fluctuations resulting from brain cells communications via electrical impulse. In Particular when neurons are activated, ions such as Na+ , K+ and CI– are pushed through the neuron membrane. EEG is a weak signal and needs to be amplified in order to be displayed or stored on a computer. Two approaches to recording the EEG signal are invasive and non-invasive. In the invasive approach, the electrode is implanted inside the human brain, which requires surgery. In a non-invasive approach, electrodes are placed on the surface of the skull, which have many benefits such as risk free, easy setting, and repeating measurement. In addition, it is more favorable in developing and designing application for normal people. The focus of this thesis will be based on this non-invasive EEG technique .The first question of this thesis is how to integrate the low-quality cheap EEG headset, which has only one electrode located in forehead, with an Internet of Things framework. In order to do this we have to first build and design the EEG Server which is able to translate EEG signals into commands. Then, we must build an algorithm that construct different patterns from these commands, and these patterns will be used to control different Internet of Things devices. The expected outcome after integration and build, the different EEG pattern is the ability to control Internet of Thing devices such as Light turning on/off, music playing and etc.The second question of this thesis is how to build the EEG Edge which is able to classify between eye close and eye open states. In order to answer this question we will use an extension of the Internet of things framework that supports intelligent edge, which is presented in [38].

So, in order to build the EEG Edge we need to extract EEG features for different subjects and build the model that is able to classify between open and close eye states. There are different types of features that could be extracted from EEG raw signal. However, for this application we need only to extract the power spectrum density features. Lastly, we need to define the feature extraction extension which will contain the EEG features and define the execution extension which will contain the classifier model. The expected outcome after the integration will be the ability to classify eye states on the edge.The third question of this thesis is how to build a model that is able to detect and classify positive and negative emotions. In order to classify the emotions, different factors must be considered, which include participants, stimuli, the temporal window, and EEG features. Different EEG features will be extracted from EEG raw signal and these features include time domain features, frequency domain features, and nonlinear features. Different video clips will be used as stimuli in order to trigger different emotions. The expected outcome will be the ability to classify two different types of emotions, positive and negative emotions.EEG is the electrical activity measurement in the brain. The first measure for EEG was recored by Has Berger in 1924 using galvanometer. Based on the internal brain behavior or external stimulus, EEG varies in amplitude and frequency of the wave. The system contains a EEG headset, and this thesis used ”NeuroSky Mindwave Mobile” which is using Bluetooth connection to transfer the EEG signal. The EEG receiver records and receives the EEG signal coming from a EEG headset which is written in Python. I used the Wukong framework to deploy WuClass for the EEG, and WuClass for a controller on Intel Edison and Raspberry Pi. Figure shows the system architecture will be used in this thesis.There are a lot of commercial EEG headsets from the simplest ones to the more sophisticated one. Table compares different EEG headset. These different EEG headsets are able to capture different mental states, and different facial expressions. Both emotiv headsets, the EPOC+ and INSIGHT, capture excitement, frustration, engagement, meditation and affinity. Also, Emotiv headsets capture EEG bands which are Delta, Theta, Alpha, Beta, and Gamma. In addition, Emotiv headsets capture some facial expressions such as blinking, smiling, clenching teeth, growing blueberries in containers laughing and smirking. On the other hand, the Neurosky Mindwave Mobile is limited to capturing only two mental states which are meditation and relaxation. Finally, the Muse headset can capture positive and negative emotions. Also, the Muse headset captures EEG bands which are Delta, Theta, Alpha, Beta, and Gamma. In addition, the Muse headset also captures some facial expressions such as jaw clenching andeye blinking. Emotive EPOC+ sensors use saline soaked felt pad technology, and emotive INSIGHT sensors use long-life semi-dry polymer technology. Neurosky Mindwave Mobile and Muse sensors use long life dry technology.NeuroSky Mindwave Mobile consists of eight parts which are ear clip, ear arm, battery area, power switch, adjustable head band, sensor tip, sensor arm, and think gear chip. The operation of this device is based on two sensors to detect and filter EEG signals. The sensor tip on the forehead detects the electrical signal from the frontal lobe of the brain. The second sensor is an ear clip which is used as ground to filter out the electrical noise. Figure shows NeuroSky Mindwave Mobile and Figure shows the electrode position of NeuroSky Mindwave Mobile.

This thesis uses NeuroSky Mindwave Mobile for many reasons. First, this project aims to offer a low-cost system, which can be used by everyone. Second, NeuroSky Mindwave Mobile is highly resistant to noise and its signal is digitized before it is transmitted throughBluetooth. Third, NeuroSky Mindwave Mobile offers unencrypted EEG signal compared to the Emotive and Muse which are encrypted.It is NeuroSky algorithm to characterize mental states. This algorithm applied on the remaining signal that is acquired from removing the noise and the muscle movements of the raw brain wave signals. Two eSense signal are produced as a result of this algorithm: attention and meditation signals. These signals detect the concentration and relaxation of subject. The values of these signals range from 0-100 in which zero indicates low in concentration or in low in relaxation, and 100 indicates high in concentration or high in relaxation.One major limitation is the accuracy of the EEG signal captured by NeuroSky Mindwave Mobile, because the NeuroSky Mindwave Mobile Mobile has only one electrode which is FP1. The problem with FP1 is its susceptibility to a lot of noise coming from eye movement and muscle movement. Another possible issue is comfort. One subject claims it is uncomfortable to wear. This is most likely due to the rigid headband design as well as the need for the ear clip as the reference sensor.In order to detect eye blinking, OpenCV library is used which is an open-source library of programming functions aimed to offer real-time computer vision. This library is used to detect eye blinking to trigger the system, hold on to a certain state, and change between different states. The algorithm used to classify blinking is Haar Cascades classifier which is a machine learning approach in which the cascade function is trained by a lot of negative and positive pictures, then used to detect objects in other image. In order to obtain EEG signal from Neurosky Mindwave Mobile I used an open-source API written in Python suggested by the Neurosky company. Two major libraries used are bluetooth headset.py and parser.py. The bluetooth headset library contains methods to connect the Mindwave Mobile to the computer via Bluetooth either by specifying a MAC address or by automatically searching for a device named ”Mindwave Mobile”. If it does not find a device automatically or with the MAC address specified, it will raise an error. The other library parser.py is specific to NeuroSky Mindwave Mobile device. There are two major classes in this library: ThinkGearParser and Time Series Recorder. It must first create a new Time Series Recorder object and then include this object into a Think Gear Parser object, which will package the EEG information and be able to display the data on the computer. The other important library that is necessary for WuKong integration is a socket which creates a new socket that points to the exact IP address of the peripheral device . In order to control a device Wuclass has to be designed. For this project, EEG Server is Wuclass that is able to receive the EEG signal and transform to different actions. Triggering an action will take around 10 seconds as shown in Figure . Figure explains and shows the flow of data and control of the system in Flow based program .

Trees on dwarfing rootstocks are smaller than those on standard or seedling rootstocks

It is noteworthy that at harvest only two transcription factors were differentially expressed, both showing higher expressions in T fruits and in the case of the ortholog of PAP2/IAA27, also at 1 week of cold storage . SlIAA27 silencing results in greater auxin sensitivity in tomato. Moreover, a gain-of-function mutation in IAA16 confers poorer responses to auxins and ABA in Arabidopsis. Thus, it is likely that high levels of these genes at harvest contribute to delay the ripening program or protect fruits LS during cold storage, at least at the beginning of cold storage. The analysis of the expression profiles during cold of the genes differentially expressed in M fruits resulted in important and unexpected expression characteristics. In fruits LS, these genes behaved like ripening genes and were able to continue with the ripening program in the cold in fruits LS, while the ripening expression of other ripening genes was normally halted , which is not the case of high sensitive fruits. The ability of cold to stop fruit ripening has been previously reported, even if no details of how this happens at the molecular level have yet been provided. Although we have no hypothesis about why these genes continued with the ripening program in the cold , we believe that this may be because these genes are part of the adaptation mechanism or simply reflected that LS fruits perform better in the cold than S fruits. In apples the ability to set up ripening during cold seems to be an adaptative mechanism to shorten ripening time in colder autumns. On the other hand, 25 liter square pot this unexpected behavior of some of the genes differentially expressed at harvest indicates that they not only can form part of a mechanism for the interaction between endogenous and exogenous signals, they could also be able to contribute to mealiness in response to cold stress.

In light of this, it is interesting to remember that environmental/ripening stage/cultural preharvest practices have a strong effect on CI sensitivity during the post harvest which, together with the genetic background, may be responsible for the differences noted in the M stage that condition the cold response.Fruit trees differ from landscape trees in that they are best kept relatively small to facilitate routine pruning, fruit thinning, managing pests, and harvesting fruit from the ground or a ladder. Fruit trees that are allowed to grow above a manageable height produce excessive fruit, leading to branch breakage, smaller-size fruit, and, in some cases, pest problems . Most fruit trees are trained to the open center system and are topped annually to reduce limb breakage. However,some fruit trees lend themselves to central leader training, so these considerations are less important. The major problem with those that do grow very tall, however, is that the fruit are borne higher in the tree and the lower branches become shaded. This results in the decline of these branches and ultimately renders them fruitless. Unlike fruits, nuts are knocked, shaken, or allowed to fall. They are not usually picked by hand, so tree height is less important. Some size control is necessary for preventing branch failures and maximizing nut production, because large trees are more difficult to knock. Pruning of nut trees generally consists of thinning or cutting back selected branches to suitable lateral branches. Walnut and pistachio trees should be trained to a modified central leader to maximize fruit production and maintain a branch structure that can support the nut crop .The best strategy for keeping trees relatively small is to use a dwarfing rootstock when available. However, semidwarf rootstocks differ in their ability to cause dwarfing, and many semidwarf trees sold in retail nurseries are only slightly dwarfing. For example, apple rootstocks can range from about 80 percent of the size of a standard tree to about 60 percent to about 30 percent . Therefore, some semidwarfs are still, practically speaking, full-size trees. Other fruit species do not have this range of dwarfing rootstocks available, and most are only slightly dwarfing. For the stone fruits, such as peaches and nectarines, the dwarfing rootstock most commonly available is Citation, which produces a tree that is somewhat smaller.

Citrus can be dwarfed to approximately 50 percent by growing trees on the ‘Flying Dragon’ trifoliate orange rootstock. However, availability of trees grafted to ‘Flying Dragon’ is limited due to the very slow growth of grafted trees. Genetic dwarf trees are very easy to manage and are aesthetically pleasing in the landscape . They naturally produce short internodes and are usually planted on standard rootstocks. A limited number of genetic dwarf varieties are available for almond, apple, apricot, nectarine, and peach.Deciduous fruit and nut trees are ideally planted bare root, but containerized trees can also be used. All bare-root trees intended solely for fruit or almond production should be headed 18 to 24 inches above the ground at planting to force low branching; walnuts and pecans should be headed higher. If this were not done, the first laterals would typically form around 5 to 6 feet above the ground, growth would be weak, and much of the fruit would be out of reach from the ground. It is important to develop a new leader in headed trees if the central leader method is to be maintained. Select one of the shoots that grow near the heading cut, and tie it to a stake in an upright position if it is not growing upright naturally. Additional pruning may be needed to eliminate branch crowding or prevent codominant trunks from forming. Higher branching may be desirable for fruit trees in some urban settings to allow for maintenance of vegetation under and around the tree. In areas where deer are a problem, lower branching may not be practical without proper protection. Containerized fruit trees are often planted in spring or summer, so they cannot be headed without removing all the foliage. Either leave the tree as it was headed in the nursery, or make the lower heading cut in the next dormant season. In hot regions where afternoon sun hits the trunk, apply a white interior latex paint diluted 50-50 with water to the trunk to prevent sunburn injury.Each fruit and nut species has a preferred training method based on the species’ growth habits and fruiting characteristics. Ideally, most of the fruit should be produced low in the tree to facilitate fruit thinning, pest management, and harvest. Because fruit is produced on spurs or 1-year-old branches that require sunlight for flower development, direct sun must penetrate into the lower portions of the canopy for fruit production low in the tree. Nearly all fruits and nuts are borne mainly on spurs, but peach and nectarine fruit are borne only on 1-year-old shoots that grew the previous summer. Most stone fruits and almonds are best trained to the “open center” or “open vase” method, where the center of the tree is routinely kept free of vigorous shoots. In this manner, lower fruiting branches receive sufficient light through the tree’s center. Apples and pears can also be trained to the open center system but are better adapted to central leader training, where lateral branches are trained outward from a vertical leader allowing sunlight to penetrate from the sides.

Persimmons are also well adapted to central leader training. For apples and pears it may be prudent to develop two or three leaders in case fire blight kills one of the leaders; lateral branches are directed to the outside of the tree. Apples and pears can also be espalier trained. This method involves pruning the tree to form a narrow, gallon pot flat plane on a trellis or against a wall or fence. Permanent, horizontal branches that resemble cordons on grapevines are selected to produce the fruiting spurs. Walnuts, pistachios, and persimmons can also be initially trained with a central leader, but then trees are allowed to develop a natural rounded crown; this method is referred to as “modified central leader” training. Fig trees can also be trained using this method or using open center training, and they can be kept fairly short or allowed to grow tall, or they can even be espalier trained. Modified central leader training can be used for pomegranates, but their rangy growth and constant root suckering make them better adapted to a system that allows them to grow into a large multit-runked bush. Pruning them typically involves heading back and thinning vigorous upright branches and removing old trunks or scaffold branches to rejuvenate trees. Citrus trees can be allowed to grow with little training except to eliminate scaffold branches with narrow crotch angles. Manage water sprouts by heading, shortening to a lower lateral, or, in some cases, completely removing. Save water sprouts that bend over, ultimately contributing to the typical mounding citrus canopy. Remove all root stock suckers; they often grow up the center and are difficult to see. Painting the trunk white not only helps to prevent sunburn but will also make root stock suckers easier to spot. Over time, the shaded inner fruiting branches of citrus trees die, and fruit production moves to the top and sides. This characteristic is considered acceptable for citrus. Citrus trees can also be hedged, and they are very adaptable to espalier training. Nearly all species can be trained as fruit bushes, using a method in which trees are trained in the first and second year by heading shoots when they reach about 2 feet in length. The resulting new shoots are headed again, and this is followed each time by some thinning of shoots as well. Once the desired tree height is achieved , pruning consists of removing shoots above the desired height and thinning remaining shoots and branches about twice a year. Nearly all pruning on fruit bushes should be done in the growing season to reduce their vigor, but touch-up pruning is useful in winter, when branch structure is more visible. Apricots and cherries should be pruned only in late summer, when dry weather is predicted for an extended period. These species are susceptible to branch canker diseases, caused mainly by Eutypa and Botryosphaeria fungi, which infect branch injuries made before or during wet weather or periods of very high humidity. However, most other fruit and nut trees can be pruned any time enough leaves have fallen to make the tree structure visible.After planting and heading bare-root stone fruit trees, two options are possible for the resulting shoots. Either they are allowed to grow through the summer, or training can begin during the first growing season by selecting three or four well-placed shoots when they are 1 to 3 feet long and heading back all other shoots to 4 to 6 inches. By winter, the primary scaffold branches are selected and headed, and all other upright branches are removed. Continue to develop the tree to a vase shape over the next 2 years. Ideally, each primary scaffold should branch into two secondary branches, which, in turn, branch into two tertiary branches. Prune out vigorous upright shoots in the tree’s center in winter, and maintain the open center by removing vigorous upright shoots once during summer. In winter, thin fruiting branches to reduce fruit load and minimize the need for fruit thinning. Head the trees back to about the same height every year, preferably to a height that can be reached using only a short ladder . With almond trees, scaffold branches are selected as with stone fruit trees. After that, however, the opening in the center of the tree can be somewhat narrower than stone fruit trees. Annual pruning involves thinning branches to avoid overcrowding.Apples, pears, persimmons, and pecan trees are best trained to develop a central leader similar to those in many shade and ornamental trees . Lateral branches grow outward from the leader, either in tiers of approximately four branches each or spaced fairly uniformly up and around the trunk. Rather than simply allowing the trunk to continue growing naturally after planting, the trunk is headed at about 18 to 24 inches above the ground, and the most vigorous and upright shoot that develops is selected to become the new leader. This practice is done to force the first tier of four lateral branches below the heading cut. When the new leader has grown about 2.5 feet, it is then headed back about 6 inches.

Fleshy fruit are a relatively recent evolutionary innovation

This family has been delimited into four ‘‘Classes’’, and 4 of the 17 members of the Class I AtHBs have been shown to be involved in ABA responses across diverse tissues . In addition, the expression of three AtHB6, 7, and 12, have been show to be up-regulated by ABA . An examination of the grape genome identified 10 orthologs that cluster with the Class 1 HBs . The PP2C protein phosphatases represent another large gene family being made up of 80 genes in Arabidopsis . Within this family group a cluster of genes containing many genes that have been characterized to function in the ABA-signaling pathway; most notably, the ABA-insensitive mutants, abi1 and abi2 . In addition, AtPP2C-A, AtHAB1, AtHAB2, and AtAHG1, also members of Group A, function in ABA signaling across diverse tissues . In Arabidopsis, all the members of this group are induced by ABA treatment. In grape, nine VvPP2Cs clustered in group A . The WRKY transcription factors are a large gene family, consisting of approximately 70 members making up three groups in Arabidopsis . The barley HvSUSIBA2, and AtWRKY2, and AtWRKY34, all fall within the same group consisting of 14 members in Arabidopsis . HvSUSIBA2 modulates the expression of a barley isoamylase gene during seed development via the binding of SURE elements . In addition, Sun et al. demonstrated that expression of HvSUSIBA2 is induced by exogenous sugar and its native expression profile during seed development correlates strongly with endogenous sucrose levels.

Hammargren et al. found that the sugar responsiveness of a nucleoside diphosphate kinase is altered in Atwrky2 and Atwrky34 mutant backgrounds. Grape contained 13 putative orthologs that fall within this group .Expression profiling was carried out in berry skins of fieldgrown Cabernet Sauvignon in order to identify those orthologs expressed during ripening. In addition, black plastic planting pots expression profiles under both control- and deficit-irrigated conditions were compared in order to identify orthologs whose expression pattern reflected the advancement of ripening under ED. Water deficits were applied continuously from fruit set until to the onset of ripening, resulting in an average difference in midday leaf water potential of 0.36 MPa before the onset of ripening and no difference during ripening . Of the 67 orthologous genes identified, 38 were expressed in grape berries. A summary of the expression profiles of all the genes examined in this study can be found in Figs. 4 and 5 while more detailed expression data is contained in Suppl. File 1. The majority of these genes were differentially regulated during berry development, with 26 exhibiting statistically significant changes with time in control and/or ED . There were few statistically significant differences in the magnitude of expression between control and ED . Six genes exhibited statistically significant differences between control and ED. Four of these instances, VvHB8, VvSnRK5, VvPP2C-3, and VvPP2C-7, all exhibit elevated levels of expression in ED at or just prior to the onset of ripening . Eight of the VvWRKYs selected for analysis were expressed in ripening grape. VvWRKY3, 5, and 6, were all differentially expressed during ripening and exhibited similar patterns of expression . They were up-regulated ranging from 4- to 16-fold at the onset of ripening. There were no significant differences in the expression of VvWRKY1, 2, 16, 18, and 19 across development or during water deficit . Of the 10 Class I, VvHB orthologs only four are expressed in fruit during ripening . Both VvHB4 and 8 were strongly up-regulated at the onset of ripening, exhibiting increases of [16-fold. VvHB8 is up-regulated much earlier under water deficit and high levels persist until late in ripening. VvHB4 is down-regulated early in development under ED. VvHB2 expression in controls generally decreased during development with a small up-regulation at the onset of ripening, although these changes were not statistically significant.

Under ED, however, this pattern of expression is more pronounced with a sharp eight fold decrease in expression at 81 DAA . VvHB3 was constitutively expressed during ripening with no significant changes over time or under ED. Six of the VvPP2Cs were expressed in grape berries. VvPP2C-3, 6, 7, and 9 were all differentially expressed during ripening, while VvPP2C-1 and 5 expression did not change significantly . VvPP2C-3, 6, 7, and 9 were all up-regulated strongly at the onset of ripening increasing as much as 16-fold. VvPP2C-3 and VvPP2C-7 expression were clearly induced earlier and to higher levels in ED. Both exhibited statistically significant two to fourfold greater levels of expression in green berries at the onset of ripening.In order to more directly test the effects of sugar and ABA on the onset of ripening, immature Cabernet Sauvignon berries were harvested from the field at 61 DAA and cultured in the presence of various combinations of sucrose and ABA until 84 DAA, a period of 23 days. The onset of ripening in the clusters from which these berries were collected in the field occurred at approximately 73 DAA, therefore the cultured berries were collected approximately 12 days prior to the onset of ripening. Ripening phenomena were induced in berry culture when treated with sucrose and ABA as evidenced by changes in color, softening, and gene expression. Berries treated with 10% sucrose and various ABA concentrations changed color while those treated with 2 or 10% sucrose alone remained green . 200 lM ABA and 2% sucrose ? 200 lM ABA treatments were included in our analyses but yielded no results because of a phenomenon where the berries exploded reproducibly . On average, cultured berries gained weight over the culturing period . Sucrose treatments of 2 and 10% showed the greatest weight gains corresponding to gains of 21 and 8%, respectively. Berries cultured in the presence of 10% sucrose with the addition of various ABA concentrations showed average weight gains of approximately 4%. Previous studies have found that a precipitous drop in grape berry elasticity occurs just prior to the onset of ripening in grape . In our current culture experiments, berry elasticity remained equal to that at T0 in the 2 and 10% sucrose treatments, while sharply decreasing with ABA treatment . Finally, the grape Myb transcription factor VvMybA1 was utilized as a molecular marker for the onset of ripening. VvMybA1 is responsible for activating anthocyanin biosynthesis and has a distinct pattern of expression; being absent prior to the onset of ripening at which time it is strongly up-regulated . In berry skins, VvMybA1 expression was completely absent from the 2 and 10% sucrose treatments and was strongly up-regulated in the 10% sucrose ? ABA treatments as expected . We hypothesized that orthologs of gene families regulated by sugar and ABA, whose expression was strongly up-regulated at the onset of ripening and advanced under ED, would be regulated similarly by sugar and/or ABA in cultured berries. To test this, changes in the expression of VvHB4, VvHB8, VvPP2C-3, and VvPP2C-6 were investigated in skins of cultured berries. In the field, all genes were strongly up-regulated at the onset of ripening and advanced under ED and in berry culture, expression was strongly induced in the presence of 10% sucrose ? ABA when compared with treatments of 2 and 10% sucrose alone . Among those genes analyzed, the magnitude of induction in the field versus in berry culture was variable. For example, when data from Fig. 7 was expressed as fold change both VvHB4 and VvPP2C-3 were induced 10-fold from 57 to 74 DAA in the field compared to 6- and 40-fold in culture, respectively.The transcriptional data in this study demonstrate that numerous sugar and ABA-signaling orthologs are expressed during ripening in grape, black plastic pots for plants and identify novel candidates in the control of non-climacteric fruit ripening. Several genes exhibited patterns of expression correlating with sugar and ABA accumulation at the onset of ripening in field-grown fruit. Changes in color, softening, and gene expression analogous to the onset of ripening in the field were induced in berry culture when treated with sucrose and ABA, demonstrating their role in controlling the onset of ripening. This study shows that many orthologous sugar and ABA-signaling components are regulated in fleshy fruit similar to their regulation originally characterized in model systems across diverse processes.These genes are easily delimited through nesting the currently available grape sequences comprising a gene family within their corresponding family in Arabidopsis. However, the current grape genome assembly, and its annotated proteome, certainly does not identify all the genes present in the grape genome so our analyses most likely failed to identify some orthologs. In the current study, we chose to use QPCR in our expression analyses instead of the current grape microarray for several reasons. First, the majority of the genes analyzed here are not present on the current Affymetrix Vitis vinifera gene chip since the chip was derived from ESTsand the present study utilizes the complete genome. Second, microarrays suffer from several limitations one of which is that they are particularly insensitive in quantifying low abundant transcripts .

Transcription factors are a largely low abundance transcripts, and a study in Arabidopsis comparing QPCR and Affymetrix microarray approaches found that the microarray could detect \55% of 1,400 transcription factors tested, compared to [85% via QPCR . Furthermore, cross-hybridization is common on microarrays, which is especially problematic when considering conserved gene families like those examined here. However, those genes represented on the chip were identified and expression profiles were compared with those determined via microarray analyses in two recent studies. Notably, Koyama et al. demonstrated via microarray that a VvHB transcription factor , identical to VvHB8 in this study, is also up-regulated at the onset of ripening and in response to exogenous sugar and ABA. For several other genes, expression profiles in the current study are nearly identical to those found by Deluc et al. . Many grape orthologs of genes shown previously to be modulated by sugar and/or ABA in model systems exhibited expression patterns during ripening, and in response to water deficit, consistent with their modulation by sugar and/or ABA in grape. More specifically, these genes are induced at the onset of ripening and induced earlier, and to higher levels, under water deficit. This characteristic pattern of expression is shared with many flavonoid pathway genes , and these correlations in expression throughout ripening suggest common regulatory mechanisms. Experiments in berry culture demonstrated that several sugar and ABA-signaling orthologs, and VvMybA1, a transcriptional activator of anthocyanin biosynthesis, are up-regulated by exogenous ABA in the presence of high sucrose. These data suggest that sugar and ABA play a predominant role in regulating the expression of a suite of genes at the onset of ripening, including those responsible for anthocyanin biosynthesis and components of their own signaling pathways. These results have interesting evolutionary implications demonstrating that some orthologs are consistently regulated by sugar and ABA across diverse developmental processes. Land plants, in general, have undergone abundant gene duplication through their evolutionary history , although there is debate on the exact nature and timing of these events across angiosperms and in grape specifically . Gene duplication is considered cornerstone to providing the raw material for evolution. A duplicate gene, now no longer essential, can undergo changes in its structure and/or regulation allowing for it to take up a novel role. The results of this study show that some of the Group A PP2Cs and Class I HBs have maintained their ABA responsiveness during fruit ripening in grape. At least with regard to ABA responsiveness, the nature of regulation has been conserved, but co-opted into a completely different developmental context. This may provide for the discovery of novel cis-regulatory elements through promoter sequence comparisons across species. In grape, advances in elucidating molecular mechanisms suffer from a lack of transgenic and related technologies on which most reverse genetic studies are based. This study demonstrates that model systems can provide fundamental knowledge and insight into function of other agronomically important plant species even in extremely divergent developmental processes. Equally, this suggests that Arabidopsis, with its wealth of tools available for facilitating reverse genetic studies, may provide a valuable system to characterize genes of interest from grape or other crop species. Already there are several examples of the successful characterization of grape genes in Arabidopsis and tobacco . This could prove especially useful considering the limitations of functional genetic analyses in perennial fruit crops where long propagation times areprohibitive. Future studies should include attempts to complement ABA and sugar signaling mutants in Arabidopsis and other model species with their orthologs implicated in ripening.

Connect the coarse positioner control cable to the cryostat

At high tuning fork amplitudes, interactions between the nanoSQUID tip and the surface can produce local variations in oscillation amplitude and appear as parasitic signals at the tuning fork frequency. Of course, the nanoSQUID is highly sensitive to local temperature, so systems with thermal gradients will generally have backgrounds associated with that. But by far the most important parasitic contrast mechanism in the nanoSQUID campaigns discussed here is electric field contrast through parasitic Coulomb blockade.Below I have included a set of instructions for execution of a nanoSQUID magnetic imaging campaign using the instruments in Andrea Young’s lab. It may be useful if you are operating or building a nanoSQUID microscope in a different lab, but I would like to emphasize that the instructions below are merely sufficient for getting the nanoSQUID sensor to a sample, they are almost certainly not optimized for expediency. I’m sure that as the technology matures many steps will be rendered superfluous. Of course, if you’re using the nanoSQUID to study a bulk material and not a microscopic heterostructure, navigation is not necessary and you will be able to skip most of these instructions. A nanoSQUID imaging campaign can begin when the microscope is cooled down and all of the necessary systems are operational. One must check that: The tuning fork has a good resonance with Q 1000, plant pot with drainage the phase-locked loop inside the Zurich lock-in amplifier locks, and you can find an AC tuning fork excitation at which clicking “Set PLL Threshold” produces a 0.25 Hz standard deviation.

If you are in Andrea Young’s lab and not some other institution running a nanoSQUID microscope, remember that this custom tuning fork amplifier needs 5 V, not 15 V, unlike most of our custom electronics. The tip has been characterized and is a SQUID. Sensitivity is good enough that magnetic field noise is 25 nT/rtHz . The SQUID interference pattern looks reasonably healthy and corresponds to a diameter that is close to the SEM diameter . It is important to remember that it is possible for the Josephson junctions producing nanoSQUIDs to end up higher on the sensor. These might produce healthy SQUIDs but will not be useful for scanning, and discovery of this failure mode comes dangerously late in the campaign, so SQUIDs high up on the pipette are very destructive failure modes. This failure mode is uncommon but worth remembering. If you have access to a vector magnet, such SQUIDs also usually have large cross sections to in-plane magnetic flux, and this can be useful for identifying them and filtering them out. The capacitances of the Attocube fine positioners are = µF. These scanners have a range of µm. They creep significantly more than the piezoelectric scanners used in most commercial STM systems, but their large range is quite useful. Damage to the scanners or the associated wiring will appear as deviations from these capacitances. Small variations around these values are fine. After you are done testing these capacitances, reconnect them. Make sure you’re testing the scanner/cryostat side of the wiring, not the outputs of the box- this is a common silly mistake that can lead to unwarranted panic. If you’re working in Andrea Young’s lab, make sure the Z piezo is ungrounded . If for whatever reason current can flow through the circuit while you’re probing the capacitance, you will see the capacitance rise and then saturate above the range of the multimeter.

Because the nanoSQUID is a sharp piece of metal that will be in close contact with other pieces of metal, it sometimes makes sense to ground the nanoSQUID circuit to the top gate of a device, or metallic contacts to a crystal, to prevent electrostatic discharge while scanning or upon touchdown. If you have decided to set up such a circuit, make sure that the sample, the gates, and the nanoSQUID circuit are all simultaneously grounded. If you forget to float one of these circuits and bias the SQUID or gate the device, you can accidentally pump destructive amounts of current through the nanoSQUID or device. However, you must make sure that the z piezoelectric scanner is not grounded. You can now begin your approach to the surface. You should ground the nanoSQUID and the device. If you are in Andrea’s lab, verify that the three high current DB-9 cables going from the coarse positioner controller box to the box-to-cable adapter are plugged in in the correct positions. The cables for each channel all have the same connectors, so it is possible to mix up the x, y, and z axes of the coarse positioners. This is a very destructive mistake, because you will not be advancing to the surface and will likely crash the nanoSQUID into a wirebond, or some other feature away from the device. The remaining instructions assume you are using the nanoSQUID control software developed in Andrea’s lab, primarily by Marec Serlin and Trevor Arp. The software is a complete and self-contained scanning probe microscopy control system and user interface based on Python 3and PyQT. Open the coarse positioner control module. Click the small capacitor symbol. You should hear a little click and see 200 nF next to the symbol . The system has sent a pulse of AC voltage to the coarse positioners; the click comes from the piezoelectric crystal moving in response. Check that you see a number around 1000 µm in the resistive encoder window for axis 3 . Note whether you see a number around 2000-3000 µm in the windows for axis 1 and axis 2. If you are in Andrea’s lab, it is possible that you will not for axis 2.

Axis 2 has had problems with its resistive encoder calibration curve at low temperature. The issue seems to be an inaccurate LUT file in the firmware; new firmware can be uploaded using Attocube’s Daisy software. It is not a significant issue if you cannot use the axis 1 and 2 resistive encoders; however, it is critical that there be an accurate number for axis 3. Set the output voltage frequency to be somewhere in the range 5-25 Hz . Set the output voltage to 50 V to start . Make sure that the check box next to Output is checked. Move 10 µm toward the sample . If Axis 3 doesn’t move, don’t panic! It’s usually the case that the coarse positioners are sticky after cooling down the probe before they’ve been used. Try moving backwards and forwards, then increase the voltage to 55 V, then 60 V. Once they’re moving, decrease the voltage back to 50 V. Note the PLL behavior- if there’s a software issue and pulses aren’t being sent, you won’t see activity in the PLL associated with the coarse positioners. Under normal circumstances you should see considerable crosstalk between the PLL and the coarse positioners while the coarse positioners are firing. There are significant transients in the resistive encoder readings after firing the coarse positioners; this is likely a result of heating, but could also have a contribution from mechanical settling and creep. We have observed that the decay times of transients are significantly longer in the 300 mK system than in the 1.5 K or 4 K systems, likely indicating that these transients are largely limited by heat dissipation, pot with drainage holes at least at very low temperatures. Go into the General Approach Settings of the Approach Control window. There’s a setting in there for coarse positioner step size- set that to 4 µm or so. This is the amount the coarse positioners will attempt to move between fine scanner extensions. They always overshoot this number . Overshooting is of course dangerous because it can produce crashes if it is too egregious. In the Approach Control window, click Set PLL Threshold, verify that standard deviation of frequency is 0.25 Hz. Enter 5 µm into the height window. Verify that Z is ungrounded . Click Constant Height. Check that the PID is producing an approach speed of 100 nm/s. It is important that you sit and watch the first few rounds of coarse positioner approach. This is boring, but it is important the first few coarse positioning steps often cause the tuning fork to settle and change, which can cause the approach to accelerate or fail. Also by observing this part of the process you can often find simple, obvious issues that you’ve overlooked while setting up the approach. Getting to the surface will take several hours. Typically you’ll want to leave during this time. When you return, the tip should be at constant height. I’d recommend clicking constant height again and approaching to contact again to verify that you’re at the surface. You should be between 10 µm and 20 µm from the surface. It may be necessary to withdraw, approach with the coarse positioners a few µm, and then approach again to ensure you have enough scanner range in the z direction. Click withdraw until you’re fully withdrawn. Click Frustrate Feedback to enable scanning with tip withdrawn. I will present instructions as if you are attempting to navigate to a device through which you can flow current. This will generate gradients in temperature from dissipation and ambient magnetic fields through the Biot-Savart law, both of which the nanoSQUID sensor can detect. I strongly recommend that you navigate with thermal gradients if at all possible. The magnetic field is a signed quantity, so you need to have a pretty strong model and a clear picture of your starting location to successfully use it to navigate.

Thermal gradients can be handled with simple gradient ascent; this will almost always lead you to the region of your circuit with the greatest resistance, which is typically an exfoliated heterostructure if that is what you’re studying. You will likely need to have a helium atmosphere inside the microscope to pursue thermal navigation. A pressure of a few mBar is plenty, but be advised that this may require that you operate at elevated temperatures.Helium 4 has plenty of vapor pressure at 1.5 K, but this is not really an option at 300 mK, and many 300 mK systems struggle with stable operation at any temperature between 300 mK and 4 K. You should run an AC current through your device at finite frequency. Higher frequencies will generally improve the sensitivity of the nanoSQUID, but if the heterostructure has finite resistance the impedance of the device might prevent operation at very high frequency. It’s worth mentioning that the ‘circuit’ you have made has some extremely nonstandard ‘circuit elements’ in it, because it relies on heat conduction and convection from the device through the helium atmosphere to the nanoSQUID. If you don’t know how to compute the frequency-dependent impedance of heat flow through gaseous helium at 1.5K, then that’s fine, because I don’t either! I only mention it because it’s important to keep in mind that just because your electrical circuit isn’t encountering large phase shifts and high impedance, doesn’t mean the thermal signal is getting to your nanoSQUID without significant impedance. I recommend operating at a relatively low frequency for these reasons, as long as the noise floor is tolerable. In practice this generally means a few kHz. I’d also like to point out that if you are applying a current to your device at a frequency ω, then generally the dominant component of the thermal signal detected by the nanoSQUID will be at 2 · ω, because dissipation is symmetric in current direction . Next you will perform your first thermal scan, 10-20 µm above the surface near your first touchdown point. If you have performed a thermal characterization, then pick a region with high thermal sensitivity, but generally this is unnecessary- I usually simply attempt to thermally navigate with a point that has good magnetic sensitivity. Bias the SQUID to a region with good sensitivity. Check the transfer function. Set the second oscillator on the Zurich to a frequency that is low noise . Connect the second output of the Zurich to the trigger of one of the transport lock-ins and trigger the transport lock-in off of it.

Scarlet Royal table grape is one of the major red varieties in California

For RNA-seq analysis, a total of 8 RNA-seq libraries were generated, comprising four biological replicates from each of the two vineyards . The libraries were constructed as previously described using the NEBNext Ultra II RNA Library Prep Kit for Illumina . Subsequently, these libraries were pooled in equal amounts and subjected to paired-end 150-base sequencing on two lanes of the NovaSeq 6000 platform at the Novogene Co., Ltd .Illumina sequencing of the multiplexed RNA-seq libraries resulted in 8 FASTQ files containing sequences, and the dataprocessing followed the methods described in our previous work . In summary, the quality of reads was assessed using FASTQ before and after trimming with Trimmomatic v0.39 . Subsequently, the trimmed reads were quantified using Salmon in non-alignment based mode to estimate transcript abundance . The transcripts were mapped to the Vitis transcriptome file “Vvinifera_457_v2.1.transcript_primaryTranscript Only.fa” extracted from Phytozome database , resulting in a mapping rate higher than 61.9% . To identify differentially expressed genes between V7 and V9 at the sampling point, we utilized the DESeq2 and EdgeR packages with default parameters . For convenience, the DEGs generated by both DESeq2 and EdgeR pipelines, large pot with drainage with a threshold of PFDR<0.05 and log2fold change > 1.5 or < –1.5, were considered as being expressed . For the analysis of Gene Ontology terms, we employed the g:Profiler website with the g:SCS multiple testing correction method, using a significance threshold of 0.05 .

Finally, to visualize the consensus result, the Web-based tool Venny was used .Co-expression network modules were constructed using the variance stabilizing transformation values and the R package WGCNA . Before analyzing the data, lowly expressed genes among all sample types were removed by DESeq2, and the remaining non-lowly expressed genes of the 8 samples were used in module construction. The co-expression modules were obtained using the default settings, except that the soft threshold power was set to 9, TOMType was set to signed, minModuleSize was set to 30, mergeCutHeight was set to 0.25, and scale-free topology fit index was set to 0.8 . A module eigengene value, which summarizes the expression profile of a given module as the first principal component, was calculated and used to evaluate the association of modules with berry biochemical characteristics of V7-berries and V9-berries at the fifth sampling time . The resultant final WGCNA matrix had 42 modules with 17,553 genes. The module membership and gene significance values were calculated, subsequently the intramodular hub genes were identified . Despite the premium fruit quality of the variety, in some cases, an undesirable taste was observed under certain unknown circumstances. To gain comprehensive insights into the development of the occasional berry astringency of Scarlet Royal and understand the underlying mechanism of this phenomenon, berries were investigated at two contrasting vineyards , both following the same commercial cultural practices. However,leaf petioles analysis of grapes from both vineyards showed considerable differences in nutrient levels, especially in the primary macronutrients .

During both seasons, the amount of nitrogen in the form of nitrate in LP-V9 was roughly 2 to 3 times higher than the normal levels, in contrast to its counterpart in LP-V7, which slightly accumulated more or less N. Similarly, LP-V9 contained higher percentages of phosphorus and potassium compared to LP-V7 . Conversely, the amounts of secondary macronutrients, calcium and magnesium , in LP-V7 were within the normal range but greater than LP-V9, which showed Mg deficiency in the first year only. Regarding the micronutrients, their levels were mainly within or around the normal range at both vineyards and during both seasons, with some differences . For example, zinc was slightly higher in LP-V9, especially in the first year. On the contrary, manganese and chlorine were roughly 2 times higher in V7 . Similarly, soil analysis shoed a higher level of nitrogen, potassium and magnesium . However, no significant difference was observed in all other soil macro and micronutrients. During the two seasons of the study, we determined the total marketable yield and the number of clusters in both vineyards. Our data revealed a higher yield in V7 compared to V9 in 2016 and 2017, respectively. The lower yield in V9 can likely be attributed to the smaller number of clusters in V9 compared to V7 during 2016 and 2017. To monitor the changes in the biochemical composition of Scarlet Royal berries, V7 and V9 berries were periodically sampled at six time points from veraison until the end of the season . The obtained data showed that berry polyphenols exhibited discernible patterns in both vineyards, most notably during the ripening stage . Of special interest were the tannin compounds, which widely affect organoleptic properties such as astringency and bitterness . Our data showed that berries from both V7 and V9 vineyards maintained lower levels of tannin from veraison up to the middle of August . Subsequently, a significant gradual increment of tannin took place. However, only V9-berries showed consistent accumulation of tannin over the two studied seasons compared to V7-berries, where the significant induction occurred only during the first season.

It is worth noting that the levels of tannin were lower in both vineyards during the second year compared to the first season. Nevertheless, they were more pronounced in V9-berries compared to V7-berries, with roughly 2- to 4.5-fold increases by the end of the harvesting time during the two seasons, respectively . The patterns of catechin and quercetin glycosides were inconsistent during both seasons, particularly within V7-berries . During the first year, for instance, the levels of catechin were similar in both vineyards, showing a dramatic increase only by the end of the season . In contrast, during the second year, such induction of catechin was exclusively restricted to V9-berries, starting from time S3 . For quercetin glycosides, V7-berries exhibited significantly higher amounts at early stages during both seasons relative to V9-berries . However, subsequent amounts were comparable in both vineyards during the first season only , but not in the second one, where V7-berries showed a significant drop at the last sample S6 . Interestingly, the levels of quercetin glycosides were roughly equal at the last V9-berries sample between both seasons despite such inconsistency. For total anthocyanins , the levels in early samples were comparable in both vineyards and seasons . Afterwards, their pattern started to vary between V7 and V9 within the same season, as well as from the first season to the second, as the nutrient amounts fluctuated as well . Nevertheless, TAC accumulation was positively correlated with the progress of ripening in V7-berries, but not V9-berries. To further confirm our data, we measured these phenolic compounds for the third time in mid-September of the next year . Overall, the results showed that the patterns of tannins and TAC were reciprocally inverted between V7-berries and V9-berries as ripening advanced. In addition, both catechin and quercetin glycosides most likely followed the pattern of tannins despite their seasonal fluctuations. To further distinguish V7-berries and V9-berries and assess their astringency development, a panel test was performed using samples at three commercial harvest times . A group of 12 nontechnical panelists scored berry astringency on a scale from 1 to 7, square pot where 1 is extremely low and 7 is extremely high. The panelists were trained using samples from contrasting standard varieties, including Flame Seedless and Crimson as non-astringent and Vintage Red known for its astringent taste . The results showed that V7-berries exhibited lower intensity of astringency compared to V9-berries . As ripening proceeded, astringency levels increased in V9-berries, but decreased in V7-berries. Moreover, we collected samples from clusters with various astringent taste and measured its tannins content. We were able to determine that the threshold level of tannins that causes the Scarlet Royal astringency taste is around 400 mg/L . Taking into account the levels of polyphenol compounds and the taste panel data together , it is evident that astringency development is positively associated with tannins’ accumulation throughout the ripening process of V9-berries. Nevertheless, organoleptic analysis revealed a significant difference in the berries of the two vineyards, particularly in terms of total soluble solids and titratable acidity . Notably, V9 berries exhibited higher titratable acidity and lower total soluble solids, especially in the later stages . It’s worth noting that the weight of V9 berries is also higher than that of V7 .To better understand the molecular events associated with the induction of tannins and astringency upon ripening, the berry transcriptome profile was analyzed in both V7-berries and V9- berries at the late commercial harvest date . Following the quality and quantity check, extracted RNA from quadruplicate samples was deeply sequenced . Of the 19.7 to 24.4 million high-quality clean reads per replicate, 61.9% to 66.1% were mapped against the V. vinifera transcriptome . Hierarchical clustering of the RNAseq data showed explicit changes in the berry transcriptome profile between V7- berries and V9-berries . The Principal Component Analysis showed high consistency among biological replicates . Samples were mainly separated along the first component , which was responsible for 97% of the variance, and was definitely associated with the site ofcultivation; V7 and V9. In contrast, the second component was trivial, accounting for only 1% of the variance and was probably attributed to experimental error. Such results were expected, as berry samples came from the same cultivar, Scarlet Royal , and the only difference between them was the vineyard locations.

To identify the differentially expressed genes in V7- berries and V9-berries at this specific time within the ripening window, the RNAseq data were analyzed using two different Bioconductor packages, DESeq2, and EdgeR . Subsequently, the DEGs with FDR < 0.05 and log2fold change > 1.5 or < –1.5 generated by both pipelines were considered . The pairwise comparison between berry transcriptomes resulted in 2134 DEGs, with 1514 up-regulated and 620 down-regulated . The data manifested the impact of the cultivation site on the transcriptional reprogramming of a large number of genes that ultimately affect berry quality. Most apparently, at the V9 vineyard, where roughly 2.5-fold higher number of berry transcripts were upregulated compared to V7 . Subsequently, the enrichment of Gene Ontology terms and Kyoto Encyclopedia of Genes and Genomes pathways were analyzed among the up- and down-regulated DEGs using the Vitis vinifera Ensembl GeneID . Among the significantly enriched GO terms, the up-regulated transcripts in V9-berries exhibited high enrichment in the molecular function GO terms for quercetin 3-O-glucosyltransferase activity and quercetin 7-Oglucosyltransferase activity . Additionally, the V9-berries induced DEGs were highly enriched in the biological process GO terms for the jasmonic acid signaling pathway and cellular response , Lphenylalanine metabolic process , L-phenylalanine biosynthetic process , and nitrogen compound metabolic process . Similarly, these DEGs were highly enriched in the KEGG pathways for the biosynthesis of secondary metabolites and phenylpropanoid biosynthesis . On the other hand, the down-regulated transcripts in V9-berries showed substantial augmentation in the MF GO terms for hormone binding , abscisic acid binding , and potassium ion transmembrane transporter activity . Correspondingly, the BP GO terms for hormone-mediated signaling pathway and response , auxin-activated signaling, cellular response, and homeostasis , abscisic acid-activated signaling, response, and cellular response , response to strigolactone , potassium ion transmembrane transport , and potassium ion transport , as well as the KEGG pathways for plant hormone signal transduction , brassinosteroid biosynthesis , and carotenoid biosynthesis were highly enriched in the down-regulated genes of V9-berries . Overall, the transcriptome analysis pointed out the substantial changes in transcript abundance that coordinate and reflect the observed induction of tannins/astringency during the maturation and ripening of V9-berries compared to the V7-berries .To elucidate which fundamental processes were altered during tannins/astringency induction within berries, the Weighted Gene CoExpression Network Analysis was applied to construct coexpression networks. Forty-two modules were identified based on pairwise correlations among the 17553 non-lowly expressed genes . Subsequently, the biochemical data from both V7-berries and V9-berries were correlated to the WGCNA modules, and only 2 modules, M21 and M30, displayed substantial correlations with berry polyphenols, containing 5349 and 4559 genes, respectively . The M21 module was positively linked with TAC , but negatively associated with tannins, catechin, and quercetin glycosides . On the contrary, the M30 module exhibited a positive correlation with tannins, catechin, and quercetin glycosides , but was negatively linked with TAC .

Tunnel culture is now a common practice in raspberry production

In years with adverse production conditions or low prices, a higher percentage of the crop may be diverted to the freezer or processed products market. Fresh market fruit is handled and sold primarily through local grower-shippers; a much smaller share is sold directly to consumers through farmers markets, community supported agriculture operations, farm stands and other direct and intermediated market channels such as restaurants, independent grocers and schools.Arguably the most momentous shift in cultural practices for strawberries was the introduction of preplant soil fumigants, beginning with chloropicrin in the 1950s and methyl bromide in the 1960s. Fumigation is a soil disinfestation practice that improves plant productivity and helps with the management of arthropods, nematodes, weeds, soilborne fungi and other plant pathogens. Some of the most difficult to control pathogens include Verticillium dahliae, Fusarium spp. and Macrophomina phaesolina. Without soil fumigation, these pathogens have the potential to completely destroy strawberry plantings. Early on, when CP and MB were mixed and applied together, the synergistic effects allowed strawberries to be produced as an annual rather than a biennial crop, and to be grown continuously on the same land without rotation to another crop, nft vertical farming resulting in an increase in annual strawberry acreage. The use of fumigants also led to higher and more predictable yields and fruit quality, and further enabled the development of more stable markets for strawberries .

Yields for strawberries statewide increased from a range of 2 to 4 tons per acre prior to the introduction of soil fumigants to 16 tons per acre by 1969 . Additional cultural improvements included the development of both UC and proprietary strawberry varieties uniquely adapted to coastal production conditions. Varieties were bred, for example, for disease resistance, yield and market potential. Notable UC-bred strawberry varieties include Tufts , Pajaro, Douglas, Chandler, and Selva , Camarosa and Seascape , and Aromas, Albion and Monterey . Irrigation practices also evolved, shifting from furrow irrigation in the 1960s to drip irrigation in the 1980s, which led to further improvements in plant disease management and greater water use efficiency. These and other enhancements meant that by 2012, yields could exceed 35 tons per acre . More recently, the strawberry industry has focused on “fine-tuning” fertility and water management for more efficient resource use, along with additional yield and fruit quality improvements . The Santa Cruz–Monterey area is also recognized for its early experience with conversion of conventional strawberry production to organic management . Organic strawberry production was shown to result in lower yields, which, when offset by premium prices could potentially offer higher net returns to growers. The importance of crop rotation for disease management was not addressed in the initial study by Gliessman et al. but has since been the focus of additional research, as have more complete analyses of the economics of organic strawberry production . Growers and area researchers continue to collaborate and advance organic strawberry production techniques. Most notably, a long-term research commitment has been made to determine organically acceptable disease management practices such as anaerobic soil disinfestation , the use of commercially available soil-applied biological organisms and the incorporation of soil amendments such as mustard seed and its derivatives.

The area is now seen as a global leader in organic strawberry research, and in 2012 the first organic strawberry production manual was published by UC Agriculture and Natural Resources . Statistics documenting expansion of the organic strawberry industry over time are not available on a county-by-county basis, but statistics for California show prodigious growth in acreage and value of production: from $9.7 million in 2000 to $93.6 million in 2012, a 621% increase in real dollars .Like strawberries, raspberries and blackberries have benefitted from enhancements in cultural practices. When well-managed, both types of caneberries can produce crops for up to 20 years. However, to maintain acceptable quality and yield Central Coast growers typically manage raspberries and blackberries so that they produce two and five crops, respectively, prior to removal and replanting. In Santa Cruz County, raspberry production was relatively flat in the 1960s and 1970s, but began to increase substantially in the 1980s . This can be explained by a shift from floricane, or spring-bearing varieties, to the then newly developed proprietary primocane, or fall-bearing varieties, that do not carry the productivity constraints associated with the inadequate chill requirements along the Central Coast. Primocane-bearing varieties allow growers to successfully produce a high quality raspberry crop in low- or no-chill coastal locations, and further manipulate time to harvest and yield with pruning and other management practices . Between 1990 and 2014, the number of acres planted to Santa Cruz and Monterey area raspberries almost tripled, tons produced increased by about 350% and the value of production was up by over 400% in real dollars . Santa Cruz County raspberry growers began to experiment with and adopt field-scale semi-permanent protective structures or tunnels in the 1990s and 2000s . Initially developed in Europe, field-scale tunnels allow growers to extend their production seasons, enhance yield and fruit quality, and capture high off-season prices for fresh market fruit . The controlled environment, and resulting security of production, also allows for greater market stability. This shift away from open-field production to protected cropping, along with breeding improvements, has had lasting impacts on the raspberry industry and its expansion. Cultural improvements geared towards fresh market blackberry production are more recent and include advances in breeding for thornless varieties and quality attributes . In 2011, a public primocane-bearing blackberry variety became commercially available for the first time and is now being planted in the area. Since that time, additional public and proprietary primocane-bearing varieties have been in development; some have already become available.

Open-field production was the norm until recently, but to ensure marketable fruit of high quality, and as growers have shifted additional acreage to primocane-bearing varieties, tunnel culture has been more widely adopted and, based on discussions with growers, is now estimated at roughly 80% of the acreage. Like organic strawberries, remarkable growth in the statewide production of organic raspberries and blackberries was documented between 2000 and 2012 . Acreage climbed by over 500% in both organic berry categories. Value of production was up over 3,000% in real dollars for organic raspberries and up by almost the same percentage for organic blackberries. It is important to note that although the organic raspberry and blackberry categories have demonstrated extraordinary growth, they still represent a relatively small percentage of all berry production in the area.Research points to several factors that have spurred consumer demand for all berries. Berries contain bioactive compounds, including essential vitamins, minerals, fiber and antioxidants that contribute to healthy diets, and that help to reduce the risks associated with some chronic diseases and cancers . This information has been widely shared with consumers through, for example, government programs promoting healthy eating , and more generic berry promotion programs . Per capita consumption of fresh strawberries in the United States almost doubled from 1994 to 2014, increasing from 4.1 to 8.0 pounds . U.S. per capita consumption of fresh raspberries was small by comparison, at just 0.5 pounds in 2014. Similar consumption data are not available for blackberries, but Cook notes that consumers generally view berries as complementary, and that sales for all berries have increased. Indeed, in 2014, berry sales increased 5.8% over 2013; berries were the number one produce category for U.S. grocery retailers, at $5.7 billion in annual sales . Some berry operations also benefit from their proximity to the area’s urban centers, vertical tower for strawberries which have sizeable cohorts of educated, high-income consumers who generally demonstrate an interest in health and wellness, local agriculture and fresh and organic products. In addition to the more traditional grower shipper and direct marketing channels, new technology-driven food marketing companies — virtual food hubs — have evolved to cater to this demographic. They promote the values of sustainable communities, local food economies and business integrity and transparency, all important attributes for new 21st century consumers . These companies form relationships with local growers, provide some technical and market support, and enhance sales and engagement with consumers. It is not yet clear what impacts these still-niche marketing businesses may have on the industry in total. However, growers have responded to the various health and market signals by ramping up production of both conventional and organic products, berries included.Specialists and farm advisors with UC Cooperative Extension have performed economic analyses for Santa Cruz and Monterey county fresh market berry crops for decades . The studies estimate production costs for a representative enterprise based on characteristics common to the area’s farms. Data are collected from established growers, input suppliers and other industry experts so that a diversity of operations and practices are taken into account. Since 1990, UCCE researchers have used a farm budget software program to analyze the data and present results in several formats detailing costs for cultural and harvest practices, monthly cash costs and business and investment overhead costs. The studies also include an analysis estimating net returns to growers for several yield and price scenarios. Representative costs for food safety and environmental quality programs have been incorporated into more recent studies as they have evolved to become standard business practices. The resulting production and economic information is specifically designed to assist growers, bankers, researchers and government agencies with business and policy decisions.The first economic analysis of fresh market strawberry production for Santa Cruz and Monterey counties was performed in 1969; at least one subsequent analysis has been conducted every decade since then. Though the level of detail and data included in each study has changed over time, some interesting trends can be noted. Annual land rent climbed from $150 per acre in 1969 to $2,700 in 2014, representing 2.5% and 5.5% of total production costs, respectively. The cost of soil fumigation for conventional strawberry production increased from $350 per acre in 1969 to $3,302 in 2010, representing 5.5% and 6.9% of total production costs, respectively. Production year water use gradually decreased from 80 acre-inches per acre in 1969 to 36 acre-inches by 1996 as drip irrigation became the standard. The amount of water used to bring a crop to harvest has remained roughly the same since that time; however, growers and researchers continue to investigate methods to increase water use efficiency even further. In some areas, soil types and fields, growers have been able to reduce per acre water use by several acre-inches more . When the above costs and water usage are assessed on a per ton rather than a per acre basis, production practice cost increases are less notable, and water savings even greater. Labor-intensive practices such as hand weeding and harvest are consistently shown as costly line items relative to other operations. Representative yields for conventionally produced fresh market strawberries rose from 20 tons per acre in the 1969 study to 30 tons in 2010, an increase of 50%. Even higher yields are discussed for some varieties and production conditions; county production statistics confirm that higher yields are indeed possible . Representative yields for organic strawberries, studied over a much shorter time period, rose from 15 tons per acre in 2006 to 17 tons in 2014, an increase of 13%. As more research is directed towards organic agriculture in general and strawberries in particular, yields will likely increase even more with time. Recent efforts include improvements in cultivar breeding, cultural practices and disease management, especially soil pathogen management. The most recent economic analyses for conventional, second year conventional and organic strawberry production were performed in 2010, 2011 and 2014, respectively. Second year conventional strawberries, or those producing a crop for a second year after having produced the first without replanting, represent about 15% of the total strawberry acreage in the area. Similarities and differences in total, cultural and pest management costs for the three management approaches are shown in figures 1 to 3. Total costs for conventional strawberries were $47,882 per acre and include expenses for all practices from land preparation to harvest . For the second year conventional strawberry crop, total costs were lower at $32,798 per acre, reflecting a reduction in expenditures for land preparation and reduced harvest costs because of lower yield. For organic strawberries, total costs were $49,044 per acre, slightly higher than for conventional production, mostly due to higher soil fertility input costs.

Plants were carefully watered to eliminate risk of contamination via water splash

Furthermore, colonization of pathogens during drought may further disrupt the carbon balance of plants as it influences defense and repair, creating a feedback loop that can drive plants toward a mortality tipping point . Thus, while dehydration tolerance may be important during typical seasonal drought conditions, it may be a much riskier strategy and lead to greater mortality during global-change type drought, especially in the presence of pathogens. These frameworks are consistent with our findings and provide further evidence that A. glauca experiencing acute levels of drought stress are highly predisposed to Bot. infection particularly at lower elevations that experience heightened levels of water stress.The results of this study provide strong evidence that A. glauca in the study region are vulnerable to Bot. disease and dieback, and possibly eventual mortality, related to acute drought. This is consistent with Venturas et al. , who found that acute drought in 2014 led to reduced abundance in A. glauca and other obligate seeder chaparral species and even type-conversion in the Santa Monica Mountains of southern California, USA. A review by Jacobsen and Pratt found similar consistencies among shallow-rooted, low round pots obligated seeding shrubs. Clearly, there is strong support that A. glauca populations are at risk for future dieback, and thus should be the focus of more intense studies aimed at understanding the possible mechanisms driving such events.

Manzanita are important members of the chaparral ecosystem and large-scale dieback and mortality of this species could reduce resource availability for wildlife , as well as increase the risk of more intense, fires in an ecosystem already associated with increasingly frequent fire activity. Additionally, our study provides valuable insight into areas of greatest risk for dieback and mortality, which are predominantly in lower elevations. These are important factors to consider when predicting vulnerabilities and potential impacts of future extreme drought events . Mediterranean shrublands like those in southern California already considered high risk for global-change type drought, , and research suggests a general trend of upwards-shifting ranges in southern California chaparral species driven by changes in climate .Therefore, populations of A. glauca occurring at the lower edge of their natural range are at high risk for dieback and mortality, and should be the focus of management efforts. Lastly, while studies on the various physiological mechanisms for plant survival during drought are critical for predicting differential responses to stress, there is an increased emphasis on the importance of understanding the diverse role of pathogens in order to accurately model species vulnerabilities to climate change . Studies that incorporate the impact of pathogens help inform new integrative approaches to protecting plants against drought and biotic infection, rather than treating these influences separately.

Examples include Jactel et al., , whose meta-analysis showed the significant effects of water stress on symptom severity in plants infected with latent pathogens like Bots, and experiments like Drake-Schultheis et al. , who found interactive effects between drought stress and infection from N. australe in driving symptoms of stress and increasing mortality rates in A. glauca. The results of our study align with these frameworks, and provide additional evidence that as climate change models are predicting more intense and frequent drought events, our need to understand the role of latent pathogens in at-risk natural systems is becoming more critical.Reports of large-scale, drought-associated mortality events in forest and woodland systems have been on the rise in recent decades . These reports have spanned across biomes, including in classically drought-tolerant species across Europe , Australia , Africa , and the United States . As a result, interest has been growing in understanding how species that are typically capable of withstanding periodic drought stress may become susceptible to drought and experience significant dieback and even large-scale mortality when exposed to acute or prolonged chronic drought . These droughts of unusual extremes are referred to as “globalchange-type drought” and are becoming more common as the climate warms . While the exact physiological mechanisms leading to dieback and mortality during such events are variable across species and conditions, drought is generally hypothesized to promote physiological decline either via loss of hydraulic functioning or carbon starvation or a combination of both . In the case of hydraulic failure, plants with insufficient soil water experience xylem cavitation , which can ultimately lead to cellular death. Alternatively, plants that avoid drought by closing their stomata to reduce water loss subsequently suffer insufficient carbon supply to meet other metabolic demands.

In either scenario, the stress that drought places on a plant is likely to cause measurable decreases in physiological functions that may be irreversible . An additional factor that can play a significant role in drought-related dieback and mortality is the presence or introduction of biotic agents. Indeed, introduced plant pathogens have been well documented to cause canopy dieback and dramatically alter community structure in a variety of forested systems . Some well-known examples in the United States include Dutch elm disease , chestnut blight , white pine blister rust , and sudden oak death . Significant pathogen events have also impacted the landscape in wild land shrub communities including sclerophyll shrub woodlands in Australia and salt desert scrub in the western United States . However, large-scale dieback of shrubs has been less documented than their arboreal counterparts, despite evidence of disease from fungal species being abundant in many scrubland systems including southern California chaparral , northern California foothill shrublands , and South African fynbos . Such studies, along with expectations of increasing threats from pathogens due to climate change and accelerating trade/movement of biological materials globally , have led scientists and land managers alike to anticipate introduced pathogens as important contributors to future changes in wild land communities.While both global-change-type drought and pathogens are likely important contributors to plant dieback and mortality, current research suggests that these two factors are not mutually exclusive . Rather, canopy dieback and mortality may result from the combined influences of environmental stress and biotic agents, and theoretical frameworks describing these influences have been put forth . These frameworks incorporate biotic agents into the drought-hydraulics complex described above, whereby pathogens and insects may amplify or be amplified by drought-associated hydraulic failure or carbon starvation . Amplification can occur when biotic agents damage host tissue—by defoliation or blocking transportive vessels, for example—to the extent that the effects of drought are greatly exacerbated . Alternatively, physiological responses to extreme environmental stress can have negative effects on plant defense systems, rendering them susceptible to mortality through biotic infection . In both scenarios, the effects of biotic agents and drought stress are strongly linked, and these interactions have been well documented in drought-tolerant systems such as South African fynbos , red pine forests , eucalyptus forests , and California chaparral . Latent or secondary pathogens are particularly likely to be involved with dieback and mortality events in these systems, as they are known to increase damage in hosts experiencing drought stress . Therefore, while drought events alone are expected to play an important role in reshaping ecosystems as the climate changes, in some cases, synergies between environmental stress and biotic influences might lead to shifts in plant community structure and composition, and thus ecosystems as a whole. In the Santa Ynez Mountains in Santa Barbara County, California, United States, big berry manzanita began exhibiting dramatic canopy dieback during the 2011–2018 drought . Shrubs in the genus Arctostaphylos are common in Mediterranean shrub communities extending from southwest Oregon to northern Baja California . They may occur in monospecific stands or in alliances with other important community members like chamise and Ceanothus spp. . Within these alliances, Arctostaphylos spp. frequently occupy >50% average cover , which along with their nutritious and prolific fruits, and fire-induced regeneration strategies, make them one of the most important members of the chaparral community . In the southern California chaparral ecosystem where hot, dry summers with high vapor pressure deficit are the norm , seasonal drought tolerance has long been considered a common strategy among dominant plant species, including A. glauca. However, the severity of recent canopy dieback observed suggests that this species is reaching a threshold in its drought-resistance capability. Concurrent with observations of canopy dieback, visible symptoms of fungal infection were observed including wood cankers and leaf discoloration , plastic pots 30 liters both of which progress during prolonged drought stress, suggesting that multiple driving forces contribute to manzanita dieback. Molecular sequencing identified the dominant fungal pathogen found on symptomatic A. glauca in this area to be Neofusicoccum australe, a member of the well-known pathogenic Botryosphaeriaceae family . Members of this family are most commonly associated with disease in plant species experiencing severe environmental stress , including Arctostaphylos spp. . They are also known to play a variety of functionally diverse roles, from asymptomatic endophytes to obligate pathogens . Yet, while N. australe has been described around the world , relatively few studies have been conducted on its specific interactions with host species, as it was only fairly recently described .

Historically, Bot. pathogens have most frequently been studied in agricultural host species , and little is known regarding their ecological role in wildland ecosystems , especially with regards to chaparral shrubland systems . The present study was aimed at identifying the possible role of N. australe in A. glauca dieback in Santa Barbara County, particularly in combination with extreme drought. Because this pathogen has only recently been reported on wild shrub species in California and is thought to be an introducedspecies native to Western Australia , this outbreak represents a new and undescribed threat to these wildland plant assemblages. This study addresses the following questions: How does A. glauca respond physiologically to drought and fungal infection, separately and together? Are these responses correlated with visual signs of stress, specifically leaf health? Can drought and fungal presence interact to increase or accelerate plant mortality compared to drought or fungi alone in A. glauca? To address these questions, a greenhouse experiment was conducted in November 2016 through February 2017 manipulating both drought and fungal infection and observing trends in plant stress symptoms, physiological function, and mortality. We predicted that both drought stress and fungal infection would lead to declines in physiological function compared to the control and that these declines would be strongly correlated with increases in stress severity. Furthermore, we expected that those individuals experiencing both drought stress and fungal infection would die sooner than those in all other treatment groups. This experimental study elucidated the potential of the interaction between drought stress and introduced pathogens to significantly impact chaparral shrub health and important implications for the future of these shrubs faced with increasingly frequent global change-type droughts.A completely randomized full-factorial design was used to organize the individuals into four treatment groups: droughted and inoculated with N. australe , droughted and not inoculated , watered and inoculated with N. australe and a control; watered and not inoculated . Data were collected for ~90 days to track declines in health and mortality rates among the different treatments. Drought-treated plants received 1 L of water on the day of inoculation and another 0.5 L on day 38. Those with no drought treatment received 0.5 to 1.0 L of water by hand once per week depending on soil moisture, which was monitored regularly using a TDR machine from Soil Moisture Co. . Soil moisture for non-drought plants was maintained between 15–25% moisture for the entire experiment. Cultures for inoculations were made from re-isolations of field samples that were collected in January 2016 and positively identified to be N. australe . Inoculations took place on 3 November 2016 , using methods adapted from Michailides and Swieckiand Bernhardt . Mycelial plugs were made from 8-d-old cultures growing on half strength potato dextrose agar amended with streptomycin to prevent bacterial contamination. Plants were first sprayed with 70% isopropyl alcohol to sterilize the surfaces and surrounding areas. Mycelial plugs were taken from the advancing margin of N. australe cultures and placed on strips of Parafilm using sterile petroleum jelly for adhesion. Plugs were then placed to superficial wounds made on the main stem . The Parafilm strips were then gently wrapped 2–3 times around the stem to keep the plugs in place and prevent contamination. Those plants not receiving fungal inoculation received a control inoculation with uncultured potato dextrose agar using the same techniques.

The RRB subtype found most consistently across studies has been self-injurious behaviors

Individual items, or in the case of the current study, individual questions from the RBS-R, were independently assessed for conceptual fit on the factor they most strongly loaded on to determine if it is an appropriate factor fit considering the other items that loaded strongly on the respective factor. The overall goal of EFA was to identify factors, based on a given dataset, and maximize the amount of variance explained by the model . Once a model has been theoretically and/or statistically established and hypotheses have been made, a confirmatory factor analysis can inform the likelihood of the hypothesized results.A Confirmatory Factor Analysis was conducted once the relationships among variables were established through statistical analyses and a theoretical model was evaluated . While the EFA allows for all items to load on any factor, the CFA restricts the factors on which items load. Each item was permitted to load on only one factor. Model fit was determined using recommended indices of model fit including Chi-Squared test, RMSEA, RMR, CFI and TLI. Additionally, the CFA model produced a weighted root mean square residual that is an empirically supported measure of model fit comparable to the other fit indices and is suggested to be highly useful for data that isn’t normally distributed . A WRMR value above 1.0 is considered good model fit. Factor loadings from the CFA were reported as the standardized model estimate loadings and associated standard errors.Cluster analysis provides a unique approach to examining which results in the identification of patterns that organize variables into taxonomies, grouping cases with similar patterns together . For the current study, blueberry grow pot the K-means cluster analysis was run to systematically and conceptually group participants with similar RRB patterns together.

The newly established factors from the CFA of the RBS-R were used to examine the various patterns of RRB presentation for this population. The goal of a k-means clustering is to partition individuals into clusters where every participant belongs to a cluster with others presenting with similar patterns . The optimal number of clusters must strike a balance between successfully compressing the data as a single cluster would, while maintaining maximum accuracy where every participant is assigned to its own cluster. The optimal number of clusters for the data was determined using both theoretical and empirical considerations. Previous research exploring RRBs have defined between two and six distinct types of RRBs; yet, there hasn’t been a clustering of those RRBs into distinct profiles to serve as a comparison or as an empirical rationale to test the fit of a specific number of clusters. Therefore, comparisons of three, four, five and six cluster solutions were conducted. One approach that was used to determine model fit for each cluster was to examine the number of iterations it took to satisfy the convergence criterion . There is no guarantee that data will cluster and iterate to convergence quickly, if at all. Therefore this is a reasonable justification for this approach in determining the fit between the number of clusters and the data being analyzed. Statisticians have concluded that it is acceptable to institute a maximum criterion of between 15 and 20 iterations for the data to reach convergence criterion where the clusters optimally fit the data. Cluster statistics were explored after running three, four, five and six cluster solutions; results are described below.The final research aim was to determine the ability of several behavioral and developmental characteristics to predict cluster membership. Correlation analyses among all predictors were conducted prior to running the MLR to determine presence of collinearity.

A multinomial logistic regression was run with individual cluster assignment as the outcome variable and participants’ standardized scores of ASD severity, nonverbal IQ, hyperactivity, anxiety, and coping skills as predictors. Age differences across clusters was independently examined by running a one-way ANOVA prior to running the MLR to determine if age significantly differed among clusters. The MLR provides a unique approach to determine the odds ratio of an individual being in one cluster relative to the odds of them being in the comparison cluster based on several characteristics . Therefore, it is important to choose a comparison cluster that will provide the most robust information in the analysis of these comparison solutions. Prior to exploring the individual cluster phenotypes to decide on a comparison cluster, the options were carefully considered and a conceptual decision was made. The comparison group should be the one that differs the most from the others, or the group that could be considered the “optimal outcome” group that possesses characteristics that researchers would want to test and discover what makes that group of participants different . Therefore, the cluster with the lowest levels across all RRBs was used as the baseline comparison cluster. Goodness of fit of the MLR model was assessed using the log-likelihood , which sums the probabilities of predicted outcomes and actual outcomes, analogous to the residual sum of squares in typical multiple regression. That is, the LL variable indicates how much unexplained data remains after the model is fit; where large values of the LL statistic tends to describe a poor fit for the model . Results of the multinomial logistic regression produced significance statistical values, which indicated the extent to which individual characteristics were able to significantly predict membership to one cluster over another. The individual parameter estimates for each comparison between the optimal profile group vs. the other profiles were individually examined to determine the significant and non significant results across predictor variables and interactions. The significance values were used to determine which of the characteristics were significant in predicting profile membership, with the odds ratio statistic indicating the odds of a participant being in a cluster when compared to the odds of them being a member of the optimal outcome profile group. Overall model fit statistics as well as individual parameter estimates of the multinomial logistic regression were examined.The recent changes to the DSM have created a more comprehensive list of RRB subtypes than were previously included and set a more stringent benchmark to meet criteria in the RRB domain. Such changes reflect the progression of research supporting the importance and independence of RRBs as an integral component of diagnosis, rather than a by-product of the “core” social communication impairments . From its earliest conception, ASD has been characterized by the presence of frequently and highly repetitious behaviors, with a marked desire for environmental sameness and consistency . Yet, this complex behavioral domain is historically under-represented in research efforts and falls secondary to social communication deficits in ASD research. Reviews of past studies on RRB presentation have highlighted issues including a lack of methodological consistency, with varying approaches to defining, organizing, and measuring RRBs. These discrepancies have led to splintered advancements in understanding the etiology, early behavioral manifestations and longitudinal developmental implications of RRBs . The primary aims of this study were to characterize RRB phenotypes of individuals with ASD and to determine the influence of developmental and behavioral characteristics on RRB profiles. This study revealed that there were five distinct RRB subtypes captured by the RBS-R, with five distinct phenotypic profiles generated from those subtypes. Hyperactivity, hydroponic bucket anxiety and coping skills significantly predicted participants’ RRB phenotype, while IQ and symptom severity had little effect.

The findings in this study provide a unique perspective when conceptualizing ASD symptomology and the influence of non-ASD specific traits on this core domain.The five-factor model result from the factor analyses of the RBS-R exhibits substantial consistency with previous studies examining the factor structure of the RBS-R . Comparisons between factor results of the RBS-R can be seen in Figure 3. Most notably, the current study excluded 3 items from the original compulsive scale as the item factor loadings were above .4 on more than two newly calculated factors. These results indicated that there wasn’t a single factor that accounted for the variability of each item, forcing those items to be excluded. Similarly, five items on the original ritualistic scale were excluded which included items regarding eating, sleeping, travel, play and self-care as they were highly loaded on multiple factors. These findings indicate they may not be sufficiently differentiating types of RRBs measured by each question, which leads investigators to wonder if the questions are adequately differentiating between RRB subtypes. Bishop, et al. investigated RRB data from both the ADI-R and the RBS-R and found that the ADI-R items resulted in a two-factor model, whereas the RBS-R resulted in a five-factor model as the best fit. When examined in conjunction with findings from the Lam & Aman study as well as the current results, it is evident that using a measure with a wider range of questions such as the RBS-R provides more in-depth and informative results when examining the specific types and severities of RRBs. Despite the utility of a measure dedicated to specific RRB types and severity, factor results from previous studies fail to be substantiated with each study, leading to the conclusion that a final set of RRB subtypes have yet to be established unequivocally across studies. Further, each analytic result has not been entirely consistent with the six conceptually derivedsub-scales that Bodfish, et al. originally established. Discrepancies between the subtypes and the original sub-scales, as well as between the previously proposed models can be seen below in Figure 3. As previously discussed, RRBs comprise a complex and heterogeneous set of behaviors that vary greatly depending on the population being measured; therefore, it is not a complete surprise that each factor analytic study has resulted in slightly altered structures. However, given the vast age range include in the current study and largest number of participants to date for an RBS-R factor analysis, the resulting factor structure warrants consideration as an organizational RRB factor structure to be analyzed for confirmatory analyses in future studies using the RBS-R. As seen in Figure 3, researchers who organized and defined more than two categories of RRBs had one striking consistency, the inclusion of an independent category of self-injurious behavior . Further, self- injury is arguably the most recognizable and disruptive RRB consistently found to be related to greater impairment with significantly lower IQ and higher severity of ASD symptoms . In fact, the most recent study examining RRB subtypes concluded that SI behaviors create significant difficulty in dichotomizing RRBs, as the SI items fail to load with the repetitive sensory motor category or with the insistence on sameness supporting the existence of additional subcategories . Further, SI is the only RRB subtype to consistently load identically as an entire sub-scale in every factor analytic study of the RBS-R, which was also true in the current study, indicating its distinctiveness .The RRB subtype found most consistently across studies has been self-injurious behaviors. As seen in Figure 3, researchers who organized and defined more than two categories of RRBs had one striking consistency, the inclusion of an independent category of self-injurious behavior . Further, self- injury is arguably the most recognizable and disruptive RRB consistently found to be related to greater impairment with significantly lower IQ and higher severity of ASD symptoms . In fact, the most recent study examining RRB subtypes concluded that SI behaviors create significant difficulty in dichotomizing RRBs, as the SI items fail to load with the repetitive sensory motor category or with the insistence on sameness supporting the existence of additional subcategories . Further, SI is the only RRB subtype to consistently load identically as an entiresub-scale in every factor analytic study of the RBS-R, which was also true in the current study, indicating its distinctiveness .Cluster analysis provides a novel approach to statistically explore phenotypic profiles and the co-occurrence of RRB types and severity across individuals with ASD. This is the first study of its kind to statistically generate clustered phenotypes, each consisting of multiple RRB subtypes. RRBs don’t occur in isolation; the pattern of behavior is fluid with minimal evidence to explain the variations seen across and within individuals. By studying RRBs in a way that allows for multiple RRBs to co-occur at varying levels, researchers may gain a more accurate and informative picture of how these behaviors manifest across individuals with ASD. However, when researchers rely solely on parent report measures, there is a limited scope of distinct behaviors from which combination or cluster phenotypes can be derived.

Forty-eight bar-coded small RNA libraries were constructed starting from 50 ng of small RNAs

Increasing temperature can accelerate metabolism, including sugar biosynthesis and transport, but the increase in metabolism is not uniform. For example, the increase in anthocyanin concentration during the ripening phase is not affected as much as the increase in sugar concentration. These responses vary with the cultivar, complicating this kind of analysis even further. Direct studies of temperature effects on Cabernet Sauvignon berry composition also are consistent with our data. In one study, the composition of Cabernet Sauvignon berries was altered substantially for vines grown in phytotrons at 20 or 30 °C temperatures. Cooler temperatures promoted anthocyanin development and malate concentrations and higher temperatures promoted TSS and proline concentrations. In a second study, vines were grown at 20 or 30 °C day temperatures with night temperatures 5 °C cooler than the day. In this study, higher temperatures increased berry volume and veraison started earlier by about 3 to 4 weeks. The authors concluded that warmer temperatures hastened berry development. In a third study, Cabernet Sauvignon berry composition was affected in a similar manner by soil temperatures that differed by 13 °C. TSS concentrations are also affected by light and the vine water status. Light is generally not a factor because there is usually a large enough leaf area and sufficient light levels to saturate this source to sink relationship. Sun-exposed Cabernet Sauvignon berries in the vineyard had higher TSS than shaded berries. This sunlight effect was attributed largely to an increase in berry temperature rather than an increase in the fluence rate perse.

A higher grapevine water status results in larger berry size and lower sugar concentrations and water deficit is known to increase sugar concentrations in Cabernet Sauvignon. However, nft hydroponic system temperature is thought to have the largest effect on sugar concentrations. Other transcriptomic data in the present study indicated that BOD berries were more mature at a lower sugar level than RNO berries. These included the transcript abundance profiles of genes involved in autophagy, auxin and ABA signaling, iron homeostasis and seed development. Many of these DEGs had an accelerated rate of change in BOD berries. While these transcripts are in the skins, they may be influenced by signals coming from the seed. In addition, there was a higher transcript abundance for most genes involved with the circadian clock in BOD berries. PHYB can regulate the circadian clock and PHYB activity is very sensitive to night temperatures ; PHYB reversion is accelerated to the inactive form at warmer temperatures. The inactivity of phytochrome promotes the expression of RVE1, which promotes auxin concentrations and seed dormancy. Thus, all things considered, it is likely that temperature and/or the temperature differentials between day and night significantly contributed to the differences in the rate of berry development and sugar accumulation in the two locations.Determining maturity of grapes is a difficult and error prone process. Reliable markers could aid in the decision of when to harvest the grapes. “Optimum” maturity is a judgement call and will ultimately depend on the winemaker’s or grower’s specific goals or preferences. A combination of empirical factors can be utilized including °Brix, total acidity, berry tasting in the mouth for aroma and tannins, seed color, etc. °Brix or total soluble solids by itself may not be the best marker for berry ripening as it appears to be uncoupled from berry maturity by temperature.

Phenylpropanoid metabolism, including anthocyanin metabolism, is also highly sensitive to both abiotic and biotic stresses and may not be a good indicator of full maturity. Thus, color may not be a good indicator either. Specific developmental signals from the seed or embryo, such as those involved with auxin and ABA signaling, may provide more reliable markers for berry ripening in diverse environments, but will not be useful in seedless grapes. Aromatic compounds may also be reliable markers but they will need to be generic, developmental markers that are not influenced by the environment. This study revealed many genes that are not reliable markers because they were expressed differently in different environments. One candidate marker that is noteworthy is ATG18G . Its transcript abundance increased and was relatively linear with increasing °Brix and these trends were offset at the two locations relative to their level of putative fruit maturity . ATG18G is required for the autophagy process and maybe important during the fruit ripening phase. It was found to be a hub gene in a gene subnetwork associated with fruit ripening and chloroplast degradation. Further testing will be required to know if it is essential for fruit ripening and whether its transcript abundance is influenced by abiotic and biotic stresses in grape berry skins.The ultimate function of a fruit is to produce fully mature seeds in order to reproduce another generation of plants. The ripe berry exhibits multiple traits that signal to other organisms when the fruit is ready for consumption and seed dispersal. In this study, we show that there were large differences in transcript abundance in grape skins in two different locations with different environments, confirming our original hypothesis. We also identified a set of DEGs with common profiles in the two locations.

The observations made in this study provide lists of such genes and generated a large number of hypotheses to be tested in the future. WGCNA was particularly powerful and enhanced our analyses. The transcript abundance during the late stages of berry ripening was very dynamic and may respond to many of the environmental and developmental factors identified in this study. Functional analysis of the genes and GO enrichment analysis were very useful tools to elucidate these factors. Some of the factors identified were temperature, moisture, light and biotic stress. The results of this study indicated that berries still have a “sense of place” during the late stages of berry ripening. Future studies are required to follow up on these observations. It appears that fruit ripening is very malleable. Manipulation of the canopy may offer a powerful lever to adjust gene expression and berry composition, since these parameters are strongly affected by light and temperature.The ability of a genotype to produce different phenotypes as a function of environmental cues is known as phenotypic plasticity . Phenotypic plasticity is considered one of the main processes by which plants, as sessile organisms, can face and adapt to the spatio-temporal variation of environmental factors . Grapevine berries are characterized by high phenotypic plasticity and a genotype can present variability within berries, among berriesin a cluster, and among vines . Berry phenotypic traits, such as the content of sugars, acids, phenolic, anthocyanins, and flavor compounds, are the result of cultivar and environmental influences , and often strong G × E interactions . Although grapevine plasticity in response to environmental conditions and viticulture practices may provide advantages related to the adaptation of a cultivar to specific growing conditions, it may also cause irregular ripening and large inter-seasonal fluctuations , which are undesirable characteristics for wine making . Due to its complex nature, hydroponic nft system the study of phenotypic plasticity is challenging and the mechanisms by which the genes affecting plastic responses operate are poorly characterized . In fact it is often difficult to assess the performance of different phenotypes in different environments . It has been suggested that genetic and epigenetic regulation of gene expression might be at the basis of phenotypic plasticity through the activation of alternative gene pathways or multiple genes . Epigenetics has been proposed as crucial in shaping plant phenotypic plasticity, putatively explaining the rapid and reversible alterations in gene expression in response to environmental changes. This fine-tuning of gene expression can be achieved through DNA methylation, histone modifications and chromatin remodeling . Small non-coding RNAs are ubiquitous and adjustable repressors of gene expression across a broad group of eukaryotic species and are directly involved in controlling, in a sequence specific manner, multiple epigenetic phenomena such as RNA-directed DNA methylation and chromatin remodeling and might play a role in genotype by environment interactions. In plants, small ncRNAs are typically 20–24 nt long RNA molecules and participate in a wide series of biological processes controlling gene expression via transcriptional and post-transcriptional regulation . Moreover, small RNAs have been recently shown to play an important role in plants environmental plasticity . Fruit maturation, the process that starts with fruit-set and ends with fruit ripening , has been largely investigated in fleshy fruits such as tomato and grapevine. These studies highlighted, among others, the vast transcriptomic reprogramming underlying the berry ripening process , the extensive plasticity of berry maturation in the context of a changing environment , and the epigenetic regulatory network which contributes to adjust gene expression to internal and external stimuli . In particular, small RNAs, and especially microRNAs , are involved, among others, in those biological processes governing fruit ripening . In this work, we assessed the role of small ncRNAs in the plasticity of grapevine berries development, by employing next-generation sequencing.

We focused on two cultivars of Vitis vinifera, Cabernet Sauvignon, and Sangiovese, collecting berries at four different developmental stages in three Italian vineyards, diversely located. First, we described the general landscape of small RNAs originated from hotspots present along the genome, examining their accumulation according to cultivars, environments and developmental stages. Subsequently, we analyzed miRNAs, identifying known and novel miRNA candidates and their distribution profiles in the various samples. Based on the in silico prediction of their targets, we suggest the potential involvement of this class of small RNAs in GxE interactions. The results obtained provide insights into the complex molecular machinery that connects the genotype and the environment.Two V. vinifera varieties Sangiovese , a red Italian grape variety, and Cabernet Sauvignon , an international variety, were grown side by side in three different Italian locations, representing traditional areas of Sangiovese cultivation in Italy with a long-standing wine making tradition. In order to reduce factors of variation, the age of the plants , the clone type , the rootstock , the cultivation techniques and the health status were the same among all the locations. Further details on the environmental conditions of the vineyards are provided in Supplementary Figure 1. Berries from four developmental stages were collected in two biological replicates, during the 2011 growing season, for a total of 48 samples . The four sampled stages corresponded to pea size , representing the first stage of berry development in this experimental plan, bunch closure also known as Lag Phase, 19–20 ◦Brix , which corresponds to 50% of sugar accumulation in berries, and harvest , when the berries are fully ripened and the onset of sugar accumulation is over. About 200 berries per each developmental berry stage were sampled from upper, central and lower part of cluster, both from sun exposed and shaded side and split in two biological replicates. Per each vineyard, the berries were collected from about 20 vines selected in a single uniform row and immediately frozen in liquid nitrogen and stored at −80◦C prior to analysis. The libraries were named using the initials of the vineyard where the berries were collected, followed by the initial of the cultivar and the developmental stage. For example, the sample containing berries of Sangiovese, collected in Montalcino at pea size, was named Mont_SG_ps.RNA extraction was performed as described in Kullan et al. . Briefly, total RNA was extracted from 200 mg of ground berries pericarp tissue using 1 ml of Plant RNA Isolation Reagent following manufacturer’s recommendations. The small RNA fraction was isolated from the total RNA using the mirPremier R microRNA Isolation kit and dissolved in DEPC water. All the steps suggested in the technical bulletin for small RNA isolation of plant tissues were followed except the “Filter Lysate” step, which was omitted. The quality and quantity of small RNAs were evaluated by a NanoDrop 1000 spectrometer , and their integrity assessed by an Agilent 2100 Bioanalyzer using a small RNA chip according to the manufacturer’s instructions. Small RNA libraries were prepared using the TruSeq Small RNA Sample Preparation Kit , following all manufacturers’ instructions. The quality of each library was assessed using an Agilent DNA 1000 chip for the Agilent 2100 Bioanalyzer. Libraries were grouped in pools with six libraries each . The pools of libraries were sequenced on an Illumina Hiseq 2000 at IGA Technology Services . The sequencing data were submitted to GEO–NCBI under the accession number GSE85611.