How chromosomal rearrangements shape adaptation and speciation: Case studies in Drosophila pseudoobscura and its sibling species Drosophila persimilis

The gene arrangements of Drosophila have played a prominent role in the history of evolutionary biology from the original quantification of genetic diversity to current studies of the mechanisms for the origin and establishment of new inversion mutations within populations and their subsequent fixation between species supporting reproductive barriers. This review examines the genetic causes and consequences of inversions as recombination suppressors and the role that recombination suppression plays in establishing inversions in populations as they are involved in adaptation within heterogeneous environments. This often results in the formation of clines of gene arrangement frequencies among populations. Recombination suppression leads to the differentiation of the gene arrangements which may accelerate the accumulation of fixed genetic differences among populations. If these fixed mutations cause incompatibilities, then inversions pose important reproductive barriers between species. This review uses the evolution of inversions in Drosophila pseudoobscura and D. persimilis as a case study for how inversions originate, establish and contribute to the evolution of reproductive isolation.     

Ready on arrival: Standing variation at a chromosomal inversion contributes to rapid adaptation in an invasive marine crab

In this issue of Molecular Ecology, Tepolt et al. (2021) illustrate how the genetic architecture of adaptation and life history influence invasive success. A marvel of many invasive species is that they are incredibly successful despite evolutionary expectations that they will have low adaptive potential and suffer inbreeding depression due to initially small founding population sizes. Determining the combinations of ecoevolutionary factors that permit this apparent “genetic paradox of invasions” is an ongoing endeavour of invasive species research. Tepolt et al. (2021) study the European green crab in its invasive range on the North American west coast. Following a single introduction into California, this crab quickly spread across a wide latitude gradient, despite low diversity in the original founding population. Adaptation of this crab to clinal variation in temperature appeared largely driven by an inferred chromosomal inversion. This inversion exists as a balanced polymorphism in the European home range of green crabs and is associated with thermal tolerance. Tepolt et al. (2021) therefore demonstrate that adaptive evolution post introduction need not be impeded by bottlenecks if variation at key parts of the genome is available and can be maintained in introduced populations. Moreover, Tepolt et al. (2021) show how chromosomal inversions acting as large-effect loci might facilitate rapid responses to selection in introduced populations.     

Development and preliminary evaluation of a genomewide single nucleotide polymorphisms resource generated by RAD‐seq for the small yellow croaker (Larimichthys polyactis)

Recent advances in high‐throughput sequencing technologies have offered the possibility to generate genomewide sequence data to delineate previously unidentified genetic structure, obtain more accurate estimates of demographic parameters and to evaluate potential adaptive divergence. Here, we identified 27 556 single nucleotide polymorphisms for the small yellow croaker (Larimichthys polyactis) using restriction‐site‐associated DNA (RAD) sequencing of 24 individuals from two populations. Significant sources of genetic variation were identified, with an average nucleotide diversity (π) of 0.00105 ± 0.000425 across individuals, and long‐term effective population size was thus estimated to range between 26 172 and 261 716. According to the results, no differentiation between the two populations was detected based on the SNP data set of top quality score per contig or neutral loci. However, the two analysed populations were highly differentiated based on SNP data set of both top F_ST value per contig and the outlier SNPs. Moreover, local adaptation was highlighted by an F_ST‐based outlier tests implemented in LOSITAN and a total of 538 potentially locally selected SNPs were identified. blast2go annotation of contigs containing the outlier SNPs yielded hits for 37 (66%) of 56 significant blastx matches. Candidate genes for local adaptation constituted a wide array of biological functions, including cellular response to oxidative stress, actin filament binding, ion transmembrane transport and synapse assembly. The generated SNP resources in this study provided a valuable tool for future population genetics and genomics studies of L. polyactis.     

Linkage disequilibrium network analysis (LDna) gives a global view of chromosomal inversions,local adaptation and geographic structure

