Horticulture

Many crystalline compounds have cleavage planes; that is, planes along which cracks propagate most readily. When such compounds are stressed beyond their yield strength, they tend to break up into pieces with characteristic shapes that inherit the anisotropy of the chemical bonds forming the crystal out of which they are composed. Indeed, this observation was a compelling piece of early evidence for the existence of crystallinity, and even atoms themselves. There exists a class of materials withcovalent bonds between unit cells in a two dimensional plane and much weaker van der Waals bonds in the out-of-plane direction, producing extraordinarily strong chemical bond anisotropy. In these materials, known as ‘van der Waals’ or ‘two dimensional’ materials, this anisotropy produces cleavage planes that tend to break bulk crystals up into two dimensional planar pieces. Exfoliation is the process of preparing a thin piece of such a crystal through mechanical means. In some of these materials, growing strawberries hydroponically the chemical bond anistropy is so strong that it is possible to prepare large flakes that are atomically thin . These two dimensional crystals have properties quite different from their bulk counterparts.

They do have a set of discrete translation symmetries, which makes them crystals, but they only have these symmetries along two axes- there is no sense in which a one-atom-thick crystal has any out-of-plane translation symmetries. For this reason they have band structures that differ markedly from their three dimensional counterparts. A variety of techniques have been developed for preparing atomically thin flakes from van der Waals materials, but by far the most successful has been scotch tape exfoliation. In this process, a chunk of a van der Waals crystal is placed on a piece of scotch tape. This separates the chunk of van der Waals crystal along its cleavage planes into two pieces of comparable size on opposite sides of the piece of tape. This process is repeated several times, further dividing the number of atomic layers in each chunk with each successive repetition. If we assume we are dividing the number of atomic layers in each piece roughly in half with each round, after N repetitions the resulting crystals should have thicknesses reduced by a factor of 2 1 N . This is enough to reduce each flake to atomic dimensions after a small number of repetitions of the process. This process is, of course, self-limiting; once the flakes reach atomic dimensions, they cannot be further subdivided. Graphite can be exfoliated through this process into flakes one or a few atoms thick and dozens of microns in width, a width-to-thickness aspect ratio of 105 .

As previously discussed, this process cannot be executed on every material. It depends critically on scotch tape bonding more strongly to a layer of the crystal than that layer bonds to other layers within the crystal. It also depends on very strong in-plane bonds within the material, which must support the large stresses associated with reaching such high aspect ratios; materials with weaker in-plane bonds will rip or crumble. In practice these materials are almost always processed further after they have been mechanically exfoliated, and the preparation process typically begins when they are pressed onto a silicon wafer to facilitate easy handling. Samples prepared in this way are called ‘exfoliated heterostructures.’ It is of course interesting that this process allows us to prepare atomically thin crystals, but another important advantage it provides is a way to produce monocrystalline samples without investing much effort in cleanly crystallizing the material; mechanical separation functions in these materials as a way to separate the domains of polycrystalline materials. Graphene was the first material to be more or less mastered in the context of mechanical exfoliation, but a variety of other van der Waals materials followed, adding substantial diversity to the kinds of material properties that can be integrated into devices composed of exfoliated heterostructures. Monolayer graphene is metallic at all available electron densities and displacement fields, but hexagonal boron nitride, or hBN, is a large bandgap insulator, making it useful as a dielectric in electronic devices.

Exfoliatable semiconductors exist as well, in the form of a large class of materials known as transition metal dichalcogenides, or TMDs, including WSe2, WS2, WTe2, MoSe2, MoS2, and MoTe2. Exfoliatable superconductors, magnets, and other exotic phases are all now known, and the preparation and mechanical exfoliation of new classes of van der Waals materials remains an area of active research. Once two dimensional crystals have been placed onto a silicon substrate, they can be picked up and manipulated by soft, sticky plastic stamps under an optical microscope. This allows researchers to prepare entire electronic devices composed only of two dimensional crystals; these are known as ‘stacks.’ These structures have projections onto the silicon surface that are reasonably large, but remain atomically thin- capacitors have been demonstrated with gates a single atom thick, and dielectrics a few atoms thick. Researchers have developed fabrication recipes for executing many of the operations with which an electrical engineer working with silicon integrated circuits would be familiar, including photolithography, etching, and metallization. I think it is important to be clear about what the process of exfoliation is and what it isn’t. It is true that mechanical exfoliation makes it possible to fabricate devices that are smaller than the current state of the art of silicon lithography in the out-of-plane direction. However, these techniques hold few advantages for reducing the planar footprint of electronic devices, so there is no meaningful sense in which they themselves represent an important technological breakthrough in the process of miniaturization of commercial electronic devices. Furthermore, and perhaps more importantly, it has not yet been demonstrated that these techniques can be scaled to produce large numbers of devices, and there are plenty of reasons to believe that this will be uniquely challenging. What they do provide is a convenient way for us to produce two dimensional monocrystalline devices with exceptionally low disorder for which electron density and band structure can be conveniently accessed as independent variables. That is valuable for furthering our understanding of condensed matter phenomena, independent of whether the fabrication procedures for making these material systems can ever be scaled up enough to be viable for use in technologies.The properties of crystals differ from the properties of atoms floating in free space because the atomic orbitals of the atoms in a crystal are close enough to those of adjacent atoms for electrons to hop between atoms. The resulting hybridization of atomic orbitals produces quantum states delocalized over the entire crystal with the capacity to carry momentum. This situation is shown in schematic form in Fig. 1.9A. For quantum states delocalized over the entire crystal, growing tomatoes hydroponically position ceases to be a useful basis. Instead, under these conditions we label electronic wave functions by their momenta, kx and ky. The atomic orbitals that prior to hybridization had discrete energy spectra now have energy spectra given by discrete functions of momentum, f. We call these functions electronic bands.Electrons loaded into the electronic bands of a two dimensional crystal will occupy the quantum states with the lowest available energies, so we can specify a maximum energy at which we expect to find electrons for any given electron density. We call that energy EF , the Fermi level. We can raise the Fermi level by adding additional electrons to the crystal, as shown in Fig. 1.9B. We have already discussed how two dimensional crystals naturally allow for manipulation of the electron density, and thus the Fermi level. We have also already discussed how the application of an out-of-plane electric field to a two dimensional crystal will change the structure of the atomic orbitals supported by that crystal.

It naturally follows that atomic orbitals so modified will produce different electronic bands, as shown in Fig. 1.9C. It is relatively straightforward to compute how electronic bands will respond to the application of a displacement field. We will be using the momentum and energy basis for the rest of this document; this basis is known as momentum space. The simplest experiment we can perform to probe the electronic properties of a two dimensional crystal in this geometry is an electronic transport experiment, in which a voltage is applied to a region of the crystal with another region grounded, so that electrical current flows through the crystal. We can check if the crystal supports any electrical transport at all, and if it does we can measure the electrical resistance of the crystal this way, in close analogy to how this is done for three dimensional crystals. Crystals will only accept and thus conduct electrons if there are available quantum states at the Fermi level; we call these crystals metals , and they can be identified in band structure diagrams by the intersection of the Fermi level with an electronic band . Crystals without empty quantum states at the Fermi level will not accept and conduct electrons , and they can be identified in band structure diagrams with crystals for which the Fermi level does not intersect with an electronic band . There exists a variety of other experiments we can perform on two dimensional crystals in order to understand their properties. Two dimensional crystals can support electronic transport in the in-plane direction if they are metals, as shown in Fig. 1.10. Capacitors can also support electronic transport in the out-of-plane direction, as long as that electronic transport occurs at finite frequency. The same structure that we use to modify the electron density and ambient out-of-plane electric field of a two dimensional crystal can also be used as a capacitive AC conductor, as illustrated in Fig. 1.11A. The conductance will depend only on the frequency at which an AC voltage is applied and the geometry of the parallel plate capacitor. However, if a two dimensional crystal is added in series, the capacitance of the top gate to the bottom gate may be substantially modified. If the two dimensional crystal is an insulator, electric fields will penetrate it and the capacitance between the two gates will not change. However, if the two dimensional crystal is a metal it will accept electrons and cancel the applied electric field, dramatically reducing the capacitance between the top and bottom gates and neutralizing the AC current through the capacitor. This technique can be used to measure the electronic properties of specifically the bulk of a two dimensional crystal; it is the property that was both calculated and measured in Fig. 1.2C and D. These two techniques are the bread and butter of the experimental study of two dimensional crystals, because they require only the ability to create stacks of two dimensional crystals and access to tools common to the study of all other microelectronic systems. We will discuss a considerable amount of electronic transport and capacitance data as well. However, the primary focus of this thesis will be on systems for which the nanoSQUID microscope can provide important information that is inaccessible to these techniques, and so we will discuss a few such systems next.Consider the following procedure: we obtain a pair of identical two dimensional atomic crystals. We slightly rotate one relative to the other, and then place the rotated crystal on top of the other . The resulting pattern brings the top layer atoms in alignment with the bottom layer atoms periodically, but with a lattice constant that is different from and in practice often much larger than the lattice constant of the original two atomic lattices. We call the resulting lattice a ‘moir´e superlattice.’ The idea to do this with two dimensional materials is relatively new, but the notion of a moir´e pattern is much older, and it applies to many situations outside of condensed matter physics. Pairs of incommensurate lattices will always produce moir´e patterns, and there are many situations in daily life in which we are exposed to pairs of incommensurate lattices, like when we look out a window through two slightly misaligned screens, or try to take pictures of televisions or computer screens with our camera phones. Of course these ‘crystals’ differ pretty significantly from the vast majority of crystals with which we have practical experience, so we’ll have to tread carefully while working to understand their properties.

Our studies suggest a combination of its transmission dynamics and how they are affected by management issues, such as the quantity of shade and the density of planting, plus a variety of control from above elements represent a source of control, which sometimes fails .Understanding the general structure of ecological communities has long been a central goal of ecology, from Haeckel to us. Empiricists commonly, and probably necessarily, focus on the community of X, which is to say an assemblage of species defined by some set of criteria: the fungal community of Lake Wobegon, the community of gall-forming insects of oak trees, the microbial community of the human gut, the community of four ciliate species, and so on. Theoreticians perhaps feel less constraint. In the present article, we have defined the community as the herbivores of the coffee plant and their associates, in which top-down control is the goal of management . The framing of regulation from above from theoretical ecology translates directly into biological control from agroecology. Indeed, in agroecology regulation from above is elementary, large round plant pots in that the top-down agents are frequently obvious . However, stopping at that level of understanding may obscure more than clarify, much as the simple phrase controlled from above may indeed obscure .

Precisely how that control is affected may involve many complicated interactions and contingencies, making, we argue, the framing of complex systems a necessary one. The fungus that attacks the scale is most efficient when the scale is hyper dense at a local level, something that cannot happen unless it is under the protection of a mutualistic ant, which deters the other predator , which, however, is able to take advantage of a spatial pattern that is self-organized through a Turing-like process, and so forth. Indeed, we argue that the understanding we claim to have of this system so far comes from detailed study, both empirical and theoretical, and, most importantly is dramatically enriched through the application of some of the concepts newly developed in the distinct field of complex systems. Almost 10 years ago, some of us published a summary of this overall system , suggesting that understanding it required more than just an identification of who eats whom. This update emphasizes that point. Our narrative in the present article is perhaps a bit heterodox. We study a very complicated system , and we seek to understand it through theoretical ecology. To some, at least in the recent past, this might imply a large-scale computer model or sophisticated data manipulation. Our approach is distinct, recalling the wisdom of Levins’ paper on the strategy of model building. We seek to understand, at a deep level, how this system works, not necessarily for the purpose of predicting its future state.

We offer theoretical propositions, many of which are stimulated by mathematical arguments, but we do not seek what postmodern thinkers would have called a “totalizing discourse” with a large-scale model. Rather, we seek to use recent advances in complex systems as a way of stimulating thought, with the mathematical models that go along with them as “educating our intuition,” as Levins urged frequently. The models themselves represent approximate metaphors for this complex reality, all fitting into a hierarchy of understanding , which is mainly qualitative even though originally formulated through mathematical reasoning.Furthermore, our claim that this is a complex reality is meant to imply something deeper than the obvious claim that it is complicated. It is a complex system. For example, if the only players in the system were Azteca, C. viridis, and A. orbigera, the system wouldn’t be exceptionally complicated , but it would be a complex system, because it would have a clear emergent property . Even adding the phorid would mean two predators and two prey, but the spatial pattern that emerges and the dependence of one system on a second system operative at a completely distinct time scale is an essential structural component of the system as a whole. The emergence would defy understanding if only the separate component parts were studied, which is to say if it were approached from a purely reductionist perspective. If the only players were the ants and the coffee berry borer, but the ants did not exhibit trait-mediated indirect interactions, the system would be complicated but not necessarily complex. This distinction between complicated and complex is important for our narrative. Because it is a complex system, it requires a more holistic approach to understand and manage, and there’s more potential for surprise .

