Completion of the canonical pathway for assembly of anticancer drugs vincristine/vinblastine in Catharanthus roseus

Glycosylation is a key modification for most molecules including plant natural products, for example, flavonoids and isoflavonoids, and can enhance the bioactivity and bioavailability of the natural products. The crystal structure of plant rhamnosyltransferase UGT89C1 from Arabidopsis thaliana was determined, and the structures of UGT89C1 in complexes with UDP‐β‐l ‐rhamnose and acceptor quercetin revealed the detailed interactions between the enzyme and its substrates. Structural and mutational analysis indicated that Asp356, His357, Pro147 and Ile148 are key residues for sugar donor recognition and specificity for UDP‐β‐l ‐rhamnose. The mutant H357Q exhibited activity with both UDP‐β‐l ‐rhamnose and UDP‐glucose. Structural comparison and mutagenesis confirmed that His21 is a key residue as the catalytic base and the only catalytic residue involved in catalysis independently as UGT89C1 lacks the other catalytic Asp that is highly conserved in other reported UGTs and forms a hydrogen bond with the catalytic base His. Ser124 is located in the corresponding position of the catalytic Asp in other UGTs and is not able to form a hydrogen bond with His21. Mutagenesis further showed that Ser124 may not be important in its catalysis, suggesting that His21 and acceptor may form an acceptor‐His dyad and UGT89C1 utilizes a catalytic dyad in catalysis instead of catalytic triad. The information of structure and mutagenesis provides structural insights into rhamnosyltransferase substrate specificity and rhamnosylation mechanism.     

The assembly of (+)‐vincadifformine‐ and (−)‐tabersonine‐derived monoterpenoid indole alkaloids in Catharanthus roseus involves separate branch pathways

The biological activity of monoterpenoid indole alkaloids (MIAs) has led to their use in cancer treatment and other medical applications. Their biosynthesis has involved the formation of reactive intermediates by responsible enzymes to elaborate several different chemical scaffolds. Modification of scaffolds through different substitution reactions has produced chemically diverse MIAs and related biological activities. The present study characterizes the three‐step pathway involved in the formation of (+)‐echitovenine, the major O‐acetylated MIA of Catharanthus roseus roots, and differentiates it from a parallel pathway involved in the formation of hörhammericine. Separate hydrolases convert a common reactive MIA intermediate to aspidosperma skeletons of opposite specific rotations, that is (+)‐vincadifformine and (?)‐tabersonine, respectively. The formation of (+) minovincinine from (+) vincadifformine 19‐hydroxylase (V19H) is catalyzed by a root‐specific cytochrome P450 with high amino acid sequence similarity to the leaf‐specific tabersonine‐3‐hydroxylase involved in vindoline biosynthesis. Similarly, O‐acetylation of (+)‐minovincinine to form (+) echitovenine involves minovincinine‐O‐acetytransferase. The substrate specificity of V19H and MAT for their respective (+)‐enantiomers defines the separate enantiomer‐specific pathway involved in (+)‐echitovenine biosynthesis and differentiates it from a parallel (?)‐enantiomer‐specific pathway involved in the formation of hörhammericine from (?)‐tabersonine.     

Genome‐wide detection of signatures of selection in Korean Hanwoo cattle