Recent advances in sequencing allow population‐genomic data to be generated for virtually any species. However, approaches to analyse such data lag behind the ability to generate it, particularly in nonmodel species. Linkage disequilibrium (LD, the nonrandom association of alleles from different loci) is a highly sensitive indicator of many evolutionary phenomena including chromosomal inversions, local adaptation and geographical structure. Here, we present linkage disequilibrium network analysis (LDna), which accesses information on LD shared between multiple loci genomewide. In LD networks, vertices represent loci, and connections between vertices represent the LD between them. We analysed such networks in two test cases: a new restriction‐site‐associated DNA sequence (RAD‐seq) data set for Anopheles baimaii, a Southeast Asian malaria vector; and a well‐characterized single nucleotide polymorphism (SNP) data set from 21 three‐spined stickleback individuals. In each case, we readily identified five distinct LD network clusters (single‐outlier clusters, SOCs), each comprising many loci connected by high LD. In A. baimaii, further population‐genetic analyses supported the inference that each SOC corresponds to a large inversion, consistent with previous cytological studies. For sticklebacks, we inferred that each SOC was associated with a distinct evolutionary phenomenon: two chromosomal inversions, local adaptation, population‐demographic history and geographic structure. LDna is thus a useful exploratory tool, able to give a global overview of LD associated with diverse evolutionary phenomena and identify loci potentially involved. LDna does not require a linkage map or reference genome, so it is applicable to any population‐genomic data set, making it especially valuable for nonmodel species.     

Distance,elevation and environment as drivers of diversity and divergence in bumble bees across latitude and altitude

Identifying drivers of dispersal limitation and genetic differentiation is a key goal in biogeography. We examine patterns of population connectivity and genetic diversity using restriction site‐associated DNA sequencing (RADseq) in two bumble bee species, Bombus vosnesenskii and Bombus bifarius, across latitude and altitude in mountain ranges from California, Oregon and Washington, U.S.A. Bombus vosnesenskii, which occurs across a broader elevational range at most latitudes, exhibits little population structure while B. bifarius, which occupies a relatively narrow higher elevation niche across most latitudes, exhibits much stronger population differentiation, although gene flow in both species is best explained by isolation with environmental niche resistance. A relationship between elevational habitat breadth and genetic diversity is also apparent, with B. vosnesenskii exhibiting relatively consistent levels of genetic diversity across its range, while B. bifarius has reduced genetic diversity at low latitudes, where it is restricted to high‐elevation habitat. The results of this study highlight the importance of the intersect between elevational range and habitat suitability in influencing population connectivity and suggest that future climate warming will have a fragmenting effect even on populations that are presently well connected, as they track their thermal niches upward in montane systems.     

Carbon sequestration in soil in a semi‐natural Miscanthus sinensis grassland and Cryptomeria japonica forest plantation in Aso,Kumamoto, Japan

Although Miscanthus sinensis grasslands (Misc‐GL) and Cryptomeria japonica forest plantations (Cryp‐FP) are proposed bioenergy feedstock systems, their relative capacity to sequester C may be an important factor in determining their potential for sustainable bioenergy production. Therefore, our objective was to quantify changes in soil C sequestration 47 years after a Misc‐GL was converted to a Cryp‐FP. The study was conducted on adjacent Misc‐GL and Cryp‐FP located on Mt. Aso, Kumamoto, Japan. After Cryp‐FP establishment, only the Misc‐GL continued to be managed by annual burning every March. Mass C and N, δ¹³C, and δ¹⁵N at 0–30 cm depth were measured in 5 cm increments. Carbon and N concentrations, C:N ratio, δ¹³C, and δ¹⁵N were measured in litter and/or ash, and rhizomes or roots. Although C input in Misc‐GL by M. sinensis was approximately 36% of that in Cryp‐FP by C. japonica, mean annual soil C sequestration in Misc‐GL (503 kg C ha^?1 yr^?1) was higher than that in Cryp‐FP (284 kg C ha^?1 yr^?1). This was likely the result of larger C input from aboveground litter to soil, C‐quality (C:N ratio and lignin concentration in aboveground litter) and possibly more recalcitrant C (charcoal) inputs by annual burning. The difference in soil δ¹⁵N between sites indicated that organic C with N had greater cycling between heterotrophic microbes and soil and produces more recalcitrant humus in Misc‐GL than in Cryp‐FP. Our data indicate that in terms of soil C sequestration, maintenance of Misc‐GL may be more advantageous than conversion to Cryp‐FP in Aso, Japan.     

Trophic structure of a coastal fish community determined with diet and stable isotope analyses