A merely complicated system would not have these characteristics. That our model system is coffee is significant in several ways. First, traditional coffee management, with its characteristic shade trees, helps to create landscapes that are friendly to biodiversity conservation . It is a classic high-quality matrix for all sorts of animals and plants. Second, it involves a commodity that is of extremely high value, sometimes the main source of wealth for entire countries. Third, it is the basis of livelihood for millions of small farmers the world over. Fourth, when properly cultivated with shade, it joins other agroforestry systems in the worldwide struggle against climate change. Given all that, understanding the details of its operation would seem worthwhile, and marshaling recent insights from complex systems to anchor that narrative brings one of the classical questions of community ecology into focus as a practical issue. Consequently, besides being of potential importance for ecology, it makes ecology important for some practical aspects of this important crop. It is, for example, evident from only a qualitative understanding of the control from above system that a key element is the species of ant that nests in the shade trees and that, if those shade trees are eliminated , the whole control structure will be dramatically interrupted. Questions also arise about generality. Does this model system reflect something more general about the structure of control from above, or does it simply reflect interactions of this one particular system? First, most terrestrial systems have a spatial component involved, and framing the spatial component as one in which a subsystem operates to effectively create a spatial pattern in which other subsystems may operate is likely to apply frequently. Indeed, the idea of a predator–prey system generating a Turing pattern may be increasingly appreciated as more research programs interrogate the idea . Second, population dynamics unfolding on this space are likely to be nonlinear, and this nonlinearity will frequently be of the form that critical transitions lead to an alternative equilibrium within hysteretic zones, which may be multiple and constrain the herbivores above which control is being exerted . Third, the idea that multiple herbivores have their own suite of controlling factors is almost certainly true, but the idea that there will be connections, even if weak, with other subcomponents of the control from above, is likely to be characteristic. These three generalities encompass the complex systems topics of Turing pattern formation, critical transitions, hysteresis, chaos, basin boundary collisions, trait-mediated indirect interactions, and scale-dependent spatial processes, all of which are exemplified in our model system, plant pots round and certainly may be embedded in other systems of control from above. The message is not that these particular topics are essential but, rather, that control from above is not the one-dimensional process frequently imagined of a predator guild preying on a prey guild but, instead, a complex community of predators and parasites and diseases that interact with one another in complicated ways to eventually generate a self-organized system that exerts effective control over the herbivory. Much as one might say that the vertebrate circulatory system is responsible for bringing oxygen to each cell in the body, one might similarly simplify and say that natural enemies in the coffee agroecosystem are responsible for the regulation of potential pests. However, it is the heart, the veins, the arteries, exchanges across membranes, and so forth that tell the real story of how the delivery of oxygento our tissues actually happens. It is a complex system, the details of which are certainly of interest to health and healing.

Similarly, in our agroecosystem example, the subsystem that creates large-scale pattern sets the stage for a subsystem involving a predator and a disease that affect regulation of one pest, whereas the community structure of ants determines the efficiency of their predacious activities on a second pest and the disease that helps regulate the first pest is an antagonist to the third pest. This is all to say that yes, it is control from above, but that control is delivered through the ecological complexity of the community of natural enemies. It is misleading to suggest that listing the natural enemies and merely identifying them as such is sufficient. It is only through the lens of the reality of its state as a complex system that we may gain full appreciation of the ecological principle of top-down control, which then can be fully exploited in attempts to aid the management of this important agroecosystem. There is something of a conundrum in this narrative. Although it is clear that knowledge of all the ecological complexity could inform practical decisions that producers might want to make, is such detailed knowledge really necessary to provide useful advice to the farmer? If ecological knowledge of the particular system is primitive, could well meaning agroecological advisors give advice that will have unintended negative consequences? Post-WWII industrial agriculture enthusiasts embraced DDT and other pesticides creating the well-known pesticide treadmill that haunts us still today. Indeed, that is one of the issues that caused many environmentally conscious analysts to call for the science of ecology to be more actively embraced by agricultural planners. However, ecology is complicated. Secondary consequences cannot necessarily be predicted short of detailed study and the normal rules of thumb extrapolated from a few experiments or extralocal traditions could backfire. Perhaps the famous medical practitioner’s oath primum non nocere makes sense in agriculture as well. As farmers seek solutions to perceived problems on their farms, agroecologists rightly wish to use the science of ecology to help. However, frequently , ecological knowledge of the particular system is not very well understood because it is only recently that agroecological advocates have begun to break into the mainstream, and the basic research required to understand some of the vexing problems the farmers face has yet to be done. It is therefore common to use a few rules of thumb: avoid monocultures, don’t poison your natural enemies, maintain healthy soil, and so on. Such rules of thumb, on the basis of perceived ecological rules, for the most part make sense and probably conform well to the admonition primum non nocere. However, it is worth remembering the dust bowl, pest resurgence following pesticides, ocean dead zones, and other consequences that we live with today because a previous generation of farm advocates, equally sincere in their desires to help farmers, were prematurely confident in the ability of their tools to help the farmer.In recent years, a variety of novel two-dimensional van der Waals magnets have been discovered, founding the active field of 2D magnetism. Among these prospective compounds, binary chromium tellurides Cr1–δTe are attractive owing to their rich magnetic properties, as well as inherent chemical and structural compatibility when forming heterostructures with other topological systems, such as tetradymite-type topological insulators or chalcogenide-based Dirac/ Weyl semimetals. Furthermore, the broken time-reversal symmetry and spin-orbit coupling offer unique opportunities for the interplay between spin configurations and reciprocal-space topology. In this regard, ferromagnetic Cr2Te3 with strong perpendicular magnetic anisotropy is an intriguing platform to host non-trivial topological physics, particularly for the high-quality thin films grown by molecular beam epitaxy . An important consequence of the band topology in Cr2Te3 is the Berry curvature underlying the anomalous Hall effect. The intrinsic AHE is topological in nature and a hallmark of itinerant ferromagnets, which has also been observed in more exotic systems even without a net magnetization, such as spin liquids, antiferromagnets, and Weyl semimetals. When SOC coexists with long-range magnetic order, the Berry curvature can be significantly influenced near avoided band crossings, rendering the system an incredibly rich playground combining topology and magnetism.

The size of this group is larger than the number of new residents who moved from any state in the union. Among current U.S. residents who moved to Colorado, the largest shares arrived from California , Texas , Florida , Arizona , Illinois , New York , Virginia , Missouri , Georgia , and North Carolina . The most popular destinations for Coloradoans moving out of state were California, Texas, Washington, Arizona, and Florida. Because most of these population gains have occurred along the Front Range and in the Denver metro area, the housing market in the state’s largest city remains one of the least affordable in the nation. Among U.S. metropolitan areas, the Denver housing market is the most expensive of any city not located in a coastal state. Median home prices in Denver exceeded the national average by more than $100,000 in 2017, and the average price of a single-family homes sold in Denver exceeded $500,000 for the first time in February 2018 . Many residents priced out of the home-buying market have also encountered affordability challenges in the rental market. According to one real estate research firm, plant plastic pots apartment rental rates in Denver have increased by nearly 50 percent since 2010. This is the largest increase in rent for any city outside of the San Francisco Bay area. Increased demand for housing units has spurred development in the Denver metro area as evidenced by the construction of 23,000 new apartments between 2016 and 2017 .

A growth in the number of available housing units has increased the vacancy rate in the metro area to 6.4 percent in 2018. This figure is the highest seen in Denver since the 2009 recession. It has also resulted in a modest reversal in the long-term trend of rent increases, as the median monthly rental price in Denver decreased to $1,353 from $1,370 at the end of the third quarter in 2017. Table 1 reports selected data from the Census Bureau and Bureau of Labor Statistics  for Denver County, Colorado, and the United States. As seen in Table 1, per capita and household income in Denver and Colorado remain above the national average. The OSBP projects personal income in Colorado to outpace the national average of 3.1 percent in 2017 with a 5.4 percent growth rate .The 2017 population estimate for Colorado is slightly greater than 5.6 million people. Colorado remains one of the fastest growing states in the nation, although population growth has slowed during the past year. After reaching nearly 2 percent in 2015, population growth fell to 1.4 percent in 2017. This figure remains twice the national growth rate of 0.7 percent. The state projects the number of residents to increase to 5.8 million by the 2020 Census. For most racial minority groups, the state is less diverse than the nation as a whole. Native Americans are one exception to this, as the share of those with American Indian heritage is slightly greater in Denver County and Colorado compared to their percentage nationally . The proportion of residents who identify as Hispanic or Latino is also greater in Colorado and Denver County than in the entire country .

The Census estimates that 13.2 percent of current U.S. residents were born abroad. Among all Coloradoans, the percentage of foreign-born residents is about 3 percentage points lower, while the percent foreign born residing in Denver County is nearly 3 percentage points greater.The BLS estimates that Colorado’s unemployment rate in December 2017 was 3.0 percent, while Denver’s unemployment rate was 2.9 percent . Earlier in 2017, the state unemployment reached a record low of 2.3 percent. The other six metropolitan statistical areas tracked by the BLS averaged an unemployment rate of 3.4 percent. Fort Collins had the lowest unemployment rate at 2.5 percent. Only two cities had an unemployment rate greater than the national unemployment rate of 4.1 percent—Grand Junction and Pueblo . Personal income growth among state residents reached nearly 8 percent in 2014, but growth slowed over the next two years. This trend reversed in 2017 as personal income growth increased to 5.4 percent. This growth rate exceeds the national rate by more than 2 percentage points . According to the OSPB, per capita income and wage growth in Colorado over the past year also outpaced the national figures. With regard to party registration in Colorado, voters were nearly evenly divided among Democratic, Republican, and unaffiliated categories at the time of the 2016 election. According to voter registration data from the Secretary of State’s office, the state had nearly 3.3 million active voters in November 2016. A plurality of these voters registered as unaffiliated , while the share of Democratic and Republican were nearly equivalent. In February 2018, the number of active voters decreased relative to November 2016 by about 1.7 percent.

This may be partially attributable to controversy surrounding President Donald Trump’s Commission on Voter Fraud, which made data privacy a concern after the commission requested, “voluminous information on voters, including names, addresses, dates of birth, political affiliations and the last four digits of Social Security numbers, along with voting history” . Fifteen months after the 2016 election, the proportion of unaffiliated voters in the state increased to 36.3 percent, while the share of Democrats and Republicans decreased to 31.1 percent and 30.8 percent, respectively . This change is likely driven by primary election reforms approved by voters in the 2016 election. Voters overwhelmingly approved Proposition 107, which adopted a presidential primary in lieu of the existing caucus system, and Proposition 108, which allowed unaffiliated voters to participate in the party primary of their choice. Previously, unaffiliated voters were prohibited from participating in any primary elections or caucus meetings. Because this reform allows unaffiliated voters to participate in the primary of their choosing, it appears that a substantial number of Coloradans changed their party registration status to take advantage of this new opportunity.Colorado’s economic trajectory remains generally positive. In late March, the Governor’s Office of State Planning and Budgeting released its revised economic forecast. The report summarized the condition of Colorado’s economy by stating, “Colorado’s economy is on solid footing with strong employment growth and expectations of an ongoing expansion. New business formation continues to grow, while Colorado oil production is at record levels. Although much of the state’s economic growth has occurred along the Front Range, stabilizing farmland values and increases in energy prices and production have recently supported rural areas as well. Looking forward, higher costs of living and tight labor market conditions are expected to constrain further growth through the forecast period” . The OSPB characterized the 3.1 percent increase in General Fund revenue in the 2016–2017 fiscal year as decreases in oil, gas, and other commodity prices as the catalysts for the three-year decline in corporate tax revenue beginning in 2014. After falling 21.9 percent in the third year of this recent decline, corporate tax revenue is estimated to increase by 38.6 percent in the current fiscal year. In addressing how federal tax reform may affect state corporate tax revenue, the OSPB projects continued growth in corporate tax revenue, but cautions that “future increases will be constrained by higher business costs, especially for employee compensation, black plastic pots which will reduce profit margins and result in lower tax liabilities” . Taxes from the state’s legal marijuana market continue to grow. According to data released by the Colorado Department of Revenue, total marijuana sales surpassed $1.5 billion in 2017 . Of this total, $1.09 billion in sales came from the retail market, while medical marijuana sales were approximately $0.42 billion. Table 2 reports annual marijuana sales and tax revenue data. Sales have increased each calendar year since the retail market began operation on January 1, 2014, but the growth rate has gradually declined each year. In 2014, marijuana sales totaled nearly $680 million. This figure increased by 45.7 percent in 2015 to $990 million. Sales increased by 31.3 percent and surpassed the $1 billion mark for the first time in 2016 to reach a total of $1.3 billion. Total sales in 2017 amounted to $1.5 billion, which represents a 15.3 percent increase.