The Korean Hanwoo cattle have been intensively selected for production traits, especially high intramuscular fat content. It is believed that ancient crossings between different breeds contributed to forming the Hanwoo, but little is known about the genomic differences and similarities between other cattle breeds and the Hanwoo. In this work, cattle breeds were grouped by origin into four types and used for comparisons: the Europeans (represented by six breeds), zebu (Nelore), African taurine (N'Dama) and Hanwoo. All animals had genotypes for around 680 000 SNPs after quality control of genotypes. Average heterozygosity was lower in Nelore and N'Dama (0.22 and 0.21 respectively) than in Europeans (0.26–0.31, with Shorthorn as outlier at 0.24) and Hanwoo (0.29). Pairwise F_ST analyses demonstrated that Hanwoo are more related to European cattle than to Nelore, with N'Dama in an intermediate position. This finding was corroborated by principal components and unsupervised hierarchical clustering. Using genome‐wide smoothed F_ST, 55 genomic regions potentially under positive selection in Hanwoo were identified. Among these, 29 were regions also detected in previous studies. Twenty‐four regions were exclusive to Hanwoo, and a number of other regions were shared with one or two of the other groups. These regions overlap a number of genes that are related to immune, reproduction and fatty acid metabolism pathways. Further analyses are needed to better characterize the ancestry of the Hanwoo cattle and to define the genes responsible to the identified selection peaks.     

Genome‐wide prediction of age at puberty and reproductive longevity in sows

Traditional selection for sow reproductive longevity is ineffective due to low heritability and late expression of the trait. Incorporation of DNA markers into selection programs is potentially a more practical approach for improving sow lifetime productivity. Using a resource population of crossbred gilts, we explored pleiotropic sources of variation that influence age at puberty and reproductive longevity. Of the traits recorded before breeding, only age at puberty significantly affected the probability that females would produce a first parity litter. The genetic variance explained by 1‐Mb windows of the sow genome, compared across traits, uncovered regions that influence both age at puberty and lifetime number of parities. Allelic variants of SNPs located on SSC5 (27–28 Mb), SSC8 (36–37 Mb) and SSC12 (1.2–2 Mb) exhibited additive effects and were associated with both early expression of puberty and a greater than average number of lifetime parities. Combined analysis of these SNPs showed that an increase in the number of favorable alleles had positive impact on reproductive longevity, increasing number of parities by up to 1.36. The region located on SSC5 harbors non‐synonymous alleles in the arginine vasopressin receptor 1A (AVPR1A) gene, a G‐protein‐coupled receptor associated with social and reproductive behaviors in voles and humans and a candidate for the observed effects. This region is characterized by high levels of linkage disequilibrium in different lines and could be exploited in marker‐assisted selection programs across populations to increase sow reproductive longevity.     

Genome‐wide association reveals the locus responsible for four‐horned ruminant

Phenotypic variability in horn characteristics, such as their size, number and shape, offers the opportunity to elucidate the molecular basis of horn development. The objective of this study was to map the genetic determinant controlling the production of four horns in two breeds, Jacob sheep and Navajo‐Churro, and examine whether an eyelid abnormality occurring in the same populations is related. Genome‐wide association mapping was performed using 125 animals from the two breeds that contain two‐ and four‐horned individuals. A case–control design analysis of 570 712 SNPs genotyped with the ovine HD SNP Beadchip revealed a strong association signal on sheep chromosome 2. The 10 most strongly associated SNPs were all located in a region spanning Mb positions 131.9–132.6, indicating the genetic architecture underpinning the production of four horns is likely to involve a single gene. The closest genes to the most strongly associated marker (OAR2_132568092) were MTX2 and the HOXD cluster, located approximately 93 Kb and 251 Kb upstream respectively. The occurrence of an eyelid malformation across both breeds was restricted to polled animals and those carrying more than two horns. This suggests the eyelid abnormality may be associated with departures from the normal developmental production of two‐horned animals and that the two conditions are developmentally linked. This study demonstrated the presence of separate loci responsible for the polled and four‐horned phenotypes in sheep and advanced our understanding of the complexity that underpins horn morphology in ruminants.     

Camptotheca acuminata 10‐hydroxycamptothecin O‐methyltransferase: an alkaloid biosynthetic enzyme co‐opted from flavonoid metabolism

The medicinal plant Camptotheca acuminata accumulates camptothecin, 10‐hydroxycamptothecin, and 10‐methoxycamptothecin as its major bioactive monoterpene indole alkaloids. Here, we describe identification and functional characterization of 10‐hydroxycamptothecin O‐methyltransferase (Ca10OMT), a member of the Diverse subclade of class II OMTs. Ca10OMT is highly active toward both its alkaloid substrate and a wide range of flavonoids in vitro and in this way contrasts with other alkaloid OMTs in the subclade that only utilize alkaloid substrates. Ca10OMT shows a strong preference for the A‐ring 7‐OH of flavonoids, which is structurally equivalent to the 10‐OH of 10‐hydroxycamptothecin. The substrates of other alkaloid OMTs in the subclade bear little similarity to flavonoids, but the 3‐D positioning of the 7‐OH, A‐ and C‐rings of flavonoids is nearly identical to the 10‐OH, A‐ and B‐rings of 10‐hydroxycamptothecin. This structural similarity likely explains the retention of flavonoid OMT activity by Ca10OMT and also why kaempferol and quercetin aglycones are potent inhibitors of its 10‐hydroxycamptothecin activity. The catalytic promiscuity and strong inhibition of Ca10OMT by flavonoid aglycones in vitro prompted us to investigate the potential physiological roles of the enzyme in vivo. Based on its regioselectivity, kinetic parameters and absence of 7‐OMT flavonoids in vivo, we conclude that the major and likely only substrate of Ca10OMTin vivo is 10‐hydroxycamptothecin. This is likely accomplished by Ca10OMT being kept spatially separated at the tissue levels from potentially inhibitory flavonoid aglycones, and flavonoid aglycones being rapidly glycosylated to non‐inhibitory flavonoid glycosides.     

Homology‐dependent repair is involved in 45S rDNA loss in plant CAF‐1 mutants

Arabidopsis thaliana mutants in FAS1 and FAS2 subunits of chromatin assembly factor 1 (CAF1) show progressive loss of 45S rDNA copies and telomeres. We hypothesized that homology‐dependent DNA damage repair (HDR) may contribute to the loss of these repeats in fas mutants. To test this, we generated double mutants by crossing fas mutants with knock‐out mutants in RAD51B, one of the Rad51 paralogs of A. thaliana. Our results show that the absence of RAD51B decreases the rate of rDNA loss, confirming the implication of RAD51B‐dependent recombination in rDNA loss in the CAF1 mutants. Interestingly, this effect is not observed for telomeric repeat loss, which thus differs from that acting in rDNA loss. Involvement of DNA damage repair in rDNA dynamics in fas mutants is further supported by accumulation of double‐stranded breaks (measured as γ‐H2AX foci) in 45S rDNA. Occurrence of the foci is not specific for S‐phase, and is ATM‐independent. While the foci in fas mutants occur both in the transcribed (intranucleolar) and non‐transcribed (nucleoplasmic) fraction of rDNA, double fas rad51b mutants show a specific increase in the number of the intranucleolar foci. These results suggest that the repair of double‐stranded breaks present in the transcribed rDNA region is RAD51B dependent and that this contributes to rDNA repeat loss in fas mutants, presumably via the single‐stranded annealing recombination pathway. Our results also highlight the importance of proper chromatin assembly in the maintenance of genome stability.     

Sex‐specific responses of Asian citrus psyllid to volatiles of conspecific and host‐plant origin

Diaphorina citri Kuwayama (Hemiptera: Psyllidae) is the primary vector of Candidatus Liberibacter spp. bacteria that cause citrus greening, a disease of worldwide importance. Olfactometry was employed to test responses of D. citri to odours from intact citrus plants (Mexican lime, Citrus aurantifolia, sour orange, Citrus aurantium, Marsh grapefruit, Citrus paradisi and Valencia orange, Citrus sinensis), citrus plants previously infested with D. citri, and odours of conspecifics including nymphs, adult insects of same and opposite sex, and their products (honeydew), both alone and in combination. In contrast to other studies, psyllids of both sexes were attracted to volatiles of undamaged Mexican lime leaves, whereas undamaged grapefruit attracted only females, and leaves of Valencia and sour orange did not attract either sex. All four plant species attracted female psyllids when previously infested, but only Mexican lime and sour orange‐attracted males. Thus, Citrus species appear to vary in the production of both constituitive and induced volatiles that attract adult psyllids. Volatiles emitted by nymphs did not attract either sex, but psyllid honeydew was attractive to males, likely due to female pheromone residues. Males oriented to the odour of females, whereas the reverse was not true, and neither males nor females oriented to same‐sex volatiles. The addition of conspecific cues (adults, nymphs or honeydew) did not increase female attraction to previously infested leaves, but male response was increased by the presence of adults and honeydew, regardless of plant species. Thus, female psyllids appear to orient more strongly to volatiles of plant origin, whereas males respond more strongly to cues emanating from females and conspecific excretions. These results suggest that female psyllids drive the initial colonization of host plants, whereas males orient to females and infested plants. Identification of the specific volatiles involved may permit their use in monitoring and management of this pest.     

Genome‐wide association study for body weight in cattle populations from Siberia

Body weight is a complex trait in cattle associated with commonly used commercial breeding measurements related to growth. Although many quantitative trait loci (QTL) for body weight have been identified in cattle so far, searching for genetic determinants in different breeds or environments is promising. Therefore, we carried out a genome‐wide association study (GWAS) in two cattle populations from the Russian Federation (Siberian region) using the GGP HD150K array containing 139 376 single nucleotide polymorphism (SNP) markers. Association tests for 107 550 SNPs left after filtering revealed five statistically significant SNPs on BTA5, considering a false discovery rate of less than 0.05. The chromosomal region containing these five SNPs contains the CCND2 gene, which was previously associated with average daily weight gain and body mass index in US beef cattle populations and in humans respectively. Our study is the first GWAS for body weight in beef cattle populations from the Russian Federation. The results provided here suggest that, despite the existence of breed‐ and species‐specific QTL, the genetic architecture of body weight could be evolutionarily conserved in mammals.     

Genome‐wide screening for DNA variants associated with reading and language traits

Reading and language abilities are heritable traits that are likely to share some genetic influences with each other. To identify pleiotropic genetic variants affecting these traits, we first performed a genome‐wide association scan (GWAS) meta‐analysis using three richly characterized datasets comprising individuals with histories of reading or language problems, and their siblings. GWAS was performed in a total of 1862 participants using the first principal component computed from several quantitative measures of reading‐ and language‐related abilities, both before and after adjustment for performance IQ. We identified novel suggestive associations at the SNPs rs59197085 and rs5995177 (uncorrected P ≈ 10^–7 for each SNP), located respectively at the CCDC136/FLNC and RBFOX2 genes. Each of these SNPs then showed evidence for effects across multiple reading and language traits in univariate association testing against the individual traits. FLNC encodes a structural protein involved in cytoskeleton remodelling, while RBFOX2 is an important regulator of alternative splicing in neurons. The CCDC136/FLNC locus showed association with a comparable reading/language measure in an independent sample of 6434 participants from the general population, although involving distinct alleles of the associated SNP. Our datasets will form an important part of on‐going international efforts to identify genes contributing to reading and language skills.     

Genome‐wide association study identifies variants in the CAPN9 gene associated with umbilical hernia in pigs

Pig umbilical hernia (UH) affects pig welfare and brings considerable economic loss to the pig industry. To date, the molecular mechanisms underlying pig UH are still poorly understood. To identify potential loci for susceptibility to this disease, we performed a genome‐wide association study in an Erhualian × Shaziling F₂ intercross population. A total of 45 animals were genotyped using Illumina Porcine SNP60 BeadChips. We observed a SNP (rs80993347) located in the calpain‐9 (CAPN9) gene on Sus scrofa chromosome 14 that was significantly associated with UH (P = 1.97 × 10^?10). Then, we identified a synonymous mutation rs321865883 (g.20164T>C) in exon 10 of the CAPN9 gene that distinguished two affected individuals (CC) from their normal full‐sibs (TC). Finally, quantitative polymerase chain reaction was explored to investigate the mRNA expression profile of the CAPN9 gene in 12 tissues in Yorkshire pigs at different developmental stages (3, 90 and 180 days). CAPN9 showed high expression levels in the gastrointestinal tract at these three growth stages. The results of this study indicate that the CAPN9 gene might be implicated in UH. Further studies are required to establish a role of CAPN9 in pig UH.     

Genome‐wide association study identifies an NLR gene that confers partial resistance to Magnaporthe oryzae in rice

Because of the frequent breakdown of major resistance (R) genes, identification of new partial R genes against rice blast disease is an important goal of rice breeding. In this study, we used a core collection of the Rice Diversity Panel II (C‐RDP‐II), which contains 584 rice accessions and are genotyped with 700 000 single‐nucleotide polymorphism (SNP) markers. The C‐RDP‐II accessions were inoculated with three blast strains collected from different rice‐growing regions in China. Genome‐wide association study identified 27 loci associated with rice blast resistance (LABRs). Among them, 22 LABRs were not associated with any known blast R genes or QTLs. Interestingly, a nucleotide‐binding site leucine‐rich repeat (NLR) gene cluster exists in the LABR12 region on chromosome 4. One of the NLR genes is highly conserved in multiple partially resistant rice cultivars, and its expression is significantly up‐regulated at the early stages of rice blast infection. Knockout of this gene via CRISPR‐Cas9 in transgenic plants partially reduced blast resistance to four blast strains. The identification of this new non‐strain specific partial R gene, tentatively named rice blast Partial Resistance gene 1 (PiPR1), provides genetic material that will be useful for understanding the partial resistance mechanism and for breeding durably resistant cultivars against blast disease of rice.     

Genome‐wide sequence information reveals recurrent hybridization among diploid wheat wild relatives

Many conflicting hypotheses regarding the relationships among crops and wild species closely related to wheat (the genera Aegilops, Amblyopyrum, and Triticum) have been postulated. The contribution of hybridization to the evolution of these taxa is intensely discussed. To determine possible causes for this, and provide a phylogeny of the diploid taxa based on genome‐wide sequence information, independent data were obtained from genotyping‐by‐sequencing and a target‐enrichment experiment that returned 244 low‐copy nuclear loci. The data were analyzed using Bayesian, likelihood and coalescent‐based methods. D statistics were used to test if incomplete lineage sorting alone or together with hybridization is the source for incongruent gene trees. Here we present the phylogeny of all diploid species of the wheat wild relatives. We hypothesize that most of the wheat‐group species were shaped by a primordial homoploid hybrid speciation event involving the ancestral Triticum and Am. muticum lineages to form all other species except Ae. speltoides. This hybridization event was followed by multiple introgressions affecting all taxa except Triticum. Mostly progenitors of the extant species were involved in these processes, while recent interspecific gene flow seems insignificant. The composite nature of many genomes of wheat‐group taxa results in complicated patterns of diploid contributions when these lineages are involved in polyploid formation, which is, for example, the case for tetraploid and hexaploid wheats. Our analysis provides phylogenetic relationships and a testable hypothesis for the genome compositions in the basic evolutionary units within the wheat group of Triticeae.     

Genome‐wide association studies reveal genetic loci associated with plasma cholinesterase activity in ducks

Plasma cholinesterase (PCHE) activity is an important auxiliary test in human clinical medicine. It can distinguish liver diseases from non‐liver diseases and help detect organophosphorus poisoning. Animal experiments have confirmed that PCHE activity is associated with obesity and hypertension and changes with physiological changes in an animal's body. The objective of this study was to locate the genetic loci responsible for PCHE activity variation in ducks. PCHE activity of Pekin duck × mallard F2 ducks at 3 and 8 weeks of age were analyzed, and genome‐wide association studies were conducted. A region of about 1.5 Mb (21.8–23.3 Mb) on duck chromosome 9 was found to be associated with PCHE activity at both 3 and 8 weeks of age. The top SNP, g.22643979C>T in the butyrylcholinesterase (BCHE) gene, was most highly associated with PCHE activity at 3 weeks (?logP = 21.45) and 8 weeks (?logP = 27.60) of age. For the top SNP, the strong associations of CC and CT genotypes with low PCHE activity and the TT genotype with high PCHE activity indicates the dominant inheritance of low PCHE activity. Problems with block inheritance or linkage exist in this region. This study supports that BCHE is a functional gene for determining PCHE levels in ducks and that the genetic variations around this gene can cause phenotypic variations of PCHE activity.     

Genome‐wide scans detect adaptation to aridity in a widespread forest tree species

Patterns of adaptive variation within plant species are best studied through common garden experiments, but these are costly and time‐consuming, especially for trees that have long generation times. We explored whether genome‐wide scanning technology combined with outlier marker detection could be used to detect adaptation to climate and provide an alternative to common garden experiments. As a case study, we sampled nine provenances of the widespread forest tree species, Eucalyptus tricarpa, across an aridity gradient in southeastern Australia. Using a Bayesian analysis, we identified a suite of 94 putatively adaptive (outlying) sequence‐tagged markers across the genome. Population‐level allele frequencies of these outlier markers were strongly correlated with temperature and moisture availability at the site of origin, and with population differences in functional traits measured in two common gardens. Using the output from a canonical analysis of principal coordinates, we devised a metric that provides a holistic measure of genomic adaptation to aridity that could be used to guide assisted migration or genetic augmentation.     

A genomic investigation of ecological differentiation between free‐living and Drosophila‐associated bacteria

Various bacterial taxa have been identified both in association with animals and in the external environment, but the extent to which related bacteria from the two habitat types are ecologically and evolutionarily distinct is largely unknown. This study investigated the scale and pattern of genetic differentiation between bacteria of the family Acetobacteraceae isolated from the guts of Drosophila fruit flies, plant material and industrial fermentations. Genome‐scale analysis of the phylogenetic relationships and predicted functions was conducted on 44 Acetobacteraceae isolates, including newly sequenced genomes from 18 isolates from wild and laboratory Drosophila. Isolates from the external environment and Drosophila could not be assigned to distinct phylogenetic groups, nor are their genomes enriched for any different sets of genes or category of predicted gene functions. In contrast, analysis of bacteria from laboratory Drosophila showed they were genetically distinct in their universal capacity to degrade uric acid (a major nitrogenous waste product of Drosophila) and absence of flagellar motility, while these traits vary among wild Drosophila isolates. Analysis of the competitive fitness of Acetobacter discordant for these traits revealed a significant fitness deficit for bacteria that cannot degrade uric acid in culture with Drosophila. We propose that, for wild populations, frequent cycling of Acetobacter between Drosophila and the external environment prevents genetic differentiation by maintaining selection for traits adaptive in both the gut and external habitats. However, laboratory isolates bear the signs of adaptation to persistent association with the Drosophila host under tightly defined environmental conditions.