A combination of dietary guild analysis and nitrogen (δ¹⁵N) and carbon (δ¹³C) stable‐isotope analysis was used to assess the trophic structure of the fish community in Rhode Island and Block Island Sounds, an area off southern New England identified for offshore wind energy development. In the autumn of 2009, 2010 and 2011, stomach and tissue samples were taken from 20 fish and invertebrate species for analysis of diet composition and δ¹⁵N and δ¹³C signatures. The food chain in Rhode Island and Block Island Sounds comprises approximately four trophic levels within which the fish community is divided into distinct dietary guilds, including planktivores, benthivores, crustacivores and piscivores. Within these guilds, inter‐species isotopic and dietary overlap is high, suggesting that resource partitioning or competitive interactions play a major role in structuring the fish community. Carbon isotopes indicate that most fishes are supported by pelagic phytoplankton, although there is evidence that benthic production also plays a role, particularly for obligate benthivores such as skates Leucoraja spp. This type of analysis is useful for developing an ecosystem‐based approach to management, as it identifies species that act as direct links to basal resources as well as species groups that share trophic roles.     

Water temperature and fish growth: otoliths predict growth patterns of a marine fish in a changing climate

Ecological modeling shows that even small, gradual changes in body size in a fish population can have large effects on natural mortality, biomass, and catch. However, efforts to model the impact of climate change on fish growth have been hampered by a lack of long‐term (multidecadal) data needed to understand the effects of temperature on growth rates in natural environments. We used a combination of dendrochronology techniques and additive mixed‐effects modeling to examine the sensitivity of growth in a long‐lived (up to 70 years), endemic marine fish, the western blue groper (Achoerodus gouldii), to changes in water temperature. A multi‐decadal biochronology (1952–2003) of growth was constructed from the otoliths of 56 fish collected off the southwestern coast of Western Australia, and we tested for correlations between the mean index chronology and a range of potential environmental drivers. The chronology was significantly correlated with sea surface temperature in the region, but common variance among individuals was low. This suggests that this species has been relatively insensitive to past variations in climate. Growth increment and age data were also used in an additive mixed model to predict otolith growth and body size later this century. Although growth was relatively insensitive to changes in temperature, the model results suggested that a fish aged 20 in 2099 would have an otolith about 10% larger and a body size about 5% larger than a fish aged 20 in 1977. Our study shows that species or populations regarded as relatively insensitive to climate change could still undergo significant changes in growth rate and body size that are likely to have important effects on the productivity and yield of fisheries.     

Design, synthesis, and evaluation of a novel series of alpha-substituted phenylpropanoic acid derivatives as human peroxisome proliferator-activated receptor (PPAR) alpha/delta dual agonists for the treatment of metabolic syndrome

A series of alpha-alkyl-substituted phenylpropanoic acids was prepared as dual agonists of peroxisome proliferator-activated receptors alpha and delta (PPARalpha/delta). Structure-activity relationship studies indicated that the shape of the linking group and the shape of the substituent at the distal benzene ring play key roles in determining the potency and the selectivity of PPAR subtype transactivation. Structure-activity relationships among the amide series (10) and the reversed amide series (13) are similar, but not identical, especially in the case of the compounds bearing a bulky hydrophobic substituent at the distal benzene ring, indicating that the hydrophobic tail part of the molecules in these two series binds at somewhat different positions in the large binding pocket of PPAR. alpha-Alkyl-substituted phenylpropanoic acids of (S)-configuration were identified as potent human PPARalpha/delta dual agonists. Representative compounds exhibited marked nuclear receptor selectivity for PPARalpha and PPARdelta. Subtype-selective PPAR activation was also examined by analysis of the mRNA expression of PPAR-regulated genes.     

Trophic level and basal resource use of soil animals are hardly affected by local plant associations in abandoned arable land