If the trend continues, Colorado could expect to see sales numbers plateau since the growth rate has decreased by an average of 15 percentage points each year. Sales data through September 2018 indicate sales of nearly $1.16 billion. The monthly sales average for 2018 puts the state on track for an annual marijuana sales total of about $1.54 billion, which would be the largest sales amount to date and constitute a slight increase of 2.4 percent from 2017.State revenue from the medical and retail marijuana markets in the form of taxes and fees reached nearly $250 million in 2017. This is a 27.8 percent increase from 2016 following a 48.5 percent increase in tax revenue from 2015 to 2016. During the first 10 months of 2018, the state reported tax and fee revenues of $223 million. A linear projection of tax revenue for the remainder of the year suggests that the annual tax revenue would reach nearly $268 million, an increase of about 8 percent from the prior year. The largest annual increase in marijuana taxGovernor Hickenlooper’s November letter to the JBC highlighted education and the criminal justice system as two areas where additional investments were necessary . The proposed and enacted funding levels for state departments and changes from the prior year’s budget are reported in Table 3. In total, the General Fund budget request is $11.5 billion, which is an increase of $292.1 million from the prior year. This represents a 2.1 percent increase, while total fund spending is proposed to increase by 3.7 percent to $30.5 billion. The budget signed by the governor into law on April 30 authorized a General Fund budget of $11.42 billion and $30.63 billion in total fund spending . As the data reported in Table 3 indicate, the governor’s budget proposes increasing the funding levels of all state department relative to last year, exempting the Department of Local Affairs and Department of Labor and Employment. The final budget ultimately increased the funding levels for all state departments except for Labor and Employment. The median proposed increase is 5.1 percent, while the enacted budget imposes a median departmental spending increase of 5.7 percent. In addition to funding existing government programs, the governor’s budget proposed new spending initiatives for the Department of Health Care Policy and Financing , Higher Education , K-12 Education , Department of Corrections , and Department of Human Services , in addition to an increase in the General Funds reserve .Relative to last year’s budget, Governor Hickenlooper’s budget request for the 2018– 2019 fiscal year is much more optimistic. This positive tone is a product of a good economic climate and the successful passage of a major budgetary reform in the prior session. Transitioning hospital provider fees into a government enterprise fund made these funds TABOR exempt. In the 2016–2017 fiscal year, hospital provider fee revenue was $654.4 million. In future years, these funds do not count toward the TABOR revenue cap. The governor emphasized the importance of this reform by noting that “The passage of S.B. 17-267 has materially and positively changed the State’s financial outlook compared with one year ago, when the request had to close a $500 million funding gap in the General Fund” . Beyond this important budgetary reform, Governor Hickenlooper also lauded the strong upward trajectory of the state’s economy by claiming, “Colorado’s economy continues to outperform nearly every state and the national economy overall” . Statewide unemployment remains low and job creation numbers are strong—approximately 53,000 new jobs are projected for 2018. As the state’s population continues to grow, Hickenlooper’s budget request reflects spending priorities to address increased demand for certain state services. Increases in education spending are one of the most notable aspects of the budget request. K-12 spending is proposed to increase by $84.6 million, which represents an increase of $343.38 per student. According to state estimates, the 178 school districts in Colorado currently serve the educational needs of more than 865,000 students. The Department of Education projects that enrollments will grow by 0.7 percent in the current fiscal year. The education budget request remarks that “additional funding proposed by the Governor will allow local districts and charter schools to decide how to best improve the education opportunities of their students.” As required by Amendment 23, the inflationary increase to the State Education Fund amounts to $8.9 million.

In addition to the worst public health crisis in a century and a major economic recession, the year 2020 also witnessed a highly charged presidential election unlike any the nation has ever seen.The global COVID pandemic had less of an effect on election administration in Colorado— relative to other states—because Colorado has employed a vote by mail system since 2013. Most of the nearly 3.3 million ballots in the 2020 election were submitted by mail or returned to ballot drop-off sites. Turnout among active voters was 86.9 percent, which corresponds to a voting eligible turnout rate of 76.4 percent. The results of the 2020 elections in Colorado produced few legitimate surprises. The Joe Biden and Kamala Harris ticket carried the state with 55.4 percent of the vote to 41.9 percent for Donald Trump and Mike Pence. The party composition of the state’s House delegation remained unchanged, with four Democrats and three Republicans elected to 2-year terms. Political novice and business owner Lauren Boebert defeated Republican incumbent Scott Tipton in a primary election for the state’s third congressional district.

Boebert’s firebrand political style, characterized by incendiary rhetoric and enthusiastic allegiance for President Trump, large plastic growing pots catapulted her to a surprising 9-point primary win over Tipton, who had served five terms in the House. Boebert’s brash campaign included vociferous support for the second amendment and tacit backing of various QAnon conspiracy theories. Despite numerous campaign controversies, Boebert emerged victorious in the general election where she defeated Democrat Diane Mitsch Bush by 6 percentage points in Colorado’s most competitive congressional election of 2020. Party control of the General Assembly remained with the Democrats, who maintained a 41-24 advantage in the House and increased the size of their majority by an additional seat in the Senate where party division stands at 20 Democrats and 15 Republicans. There was not a governor’s election in the 2020 election cycle following the election of Democratic Governor Jared Polis in 2018. The contentious presidential election took center stage as a result, although neither Trump nor Biden held a campaign event in Colorado, and the Democrats ultimately won the presidential popular vote in the state by a substantial margin. As Colorado solidified its status as a blue state in the Electoral College, many polls suggested a more competitive race for U.S.

Senate between Republican incumbent Cory Gardner and former governor John Hickenlooper. Following a short-lived candidacy for the Democratic presidential nomination, Hickenlooper launched a bid to unseat the one-term incumbent Gardner, viewed by many as the most vulnerable Senate Republican running for reelection. Gardner was widely considered the second most vulnerable incumbent Senator in this cycle after Democratic Senator Doug Jones from Alabama. Hickenlooper emerged victorious as expected on election night, defeating Gardner by nearly 10 percentage points receiving 53.5 percent of the vote. Several economic and budget related ballot measures appeared on the statewide ballot in 2020, including a proposal to reduce the state’s income tax rate from 4.63 percent to 4.55 percent, which passed with robust support. Voters also approved Amendment B, which repealed the Gallagher Amendment that established residential and non-residential property tax rates in the Colorado Constitution. Although tax increases have historically been unpopular in the Centennial State, the passage of proposition EE raised taxes on cigarettes, tobacco, and other nicotine products such as e-cigarettes and vaping cartridges. Stemming from recent controversies surrounding the use of state fees to circumvent constitutional requirements imposed by the Taxpayer Bill of Rights , Proposition 117 will require voter approval of future state enterprise funds with projected or actual fee revenue exceeding $100 million within its first five years. The passage of Proposition 118 created a new family and medical leave program funded by a payroll tax split by employers and employees that will allow workers to take up to 12 weeks of paid leave. Colorado’s historical passage rate is about 35 percent, but in 2020 a total of nine of the eleven statewide measures passed.

The only two that failed were Proposition 115, which sought to ban abortion 22 weeks following gestation, and amendment C, which would have reformed lottery regulations for charitable organizations. In sum, however, statewide election outcomes resulted in an uncharacteristic number of policy changes, many with important economic consequences. Unified Democratic government remained in place in Denver, which set the stage for the 2021–2022 budget cycle. The three sections that follow discuss the state’s demographics, revenue, and spending.The Census Bureau estimates the state’s population in 2020 to be 5.7 million residents, which is an increase of 14.8 percent from 2010. This population increase has outpaced growth in new housing units over the past decade. Only Utah and Texas have experienced larger growth rates since 2010. Colorado joined Oregon, Montana, Texas, Florida, and North Carolina as the six states that gained at least one seat in the House of Representatives following reapportionment. Because of the redistricting reforms approved by voters in 2018, boundaries for all eight congressional districts will be drawn by an independent redistricting commission. Table 1 presents further demographic and economic data from the Census Bureau and Bureau of Labor Statistics  for Denver County, Colorado, and the United States. The pre-pandemic unemployment rate in Colorado commonly fluctuated between 2 and 3 percent. At the beginning of 2020, the state’s unemployment rate stood at just 2.7 percent. Job losses from the COVID recession produced double-digit unemployment figures for three consecutive months in the spring, with a maximum of 12.1 percent unemployment in April . Industries disproportionately affected included retail trade, accommodation, hospitality, food service, arts, entertainment, and recreation. Similar to patterns in other states, job losses in Colorado were sharper and longer lasting among lower-income earners. State economists reported a change in employment among lower-income individuals of about 35 percent from January to May 2020. After reaching an apex in April, unemployment figures have improved or remained steady in each month thereafter. BLS data indicate that Colorado unemployment returned to single digits in July 2020 and fell below 7 percent in October. Unemployment levels in the state steadily decreased by a monthly average of 0.14 percent, reaching a pandemic low of 5.1 percent in November. According to one estimate, the state had recovered about 85 percent of the COVID job losses, which is better than the national average, but only one metropolitan area—Colorado Springs—had so far exceeded pre-pandemic employment levels . The civilian labor force participation rates of 72.2 percent in Denver County and 67.6 percent in Colorado are both better than the national figure of 63 percent.Per capita and household income levels in Colorado and Denver County remained greater than the national average. Median home prices in the state and county were likewise substantially above the median home price in the United States . The influx of new residents coupled with high real estate prices has kept affordable housing a persistent problem, particularly in the Denver metro area and many mountain communities across the state. For example, large plastic pots the Crested Butte city council recently declared a housing emergency while debating policy reforms to address the affordable housing crisis, including a one-year moratorium on short-term rental licenses. The city was also considering other measures, including easing restrictions on camping and imposing greater taxes on rental or vacation properties. Other innovative efforts to increase the availability of affordable housing in places such as Ouray, Ridgway, and Telluride included revising zoning ordinances to permit greater resident density. Efforts also included designating housing as only available to the local workforce by attaching income and deed restrictions to certain property developments. The rapid growth of the state’s population has also had effects on the composition of the Colorado electorate. According to the state’s active voter rolls, a plurality of voters in Colorado are registered as unaffiliated . Thirty percent of the electorate is registered Democratic, and Republican registration places third at 27.3 percent .

