Abstract Harpagoside, a high-value anti-inflammatory iridoid compound, is traditionally extracted from the roots of Harpagophytum procumbens (Pedaliaceae, Lamiales), a Southern African desert plant widely used in traditional medicine but currently threatened by overexploitation. Scrophularia nodosa (Scrophulariaceae, Lamiales) is a perennial annual plant widely distributed in Western Europe and accumulates several iridoid compounds with known biological activities, such as catalpol, aucubin, harpagide and particularly harpagoside. We gathered extensive genomic and transcriptomics resources for this species and aimed at deciphering the biosynthetic pathway leading to the most abundant iridoid in S. nodosa , harpagoside. We found that the early iridoid pathway is well conserved with other iridoid-producing plants and validated the enzyme activities by transient co-expression in N. benthamiana . Investigation into the large BAHD family showed subclade 6i expanding in Scrophulariaceae, with an atypical VYPWG motif instead of the canonical DFGWG. In this branch, we discovered and characterized harpagoside synthase, a BAHD-type cinnamoyl transferase enzyme showing unique high specificity to the uncommon cinnamoyl-CoA acyl donor and catalyzing the final step of harpagoside biosynthesis. These results establish S. nodosa as a new model to investigate unexplored branches of the iridoids metabolism, and are a first step towards sustainable harpagoside and high-value cinnamoyl-containing conjugates production.

BACKGROUND:Plant microbiota has received increasing attention in recent years. In particular, the microbiota associated with cereals is being extensively studied to identify bacterial strains that can promote plant health and growth. Barley is the fourth most important cereal worldwide in terms of agricultural production. Intensive barley agriculture requires the use of chemical fertilizers to compensate for nutrient deficiencies in soils and limit pathogen development. The isolation and use of bacteria that can enhance the bioavailability of soil nutrients and inhibit the development of plant pathogens could ultimately limit the use of these chemicals. In this study, we have isolated from a barley microbiota three bacterial strains belonging to the genus Streptomyces. These strains were characterized and named GPA1, GPAT2, and GPN2. RESULTS:These three closely related isolates were from the same bacterial genus Streptomyces. Based on a phylogenetic analysis, the strains GPAT2 and GPN2 were classified as Streptomyces murinus, while GPA1 was identified as a new species. All strains showed antagonistic activity against two microorganisms that inhibit barley germination: Pseudomonas sp. MRN1 and Fusarium sp. CK. In addition, these strains exhibited different effects on the growth of barley cultivated under hydroponic and axenic conditions. In fact, GPN2 appeared to have no effect whereas the inoculation of barley seedlings with GPAT2 and GPA1 resulted in a reduction and an increase in root length after two weeks of growth, respectively. GPA1 had various Plant Growth-Promoting (PGP) abilities, including phosphate and zinc solubilization and siderophore production. A metabolite profiling of the GPA1 bacterial culture also showed its production and excretion of indole-3-acetic acid (IAA). CONCLUSION:In this study, we have characterized three closely related bacteria, which display different effects on barley seedlings growth. These results revealed that the type of interactions of Streptomyces with barley is strain-dependent, suggesting that these interactions may arise from specific molecular mechanisms acquired through coevolutionary processes.

Streptomycetaceae are found ubiquitously within plant microbiota. Several species belonging to this family are plant growth-promoting bacteria or may inhibit phytopathogens. Such bacteria therefore exert crucial functions in host development and resistance to stresses. Recent studies have shown that plants select beneficial bacteria into their microbiota. However, the selection process and the molecular mechanisms by which selected bacteria modulate the physiology of their host are not yet fully understood. Previous work revealed that the metabolic status of Arabidopsis thaliana was crucial for the selection of Streptomycetaceae into the microbiota, in particular bacteria phylogenetically related to Streptomyces cocklensis or Actinacidiphila bryophytorum (previously named Streptomyces bryophytorum). Here, the Arabidopsis-Streptomycetaceae interaction was further depicted by inoculating axenic A. thaliana with S. cocklensis DSM 42063 or A. bryophytorum DSM 42138. We showed that these two bacteria colonize A. thaliana ecotype Columbia-0 plants, but the colonization efficiency is reduced in a chs5 mutant of the same ecotype, being altered in isoprenoid, phenylpropanoid, and lipid profiles. We observed that A. bryophytorum inhibits growth of the chs5 mutant but not of the wild type, suggesting that the Arabidopsis-Actinacidiphila interaction depends on the metabolic status of the host. Using a mass spectrometry-based proteomic approach, we showed that S. cocklensis and A. bryophytorum modulate the A. thaliana proteome, in particular components involved in photosynthesis or phytohormone homeostasis. This study unveils specific aspects of the Arabidopsis-Streptomycetaceae interaction and highlights its complexity and diversity. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license .