Plants provide resources and shape the habitat of soil organisms thereby affecting the composition and functioning of soil communities. Effects of plants on soil communities are largely taxon‐dependent, but how different functional groups of herbaceous plants affect trophic niches of individual animal species in soil needs further investigation. Here, we studied the use of basal resources and trophic levels of dominating soil meso‐ and macrofauna using stable isotope ratios of carbon and nitrogen in arable fallow systems 3 and 14–16 years after abandonment. Animals were sampled from the rhizosphere of three plant species of different functional groups: a legume (Medicaco sativa), a nonlegume herb (Taraxacum officinale), and a grass (Bromus sterilis). We found virtually no consistent effects of plant identity on stable isotope composition of soil animals and on thirteen isotopic metrics that reflect general food‐web structure. However, in old fallows, the carbon isotope composition of some predatory macrofauna taxa had shifted closer to that of co‐occurring plants, which was particularly evident for Lasius, an aphid‐associated ant genus. Trophic levels and trophic‐chain lengths in food webs were similar across plant species and fallow ages. Overall, the results suggest that variations in local plant diversity of grassland communities may little affect the basal resources and the trophic level of prey consumed by individual species of meso‐ and macrofauna belowground. By contrast, successional changes in grassland communities are associated with shifts in the trophic niches of certain species, reflecting establishment of trophic interactions with time, which shapes the functioning and stability of soil food webs.     

Neuroprotection by quercetin via mitochondrial function adaptation in traumatic brain injury: PGC‐1α pathway as a potential mechanism

The aim of this study was to investigate the neuroprotective effects of quercetin in mouse models of traumatic brain injury (TBI) and the potential role of the PGC‐1α pathway in putative neuroprotection. Wild‐type mice were randomly assigned to four groups: the sham group, the TBI group, the TBI+vehicle group and the TBI+quercetin group. Quercetin, a dietary flavonoid used as a food supplement, significantly reduced TBI‐induced neuronal apoptosis and ameliorated mitochondrial lesions. It significantly accelerated the translocation of PGC‐1α protein from the cytoplasm to the nucleus. In addition, quercetin restored the level of cytochrome c, malondialdehyde and superoxide dismutase in mitochondria. Therefore, quercetin administration can potentially attenuate brain injury in a TBI model by increasing the activities of mitochondrial biogenesis via the mediation of the PGC‐1α pathway.     

Foliar δ13C response patterns along a moisture gradient arising from genetic variation and phenotypic plasticity in grassland species of Inner Mongolia

Plants depend upon both genetic differences and phenotypic plasticity to cope with environmental variation over different timescales. The spatial variation in foliar δ¹³C levels along a moisture gradient represents an overlay of genetic and plastic responses. We hypothesized that such a spatial variation would be more obvious than the variation arising purely from a plastic response to moisture change. Leymus chinensis and Stipa spp. were sampled from Inner Mongolia along a dry‐wet transect, and some of these species were transplanted to an area with a moisture gradient. For Stipa spp., the slope of foliar δ¹³C and mean annual precipitation along the transect was significantly steeper than that of foliar δ¹³C and mean annual precipitation after the watering treatment. For L. chinensis, there was a general decreasing trend in foliar δ¹³C under the different (increasing) watering levels; however, its populations showed an irregular relationship between foliar δ¹³C and moisture origin. Therefore, support for our hypothesis was obtained from Stipa spp., but not from L. chinensis.     

The role of immigration and local adaptation on fine‐scale genotypic and phenotypic population divergence in a less mobile passerine

Dispersal and local patterns of adaptation play a major role on the ecological and evolutionary trajectory of natural populations. In this study, we employ a combination of genetic (25 microsatellite markers) and field‐based information (seven study years) to analyse the impact of immigration and local patterns of adaptation in two nearby (< 7 km) blue tit (Cyanistes caeruleus) populations. We used genetic assignment analyses to identify immigrant individuals and found that dispersal rate is female‐biased (72%). Data on lifetime reproductive success indicated that immigrant females produced fewer local recruits than their philopatric counterparts whereas immigrant males recruited more offspring than those that remained in their natal location. In spite of the considerably higher immigration rates of females, our results indicate that, in absolute terms, their demographic and genetic impact in the receiving populations is lower than that in immigrant males. Immigrants often brought novel alleles into the studied populations and a high proportion of them were transmitted to their recruits, indicating that the genetic impact of immigrants is not ephemeral. Although only a few kilometres apart, the two study populations were genetically differentiated and showed strong divergence in different phenotypic and life‐history traits. An almost absent inter‐population dispersal, together with the fact that both populations receive immigrants from different source populations, is probably the main cause of the observed pattern of genetic differentiation. However, phenotypic differentiation (P_ST) for all the studied traits greatly exceeded neutral genetic differentiation (F_ST), indicating that divergent natural selection is the prevailing factor determining the evolutionary trajectory of these populations. Our study highlights the importance of integrating individual‐ and population‐based approaches to obtain a comprehensive view about the role of dispersal and natural selection on structuring the genotypic and phenotypic characteristics of natural populations.     