The proportion of unaffiliated voters increased by nearly 5 percentage points since the 2016 election. Contributing factors to this increase are election reforms passed in 2016, which replaced the existing caucus system with a presidential primary, and a second reform that allowed unaffiliated voters to participate in the party primary of their choice. Unaffiliated voters were previously prohibited from participating in primary elections or caucus meetings. Voters in Colorado who self-identified as unaffiliated or independent voted 55-40 in favor of the Democratic ticket in the 2020 presidential election, while the split among ideological moderates was 65-30 in favor of Biden/Harris . Regarding the racial and gender composition of the electorate, exit poll data from 2020 estimate white voters to constitute 79 percent of the electorate. Latino and Latina voters were the largest non-white voting bloc at 12 percent of voters, while the proportion of Black voters was about 3 percent . These numbers, however, do not reflect the actual size of racial minority groups in the state. According to Census data, Colorado’s Hispanic population is nearly 22 percent, the Black population is 4.6 percent, and the Native American population is 1.6 percent. A plurality of the state’s voters in 2020 were white women . An additional 35 percent of ballots were cast by white men. White women voted Democratic in the presidential election by a 2-1 margin, while Biden received 49 percent of the state’s white male vote to Trump’s 48 percent .Decreased economic activity and substantial job losses beginning in March 2020 resulted in a projected decline in state revenue for the first time since 2009. The Office of State Planning & Budgeting revenue forecast projected an overall decrease in state revenue of 4.7 percent to $12.3 billion in the 2020–21 fiscal year . This was a substantial decrease in state revenue but an improvement upon earlier projections, which previously estimated a percentage drop in revenues exceeding double digits. Revenues for the 2021–22 fiscal year, however, are predicted to return to pre-pandemic levels. Prior to the COVID pandemic, individual income tax revenues grew by nearly 5 percent. The projected decline in income tax revenue during the first year of the pandemic was 8.1 percent totaling $8.2 billion . The OPSB projects income tax revenue to rebound in 2022 and beyond—however, not reaching pre-pandemic growth levels, claiming that the “reduction is expected to be short-lived and mitigated by the fact that low income earners were disproportionately impacted by the pandemic” . Such a turnaround would be critical to Colorado’s economic recovery as income taxes represent about two-thirds of General Fund revenues. Corporate income tax collections were also deeply affected by the pandemic and fell by 20.8 percent to $728.3 million in the 2019–20 fiscal year. Economists predict that corporate taxes will continue in a downward trajectory by an additional 10 percent in 2021. Unlike income and corporate taxes, sales tax revenues increased during the pandemic on account of increased online shopping and a gradual rebound in retail sales. Overall, sales tax revenue increased 4.7 percent in 2020, with state projections for an additional 6.5 percent this year. Beginning this year, Colorado will also generate a greater amount of cigarette tax revenue, which equaled $32.5 million in 2020, following the passage of proposition EE. Cigarette taxes increased from 84 cents to $1.94 per pack in 2021 and will reach $2.64 in 2027. Tax rates on other tobacco products, including vaping products containing nicotine, will similarly increase to 62 percent in 2027.Governor Polis’ budget proposal for the 2021–2022 fiscal year proposed a total funds budget of $35.4 billion, which is a 9 percent increase from the prior year. The General Fund spending budget of $13.6 billion is a 20 percent increase of $2.25 billion. His letter to the Joint Budget Committee began by noting the extraordinary challenges and opportunities brought about by the global pandemic, an economic recession, the Black Lives Matter movement for racial justice, and historic wildfires and other natural disasters. Introducing his second budget since the pandemic began, Polis stated, “Each of these crises has exposed weaknesses in our systems of public health and health care, commerce, and education, and highlighted the urgency to take bold climate action in order to protect our Colorado way of life. To meet the moment and to build resiliency for the future, we must make critical investments to position Colorado for a strong, inclusive, and rapid recovery. Colorado will not only recover from these crises, but together we will build back stronger than before” .

The slight difference in mapping rate between the two references can be explained by either the absence of some low-expression transcripts in the ISNT or the residual error rate in isoform sequences. Gene expression analysis using the ISNT as reference showed similar results compared to the Cabernet Sauvignon genome assembly, with a very high correlation of expression level and differential gene expression, and with similar global transcriptomic changes. However, we observed differences in the number of expressed and differentially expressed features that depend on the reference used. Those differences could be explained by the diploid phasing of the Cabernet Sauvignon genome assembly and that multiple ISNT transcripts might correspond to a single gene locus. Nonetheless, similar relative amounts of Biological Process GO terms were found among the differentially expressed genes, confirming that the transcriptome obtained using Iso-Seq captured the transcriptional reprogramming underlying the main physiological and biochemical changes during grape berry development. In addition, gene expression analysis revealed that some private isoforms are significantly modulated during berry development, square pots for plants indicating that in addition to identifying the private gene space, the ISNT reference makesit possible to observe its expression.

In conclusion, this study demonstrates that Iso-Seq data can be used to create and refine a comprehensive reference transcriptome that represents most genes expressed in a tissue undergoing extensive transcriptional reprogramming during development. In grapes, this approach can aid developing transcriptome references and is particularly valuable given diverse cultivars with private transcripts and accessions that are genetically distant from available genome references, like the non-vinifera Vitis species used as rootstocks or for breeding. The pipeline described here can be useful in efforts to reconstruct the gene space in plant species with large and complex genomes still unresolved.Grapevines are susceptible to the highest number of pathogens to infect a single crop, with over 70 viruses detected. In 2008, a new virus was first observed in Napa County, California, which economically threatened grapevines: grapevine red blotch virus. This virus is the causative agent of grapevine red blotch disease, which has been identified in vineyards across the United States, Canada, Argentina, Mexico, South Korea, and India . Reports indicate GRBV primarily spreads through propagation material and secondarily through an insect vector. Spissistilus festinus was shown to successfully transmit GRBV in greenhouse settings, yet this has not been replicated in vineyards to date.

GRBV has been identified as a virus from the Geminiviridae family containing a circular single-stranded DNA genome similar to other geminiviruses. GRBD expresses symptoms of reddening of leaf blades and margins, with reddening of the primary, secondary, and tertiary veins in red grape cultivars. GRBV affects grapevines in various ways. For example, leaves on infected vines show increased levels of sugar, phenolics, particular amino acids, and enzymatic activity related to plant defense, as well as a reduction in carbon fixation. However, the most damaging are the effects on grape composition which has been shown to be detrimental to final wine quality. GRBV delays ripening by decreasing the accumulation ofsugar and anthocyanin in berries, potentially due to the impairment of translocation mechanisms. The virus has variable impacts on primary and secondary metabolites, specifically phenolic and aroma compounds. In summary, detrimental economic impacts to vineyards in the United States could reach $68,548/ha with vine removal being the only current method of alleviation. Consequently, recent research has strived to understand the effects and functioning of GRBV to establish mitigation strategies to alleviate the impact on grape composition and wine quality. The first phase is characterized by cell division and production of seeds, as well as synthesis of tannins and organic acids. The second phase is characterized with the onset of veraison, which is when the grape berry begins to soften and change color. The final and third phase is berry engustment/ripening, where berries increase in size, sugar accumulates, acidity declines, and secondary metabolites such as anthocyanins and aromatic compounds are synthesized inside the berry.

Studies have shown that volatile compounds such as terpenoids and C6 compounds, begin to accumulate in berries after veraison, and are controlled by numerous factors. The synthesis of these compounds in berries is also altered by external factors such as light exposure or pathogens. In addition, these secondary metabolites are crucial to grape growers and winemakers due to their importance in the quality of a final wine. Grape maturity has shown to be a key driver in the composition of a final wine, where later harvested fruit produces wines with lower concentrations of C6 alcohols and higher in concentration of esters. However, the impacts of GRBV on volatile compound abundance in harvested grapes has not been investigated. A plant’s genetic material may influence susceptibility to viral infections. Additionally, rootstocks can impact grapevine physiology and impact the overall composition of a grape berry. For instance, rootstock 110R causes high vigor and high drought tolerance in grapevines; whereas 420A is a rootstock of low to moderate vigor and low drought resistance. Vigor, resulting in greater shoot length and hence leaf area, may impact net carbon assimilation and the translocation of metabolites into the berry, consequently affecting the final wine composition. These hydric differences affecting carbon metabolism in rootstocks can also impact the plant-pathogen interactions. Therefore, it is plausible that severity of GRBD symptoms will be dependent on the interaction between scion cultivar-rootstock. However, this has not been fully investigated. Macro and micro climate fluctuations may also be a factor in pathogen-plant interactions, and should be considered. This study investigated the impact of GRBV on the biosynthesis and accumulation of primary and secondary metabolites in grape berries throughout ripening and at harvest. Additionally, the influence of seasonal and genotypic factors on disease expression within grapevines were studied.Figures 1–3 depict sugar accumulation, anthocyanin levels, TA, and pH through ripening. Sugar accumulation was determined by converting ◦Brix to mg of sugar per berry. Anthocyanin content was lower in RB. grapes when compared to RB grapes for both years and rootstocks during ripening . However, the degree of impact varied depending on season and rootstock. In 2016, both rootstocks were equally impacted throughout ripening regarding sugar accumulation and anthocyanins levels. However, in 2017 sugar accumulation was generally not significantly impacted by disease status. In 2017, grape anthocyanin levels were more significantly impacted for infected vines on 110R rootstock than 420A rootstock, whereas rootstock impact was less apparent in 2016. At harvest in 2016 , CS 110R and 420A rootstocks, respectively, had a 2% and 11% decrease in anthocyanin content and a 12% and 18% decrease in sugar content in RB grapes when compared to RB grapes. In 2017, atharvest , anthocyanin content was 35% and 11% lower in RB when compared to RB , and sugar content was 9% and 7% lower, for 110R and 420A, respectively. By plotting ◦Brix over ripening and fitting a linear trendline, it is possible to compare the rate of ripening for RB and RB grapevines .

As indicated by the slope of the best fit line, the rate of ripening was always higher for RB data vines when compared to RB data vines, with the exception for CS 420A in 2017. In addition, large square plant pots the rate was also lower in 2016 than 2017 across virus status and rootstocks. Interestingly, the difference in the rate of ripening between RB and RB data vines was larger in 2016 than in 2017 which correlates to the larger differences in accumulated sugar at harvest. In 2017, the rate was lower for CS 420A than CS 110R across virus status.The volatile compound profiles of RB and RB grapes were determined in both 2016 and 2017 seasons . PCA was performed to plot the variability between RB and RB grape samples . Between 80.6- 94.5% of the variance is explained by the PCA in Figures 4 and 5. For CS 110R, only the significantly different volatile compounds between RB and RB are plotted. There were, respectively ten and nine significant volatile compounds that explained the difference between treatments for CS 110R in 2016 and 2017. For CS 420A, the volatile compounds that contributed most to the variance of the PCA were plotted, due to few volatile compounds being significantly different. This selection was based on the squared cosine which shows the importance of the volatile compounds to explain the variance in the data. For CS 420A, in 2016, only cis-3-hexen-1-ol was significantly different, and in 2017, only β-linalool and β-citronellol were significantly different with a α level of 0.95.After veraison, volatile compound accumulation begins in grapes and changes through ripening. However, the impacts of GRBV on grape volatile compounds have not been investigated. We found 35 different aromatic compounds in grapes from the two rootstocks over two seasons, of which 24 were similar between the two years studied . In 2016, across rootstock, RB grapes were generally lower in volatile compound levels than RB , except for C6 compounds such as 2-hexenal, hexanal, cis-3-hexen-1-ol, trans-2-hexen-1-ol, and hexanol . These C6 volatile compounds are synthesized in the grape skin through the lipoxygenase pathway, are generally responsible for green or grassy aromas and accumulate in CS grapes up to 18 ◦Brix . With the exception of hexanol, the levels of these compounds begin to significantly decrease thereafter, with a 67% decrease in grapes at 25 ◦Brix when compared to grapes at 18 ◦Brix. These observations correlate with the common finding that GRBV causes a delay in ripening, with green aromas being present and correlated with the lower sugar accumulation. On the other hand, RB grapes were highly correlated with monoterpenes such as limonene, ß-myrcene, α-terpinene, geranial and p-cymene , which are responsible for floral and fruity aromas. These compounds have been associated with CS grapes at harvest and are known to increase through grape ripening and may decrease at over ripeness. In addition, RB grapes were also highly correlated with esters such as ethyl acetate, ethyl hexanoate, hexyl acetate, and ethyl butyrate. Although ester formation is mainly related to yeast or bacteria metabolism during winemaking , grapes are also known to synthesize esters. Anthraniloyl-coenzyme A :methanol acyltransferase is known to be responsible for the formation of methyl anthranilate in grapes and it is also classified as an ester-forming acyltransferase, which could be responsible for the formation of esters in grapes. The esters found in the current work that in general related mostly to RB grapes, are known to produce red and black fruit aromas. Collectively, these results confirm that RB grapes underwent normal ripening processes and produced more fruity aromas, while RB grapes at harvest have aroma characteristics more related to early ripening stages.Results indicated that in 2016 GRBD had a larger impact regarding sugar accumulation, pH, TA, and final sugar content than in 2017. In addition, the harvest dates were two to three weeks later in 2017 than in 2016. These observations can potentially be explained by the difference in temperature between the two seasons. In 2017, Napa County experienced a heat wave from 26 August–11 September, where nine days were over 35 ◦C, and four days were over 40 ◦C. The cumulative growing degree days for both years can be seen in Figure 1e. Extreme heat conditions during grape maturation have been shown to inhibit enzymatic activity and halt the biosynthesis of metabolites inside the grape berry. Inhibition of these processes due to heat leads to decreases in sugar accumulation and increases in acidity in healthy fruit. This is thought to be caused from a decrease in rate of translocation of sugars from leaves to fruit, through the reduction in photosynthesis at temperatures greater than 30 ◦C. The rate of ripening in 2017 was faster than 2016prior to the heat spike . However, during the heat spike in late August to harvest, sugar accumulation plateaued resulting in extended harvest times in 2017. In addition, research has shown that temperature can alter virus-induced gene silencing which is triggered with the infection of a virus as a plant-derived defense mechanism to downregulate the genes of interest. Previous work on other plant species infected with a geminivirus has shown that the extent of gene silencing is related to temperature. Specifically, Chellappan et al. showed that temperatures over 30 ◦C induced gene silencing, which interfere with gene expression, resulting in decreases in viral DNA accumulation and decreases in symptoms.