Plants have evolved sophisticated DNA repair mechanisms to cope with the deleterious effects of ultraviolet (UV)-induced DNA damage. Indeed, DNA repair pathways cooperate with epigenetic-related processes to efficiently maintain genome integrity. However, it remains to be deciphered how photodamages are recognized within different chromatin landscapes, especially in compacted genomic regions such as constitutive heterochromatin. Here we combined cytogenetics and epigenomics to identify that UV-C irradiation induces modulation of the main epigenetic mark found in constitutive heterochromatin, H3K9me2. We demonstrated that the histone demethylase, Jumonji27 (JMJ27), contributes to the UV-induced reduction of H3K9me2 content at chromocentres. In addition, we identified that JMJ27 forms a complex with the photodamage recognition factor, DNA Damage Binding protein 2 (DDB2), and that the fine-tuning of H3K9me2 contents orchestrates DDB2 dynamics on chromatin in response to UV-C exposure. Hence, this study uncovers the unexpected existence of an interplay between photodamage repair and H3K9me2 homeostasis. UV exposure modulates H3K9me2 contents at chromocentres via the histone demethylase JMJ27 which interacts with the DNA damage recognition factor, DDB2, and facilitates photodamage repair. The loading or release of DDB2 at nucleosome sites is optimized by the fine-tuning of H3K9me2.

Over the course of evolution, land plant mitochondrial genomes have lost many transfer RNA (tRNA) genes and the import of nucleus-encoded tRNAs is essential for mitochondrial protein synthesis. By contrast, plastidial genomes of photosynthetic land plants generally possess a complete set of tRNA genes and the existence of plastidial tRNA import remains a long-standing question. The early vascular plants of the Selaginella genus show an extensive loss of plastidial tRNA genes while retaining photosynthetic capacity, and represent an ideal model for answering this question. Using purification, northern blot hybridization, and high-throughput tRNA sequencing, a global analysis of total and plastidial tRNA populations was undertaken in Selaginella kraussiana. We confirmed the expression of all plastidial tRNA genes and, conversely, observed that nucleus-encoded tRNAs corresponding to these plastidial tRNAs were generally excluded from the chloroplasts. We then demonstrated a selective and differential plastidial import of around forty nucleus-encoded tRNA species, likely compensating for the insufficient coding capacity of plastidial-encoded tRNAs. In-depth analysis revealed differential import of tRNA isodecoders, leading to the identification of specific situations. This includes the expression and import of nucleus-encoded tRNAs expressed from plastidial or bacterial-like genes inserted into the nuclear genome. Overall, our results confirm the existence of molecular processes that enable tRNAs to be selectively imported not only into mitochondria, as previously described, but also into chloroplasts, when necessary.

Abstract Viral RNAs can be uridylated in eukaryotic hosts. However, our knowledge of uridylation patterns and roles remains rudimentary for phytoviruses. Here, we report global 3′ terminal RNA uridylation profiles for representatives of the main families of positive single-stranded RNA phytoviruses. We detected uridylation in all 47 viral RNAs investigated here, revealing its prevalence. Yet, uridylation levels of viral RNAs varied from 0.2% to 90%. Unexpectedly, most poly(A) tails of grapevine fanleaf virus (GFLV) RNAs, including encapsidated tails, were strictly monouridylated, which corresponds to an unidentified type of viral genomic RNA extremity. This monouridylation appears beneficial for GFLV because it became dominant when plants were infected with nonuridylated GFLV transcripts. We found that GFLV RNA monouridylation is independent of the known terminal uridylyltransferases (TUTases) HEN1 SUPPRESSOR 1 (HESO1) and UTP:RNA URIDYLYLTRANSFERASE 1 (URT1) in Arabidopsis (Arabidopsis thaliana). By contrast, both TUTases can uridylate other viral RNAs like turnip crinkle virus (TCV) and turnip mosaic virus (TuMV) RNAs. Interestingly, TCV and TuMV degradation intermediates were differentially uridylated by HESO1 and URT1. Although the lack of both TUTases did not prevent viral infection, we detected degradation intermediates of TCV RNA at higher levels in an Arabidopsis heso1 urt1 mutant, suggesting that uridylation participates in clearing viral RNA. Collectively, our work unveils an extreme diversity of uridylation patterns across phytoviruses and constitutes a valuable resource to further decipher pro- and antiviral roles of uridylation.

