Genome-wide association study of <Emphasis Type="Italic">Mycobacterium avium</Emphasis> subspecies <Emphasis Type="Italic">Paratuberculosis</Emphasis> infection in Chinese Holstein

Background

Paratuberculosis is a contagious, chronic and enteric disease in ruminants, which is caused by Mycobacterium avium subspecies paratuberculosis (MAP) infection, resulting in enormous economic losses worldwide. There is currently no effective cure for MAP infection or a vaccine, it is thus important to explore the genetic variants that contribute to host susceptibility to infection by MAP, which may provide a better understanding of the mechanisms of paratuberculosis and benefit animal genetic improvement. Herein we performed a genome-wide association study (GWAS) to identify genomic regions and candidate genes associated with susceptibility to MAP infection in dairy cattle.

Results

Using Illumina Bovine 50?K (54,609 SNPs) and GeneSeek HD (138,893 SNPs) chips, two analytical approaches were performed, GRAMMAR-GC and ROADTRIPS in 937 Chinese Holstein cows, among which individuals genotyped by the 50?K chip were imputed to HD SNPs with Beagle software. Consequently, 15 and 11 significant SNPs (P?<?5?×?10^??5) were identified with GRAMMAR-GC and ROADTDRIPS, respectively. A total of 10 functional genes were in proximity to (i.e., within 1?Mb) these SNPs, including IL4, IL5, IL13, IRF1, MyD88, PACSIN1, DEF6, TDP2, ZAP70 and CSF2. Functional enrichment analysis showed that these genes were involved in immune related pathways, such as interleukin, T cell receptor signaling pathways and inflammatory bowel disease (IBD), implying their potential associations with susceptibility to MAP infection. In addition, by examining the publicly available cattle QTLdb, a previous QTL for MAP was found to be overlapped with one of regions detected currently at 32.5?Mb on BTA23, where the TDP2 gene was anchored.

Conclusions

In conclusion, we identified 26 SNPs located on 15 chromosomes in the Chinese Holstein population using two GWAS strategies with high density SNPs. Integrated analysis of GWAS, biological functions and the reported QTL information helps to detect positional candidate genes and the identification of regions associated with susceptibility to MAP traits in dairy cattle.

     

Ectopic expression of cucumber (<Emphasis Type="Italic">Cucumis sativus</Emphasis> L.) <Emphasis Type="Italic">CsTIR</Emphasis>/<Emphasis Type="Italic">AFB</Emphasis> genes enhance salt tolerance in transgenic <Emphasis Type="Italic">Arabidopsis</Emphasis>

Auxin receptors TIR1/AFBs play an essential role in a series of signaling network cascades. These F-box proteins have also been identified to participate in different stress responses via the auxin signaling pathway in Arabidopsis. Cucumber (Cucumis sativus L.) is one of the most important crops worldwide, which is also a model plant for research. In the study herein, two cucumber homologous auxin receptor F-box genes CsTIR and CsAFB were cloned and studied for the first time. The deduced amino acid sequences showed a 78% identity between CsTIR and AtTIR1 and 76% between CsAFB and AtAFB2. All these proteins share similar characteristics of an F-box domain near the N-terminus, and several Leucine-rich repeat regions in the middle. Arabidopsis plants ectopically overexpressing CsTIR or CsAFB were obtained and verified. Shorter primary roots and more lateral roots were found in these transgenic lines with auxin signaling amplified. Results showed that expression of CsTIR/AFB genes in Arabidopsis could lead to higher seeds germination rates and plant survival rates than wild-type under salt stress. The enhanced salt tolerance in transgenic plants is probably caused by maintaining root growth and controlling water loss in seedlings, and by stabilizing life-sustaining substances as well as accumulating endogenous osmoregulation substances. We proposed that CsTIR/AFB-involved auxin signal regulation might trigger auxin mediated stress adaptation response and enhance the plant salt stress resistance by osmoregulation.     

Role of Arabidopsis <Emphasis Type="Italic">ABF1</Emphasis>/<Emphasis Type="Italic">3</Emphasis>/<Emphasis Type="Italic">4</Emphasis> during <Emphasis Type="Italic">det1</Emphasis> germination in salt and osmotic stress conditions

Key message