Target and non‐target effects of a spider venom toxin produced in transgenic cotton and tobacco plants

The peptide ω‐Hexatoxin‐Hv1a (Hvt) is one of the most studied spider toxins. Its insecticidal potential has been reported against species belonging to the arthropod orders Lepidoptera, Diptera and Orthoptera. The gene encoding Hvt has been transformed into cotton and tobacco to protect the plants from damage by lepidopteran pests. This study evaluated the expression of the ω‐HXTX‐Hv1a gene in transgenic plants, and the toxicity of plant‐expressed and purified Hvt on target lepidopteran insects and on several non‐target species. Transgenic Bollgard II cotton plants, which produce Cry1Ac and Cry2Ab2 and purified Cry2Ab2 protein were included in the study as comparators. LC₉₅ values of purified Hvt against Spodoptera littoralis and Heliothis virescens were 28.31 and 27.57 μg/ml of artificial diet, respectively. Larval mortality was 100% on Hvt‐transgenic tobacco plants but not on Hvt‐transgenic cotton, probably because of the significantly lower toxin expression level in the transgenic cotton line. Non‐target studies were conducted with larvae of the predators Chrysoperla carnea and Coccinella septempunctata, adults of the aphid parasitoid Aphidius colemani, and adult workers of the honey bee, Apis mellifera. Even at 40 μg/ml, Hvt did not adversely affect the four non‐target species. Purified Cry2Ab2 at 10 μg/ml also did not adversely affect any of the non‐target species. Our results show that Hvt might be useful for developing insecticidal plant varieties to control pest Lepidoptera.     

Identification of C‐terminal Hsp70‐interacting protein as a mediator of tumour necrosis factor action in osteoblast differentiation by targeting osterix for degradation

In patients with inflammatory arthritis, tumour necrosis factor (TNF)‐α are overproduced in inflamed joints. This leads to local erosion of cartilage and bone, periarticular osteopenia, as well as osteoporosis. But less is known regarding the molecular mechanisms that mediate the effect of TNF‐α on osteoblast function. The purpose of this study was to test that C terminus of Hsc70‐interacting protein (CHIP) has a specific role in suppressing the osteogenic activity of osteoblasts under inflammatory conditions. C2C12, MC3T3‐E1 and HEK293T cell lines were cultured and cotransfected with related plasmids. After transfection, the cells were cultured further in the presence or absence of murine TNF‐α and subjected to real time RT‐PCR, Western blot, Ubiquitination assay, Co‐immunoprecipitation, Luciferase reporter assay, Small interfering RNAs and Mineralization assay. The expression levels of TNF‐α‐induced CHIP and Osx were examined by RT‐PCR and Western blot analysis. Co‐immunoprecipitation and ubiquitination assays revealed ubiquitinated Osx, confirmed that CHIP indeed interacted with Osx and identified K55 and K386 residues as the ubiquitination sites in Osx, Luciferase reporter assay and Small interfering RNAs examined whether TNF‐α target the bone morphogenetic protein signalling through CHIP. We established stable cell lines with the overexpression of HA‐CHIP, Mineralization assay and CHIP siRNA demonstrated the important roles of CHIP on osteoblast function in conditions in which TNF‐α is overexpressed. We found that the K55 and K386 residues are ubiquitination site(s) in Osx, and that TNF‐α inhibits osteoblast differentiation by promoting Osx degradation through up‐regulation of E3 ubiquitin ligase CHIP in osteoblast. Thus, CHIP targets Osx for ubiquitination and degradation in osteoblasts after chronic exposure to TNF‐α, and inhibition of CHIP expression in osteoblasts may be a new mechanism to limit inflammation‐mediated osteoporosis by promoting their differentiation into osteoblasts.     

Relationships of growth,stable carbon isotope composition and anatomical properties of leaf and xylem in seven mulberry cultivars: a hint towards drought tolerance