How specific interactions between plant and pathogenic, commensal, or mutualistic microorganisms are mediated and how bacteria are selected by a plant are important questions to address. Here, an Arabidopsis thaliana mutant called chs5 partially deficient in the biogenesis of isoprenoid precursors was shown to extend its metabolic remodeling to phenylpropanoids and lipids in addition to carotenoids, chlorophylls, and terpenoids. Such a metabolic profile was concomitant to increased colonization of the phyllosphere by the pathogenic strain Pseudomonas syringae pv. tomato DC3000. A thorough microbiome analysis by 16S sequencing revealed that Streptomyces had a reduced colonization potential in chs5. This study revealed that the bacteria–Arabidopsis interaction implies molecular processes impaired in the chs5 mutant. Interestingly, our results revealed that the metabolic status of A. thaliana was crucial for the specific recruitment of Streptomyces into the microbiota. More generally, this study highlights specific as well as complex molecular interactions that shape the plant microbiota.

Uridylation is a widespread modification destabilizing eukaryotic mRNAs. Yet, molecular mechanisms underlying TUTase-mediated mRNA degradation remain mostly unresolved. Here, we report that the Arabidopsis TUTase URT1 participates in a molecular network connecting several translational repressors/decapping activators. URT1 directly interacts with DECAPPING 5 (DCP5), the Arabidopsis ortholog of human LSM14 and yeast Scd6, and this interaction connects URT1 to additional decay factors like DDX6/Dhh1-like RNA helicases. Nanopore direct RNA sequencing reveals a global role of URT1 in shaping poly(A) tail length, notably by preventing the accumulation of excessively deadenylated mRNAs. Based on in vitro and in planta data, we propose a model that explains how URT1 could reduce the accumulation of oligo(A)-tailed mRNAs both by favoring their degradation and because 3’ terminal uridines intrinsically hinder deadenylation. Importantly, preventing the accumulation of excessively deadenylated mRNAs avoids the biogenesis of illegitimate siRNAs that silence endogenous mRNAs and perturb Arabidopsis growth and development.

The involvement of the different Lactuca species in the domestication and diversification of cultivated lettuce is not totally understood. Lactuca serriola is considered as the direct ancestor and the closest relative to Lactuca sativa, while the other wild species that can be crossed with L. sativa, Lactuca virosa, and Lactuca saligna, would have just contributed to the latter diversification of cultivated typologies. To contribute to the study of Lactuca evolution, we assembled the mtDNA genomes of nine Lactuca spp. accessions, among them three from L. virosa, whose mtDNA had not been studied so far. Our results unveiled little to no intraspecies variation among Lactuca species, with the exception of L. serriola where the accessions we sequenced diverge significantly from the mtDNA of a L. serriola accession already reported. Furthermore, we found a remarkable phylogenetic closeness between the mtDNA of L. sativa and the mtDNA of L. virosa, contrasting to the L. serriola origin of the nuclear and plastidial genomes. These results suggest that a cross between L. virosa and the ancestor of cultivated lettuce is at the origin of the actual mitochondrial genome of L. sativa.

Uridylation is a widespread modification destabilizing eukaryotic mRNAs. Yet, molecular mechanisms underlying TUTase-mediated mRNA degradation remain mostly unresolved. Here, we report that the Arabidopsis TUTase URT1 participates in a molecular network connecting several translational repressors/decapping activators including DECAPPING 5 (DCP5), the Arabidopsis ortholog of human LSM14 and yeast Scd6. A conserved Helical Leucine-rich Motif (HLM) within an intrinsically disordered region of URT1 binds to the LSm domain of DCP5. This interaction connects URT1 to additional decay factors like DDX6/Dhh1-like RNA helicases. The combination of in planta and in vitro analyses supports a model that explains how URT1 reduces the accumulation of oligo(A)-tailed mRNAs: first, by connecting decapping factors and second, because 3’ terminal uridines can intrinsically hinder deadenylation. Importantly, preventing the accumulation of excessively deadenylated mRNAs in Arabidopsis avoids the biogenesis of illegitimate siRNAs that silence endogenous mRNAs and perturb plant growth and development.