Arabidopsis det1 mutants exhibit salt and osmotic stress resistant germination. This phenotype requires HY5, ABF1, ABF3, and ABF4.

Abstract

While DE-ETIOLATED 1 (DET1) is well known as a negative regulator of light development, here we describe how det1 mutants also exhibit altered responses to salt and osmotic stress, specifically salt and mannitol resistant germination. LONG HYPOCOTYL 5 (HY5) positively regulates both light and abscisic acid (ABA) signalling. We found that hy5 suppressed the det1 salt and mannitol resistant germination phenotype, thus, det1 stress resistant germination requires HY5. We then queried publically available microarray datasets to identify genes downstream of HY5 that were differentially expressed in det1 mutants. Our analysis revealed that ABA regulated genes, including ABA RESPONSIVE ELEMENT BINDING FACTOR 3 (ABF3), are downregulated in det1 seedlings. We found that ABF3 is induced by salt in wildtype seeds, while homologues ABF4 and ABF1 are repressed, and all three genes are underexpressed in det1 seeds. We then investigated the role of ABF3, ABF4, and ABF1 in det1 phenotypes. Double mutant analysis showed that abf3, abf4, and abf1 all suppress the det1 salt/osmotic stress resistant germination phenotype. In addition, abf1 suppressed det1 rapid water loss and open stomata phenotypes. Thus interactions between ABF genes contribute to det1 salt/osmotic stress response phenotypes.

     

Horizontal transfer of the C-termini of <Emphasis Type="Italic">tccC</Emphasis> genes in <Emphasis Type="Italic">Photorhabdus</Emphasis> and <Emphasis Type="Italic">Xenorhabdus</Emphasis>

The high molecular weight insecticidal toxin complexes (Tcs), including four toxin-complex loci (tca, tcb, tcc and tcd), were first identified in Photorhabdus luminescens W14. Each member of tca, tcb or tcc is required for oral toxicity of Tcs. However, the sequence sources of the C-termini of tccC3, tccC4, tccC6 and tccC7 are unknown. Here, we performed a whole genome survey to identify the orthologs of Tc genes, and found 165 such genes in 14 bacterial genomes, including 40 genes homologous to tccC1-7 in P. luminescens TT01. The sequence sources of the C-termini of tccC2-6 were determined by sequence analysis. Further phylogenetic investigations suggested that the C-termini of 6 tccC genes experienced horizontal gene transfer events.     

Genome-wide characterization and stress-responsive expression profiling of <Emphasis Type="Italic">MCM</Emphasis> genes in <Emphasis Type="Italic">Brassica oleracea</Emphasis> and <Emphasis Type="Italic">Brassica rapa</Emphasis>

The Minichromosome maintenance protein [MCM (2-7)] complex is associated with helicase activity for replication fork formation during DNA replication. We identified and characterized each 12 putative MCM genes from Brassica oleracea and Brassica rapa. MCM genes were classified into nine groups according to their evolutionary relationships. A high number of syntenic regions were present on chromosomes C03 and A03 in B. oleracea and B. rapa, respectively, compared to the other chromosomes. Expression analysis showed that most of the MCM(2-7) helicase-subunit genes and their coregulating MCM genes were upregulated during hydroxyurea (HU) induced stress in B. oleracea. In B. rapa, MCM(2-7) helicase genes BrMCM2_2, BrMCM7_1, BrMCM7_2 and their co-regulating genes were upregulated during replication stress. During cold stress, BoMCM6 in B. oleracea and BrMCM5 in B. rapa were remarkably upregulated. During salt stress, BoMCM6_2, BoMCM7_1, BoMCM8, BoMCM9, and BoMCM10 were markedly upregulated in B. oleracea. Hence, our study identified the candidate MCM family genes those possess abiotic stress-responsive behavior and DNA replication stress tolerance. As the first genome-wide analysis of MCM genes in B. oleracea and B. rapa, this work provides a foundation to develop stress responsive plants. Further functional and molecular studies on MCM genes will be helpful to enhance stress tolerance in plants.     

<Emphasis Type="Italic">LMI1</Emphasis>-like genes involved in leaf margin development of <Emphasis Type="Italic">Brassica napus</Emphasis>