• The fast growth of mulberry depends on high water consumption, but considerable variations in drought tolerance exist across different cultivars. Physiological and anatomical mechanisms are important to plant survival under drought. However, few research efforts have been made to reveal the relationships of these two aspects in relation to drought tolerance.
• In this study, growth rates, leaf functional physiology and anatomical characteristics of leaf and xylem of 1‐year‐old saplings of seven mulberry cultivars at a common garden were compared. Their relationships were also explored.
• Growth, leaf physiology and anatomy were significantly different among the tested cultivars. Foliar stable carbon isotope composition (δ¹³C) was negatively correlated with growth rates, and closely related to several leaf and xylem anatomical traits. Particularly, leaf thickness, predicted hydraulic conductivity and vessel element length jointly contributed 77% of the variability in δ¹³C. Cultivar Wupu had small stomata, intermediate leaf thickness, the smallest hydraulically weighted vessel diameter and highest vessel number, and higher δ¹³C; Yunguo1 had high abaxial stomatal density, low specific leaf area, moderate hydraulic conductivity and δ¹³C; these are beneficial features to reduce leaf water loss and drought‐induced xylem embolism in arid areas. Cultivar Liaolu11 had contrasting physiological and anatomical traits compared with the previous two cultivars, suggesting that it might be sensitive to drought.
• Our findings indicate that growth and δ¹³C are closely associated with both leaf and xylem anatomical characteristics in mulberry, which provides fundamental information to assist evaluation of drought tolerance in mulberry cultivars and in other woody trees.

     

miR‐148a inhibits early relapsed colorectal cancers and the secretion of VEGF by indirectly targeting HIF‐1α under non‐hypoxia/hypoxia conditions

Vascular endothelial growth factor (VEGF) is correlated with angiogenesis and early relapse of colorectal cancer (CRC). This study investigated the role of miR‐148a in the regulation of VEGF/angiogenesis and early relapse of CRC. We established a stable clone with miR‐148a expression in HCT116 and HT29 cell lines and created a hypoxic condition by using CoCl₂ to determine the underlying mechanism of miR‐148a. The effects of miR‐148a on the phosphoryl‐ERK (pERK)/hypoxia‐inducible factor‐1α (HIF‐1α)/VEGF pathway were evaluated through Western blotting and the inhibitory effect of miR‐148a on angiogenesis was demonstrated through a tube formation assay. Sixty‐three CRC tissues (28 early relapse and 35 non‐early relapse) were analysed to assess the relationship between miR‐148a and HIF‐1α/VEGF. The protein expression of pERK/HIF‐1α/VEGF in HCT116 and HT29 cells was significantly decreased by miR‐148a (all P < 0.05). The protein expression of VEGF/HIF‐1α was strongly inversely associated with the expression of miR‐148a in the 63 CRC tissue samples (all P < 0.05). Tube formation assay demonstrated that miR‐148a significantly obliterated angiogenesis. miR‐148a suppresses VEGF through down‐regulation of the pERK/HIF‐1α/VEGF pathway and might lead to the inhibition of angiogenesis; miR‐148a down‐regulation increased the early relapse rate of CRC. This demonstrates that miR‐148a is a potential diagnostic and therapeutic target.     

Phylogeography of var gene repertoires reveals fine‐scale geospatial clustering of Plasmodium falciparum populations in a highly endemic area

Plasmodium falciparum malaria is a major global health problem that is being targeted for progressive elimination. Knowledge of local disease transmission patterns in endemic countries is critical to these elimination efforts. To investigate fine‐scale patterns of malaria transmission, we have compared repertoires of rapidly evolving var genes in a highly endemic area. A total of 3680 high‐quality DBLα‐sequences were obtained from 68 P. falciparum isolates from ten villages spread over two distinct catchment areas on the north coast of Papua New Guinea (PNG). Modelling of the extent of var gene diversity in the two parasite populations predicts more than twice as many var gene alleles circulating within each catchment (Mugil = 906; Wosera = 1094) than previously recognized in PNG (Amele = 369). In addition, there were limited levels of var gene sharing between populations, consistent with local parasite population structure. Phylogeographic analyses demonstrate that while neutrally evolving microsatellite markers identified population structure only at the catchment level, var gene repertoires reveal further fine‐scale geospatial clustering of parasite isolates. The clustering of parasite isolates by village in Mugil, but not in Wosera was consistent with the physical and cultural isolation of the human populations in the two catchments. The study highlights the microheterogeneity of P. falciparum transmission in highly endemic areas and demonstrates the potential of var genes as markers of local patterns of parasite population structure.