It could be hypothesized that in the 1st year bearing 100% of the clusters, while not showing a reduction in soluble solids, may have taken a toll on plant reserves. However, neither root mass nor starch content were impacted by the crop level in our work. Palliotti and Cartechini performed cluster thinning on three varieties over three seasons and found that cluster thinning did not affect must soluble solids. In years where rainfall was more abundant , results of cluster thinning were compensated with larger berries , with compensation of berry size similar to our results. Precisely, these kinds of results are those that disturbed the correlation between leaf area to fruit mass and berry total soluble solids . This suggested that larger berries may offer higher resistance to increases in berry total soluble solids regardless of leaf area to fruit mass . This hypothesis was supported by the response of grapevines submitted to water deficits that had smaller berries with higher soluble solids despite having much lower carbon assimilation rates .Development, blueberry packaging which encompassed the timing of all physiological events recorded , was delayed clearly by defoliation when treatments were in place, which excluded bud break and flowering.

The initiation of each of these pheno-phases is quite complex as it may require more than one preexisting condition. For instance, the release from dormancy is often associated to the fulfillment of a chilling/thermal time accumulation requirement , which supported the observation that all grapevines in the same site as this experiment would have a similar date of bud break. However, entering and exiting dormancy is also concomitant with major events of mobilization of soluble carbohydrates that condition the response of the latent bud . Similarly, veraison may be modeled with thermal time but ultimately requires a sucrose stimulus . Leaf senescence of deciduous plants is largely induced by shorter days and cooler temperatures, but as evidenced in our work, defoliation treatments delayed it. Other studies have suggested that leaves are able to sense source strength and delay leaf senescence accordingly . However, in our work, source strength was achieved through more leaves rather than better leaf net carbon assimilation performance. In both experimental years we witnessed the 33%L treatments assimilate more carbon compared to 100%L to compensate. However, it remains to be seen in future works if this is in fact a carbon starvation effect or an artifact of plant water status.

In fact, high sugar levels are one of the signals inducing natural leaf senescence , and this can be modulated. Interestingly, this response was not conditioned by sink strength or differences in leaf area in the final year, only by the practice of defoliation in our study.In the third year of study, no treatments were applied and therefore, all effects observed are attributable to cumulative effect of previous years’ conditions. The so-called carryover effects have been discussed in relation to indirect observations, where the treatments were applied for several years or when historical series were analyzed . In the case of defoliation, much direct evidence of carryover effects exists. For instance, Jermini et al. showed how defoliation caused by downy mildew induced severe reductions in the successive year’s yield. Bennett et al. also reported severe reductions in yield, and these were attributed to reductions in clusters per vine, to reductions in berries per cluster, but not to changes in berry mass. In that study, there were changes in inflorescence per vine and flowers per inflorescence. Therefore, promoting root and canopy growth over the years has a strong cumulative effect on yields. Alternate bearing is an issue in some tree fruit crops such as mango, avocado, olive, pistachio, citrus, etc. [reviewed by Monselise and Goldschmidt ] and fruit removal in those crops is not only performed aiming for in-season effects, but also to maintain consistent yield over the seasons. The carryover effects of crop level in grapevine have been reported less frequently than defoliation.

Our results suggested that in fact grapevine is a perennial crop not very sensitive to alternate bearing. Yield was associated with dormant season precipitation or root and shoot starch content at budburst . In our results, starch content of roots was only affected by defoliation in July and September samplings, which are coetaneous with sucrose stimulus to berry ripening. As root starch content fully recovered in all treatments, root mass was the only factor that would explain changes in yield in the successive season . In fact, in 2019, grapevines that were defoliated during the two previous seasons had lower root mass and fruit load as a carryover effect, which led to a faster recovery of starch reserves. Likewise, the carry over effects of defoliation were evident in leaf area, berries per cluster, and yield in the final year.After many efforts directed at balancing grapevine canopy by focusing on fruit removal, a renewed focus on maintaining an active leaf area with proper solar radiation exposure to clusters is needed. Our data indicated that carbon balance and translocation were more influenced by leaf area rather than crop level . The canopy leaf area and architecture determined the photosynthetic capacity which in turned initiated the sugar induced growth. Once a large enough plant reduced carbon pool was available, berry enlargement and sugar allocation determined berry size. Canopy size also dictated how fast ripening progressed as well as storage of non-structural carbohydrates. Finally, we did not measure a direct physiological benefit of fruit removal in the primary metabolism of grapevine that could suggest one unit of leaf area would equate to one unit of fruit in the grapevine.The HPA axis maintains a diurnal rhythm marked by a daily peak after waking, a subsequent decline over the course of the day, and a nadir shortly after onset of continuous sleep . The diurnal pattern of HPA activity plays important roles in a variety of metabolic, immunological, and psychological processes that support our day-to-day functioning . In studies of children, the preferred assessment method of HPA axis activity is the collection of saliva and the measurement of cortisol . Cortisol is the “end-product” hormone released into the bloodstream from the adrenal glands—the final step in a biological cascade initiated by the hypothalamus and perpetuated by the pituitary gland. In addition to supporting the orchestration of several other processes , moderate cortisol levels are thought to support effective neural transmission and optimal learning and high-order cognition . In times of acute physiological or psychological stress, the HPA axis mounts a particularly pronounced response, culminating in high levels of cortisol that reach glucocorticoid receptors throughout the body and brain. Working with the ANS, these acute HPA stress responses coordinate the physiological and psychological resources needed to overcome the stressor . Yet, given negative feedback processes, high cortisol levels also play important regulatory roles in down-regulating HPA axis activation, allowing it and other systems to return to baseline . Collectively, blueberry packaging box these complex within- and cross-system dynamics support an organism’s ability to both respond to and recover from the effects of environmental stressors . Newborn infants can mount an HPA axis response to environmental stimuli , and normative circadian rhythms tend to stabilize as infants begin to forego their afternoon naps . However, the span from infancy through early childhood is also a time of meaningful developmental change. Indeed, a growing theoretical and empirical literature indicates that children’s early experiences play a critical role in the organization of their emerging adrenocortical systems .Low-income ecologies present a confluence of distal and proximal risk factors thought to influence children’s developing physiological stress systems and undermine optimal cognitive and social development . For example, children growing up in low-income contexts are more likely to face distal stressors, such as inhospitable and dangerous neighborhoods and inadequate access to services and social capital . Such distal risks are known to have trickle-down effects that undermine parents’ abilities to effectively read, interpret, and respond to their children’s needs . In turn, a convergent literature comprising experimental work with animals as well as observational studies of young children indicates that sensitive and responsive caregiving can support adaptive HPA axis functioning .

This is evident with respect to children’s acute stress responses. For example, young children with secure attachment relationships and more sensitive caregivers tend to show better regulated HPA axis responses when faced with acute psychological stressors . Changes in the quality of children’s caregiving environments have also been linked with their baseline, or resting levels of HPA axis activity. For example, at the more extreme end, children who are moved from very high-risk households into foster care have been found to evince comparatively lower resting cortisol levels than their peers who remain in high-risk homes . It is important to note that similar relations are evident with respect to parenting within the normal range. For instance, in prior work with the same sample as used in the present study, our group showed that higher levels of maternal sensitivity in infancy are predictive lower levels of resting cortisol, after adjusting for income and a number of potential confounds . Beyond psychosocial risks, children growing in the context of economic adversity are more apt be exposed to households that are more densely populated, noisy, disorganized, and unpredictable—aspects typically discussed under the umbrella term chaos . A growing literature suggests that chaotic environments may alter children’s ANS and HPA axis functioning in early and middle childhood. Recent work by researchers in our laboratory suggests similar effects with respect to infants and toddlers , with within-child increases in chaos predictive of contemporaneous increases in resting salivary cortisol in later infancy and toddlerhood.Notably, young children growing up in low-income contexts spend substantial amounts of time in settings outside of their homes—such as non-parental child care. Indeed, in the United States approximately 43% of children in poverty attend regular non-parental care by 9 months of age . A well-developed literature indicates that young children’s early child care experiences also play a meaningful role in their HPA axis functioning. Meta-analytic findings indicate that—compared to their normal diurnal patterns experienced at home—children tend to show cortisol increases across the day on days when they attend child care . Some work suggests these patterns are particularly strong in toddlerhood and the beginning of the early childhood years and for children who attend lower quality child care . There is also some, albeit limited, evidence of long-term effects; for example, Roisman and colleagues found that spending greater proportions of time in center-based care in infancy and early childhood was predictive of children’s subsequent cortisol awakening response in adolescence. Although such effects with regard to child care type have been rather mixed, there is some reason to suspect that heightened exposure to peers—such as often found in center-based child care—may be a salient stressor for young children . For instance, at abehavioral level, the replicated link between greater exposure to center-based child care and subsequently heightened levels of aggression in childhood has been shown to be at least partially explained by the typically higher levels of peer exposure faced by these children . That said, the potential effects of peer exposure on HPA axis functioning in infancy and toddlerhood are largely unknown. Indeed, some work has suggested that heightened peer exposure may be associated with lower levels of HPA activity. For instance, in their generative study of infants and toddlers, Watamura and colleagues found that toddlers who spent more time playing with their peers tended to show lower cortisol levels than their less social classmates. Of course, the direction of this relation remains unclear; for example, perhaps less physiologically stressed children choose to play more with peers . However, it nonetheless raises the possibility that peer exposure may provide valuable opportunities for young children to learn to negotiate such complex social contexts.With rare exception, the extant literature concerning child care and children’s early adrenocortical functioning has been based on rather small, homogeneous samples of children from middle- to upper income families. Informed by findings from studies of children’s behavioral outcomes, there is increasing evidence to suggest that child care effects for children growing up in poverty may be quite different those for children from more affluent home families .

Prior work has typically concerned diurnal cortisol rhythms over the course of the day, within-child child care variation across a more limited span , and rather blunt indices of child care exposure . In contrast, we considered children’s resting cortisol levels at a single point in time, within-child variation over the course of many months, and variation in the actual number of hours children spent in child care per week. Despite these differences , the commonalities between this prior work and the present findings are noteworthy. Taken collectively, the findings suggest that within-child increases in child care exposure—on both short- and long-term developmental time scales—tend to be associated with within-child increases in salivary cortisol, for children from low-risk contexts. In particular, in the present study we found that—holding all other time-invariant confounds constant—fairly long-term shifts in child care exposure across spans of several months may underlie changes in HPA axis function, at least at 24 months of age. This finding with respect to developmental timing is consistent with meta-analytic data suggesting that the link between child care and children’s cortisol increases over the course of the child care day tends to be quite small in infancy yet pronounced toward early childhood .

The cause of these developmental differences remains largely unclear. On the one hand, plastic growing bag it is reasonable to conjecture that the types of experiences that occur in child care, such as exposure to peers, noise, and the divided attention of caregivers, might be particularly salient stressors during infancy. Developmentally, infants lack the internal self-regulatory control that is presumably needed to down-regulate their stress responses . On the other hand, children’s experiences of the same child care context may also be quite different at different points in development . For instance, whereas an infant and early-childhood classroom may maintain the same number of peers, children tend to be more actively engaged with their peers in early childhood . Such normative social differences may begin to explain the emergence of the childcare effect in toddlerhood. This is an important question for future work. Critically, our findings indicated an inversed within-child relation between child care hours and cortisol at 24 month of age, for children from high-risk contexts. For instance, for toddlers with cumulative risk levels that were 1 SD above the mean, increases in weekly child care hours were associated with contemporaneous decreases in children’s cortisol levels. Although this conditional association was somewhat modest—corresponding to a standardized regression coefficient of approximately .14—it is approximately one quarter of the size of the effect of time of day, a known biological mechanism underlying within person cortisol variation. Furthermore, the estimate is robust as a within-person effect, which adjusts for all possible time-invariant confounds.

These findings align with our prior work with respect to testing between-child interaction effect as predictive of children’s resting cortisol levels at 48 months. In that work we found similarly that the direction of the association between child care hours and children’s cortisol levels differed for children from low- versus high-risk contexts . For children experiencing low levels of cumulative risk, attending a greater number of child care hours—on average, between 7 and 35 months—was predictive of comparatively higher cortisol levels, whereas for high-risk children a negative relation was evident. Indeed, the consistency between our prior and present findings is rather remarkable, given the different substantive interpretations of between- versus within-person effects. The former represents one’s average child care exposure compared to other children’s average exposures; the latter represents time-specific increases in child care hours at 24 months, compared to one’s own prior child care hours at 7 and 15 months—irrespective of one’s average level. It is important to note that our prior interactive findings at 48 months and the present interactive findings between 7 and 24 months almost certainly tap distinct effects. As above, the interpretations differ. In addition, 24-month and 48-month cortisol levels are virtually orthogonal . The present findings also align with prior studies that have considered interactive effects between environmental risk and child care exposure in the context of children’s behavioral outcomes. For example, using the same sample as the present study, in prior work we have found that, for children experiencing high levels of household chaos across early childhood, greater weekly hours in child care were predictive comparatively fewer behavior problems. In particular, consistent with the idea of a buffering effect, greater child care exposure ameliorated the detrimental relation between household chaos and children’s social problems.