In rapeseed (Brassica napus L.), leaf margins are variable and can be entire, serrate, or lobed. In our previous study, the lobed-leaf gene (LOBED-LEAF 1, BnLL1) was mapped to a 32.1 kb section of B. napus A10. Two LMI1-like genes, BnaA10g26320D and BnaA10g26330D, were considered the potential genes that controlled the lobed-leaf trait in rapeseed. In the present study, these two genes and another homologous gene (BnaC04g00850D) were transformed into Arabidopsis thaliana (L.) Heynh. plants to identify their functions. All three LMI1-like genes of B. napus produced serrate leaf margins. The expression analysis indicated that the expression level of BnaA10g26320D determined the difference between lobed- and entire-leaved lines in rapeseed. Therefore, it is likely that BnaA10g26320D corresponds to BnLL1.     

The <Emphasis Type="Italic">angora</Emphasis> gene weakens the effect of interaction of the mutant genes <Emphasis Type="Italic">wellhaarig</Emphasis> and <Emphasis Type="Italic">waved alopecia</Emphasis> in mice

The interaction of the mutant genes wellhaarig (we) and waved alopecia (wal) in mice was earlier demonstrated in our laboratory. The we gene significantly accelerates the appearance of alopecia in double we/wewal/wal homozygotes as compared to that in single +/+wal/wal homozygotes. It has been found in this work that the mutant gene angora-Y (Fgf5 ^go-Y ) weakens the effect of interaction of the we and wal genes. The first signs of alopecia appear in mice of the we/wewal/wal genotype at the age of 14 days, in triple Fgf5 ^go-Y /Fgf5 ^go-Y we/wewal/wal homozygotes alopecia is observed seven days later, i. e., in 21-day-old animals. The progression of alopecia in triple homozygotes is expressed to a lesser degree than in double +/+we/wewal/wal homozygotes. A single dose of the Fgf5 ^go-Y gene also decreases the effect of interaction of the we and wal genes, but less than a double dose of this gene. The first signs of alopecia in mice of the +/Fgf5 ^go-Y we/wewal/wal genotype appear only three days later than in double +/+we/wewal/wal homozygotes. The data obtained demonstrate that the Fgf5 ^go-Y gene is a powerful modifier of mutant genes determining the process of alopecia.     

Association between <Emphasis Type="Italic">cagA</Emphasis>, <Emphasis Type="Italic">vacAi</Emphasis>, and <Emphasis Type="Italic">dupA</Emphasis> genes of <Emphasis Type="Italic">Helicobacter pylori</Emphasis> and gastroduodenal pathologies in Chilean patients

In addition to the already known cagA gene, novel genetic markers have been associated with Helicobacter pylori (H. pylori) virulence: the dupA and vacAi genes. These genes might play an important role as specific markers to determine the clinical outcome of the disease, especially the vacAi gene, which has been expected to be a good marker of severe pathologies like gastric adenocarcinoma. In the present study, the association of cagA, dupA, and vacAi genes with gastroduodenal pathologies in Chilean patients was studied. One hundred and thirty-two patients positive for H. pylori were divided into two groups—non-severe and severe gastric pathologies—and investigated for the presence of cagA, dupA, and vacAi H. pylori virulence genes by PCR. The cagA gene was detected in 20/132 patients (15.2%), the vacAi1 gene was detected in 54/132 patients (40.9%), the vacAi2 gene was detected in 26/132 patients (19.7%), and the dupA gene was detected in 50/132 (37.9%) patients. Logistic regression model analysis showed that the vacAi1 isoform gene in the infected strains and the severity of the diseases outcome were highly associated, causing severe gastric damage that may lead to gastric cancer (p < 0.0001; OR = 8.75; 95% CI 3.54–21.64). Conversely, cagA (p = 0.3507; OR = 1.62; 95% CI 0.59–4.45) and vacAi2 (p = 0.0114; OR = 3.09; 95% CI 1.26–7.60) genes were not associated with damage, while the dupA gene was associated significantly with non-severe clinical outcome (p = 0.0032; OR = 0.25; 95% CI 0.09–0.65). In addition, dupA gene exerts protection against severe gastric pathologies induced by vacAi1 by delaying the outcome of the disease by approximately 20 years.     

<Emphasis Type="Italic">In vivo</Emphasis> RNAi screen identifies candidate signaling genes required for collective cell migration in <Emphasis Type="Italic">Drosophila</Emphasis> ovary