Rhodococcus strain IGTS8 is the most extensively studied model bacterium for biodesulfurization of fossil fuels via the non-destructive sulfur-specific 4S pathway. This strain was initially assigned to Rhodococcus rhodochrous and later to Rhodococcus erythropolis thus making its taxonomic status debatable and reflecting the limited resolution of methods available at the time. In this study, phylogenomic analyses of the whole genome sequences of strain IGTS8 and closely related rhodococci showed that R. erythropolis and Rhodococcus qingshengii are very closely related species, that Rhodococcus strain IGTS8 is a R. qingshengii strain and that several strains identified as R. erythropolis should be re-classified as R. qingshengii. The genomes of strains assigned to these species contain potentially novel biosynthetic gene clusters showing that members of these taxa should be given greater importance in the search for new antimicrobials and other industrially important biomolecules. The plasmid-borne dsz operon encoding fossil fuel desulfurization enzymes was present in R. qingshengii IGTS8 and R. erythropolis XP suggesting that it might be transferable between members of these species.

Thiomonas bacteria are ubiquitous at acid mine drainage sites and play key roles in the remediation of water at these locations by oxidizing arsenite to arsenate, favouring the sorption of arsenic by iron oxides and their coprecipitation. Understanding the adaptive capacities of these bacteria is crucial to revealing how they persist and remain active in such extreme conditions. Interestingly, it was previously observed that after exposure to arsenite, when grown in a biofilm, some strains of Thiomonas bacteria develop variants that are more resistant to arsenic. Here, we identified the mechanisms involved in the emergence of such variants in biofilms. We found that the percentage of variants generated increased in the presence of high concentrations of arsenite (5.33 mM), especially in the detached cells after growth under biofilm-forming conditions. Analysis of gene expression in the parent strain CB2 revealed that genes involved in DNA repair were upregulated in the conditions where variants were observed. Finally, we assessed the phenotypes and genomes of the subsequent variants generated to evaluate the number of mutations compared to the parent strain. We determined that multiple point mutations accumulated after exposure to arsenite when cells were grown under biofilm conditions. Some of these mutations were found in what is referred to as ICE19, a genomic island (GI) carrying arsenic-resistance genes, also harbouring characteristics of an integrative and conjugative element (ICE). The mutations likely favoured the excision and duplication of this GI. This research aids in understanding how Thiomonas bacteria adapt to highly toxic environments, and, more generally, provides a window to bacterial genome evolution in extreme environments.

Amongst iron-oxidizing bacteria playing a key role in the natural attenuation of arsenic in acid mine drainages (AMDs), members of the Ferrovum genus were identified in mine effluent or water treatment plants, and were shown to dominate biogenic precipitates in field pilot experiments. In order to address the question of the in situ activity of the uncultivated Ferrovum sp. CARN8 strain in the Carnoulès AMD, we assembled its genome using metagenomic and metatranscriptomic sequences and we determined standardized expression values for protein-encoding genes. Our results showed that this microorganism was indeed metabolically active and allowed us to sketch out its metabolic activity in its natural environment. Expression of genes related to the respiratory chain and carbon fixation suggests aerobic energy production coupled to ferrous iron oxidation and chemolithoautotrophic growth. Notwithstanding the presence of nitrogenase genes in its genome, expression data also indicated that Ferrovum sp. CARN8 relied on ammonium import rather than nitrogen fixation. The expression of flagellum and chemotaxis genes hints that at least a proportion of this strain population was motile. Finally, apart from some genes related to metal resistance showing surprisingly low expression values, genes involved in stress response were well expressed as expected in AMDs.