Similar interactive effect have been noted by others , and work conducted with predominantly lower income samples has increasingly shown that the “beneficial” relations between child care and children’s social and academic outcomes may be stronger and more consistent than those from studies of children from more affluent families . Despite some commonalities emerging across this literature, the mechanisms underlying these interactions remain largely unclear. Our results provide some indication that center based care may play a role. We found that—for children from low-risk contexts—attending greater proportions of one’s time in child care in center-based care was predictive of higher cortisol levels than spending more time on other types of non-parental care. In contrast, no such relation was evident for children from high-risk contexts. To some degree, this is may be inconsistent with prior work. For instance, on the basis of data from their heterogeneous sample, Roisman and colleagues found that greater exposure to center-based care in early childhood was predictive of a lower awakening cortisol levels in adolescence. Of note, though, is that the cortisol awakening response is typically considered to be a process that is distinct from typical HPA axis functioning throughout the day . There is also, however, some reason to suspect that our findings align with some behavioral findings. For instance, there is a good indication that greater exposure to center-based care may be linked to heightened levels of aggression in early childhood —perhaps particularly so for children from more affluent families . Meta-analytic findings indicate a positive correlation between cortisol and aggression levels in early childhood . Indeed, considering the potential secondary effects of cortisol on behaviors and cognition is an important next step we intend to take with these data. However, there is also reason to suspect that the translation of children’s cortisol levels into behavioral effects will not be a simple one. For example, in the present study we also found that children who were exposed to greater numbers of peers, on average, wholesale grow bags between 7 and 24 months, tended to have lower cortisol levels. This is consistent with the extant data , yet it is somewhat inconsistent with the idea that the link between center-based care and heightened levels of aggression is explained by its effects on heightened HPA axis activity: Peer problems may well underlie the effects of center-based care on aggression . Suffice it to say, there is much left to clarify. Interestingly, prior work has been suggested that the negative relation between greater peer play and lower daily cortisol levels may reflect the fact children who are less physically aroused may “select into” greater amounts of peer play . Of course, this may well be the case, yet it is noteworthy that our finding with regard to the number of children in child care—presumably reflecting peer exposure that is considerably less child driven—shows similar negative relation. Although it remains a conjecture to be tested directly, ongoing practice negotiating the social world of one’s peers may, in fact, play a “steeling” or “inoculating” role in the way that children respond physiologically to subsequent peer interactions. Evidence from nonhuman primates suggests that such steeling effects may be common and play an important role in resiliency .As such, our present findings suggest that child care type and peer exposure—at least when considered cumulatively between ages 7 and 24 months—are linked with HPA axis functioning across this period.

However, it is important to note that they do not appear to explain the effects of childcare quantity. Child care quantity, type, and peer exposure were additively and independently predictive of children’s cortisol levels, and no high-order interactions were evident. Like our prior work with children’s 48-month cortisol levels, caregiver responsivity was not related with children’s cortisol levels, between or within persons. Thus, as noted by others in prior work the mechanism underlying the child care quantity effect remains largely unclear. One possibility is that, for children facing the confluence of environmental challenges of poverty, child care may be comparatively less physiologically stressing than their experiences outside of child care. Indeed, the limited available evidence suggests that— unlike their more affluent peers—children from higher risk contexts tend to show a normative diurnal cortisol decline in their cortisol levels across the child care day . If this is the case, then the “beneficial” relation between greater child care exposure and cortisol observed for children from high-risk families may partially reflect a dose–response relation, whereby greater child care exposure leads to less exposure to the pronounced environmental stressors outside of child care. This, of course, remains a question to be addressed directly. Indeed, there are several potentially important mechanisms that were not addressed in the present study—such as the quality of children’s peer relations , teacher–child relationship quality , and degree of classroom chaos —that remain critical areas inquiry, with regard to clarifying the relation between child care exposure and HPA axis functioning.Several of our findings are robust to within-child analyses that accounted for all time invariant confound or between-child analyses that adjusted for several observed control covariates. In fact, with regard to within-child effects, the conservative nature of these estimates may underlie the fact that there were several null within-person relations . On the one hand, these specifications may strengthen the internal validity of our inferences. On the other hand, to the extent that the causal relation between these variables is contingent on cumulative exposure they would be missed by our within-person specification. They allow us to pick up cumulative, traitlike variability in these child care measures—variability that is likely quite meaningful. However, our inferences are only causally unbiased to the extent to which we have accounted for all sources of potential endogeneity . This assumption is typically dubious. As such, we make no causal claims. Nontrivially, though, these potential endogeneity problems would have to explain both selection into cumulative risk and selection into child care, as well as why the direction of the child care effects vary as a function of cumulative risk.As is common with large, comprehensive epidemiological deigns, trade-offs must often be made between breadth and depth. For instance, we were unable to collect multiple saliva samples over the course of a day or across several consecutive days. This would have strengthened the reliability of our cortisol measures and made our findings more directly informative to the extant literature concerning children’s diurnal cortisol patterns . Our saliva samples were collected modally around 10:00 a.m., when the child was at rest, after he or she had had time acclimate to the RAs in the home. There was little evidence that children were particularly stressed by the RAs’ visits. The vast majority of mothers rated their children’s behavior during the visit as being typical. Of the minority who indicated that their children were behaving differently, it was typically due to external reasons . For the few who ascribed differences as being due to the RAs’ visit, they were evenly split between those showing emotions that were comparatively more negative and comparatively more positive than usual. This is consistent with the extant laboratory-based literature that suggests that exposure to novel adults typically fails to evoke an acute stress response—even when it is intended to do so . That being said, we cannot rule out the possibility that our cortisol measures may partially tap reactivity to the RAs’ presence. Finally, by using terms like beneficial to describe, for instance, the negative relation between child care hours and cortisol, we imply that lower cortisol levels may be substantively “better” than higher cortisol levels. Other interpretations are possible. For instance, high levels of cumulative risk could make these children especially vulnerable to child care effects that cause the HPA axis to be underactive . Underactive HPA axis functioning is also a known risk factor for less optimal cognitive functioning and behavior problems . We cannot rule this interpretation out; however, we find it less plausible.

Bayesian analysis—The dataset was analyzed with Bayesian inference using the program MrBayes version 3.2.7a. Sequence data were subjected to a general time reversal model including the estimation of invariant sites and assuming a discrete gamma distribution with six rate categories . The relative fit of various models of nucleotide substitution for the ITS region, chloroplast regions, and combined data set to identify the best model was examined. The best model was selected based on the Akaike Information Criterion . The Markov Chain Monte Carlo sampling, starting from random trees and priors, was run for 1,000,000 generations and every 100th tree was sampled. Four MCMC chains comprising three heated chains and a single cold chain were used in the analyses. Majority rule consensus trees and posterior probabilities for nodes were assembled from all post-burn-in sampled trees. Phylogenetic reconstructions were estimated after a couple of independent runs to confirm that they converged on similar stationary parameter estimates. For the combined data set, plastic nursery plant pot the data for each region were merged and aligned using muscle in MEGA before running on MrBayes.

Maximum Likelihood—For Maximum Likelihood estimation, different regions used different models. For the nuclear region, ITS 4 and ITS 5, sequence alignment was performed using the muscle tool included in the MEGA 10 suite, and the Kimura-2 parameter model was applied. For the chloroplast region involving trnHpsbA primer, alignment was also performed using muscle, but the best model was the Hasegawa-Kishino model. For the combined data, alignments were performed using muscle while the best model was the Tamura 3-parameter model. The bootstrap consensus tree inferred from 1000 replicates in Maximum Likelihood, was taken to represent the evolutionary history of the taxa analyzed. The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test is shown next to the branches in the trees generated.After extraction of DNA, poor bands were obtained during PCR for many of the species studied. This could be due to the old nature of the herbarium specimens used. Several primers were tried for the chloroplast region but for many of them, only very few specimens were amplified during PCR. The best result for chloroplast region was obtained from trnH-psbA. When samples were sent for sequencing, most of the results obtained for trnH-psbA were good for one direction, and rarely did we get good results for both directions. To protect the integrity of the results, only directions with clean DNA sequences were used for alignment and analysis. This problem with trnH-psbA that is reported here was also encountered by [2]. A higher number of successful PCR reactions were obtained for the nuclear ITS region than for the chloroplast regions, unlike the result reported by [2]. Sequences obtained for the ITS region ranged from 500 to 550 bp, while those for the chloroplast region were mostly below 500 bp. Although more sequences amplified for ITS, the tree generated for both nuclear and chloroplast regions were not significantly different, so the data were combined.

Clean reads that are suitable for analysis were obtained for the two primer sets used for this study, but more samples were amplified for the nuclear region. To combine the data, only samples that we got sequences for both regions were used. This is to ensure that none of the regions will influence the topology of the tree more than the other. Also, for ease of comparison of clades obtained from the combined dataset to clades from the separate regions.Tree from combined dataset is shown below, Figure 2. The values of the posterior probabilities of the branches obtained are used as measures of branch support. Values below 0.95 are considered to have very low support. All the trees obtained from Bayesian inference and Maximum likelihood show a close relationship between Englerophytum and Synsepalum. Trees from separate analyses are available from corresponding author.Neither Synsepalum nor Englerophytum were not resolved as monophyletic and some of the clades in the polytomy corresponded to the previously recognized genera by Aubréville. Three of the species previously grouped under the genus Afrosersalisia were used in this study namely S. afzelli, S. cerasifera, and S. kassneri. As seen in Figure 2, they belong to the same clade,—clade C. For the Vincentella genus, out of the four species, data were successfully obtained for three species. These include S. muelleri, S. passargei, and S. revoluta. None of the herbaria we loaned specimens from had S. brenanii. Just like the report for Afrosersalisia, the species are in the same clade, -clade A in Figure 2. The case is slightly different for the genus Pachystella. Out of the three species in Pachystella, two species are in the same clade, clade C, separating them from S. subverticillatum, which was found in clade A. For the genus Synsepalum sensu stricto, three species were previously recognized. These include S. dulcificum—the type species for the genus, S. stipulatum, and S. subcordatum. All three species were included in this study. In this group, all the species were in one clade, clade D.

It is important to note, however, that these three species were not the only species found in the clade, other species formed the same clade with them.Two principal issues are addressed in this work. One is the merging of the small genera recognized by [7], while the second issue is the divergent view of some authors on the Synsepalum–Englerophytum complex. The previously recognized small genera Afrosersalisia, Vincentella, Pachystela, Synsepalum, Tulesta. After the revision of Sapotaceae by Pennington in 1991, 19 more species have been added to the genus Synsepalum bringing the total to 35 species currently recognized. Some authors including have called for the separation of Synsepalum from the previous small genera that were combined. As seen in the combined data tree generated, some of the clades corresponded to the previous genera that were recognized by Aubréville. Data analysis using Bayesian inference and Maximum likelihood for ITS grouped all three species of Synsepalum sensu stricto used in this study into one clade. A dichotomous tree having species of Synsepalum in one clade and species of Englerophytum in another clade would have supported the monophyly of Synsepalum, however, the tree obtained from the nuclear, chloroplast region and combined dataset resulted in a phylogenetic separation of subgroups of species belonging to both genera. In the combined tree reported here, clades of Englerophytum are nested within Synsepalum. This consistent nesting of Englerophytum within Synsepalum further suggests that the genera are closely related. Leaf venation, presence or absence of stipules, inflorescence, fusion of sepals, stamens insertion, anther position, and number of ovaries tend to overlap in the previously merged genera. The overlapping morphological character and the molecular evidence seen in tree suggest very strongly that Synsepalum sensu lato should be reversed into small genera.Molecular evidence obtained from this study and all other phylogenetic studies involving the species merged by Pennington and the newly described ones for both Synsepalum and Englerophytum shows that the two genera are not phylogenetically distinct from each other. Thus, the trees obtained in this research do not agree with the circumscription of the genera as defined by Pennington. All the trees obtained using MrBayes for Bayesian inference and MEGA X for Maximum likelihood show that Englerophytum is nested within Synsepalum. It is important to state here that only a few morphological characters were used by [7] to separate Englerophytum from Synsepalum. Synsepalum has an eucamptodromous venation and Englerophytum has a brochidodromous venation and a strong fusion of thefilament into a staminal tube. Although it is taxonomically correct, in some cases, seedling starter pot to use few morphological characters to separate genera, this is not just the case with the EnglerophytumSynsepalum complex. The nesting of Englerophytum within Synsepalum in trees generated for both nuclear chloroplast regions indicates very strongly that the group do not merit distinct generic status.In this study of the phylogenetic relationships of the Synsepalum–Englerophytum complex, neither of the two genera resolved as monophyletic.