Collective migration of loosely or closely associated cell groups is prevalent in animal development, physiological events, and cancer metastasis. However, our understanding of the mechanisms of collective cell migration is incomplete. Drosophila border cells provide a powerful in vivo genetic model to study collective migration and identify essential genes for this process. Using border cell-specific RNAi-silencing in Drosophila, we knocked down 360 conserved signaling transduction genes in adult flies to identify essential pathways and genes for border cell migration. We uncovered a plethora of signaling genes, a large proportion of which had not been reported for border cells, including Rack1 (Receptor of activated C kinase) and brk (brinker), mad (mother against dpp), and sax (saxophone), which encode three components of TGF-β signaling. The RNAi knock down phenotype was validated by clonal analysis of Rack1 mutants. Our data suggest that inhibition of Src activity by Rack1 may be important for border cell migration and cluster cohesion maintenance. Lastly, results from our screen not only would shed light on signaling pathways involved in collective migration during embryogenesis and organogenesis in general, but also could help our understanding for the functions of conserved human genes involved in cancer metastasis.     

Activation of stilbene synthesis in cell cultures of <Emphasis Type="Italic">Vitis amurensis</Emphasis> by calcium-dependent protein kinases <Emphasis Type="Italic">VaCPK1</Emphasis> and <Emphasis Type="Italic">VaCPK26</Emphasis>

Stilbenes, including trans-resveratrol (3,4′,5-trihydroxy-trans-stilbene), are known to exert beneficial health effects and contribute to plant biotic stress resistance. Much remains to be discovered about the cell signaling pathways regulating stilbene biosynthesis. It has recently been shown that overexpression of the calcium-dependent protein kinase VaCPK20 gene considerably increased t-resveratrol accumulation in cell cultures of Vitis amurensis. In this study, we analyzed the involvement of other CDPK family members, VaCPK1 and VaCPK26, on stilbene synthesis and biomass production by cell cultures of V. amurensis. We showed that overexpression of the VaCPK1 and 26 genes induced production of stilbenes by 1.7–4.6-fold (for VaCPK1) and by 2.5–6.2-fold (for VaCPK26) in several independently established cell lines compared to the empty vector-transformed control. Using HPLC-UV-MS, we detected five stilbenes in the grape cells: t-resveratrol diglucoside, t-piceid, t-resveratrol, ε- and δ-viniferin. The VaCPK1- and VaCPK26-transformed calli were capable of producing 1.4–3.1 and 1.8–4.9 mg/l of t-resveratrol, respectively (up to 0.4 for and 0.6 mg/g of dry weight for VaCPK26 and VaCPK1, respectively), while the control line synthesized only 0.5 mg/l of t-resveratrol (0.07 mg/g DW). The up-regulation of t-resveratrol production in the VaCPK1- and VaCPK26-overexpressing grape calli correlated with a significant up-regulation of stilbene synthase (STS) gene expression, especially VaSTS7. The data indicate that VaCPK1 and 26 genes, which are close homologues of VaCPK20, are positive regulators of stilbene biosynthesis in grapevine.     

Induction of two cyclotide-like genes <Emphasis Type="Italic">Zmcyc1</Emphasis> and <Emphasis Type="Italic">Zmcyc5</Emphasis> by abiotic and biotic stresses in <Emphasis Type="Italic">Zea mays</Emphasis>

Cyclotides are small plant disulfide-rich and cyclic proteins with a diverse range of biological activities. Cyclotide-like genes show key sequence features of cyclotides and are present in the Poaceae. In this study the cDNA of the nine cyclotide-like genes were cloned and sequenced using 3′RACE from Zea mays. The gene expression of two of these genes (Zmcyc1 and Zmcyc5) were analyzed by real-time PCR in response to biotic (Fusarium graminearum, Ustilago maydis and Rhopalosiphum maydis) and abiotic (mechanical wounding, water deficit and salinity) stresses, as well as in response to salicylic acid and methyl jasmonate elicitors to mimic biotic stresses. All isolated genes showed significant similarity to other cyclotide-like genes and were classified in two separate clusters. Both Zmcyc1 and Zmcyc5 were expressed in all studied tissues with the highest expression in leaves and lowest expression in roots. Wounding, methyl jasmonate and salicylic acid significantly induced the expression of Zmcyc1 and Zmcyc5 genes, but the higher expression was observed for Zmcyc1 as compared with Zmcyc5. Expression levels of these two genes were also induced in inoculated leaves with F. graminearum, U. maydis and also in response to insect infestation. In addition, the 1000-base-pairs (bp) upstream of the promoter of Zmcyc1 and Zmcyc5 genes were identified and analyzed using the PlantCARE database and consequently a large number of similar biotic and abiotic cis-regulatory elements were identified for these two genes.     