The main 3'-5' exoribonucleolytic activity of eukaryotic cells is provided by the RNA exosome. The exosome is constituted by a core complex of nine subunits (Exo9), which coordinates the recruitment and the activities of distinct types of cofactors. The RNA exosome cofactors confer distributive and processive 3'-5' exoribonucleolytic, endoribonucleolytic, and RNA helicase activities. In addition, several RNA binding proteins and terminal nucleotidyltransferases also participate in the recognition of exosome RNA substrates.To fully understand the biological roles of the exosome, the respective functions of its cofactors must be deciphered. This entails the high-resolution analysis of 3' extremities of degradation or processing intermediates in different mutant backgrounds or growth conditions. Here, we describe a detailed 3' RACE-seq procedure for targeted mapping of exosome substrate 3' ends. This procedure combines a 3' RACE protocol with Illumina sequencing to enable the high-resolution mapping of 3' extremities and the identification of untemplated nucleotides for selected RNA targets.

Thymidine kinase (TK) is a key enzyme of the salvage pathway that recycles thymidine nucleosides to produce deoxythymidine triphosphate. Here, we identified the single TK of maize (Zea mays), denoted CPTK1, as necessary in the replication of the plastidial genome (cpDNA), demonstrating the essential function of the salvage pathway during chloroplast biogenesis. CPTK1 localized to both plastids and mitochondria, and its absence resulted in an albino phenotype, reduced cpDNA copy number and a severe deficiency in plastidial ribosomes. Mitochondria were not affected, indicating they are less reliant on the salvage pathway. Arabidopsis (Arabidopsis thaliana) TKs, TK1A and TK1B, apparently resulted from a gene duplication after the divergence of monocots and dicots. Similar but less-severe effects were observed for Arabidopsis tk1a tk1b double mutants in comparison to those in maize cptk1. TK1B was important for cpDNA replication and repair in conditions of replicative stress but had little impact on the mitochondrial phenotype. In the maize cptk1 mutant, the DNA from the small single-copy region of the plastidial genome was reduced to a greater extent than other regions, suggesting preferential abortion of replication in this region. This was accompanied by the accumulation of truncated genomes that resulted, at least in part, from unfaithful microhomology-mediated repair. These and other results suggest that the loss of normal cpDNA replication elicits the mobilization of new replication origins around the rpoB (beta subunit of plastid-encoded RNA polymerase) transcription unit and imply that increased transcription at rpoB is associated with the initiation of cpDNA replication.

Nicotiana otophora contains Agrobacterium-derived T-DNA sequences introduced by horizontal gene transfer (Chen et al., 2014). Sixty-nine contigs were assembled into four different cellular T-DNAs (cT-DNAs) totalling 83 kb. TC and TE result from two successive transformation events, each followed by duplication, yielding two TC and two TE inserts. TC is also found in other Nicotiana species, whereas TE is unique to N. otophora. Both cT-DNA regions are partially duplicated inverted repeats. Analysis of the cT-DNA divergence patterns allowed reconstruction of the evolution of the TC and TE regions. TC and TE carry 10 intact open reading frames. Three of these are TE-6b genes, derived from a single 6b gene carried by the Agrobacterium strain which inserted TE in the N. otophora ancestor. 6b genes have so far only been found in Agrobacterium tumefaciens or Agrobacterium vitis T-DNAs and strongly modify plant growth (Chen and Otten, 2016). The TE-6b genes were expressed in Nicotiana tabacum under the constitutive 2 × 35S promoter. TE-1-6b-R and TE-2-6b led to shorter plants, dark-green leaves, a strong increase in leaf vein development and modified petiole wings. TE-1-6b-L expression led to a similar phenotype, but in addition leaves show outgrowths at the margins, flowers were modified and plants became viviparous, i.e. embryos germinated in the capsules at an early stage of their development. Embryos could be rescued by culture in vitro. The TE-6b phenotypes are very different from the earlier described 6b phenotypes and could provide new insight into the mode of action of the 6b genes.