Although resolution of relationships among the major clades were only weakly to moderately supported, the clades of Englerophytum are nested within Synsepalum as seen in the combined tree. This shows that the morphological circumscription of the two genera does not align with the molecular evidence. We think that morphological character states that distinguish Synsepalum might be ancestral while the character states that distinguish Englerophytum may be independently derived synapomorphies for each subclade. The polytomy observed in the genus Synsepalum shows that the genus is not a single lineage. Some of the clades in the phylogenetic tree correspond to some of the small genera merged by Pennington. The need to resurrect some of the previously merged genera, like Vincentella, is undeniably evident in the result obtained in this work and that of [2]. The Synsepalum–Englerophytum complex certainly requires some formal taxonomic changes to align with molecular evidence, but such changes cannot be made until a broader sampling of all species currently recognized in the genera is collected for comprehensive morphological as well as molecular studies. The complete genome of Synsepalum has been published. If the complete genome of Englerophytum is published, it would serve as additional information relevant to making definitive and comprehensive conclusions on taxonomic circumscriptions within the Englerophytum–Synsepalum complex.After zinc, boron is the second mostimportant micronutrient deficiency problem in California vineyards. Boron deficiencies are most common in the old flood plains and alluvial fans of the Stanislaus, Merced, San Joaquin, Kings and Kaweah rivers; the Sierra Nevada foothills; and North Coast sites with basaltic soils subject to high rainfall. Vineyard boron deficiencies are mostly associated with soils derived from basaltic and granitic parent material of the Sierra Nevada and North Coast ranges. Low boron is also associated with higher rainfall areas and soils irrigated with snowmelt water originating from the Sierra Nevada. In contrast, boron levels are typically higher and can even be toxic in soils originating from marine sedimentary material, such as in the Central Coast range. Grapevine reproductive tissues are most sensitive to boron deficiency, which results in reduced fruit-set, small “shot berries” that are round to pumpkin-shaped, and flower and fruit cluster necrosis. Boron deficiency can have a drastic effect on fruit quality and yield, even when there are only mild-to-moderate foliar symptoms. At the same time, the over-application of boron can result in plant phytotoxicity. Phytotoxicity begins as a necrosis of the leaf margins that can cause a downward cupping of the young leaves. Most commonly, grape growers have applied boron to the soil by hand or as a direct soil spray, sometimes in combination with an herbicide application . However, such applications must be carefully timed to allow for winter rainfall or irrigation to move the boron into the root zone. Boron applications by foliar spray and drip irrigation are of increased interest for their convenience and the potential for faster vine response. Foliar boron application has been studied in tree crops such as pears, prunes, cherries and almonds, and the application timing was found to A boron-deficient Thompson Seedless cluster in the trial vineyard shows reduced fruit-set, the presence of numerous pumpkin-shaped “shot berries” and necrosis of some branching. Fewer than 10% of the berries are normal size and shape. influence fruit-set and development . However, there is limited research on vine uptake and response to foliar boron and the potential for toxicity. We conducted several studies on the timing of boron foliar-spray applications in an eastern Fresno County vineyard with mild-to-moderate boron deficiency symptoms.The studies were conducted in 1998 and 1999 in an own-rooted, furrow irrigated ‘Thompson Seedless’ grape vineyard on Delhi loamy sand. The vineyard was irrigated with canal water and about 25% supplementation from well water. The extremely low boron content of canal water can contribute to low boron availability in sandy vineyard soils. A preliminary study was conducted in 1998 to determine the influence of foliar sprays on boron concentrations in vegetative and reproductive parts of the vine at bloom. Spray treatment was applied at 2 1/2 weeks pre-bloom on May 6, 1998. The two comparative treat-ments were an untreated control and a foliar boron spray at 1 pound per acre applied as Solubor at 100 gallons per acre . The trial design was a randomized complete block with four-vine plots, replicated 10 times. Vine tissue samples were taken at bloom on May 23, 1998; triple-rinsed with distilled water and oven-dried; and analyzed for boron at the ANR Analytical Laboratory at UC Davis.

The consequence of such a relation in 3D is discussed. The organization of this chapter is as follows. We introduce the cubic lattice Dirac model and a transformation on the low energy Dirac fermions to bring it into a ‘normal’ form in Sec. 3.2.2, from which one can easily read off the orders that lead to an energy gap. The chiral topological insulator is identified, and a microscopic model with hopping along the body diagonals of the cube, is shown to lead to this phase. An intuitive picture in terms of a quasi 2D starting point is developed in Sec. 3.3.2, which directly demonstrates the existence of surface Dirac states. Here, we also show how a Z2 topological insulator protected by TRS with bulk Dirac nodes at TRIM can be understood within this picture. In Sec. 3.4 the magnetoelectric coefficient ‘θ’ of such an insulator is argued to be quantized , and is calculated to be θ = π for the spinless fermion model we discuss. Introducing spin, and studying gapped superconducting states, we show in Sec. 3.6 that only a pair of singlet superconductors is allowed, which, in addition to the regular onsite s-wave paired state, bucket flower includes a singlet topological superconductor, with pairing along the body diagonals. Sec. 3.7 describes an attempt to move towards a more physical realization of the cTI phase, utilizing a layered honeycomb lattice structure.

Finally, in Sec. 3.8, we explore some topological properties of the 3D Dirac fermion system by studying its physics in the presence of point topological defects in its order parameters, and deriving the Berry’s phase terms that determine quantum interference of different orders.In many ways, a chiral topological insulator can be viewed as a close cousin of the known topological states in 3D, such as a Z2 topological insulator: A hallmark of both of these states is an appearance of non-trivial surface modes when topological bulk states are terminated by a boundary. However, for a chiral topological insulator, an arbitrary number of flavors of Dirac fermions can appear at the surface and be stable. We discussed a physically transparent picture of the chiral topological insulator, which explains the appearance of surface Dirac fermion states, and their stability in the presence of chiral symmetry. A similar picture also explains the stability of TRS Z2 topological insulators whose bulk Dirac nodes are centred at TRIM. It is shown that the θ = π axion electrodynamics can also be realized in chiral topological insulators, in addition to the known realization in Z2 topological insulators. We should also stress that chiral symmetry, which is realized in the models discussed here as sublattice symmetry, is likely broken in any realistic systems. Nevertheless, as far as breaking of chiral symmetry is sufficiently weak, θ is expected to be close to π, and can It is also worth while mentioning that chiral symmetry need not to be realized only as sublattice symmetry: class AIII symmetry can be realized in the BdG Hamiltonians for Szconserving superconductors. Thus, chiral symmetry when realized in this way is much more robust than sublattice symmetry.

Furthermore, in the 3D π-flux lattice model with inclusion of spin degree of freedom, we found a spin-chiral topological insulator , and also a singlet topological superconductor . The latter is stable as long as the physical symmetries of SU spin rotation and time reversal are present. Finally, utilizing the proximity to a Dirac state, we derived an interesting correlation, or “duality”, between the singlet topological superconductor and Neel order. These order parameters are dual in the sense that a topological defect in either one of these phases carry complementary quantum numbers: e.g.. a defect in the Neel vector can carry electric charge. We also find many such 6-tuplets of order parameters, including a six component vector consists of three Neel order and three VBS order parameters. These dualities are a natural extension of those discussed in 1D and 2D quantum spin models, the latter in the context of deconfined criticality. While in this chapter we have studied the properties of these topological defects at single particle level, and hence the topological defects are static objects, we cannot resist contemplating more interesting situation where they are dynamical entities. In particular, it is interesting to ask if there is a counterpart of the non Landau-Ginzburg transition, realized in two dimensions, can exist, possibly in the presence of strong electron correlations in three dimensions. This is left for future study. P.H. would like to thank Shinsei Ryu and Ashvin Vishwanath for the collaboration in this work. We wish to thank Ying Ran and Tarun Grover for insightful discussions, and the Center for Condensed Matter Theory at University of California, Berkeley and NSFDMR-0645691 for support. This work is published in Physical Review B, 81, 045120 .In the last chapter, we saw interesting physics emerge when appropriate sets of phases were combined together to form three component order parameters, and topological defects were created in them. A natural question to ask is whether there are other topologically non-trivial ways of combining the phases which give rise to unconventional phenomena. The strong topological insulator is a good candidate for being involved in some novel physics, since rich phenomena have already been predicted to arise when this phase is combined with conventional orders such as magnetism, crystalline order and superconductivity.

The last is particularly interesting. Superconductivity induced on the surface of a TI was predicted to have vortices harboring Majorana zero modes. These are of interest to quantum information processing, since they are intrinsically robust against errors. Recently, superconductivity was discovered in a doped TI, which could be used to induce surface superconductivity. Below we discuss a new theoretical approach to studying this remarkable superconducting phase, which provides different insights and directions for experiments. In this case, the low energy description of the system is entirely in terms of bosonic coordinates , much as the Landau Ginzburg order parameter theory describes superconductors at energies below the gap. Can fermions ever emerge is such a theory? While it is easy to imagine obtaining bosons from a fermionic theory, the reverse is harder to imagine. However, it has been shown in principle that bosonic theories that contain additional Berry’s phase terms, can accomplish this transmutation of statistics. We show that this indeed occurs in the superconductor-TI system; the order parameter theory contains a Berry phase term which implies that a particular configuration of fields – the Hopf soliton – carries fermionic statistics. While such statistics transmutation is common in one dimension, it is a rare phenomenon in higher dimensions. In the condensed matter context, an physically realizable example exists in two dimensions: solitons of quantum Hall ferromagnets are fermionic and charged, and have been observed. However, the superconductor-TI system is, to our knowledge, the first explicit condensed matter realization of this phenomenon in three dimensions. The organization of this chapter is as follows. First, we introduce our simplified model of a topological insulator with surface superconductivity, cut flower bucket and review properties of the Hopftexture. We then discuss evidence from numerical calculations on a lattice model, that demonstrate the Hopfions are fermions. We also discuss the connection between fermionic Hopfions and 3D non-abelian statistics of [128]. A simplified two dimensional example, where skyrmions are fermionic, is also discussed. Next, we provide a field theoretical derivation of the same result in the continuum, and introduce the necessary theoretical tools to compute a topological term, that leads to the fermionic statistics. Finally, we mention physical consequences, for tunneling experiments as well as for Josephson junctions.In the previous chapters, the focus has largely been on bulk topological phases and especially on the bulk of the topological insulator. Although the Majorana zero mode discussed in the previous chapter was localized on the surface, its existence was governed by the bulk topological phase. In contrast, in this chapter we focus only on the surface states of the strong topological insulator . In their simplest incarnation, the SSs of TIs correspond to the dispersion of a single Dirac particle, which cannot be realized in a purely two dimensional band structure with time reversal invariance. This dispersion is endowed with the property of spin-momentum locking, i.e., for each momentum there is a unique spin direction of the electron.