<Emphasis Type="Italic">SnRK2</Emphasis> Homologs in <Emphasis Type="Italic">Gossypium</Emphasis> and <Emphasis Type="Italic">GhSnRK2.6</Emphasis> Improved Salt Tolerance in Transgenic Upland Cotton and <Emphasis Type="Italic">Arabidopsis</Emphasis>

SnRK2s are a large family of plant-specific protein kinases, which play important roles in multiple abiotic stress responses in various plant species. But the family in Gossypium has not been well studied. Here, we identified 13, 10, and 13 members of the SnRK2 family from Gossypium raimondii, Gossypium arboreum, and Gossypium hirsutum, respectively, and analyzed the locations of SnRK2 homologs in chromosomes based on genome data of cotton species. Phylogenetic tree analysis of SnRK2 proteins showed that these families were classified into three groups. All SnRK2 genes were comprised of nine exons and eight introns, and the exon distributions and the intron phase of homolog genes among different cotton species were analogous. Moreover, GhSnRK2.6 was overexpressed in Arabidopsis and upland cotton, respectively. Under salt treatment, overexpressed Arabidopsis could maintain higher biomass accumulation than wild-type plants, and GhSnRK2.6 overexpression in cotton exhibited higher germination rate than the control. So, the gene GhSnRK2.6 could be utilized in cotton breeding for salt tolerance.     

A new <Emphasis Type="Italic">Em</Emphasis>-like protein from <Emphasis Type="Italic">Lactuca sativa</Emphasis>, <Emphasis Type="Italic">LsEm1</Emphasis>, enhances drought and salt stress tolerance in <Emphasis Type="Italic">Escherichia coli</Emphasis> and rice

Late embryogenesis abundant (LEA) proteins are closely related to abiotic stress tolerance of plants. In the present study, we identified a novel Em-like gene from lettuce, termed LsEm1, which could be classified into group 1 LEA proteins, and shared high homology with Cynara cardunculus Em protein. The LsEm1 protein contained three different 20-mer conserved elements (C-element, N-element, and M-element) in the C-termini, N-termini, and middle-region, respectively. The LsEm1 mRNAs were accumulated in all examined tissues during the flowering and mature stages, with a little accumulation in the roots and leaves during the seedling stage. Furthermore, the LsEm1 gene was also expressed in response to salt, dehydration, abscisic acid (ABA), and cold stresses in young seedlings. The LsEm1 protein could effectively reduce damage to the lactate dehydrogenase (LDH) and protect LDH activity under desiccation and salt treatments. The Escherichia coli cells overexpressing the LsEm1 gene showed a growth advantage over the control under drought and salt stresses. Moreover, LsEm1-overexpressing rice seeds were relatively sensitive to exogenously applied ABA, suggesting that the LsEm1 gene might depend on an ABA signaling pathway in response to environmental stresses. The transgenic rice plants overexpressing the LsEm1 gene showed higher tolerance to drought and salt stresses than did wild-type (WT) plants on the basis of the germination performances, higher survival rates, higher chlorophyll content, more accumulation of soluble sugar, lower relative electrolyte leakage, and higher superoxide dismutase activity under stress conditions. The LsEm1-overexpressing rice lines also showed less yield loss compared with WT rice under stress conditions. Furthermore, the LsEm1 gene had a positive effect on the expression of the OsCDPK9, OsCDPK13, OsCDPK15, OsCDPK25, and rab21 (rab16a) genes in transgenic rice under drought and salt stress conditions, implying that overexpression of these genes may be involved in the enhanced drought and salt tolerance of transgenic rice. Thus, this work paves the way for improvement in tolerance of crops by genetic engineering breeding.