Bacteria of the genus Thiomonas are found ubiquitously in arsenic contaminated waters such as acid mine drainage (AMD), where they contribute to the precipitation and the natural bioremediation of arsenic. In these environments, these bacteria have developed a large range of resistance strategies among which the capacity to form particular biofilm structures. The biofilm formation is one of the most ubiquitous adaptive response observed in prokaryotes to various stresses, such as those induced in the presence of toxic compounds. This study focused on the process of biofilm formation in three Thiomonas strains (CB1, CB2 and CB3) isolated from the same AMD. The results obtained here show that these bacteria are all capable of forming biofilms, but the architecture and the kinetics of formation of these biofilms differ depending on whether arsenite is present in the environment and from one strain to another. Indeed, two strains favoured biofilm formation, whereas one favoured motility in the presence of arsenite. To identify the underlying mechanisms, the patterns of expression of some genes possibly involved in the process of biofilm formation were investigated in Thiomonas sp. CB2 in the presence and absence of arsenite, using a transcriptomic approach (RNA-seq). The findings obtained here shed interesting light on how the formation of biofilms, and the motility processes contribute to the adaptation of Thiomonas strains to extreme environments.

the presence of toxic compounds. This study focused on the process of bio fi lm formation in three Thiomonas strains (CB1, CB2 and CB3) isolated from the same AMD. The results obtained here show that these bacteria are all capable of forming bio fi lms, but the architecture and the kinetics of formation of these bio fi lms differ depending on whether arsenite is present in the environment and from one strain to another. Indeed, two strains favoured bio fi lm formation, whereas one favoured motility in the presence of arsenite. To identify the underlying mechanisms, the patterns of expression of some genes possibly involved in the process of bio fi lm formation were investigated in Thiomonas sp. CB2 in the presence and absence of arsenite, using a transcriptomic approach (RNA-seq). The fi ndings obtained here shed interesting light on how the formation of bio fi lms, and the motility processes contribute to the adaptation of Thiomonas strains to extreme environments.

Background Both the speciation and toxicity of arsenic are affected by bacterial transformations, i.e. oxidation, reduction or methylation. These transformations have a major impact on environmental contamination and more particularly on arsenic contamination of drinking water. Herminiimonas arsenicoxydans has been isolated from an arsenic- contaminated environment and has developed various mechanisms for coping with arsenic, including the oxidation of As(III) to As(V) as a detoxification mechanism. Results In the present study, a differential transcriptome analysis was used to identify genes, including arsenite oxidase encoding genes, involved in the response of H. arsenicoxydans to As(III). To get insight into the molecular mechanisms of this enzyme activity, a Tn 5 transposon mutagenesis was performed. Transposon insertions resulting in a lack of arsenite oxidase activity disrupted aoxR and aoxS genes, showing that the aox operon transcription is regulated by the AoxRS two-component system. Remarkably, transposon insertions were also identified in rpoN coding for the alternative N sigma factor (σ 54 ) of RNA polymerase and in dnaJ coding for the Hsp70 co-chaperone. Western blotting with anti-AoxB antibodies and quantitative RT-PCR experiments allowed us to demonstrate that the rpoN and dnaJ gene products are involved in the control of arsenite oxidase gene expression. Finally, the transcriptional start site of the aoxAB operon was determined using rapid amplification of cDNA ends (RACE) and a putative -12/-24 σ 54 -dependent promoter motif was identified upstream of aoxAB coding sequences. Conclusion These results reveal the existence of novel molecular regulatory processes governing arsenite oxidase expression in H. arsenicoxydans . These data are summarized in a model that functionally integrates arsenite oxidation in the adaptive response to As(III) in this microorganism.

Arsenic is a toxic metalloid known to generate an important oxidative stress in cells. In the present study, we focused our attention on an alga related to the genus Coccomyxa, exhibiting an extraordinary capacity to resist high concentrations of arsenite and arsenate. The integrated analysis of high-throughput transcriptomic data and non-targeted metabolomic approaches highlighted multiple levels of protection against arsenite. Indeed, Coccomyxa sp. Carn induced a set of transporters potentially preventing the accumulation of this metalloid in the cells and presented a distinct arsenic metabolism in comparison to another species more sensitive to that compound, i.e. Euglena gracilis, especially in regard to arsenic methylation. Interestingly, Coccomyxa sp. Carn was characterized by a remarkable accumulation of the strong antioxidant glutathione (GSH). Such observation could explain the apparent low oxidative stress in the intracellular compartment, as suggested by the transcriptomic analysis. In particular, the high amount of GSH in the cell could play an important role for the tolerance to arsenate, as suggested by its partial oxidation into oxidized glutathione in presence of this metalloid. Our results therefore reveal that this alga has acquired multiple and original defence mechanisms allowing the colonization of extreme ecosystems such as acid mine drainages.