Since several materials were theoretically predicted to be in this phase, most of the experimental focus on TIs so far has been towards trying to directly observe these exotic SSs in real or momentum space, in tunneling and photo emission experiments, respectively, and establish their special topological nature. However, there has so far been a dearth of experiments which study the response of these materials to external perturbations, such as an external electromagnetic field. In order to fill this gap, we calculate here the response of TI surfaces to circularly polarized light. Since photons in CP light have a well-defined angular momentum, CP light can couple to the spin of the surface electrons. Then, because of the spin-momentum-locking feature of the SSs, this coupling can result in dc transport which is sensitive to the helicity of the incident light. This phenomenon is known as the circular photogalvanic effect . In this work, we derive general expressions for the direct current on a TI surface as a result of the CPGE at normal incidence within a two-band model and estimate its size for the surface of Bi2Se3, an established TI, and find it to be well within measurable limits. Since bulk Bi2Se3 has inversion symmetry and the CPGE, which is a second-order non-linear effect, is forbidden for inversion symmetric systems, this current can only come from the surface.We find, remarkably, that the dominant contribution to the current is controlled by the Berry curvature of the electron bands and grows linearly with time. In practice this growth is cut-off by a scattering event which resets the current to zero. At the microscopic level, this part of the current involves the absorption of a photon to promote an electron from the valence to the conduction band. The total current contains two other terms – both time-independent – one again involving an interband transition and the other resulting from intraband dynamics of electrons. However, for clean samples at low temperatures, the scattering or relaxation time is expected to be large, and these contributions will be eclipsed by the linear-in-time one. Hence, this experiment can also be used to measure the relaxation time for TI SSs. Historically, the Berry curvature has been associated with fascinating phenomena such as the anomalous Hall effect and the integer quantum Hall effect and therefore, it is exciting that it appears in the response here. Its main implication here is that is gives us a simple rule, in addition to the requirement of the right symmetries, for identifying the perturbations that can give a linear-in-time CPGE at normal incidence: we look for perturbations that result in a non-zero Berry curvature. Put another way, we can identify perturbations that have the right symmetries but still do not give this current because the Berry curvature vanishes for these perturbations. Importantly, for TI SSs, the requirement of a non-zero Berry curvature amounts to the simple physical condition that the spin-direction of the electrons have all three components non-zero. In other words, if the electron spin in the SSs is completely in-plane, the Berry curvature is zero and no linear-in-time CPGE is expected. The spins must somehow be tipped slightly out of the plane, as shown in Figure 6.1.1a, in order to get such a response. Thus, a pure Dirac dispersion, for which the spins are planar, cannot give this response; deviations from linearity, such as the hexagonal warping on the surface of Bi2Te3 , are essential for tilting the spins out of the plane. CPGE has been observed in the past in GaAs, SiGe and HgTe/CdHgTe quantum wells – all systems with strong spin-orbit coupling. The effect in these systems can be understood within a four-band model consisting of two spin-orbit split valence bands and two spin-degenerate conduction bands. In contrast, TI SSs can be faithfully treated within a two-band model. The simplicity of the latter system makes it more convenient for studying theoretically compared to semiconductor quantum wells, and hence, enables us to determine a connection between the CPGE and the Berry curvature.

Older adults with greatest engagement of rostral anterior cingulate, a prime regulator of amygdala activity , self-reported less unpleasantness for lowarousal negative stimuli. Together, these findings provide insight into the neural basis of agerelated changes in affective experience by which older adults preferentially enhance positive information, preferentially suppress negative information, or both. Given the recent reinvigoration of research establishing roles of the amygdala in reinforcement learning and decision-making , future studies should pay specific attention to the amygdala as a possible hub where emotional affective processes and dopaminergic reward-based processes intersect. The positivity effect may also influence decision-making through its modulation of medial temporal lobe memory systems. If positive information has privileged access to memory in aging , future choice behavior may be preferentially biased by positive experiences rather than negative experiences . There is some evidence to suggest that in the context of reward-based learning paradigms, older adults show intact episodic memory performance for the visual cues that predicted reward .

Intact subsequent memory for positive cues is particularly notable because it occurred despite reduced levels of reward-based probabilistic learning in the same subjects . How the strength of episodic memory for reward-related information is related to individual differences in the positivity effect is an open question, procona London container but represents one mechanism by which age-related changes in affective attention may shape decision-making. There are important boundary conditions that affect the occurrence and presentation of the positivity effect that are relevant for understanding how it may impact decision-making in aging. The positivity effect in aging is diminished in scenarios in which older adults are given explicit instructions on how to behave, and which stimuli to attend to . Therefore, it is more likely to impact performance on tasks that permit multiple strategies. Supporting this view, in dynamic learning environments that support exploratory behavior, older adults preferentially bias decision-making on recent positive feedback . The positivity effect is also sensitive to cognitive load. Specifically, enhancement of attention and memory for positively valenced information is absent or reversed relative to young adults when cognitive load is increased .

Therefore, the positivity effect is more likely to impact decision-making performance on tasks that support multiple strategies and are low in their computational requirements. To illustrate how the age-related positivity effect may influence decision-making performance, we use PST as an example. The PST, which does not have explicit instructions guiding choice strategies and is relatively low on working memory demands, is potentially sensitive to the age-related positivity effect. For young adults in which the positivity effect is largely absent, performance may be primarily predicted from measures of dopamine function. For example, D1 BPND would be positively correlated with choice selection to approach reward . For older adults in which the positivity effect is present, performance may be shifted from predictions based solely on measures of dopamine function . Reduced D1 BPND in aging would predict a general reduction in choice selection to approach reward relative to young participants. However, socioemotional selectivity accounts of PST performance in aging would predict greater choice selection driven by motivation to approach reward, especially if memory for the reward-predicting cue is particularly salient . Here, we suggest that the positivity effect can shift choice bias in directions that are opposite to what would be predicted based on dopamine effects alone. Such shifts may rely on the strength of the positivity effect within individuals, but nonetheless likely reflects a factor contributing to mixed and null effects in studies examining decision-making in aging.

Lines of future research should empirically test the possibility that dopaminergic decline and enhancement of positive affective processing counteract one another in laboratory-based decision-making tasks. This could be achieved by generating a normative measure of memory for positive, negative and neutral images in a large group of young adults from which to measure the magnitude of biases in individual older adults. Lang, Bradley, and Cuthbert have developed a catalogue of images with affective normative ratings, which could support these efforts . Assessment of the degree to which affective bias measures are related to variability in reward-based and punishment based decision-making would be an important contribution. Future studies could also test the extent to which individual differences in amygdala reactivity and hippocampal engagement during learning positive versus negative associations are predicted by individual differences in the positivity effect.There is an increasing number of innovative new studies examining the neural mechanisms underlying age-related changes in decision-making. The integration of in vivo imaging measures such as PET and neuromelanin-sensitive MR with computational models is a promising avenue for resolving the basic neurochemical drivers of individual differences in performance. This approach will be essential for disentangling age effects of dopamine from age effects of affective attentional bias that contribute to biases in decision-making. This approach will also be useful for defining the nature of changes in reward processing in psychological and psychiatric disorders associated with altered dopamine function and dys regulated affective processing such as schizophrenia, addiction, and depression. Critical steps for future research will be to determine the degree to which findings of laboratory-based assessments of decision-making and financial risk taking can be extended to real-world scenarios. To date, the relatively few cognitive neuroscience studies that have tested these relationships in aging have shown general agreement between performance on financial reward tasks and real-world financial measures . Future studies may focus on predicting the conditions under which older adults are vulnerable to suboptimal decision making. If changes in affective focus in aging partly counteract the detrimental effects of dopamine loss, older adults may be differentially impacted in decisions that occur in conditions in which the positivity effect is suppressed. For example, the positivity effect may be diminished for decisions made in high cognitive load conditions, where there is excessive or irrelevant information presented.Grapes are profitable fruit crop that are widely grown in the state of California, with an increasing need to accomplish cultural tasks mechanically , 2020; Kurtural and Fidelibus, 2021. However, there are many factors that are currently challenging the productivity, quality, and sustainability in wine grape vineyards, one being the increasingly significant global warming trend affecting California and the whole world , where more frequent heat waves and continued warming of air temperature imposes great threats to vineyard yield, berry and wine composition . Grape berry and wine quality are determined by the composition and concentration of secondary metabolites accumulated in berries. Flavonoids are the most abundant secondary metabolites and contribute to many quality determining traits, including color, mouthfeel, cut flower transport bucket and aging potential of wine . There are generally three classes of flavonoids in wine grapes, including anthocyanins, flavonols, and proanthocyanidins. Anthocyanins are responsible for grape berry and wine color, and they are sensitive to external environmental conditions when clusters are exposed to solar radiation and heat, with overexposure resulting in anthocyanin degradation . On the other hand, flavonols tend to be positively related to solar radiation . Solar radiation, especially UV-B, can often up-regulate flavonols’ biosynthesis, resulting in more flavonols accumulated in berry skins. However, excessive solar radiation received and heat accumulated in California would accelerate the degradation of not only anthocyanins, but also flavonols, which will cause a decline in the antioxidant capacity of resultant wine and a possible reduction in wine aging potential . In viticulture, trellis system selection is a critical aspect grower needs to consider when establishing a vineyard. An ideal trellis can promote grapevines’ photosynthetic capacity through optimizing light interception by the grapevine canopy.

Most importantly, a suitable trellis can optimize canopy microclimate by providing sufficient solar penetration into canopies since solar radiation is necessary to enhance the berry composition without excessive exposure of clusters to direct sunlight to avoid flavonoid degradation . There is evidence that grape clusters over-exposed to solar radiation are prone to occur with some of the widely used trellis systems. For example, vertical shoot position , a traditional and commonly used trellis system in viticulture production, has been found to produce canopies with high porosity which increases vulnerability of clusters to over-exposure , causing overly enhanced maturity and considerable degradation in berry anthocyanins . However, there is a lack of evaluations of the performance among various trellis systems in relation to the warming climate trends, and how their specific architectures contribute to variations in berry chemical profiles. In warm climates such as California, viticulture relies on irrigation for maintaining production, and previous work in the area showed that the application of different amounts of crop evapotranspiration can significantly modify polyphenolic and aromatic profiles in wine . Due to the increasingly frequent drought condition in many wine grape growing regions, recent studies have been focusing on the grapevine physiological and berry chemical responses towards specific levels of water deficits imposed by different ETc replacements, where water deficits are affective in manipulating grapevine water status, leaf gas exchange, components of yield, and berry composition: often, more water deficits applied to the grapevines would diminish photosynthetic capacities, but promote berry maturity . In some extremely drought conditions, however, severe water deficit might lower flavonoid concentration due to encouraged chemical degradation . Moreover, these effects resulted from different irrigation regimes can be modified by the canopy architecture as functions of trellis system since trellis systems can directly determine canopy sizes, hence resulting in different water demands from grapevines accordingly . On the other hand, over extraction of ground water to irrigate permanent crops have recently been questioned and legislation has been enacted in the state of California called the ‘Sustainable Groundwater Management Act’ . As a result, in some regions such as Napa Valley of California, grape growers will only be allowed to irrigate 120 mm per year. However, there is a lack of information on how the existing vineyards will cope with this water limitation in terms of irrigation scheduling. Therefore, the objectives of this study were to evaluate and compare 6 different trellis systems in combination with 3 irrigation strategies to understand the impact of trellis system and applied water amount on canopy architecture, grapevine physiology and berry composition. We hypothesized that traditional VSP systems would not be as efficient as the other trellis systems in terms of yield production and flavonoid accumulation, leading to greater berry flavonoid degradation and overall lower flavonoid concentrations.The second subset of 20 berries was used for the determination of skin flavonoids from each individual treatment-replicate. Skins were manually removed from the subset of 20 berries and subsequently lyophilized . After lyophilization, dry skin weights were recorded and then, the dried skins were ground into finepowder with a mixing mill . 50 mg of the freeze-dried berry skin powder were collected, and the skin flavonoids were extracted with 1 mL of methanol:water:7M hydrochloric acid to simultaneously determine flavonol and anthocyanin concentration and profile as previously described by MartınezLüscher et al. . The extracts were stored overnight in a refrigerator at 4°C. In the next day, the extracts were centrifuged at 30,000 g for 15 minutes, and the supernatants were separated from the solids and transferred into HPLC vials after being filtered by PTFE membrane filters . Then, the samples were injected into HPLC for chromatographic analysis.Anthocyanin and flavonol concentrations in berry skin tissues were analyzed with a reversed-phase HPLC with the use of two mobile phases: 5.5% formic acid in water and 5.5% formic acid in acetonitrile. The specific method used for this study required a C18 reversed-phase HPLC column for the analysis . The flow rate of the mobile phase was 0.5 mL min- 1 and the flow gradient started with 91.5% A with 8.5% B, 87% A with 13% B at 25 min, 82% A with 18% B at 35 min, 62% A with 38% B at 70 min, 50% A with 50% B at 70.01 min, 30% A with 70% B at 75 min, 91.5% A with 8.5% B from 75.01 min to 90 min. The column temperature was maintained at 25°C on both left and right sides of the column. All chromatographic solvents were of high-performance liquid chromatography grade, including acetonitrile, methanol, hydrochloric acid, formic acid. These solvents were purchased from Thermo-Fisher Scientific . Detection of flavonols and anthocyanins was recorded by the diode array detector at 365 and 520 nm, respectively. Investigated anthocyanin derivatives included di-hydroxylated forms: cyanidin and peonidin, and trihydroxylated forms: delphinidin, petunidin, and malvidin; investigated flavonols included a mono-hydroxylated form: kaempferol, di-hydroxylated forms: quercetin and isorhamnetin, and tri-hydroxylated forms: myricetin, laricitin, and syrigintin.