Marta Fonseca Martins

Nosema ceranae has been present in Brazil for more than three decades infecting Africanized honey bees.

Until the mid-1990s, the only microsporidium known to infect bees of the genus Apis was Nosema apis. A second species, Nosema ceranae, was first identified in 1996 from Asian honey bees; it is postulated that this parasite was transmitted from the Asian honey bee, Apis cerana, to the European honey bee, Apis mellifera. Currently, N. ceranae is found on all continents and has often been associated with honey bee colony collapse and other reports of high bee losses. Samples of Africanized drones collected in 1979, preserved in alcohol, were analyzed by light microscopy to count spores and were subjected to DNA extraction, after which duplex PCR was conducted. All molecular analyses (triplicate) indicated that the drones were infected with both N. ceranae and N. apis. PCR products were sequenced and matched to sequences reported in the GenBank (Acc. Nos. JQ639316.1 and JQ639301.1). The venation pattern of the wings of these males was compared to those of the current population living in the same area and with the pattern of drones collected in 1968 from Ribeirão Preto, SP, Brazil, from a location close to where African swarms first escaped in 1956. The morphometric results indicated that the population collected in 1979 was significantly different from the current living population, confirming its antiquity. Considering that the use of molecular tools for identifying Nosema species is relatively recent, it is possible that previous reports of infections (which used only light microscopy, without ultrastructural analysis) wrongly identified N. ceranae as N. apis. Although we can conclude that N. ceranae has been affecting Africanized honeybees in Brazil for at least 34 years, the impact of this pathogen remains unclear.

A new approach to characterization of the resistance of populations of Rhipicephalus microplus (Acari: Ixodidae) to organophosphate and pyrethroid in the state of Minas Gerais, Brazil

The monitoring of resistance of cattle tick populations in Brazil to the chemical bases in use is largely limited to investigation of the phenotypic profile. There are few studies investigating the role played by the genotypic profile in acaricide resistance in the country. Therefore, the aim of the present study was to carry out molecular characterization and trace out the genetic profile of populations of Rhipicephalus microplus with respect to resistance to the organophosphate and pyrethroid chemical groups. For that purpose, larvae were genotyped belonging to 587 populations for pyrethroids and 306 for organophosphates, using the polymerase chain reaction technique. It was found that 75.49% and 97.44% of the larvae studied showed resistance to the organophosphates and pyrethroids, respectively. Among the populations resistant to pyrethroids, 91.9% were heterozygotes, showing that most of the resistant populations have only one allele responsible for resistance. Therefore, it is possible to conclude that the genotyped populations have high resistance to organophosphates, and even more so to pyrethroids. This information is fundamental for understanding the mechanisms of resistance of R. microplus to acaricides, to enable improvement of control techniques.

Honey bee colonies act as reservoirs for two Spiroplasma facultative symbionts and incur complex, multiyear infection dynamics

Two species of Spiroplasma (Mollicutes) bacteria were isolated from and described as pathogens of the European honey bee, Apis mellifera, ~30 years ago but recent information on them is lacking despite global concern to understand bee population declines. Here we provide a comprehensive survey for the prevalence of these two Spiroplasma species in current populations of honey bees using improved molecular diagnostic techniques to assay multiyear colony samples from North America (U.S.A.) and South America (Brazil). Significant annual and seasonal fluctuations of Spiroplasma apis and Spiroplasma melliferum prevalence in colonies from the U.S.A. (n = 616) and Brazil (n = 139) occurred during surveys from 2011 through 2013. Overall, 33% of U.S.A. colonies and 54% of Brazil colonies were infected by Spiroplasma spp., where S. melliferum predominated over S. apis in both countries (25% vs. 14% and 44% vs. 38% frequency, respectively). Colonies were co‐infected by both species more frequently than expected in both countries and at a much higher rate in Brazil (52%) compared to the U.S.A. (16.5%). U.S.A. samples showed that both species were prevalent not only during spring, as expected from prior research, but also during other seasons. These findings demonstrate that the model of honey bee spiroplasmas as springtime‐restricted pathogens needs to be broadened and their role as occasional pathogens considered in current contexts.

Differential expression of genes during mastitis in Holstein-Zebu crossbreed dairy cows.

Among the potential public health problems of animal production, infectious-contagious diseases stand out. Mastitis is among the main diseases affecting dairy cattle. One of the most promising options to reduce the problems caused by this disease, besides proper sanitary and management practices, is selective breeding of resistant animals. To shed light on the immune response mechanisms involved in the resistance/susceptibility phenotype to this disease, we quantified the relative expression of the genes IL-2, IL-6, IL-8, IL-12, IFN-γ, TNF-α, TLR-2, SEMA5A, and FEZL in cells of crossbreed dairy cows, divided into two groups, one healthy and the other suffering from clinical mastitis. Total RNA was extracted from the cells in the milk from the animals in each group (with and without clinical mastitis). Gene expression was determined using the real-time PCR method. The levels of gene expression were compared, and the cows with mastitis were found to express 2.5 times more TLR-2 than those free of mastitis (P < 0.05). There were no significant differences in the expression of the other genes.

Microfluidic chip electrophoresis investigation of major milk proteins: study of buffer effects and quantitative approaching

The separation and quantification of major milk proteins are fundamental in dairy research. Therefore, accurate and rapid methods are profoundly important. The microfluidic chip technique is faster, and uses considerably fewer chemicals and materials than traditional techniques. The objective of this study was to improve experimental methods for separating and quantifying major milk proteins using the microfluidic chip technique. Deionized water, a total protein solubilization buffer (TPS buffer) and a separating milk protein buffer (SEP buffer) were used for the treatment of milk samples and their effects were evaluated. The results showed an excellent separation for whey proteins with α-lactalbumin migrating first, followed by β-lactoglobulin in the presence of both buffers. However, better results for major casein separation were obtained when SEP buffer was added. The order of the migration time was: β-casein first, followed by αs-casein and κ-casein. The quantitative analysis showed significant differences among the percentages of protein fractions from both buffers. The results of microfluidic chip technology using the SEP buffer solution were comparable to those obtained by SDS-PAGE for these proteins and to the data reported in the literature.

Transcript profiling of expressed sequence tags from semimembranosus muscle of commercial and naturalized pig breeds.

In general, genetic differences across different breeds of pig lead to variation in mature body size and slaughter age. The Commercial breeds Duroc and Large White and the local Brazilian breed Piau are ostensibly distinct in terms of growth and muscularity, commercial breeds are much leaner while local breeds grow much slower and are fat type pigs. However, the genetic factors that underlie such distinctions remain unclear. We used expressed sequence tags (ESTs) to characterize and compare transcript profiles in the semimembranosus muscle of these pig breeds. Our aim was to identify differences in breed-related gene expression that might influence growth performance and meat quality. We constructed three non-normalized cDNA libraries from semimembranosus muscle, using two samples from each one, of these three breeds; 6902 high-quality ESTs were obtained. Cluster analysis was performed and these sequences were clustered into 3670 unique sequences; 24.7% of the sequences were categorized as contigs and 75.3% of the sequences were singletons. Based on homology searches against the SwissProt protein database, we were able to assign a putative protein identity to only 1050 unique sequences. Among these, 58.5% were full-length protein sequences and 17.2% were pig-specific sequences. Muscle structural and cytoskeletal proteins, such as actin, and myosin, were the most abundant transcripts (16.7%) followed by those related to mitochondrial function (12.9%), and ribosomal proteins (12.4%). Furthermore, ESTs generated in this study provide a rich source for identification of novel genes and for the comparative analysis of gene expression patterns in divergent pig breeds.

Single nucleotide variants and InDels identified from whole-genome re-sequencing of Guzerat, Gyr, Girolando and Holstein cattle breeds.

Whole-genome re-sequencing, alignment and annotation analyses were undertaken for 12 sires representing four important cattle breeds in Brazil: Guzerat (multi-purpose), Gyr, Girolando and Holstein (dairy production). A total of approximately 4.3 billion reads from an Illumina HiSeq 2000 sequencer generated for each animal 10.7 to 16.4-fold genome coverage. A total of 27,441,279 single nucleotide variations (SNVs) and 3,828,041 insertions/deletions (InDels) were detected in the samples, of which 2,557,670 SNVs and 883,219 InDels were novel. The submission of these genetic variants to the dbSNP database significantly increased the number of known variants, particularly for the indicine genome. The concordance rate between genotypes obtained using the Bovine HD BeadChip array and the same variants identified by sequencing was about 99.05%. The annotation of variants identified numerous non-synonymous SNVs and frameshift InDels which could affect phenotypic variation. Functional enrichment analysis was performed and revealed that variants in the olfactory transduction pathway was over represented in all four cattle breeds, while the ECM-receptor interaction pathway was over represented in Girolando and Guzerat breeds, the ABC transporters pathway was over represented only in Holstein breed, and the metabolic pathways was over represented only in Gyr breed. The genetic variants discovered here provide a rich resource to help identify potential genomic markers and their associated molecular mechanisms that impact economically important traits for Gyr, Girolando, Guzerat and Holstein breeding programs.

Simultaneous detection of Nosema spp., Ascosphaera apis and Paenibacillus larvae in honey bee products

Honey bees are responsible for pollinating many native and cultivated plant species. These insects can be affected by many pathogens, including fungi and bacteria, both of which can form spores that are easily dispersed within the colony by means of the stored products, among other routes. The objective of this study was to develop a method to detect spores of the honey bee pathogens Nosema apis, Nosema ceranae, Ascosphaera apis and Paenibacillus larvae in samples of honey, bee pollen and royal jelly. The method was standardized for each product individually, and then analyzed by monoplex and multiplex PCR, which showed the same detection thresholds: 1.25 spores/mL of honey for N. ceranae; 7.5 spores/mL of honey for A. apis; and 0.4 spore/mL of honey for P. larvae, respectively. The standardized technique was effective and rapid for the detection of these pathogens in bee products and can be used for the establishment of official methods of sanitary control of bee products, considering the growing national and international trade of these products and the movement or migration of colonies between regions.

Expressed sequenced tags profiling of resistant and susceptible Gyr x Holstein cattle infested with the tick Rhipicephalus (Boophilus) microplus

Tick resistance in cattle is mainly found in zebu (Bos indicus) animals, although it is also present in some taurine (B. taurus) breeds. In order to characterize functional genes involved in tick resistance/susceptibility in cattle, two cDNA libraries were generated using skin tissues of selected Holstein x Gyr animals. A total of 2700 high-quality reads from both resistant and susceptible cDNA were assembled into 458 sequences (contigs) and 834 singletons, with a mean size of 447.7 nucleotides. Assignment of homologous proteins by BLASTX revealed 790 (61.1%) and 300 (23.2%) hits in resistant and susceptible cDNA, respectively; 121 of these hits matched bovine proteins. A total of 502 (38.9%) unique sequences were found to have no significant homology with known sequences and were classified as novel sequences. In general, the most abundant sequences consisted of those coding for hypothetical proteins whose function had not yet been determined, in addition to ribosomal proteins, binding proteins and structural proteins, such as keratin and collagen. The most abundant protein found was collagen type III alpha, although ribosomal proteins accounted for half of the 40 most frequent hits. In addition, five matches within the top 40 best hits corresponded to immune response proteins. These sequences could be used for future studies on functional genomics of cattle tick resistance as well as for genomic sequencing projects.

Whole genome sequencing of Guzera´ cattle reveals genetic variants in candidate genes for production, disease resistance, and heat tolerance

In bovines, artificial selection has produced a large number of breeds which differ in production, environmental adaptation, and health characteristics. To investigate the genetic basis of these phenotypical differences, several bovine breeds have been sequenced. Millions of new SNVs were described at every new breed sequenced, suggesting that every breed should be sequenced. Guzerat or Guzerá is an indicine breed resistant to drought and parasites that has been the base for some important breeds such as Brahman. Here, we describe the sequence of the Guzerá genome and the in silico functional analyses of intragenic breed-specific variations. Mate-paired libraries were generated using the ABI SOLiD system. Sequences were mapped to the Bos taurus reference genome (UMD 3.1) and 87% of the reference genome was covered at a 26X. Among the variants identified, 2,676,067 SNVs and 463,158 INDELs were homozygous, not found in any database searched, and may represent true differences between Guzerá and B. taurus. Functional analyses investigated with the NGS-SNP package focused on 1069 new, non-synonymous SNVs, splice-site variants (including acceptor and donor sites, and the conserved regions at both intron borders, referred to here as splice regions) and coding INDELs (NS/SS/I). These NS/SS/I map to 935 genes belonging to cell communication, environmental adaptation, signal transduction, sensory, and immune systems pathways. These pathways have been involved in phenotypes related to health, adaptation to the environment and behavior, and particularly, disease resistance and heat tolerance. Indeed, 105 of these genes are known QTLs for milk, meat and carcass, production, reproduction, and health traits. Therefore, in addition to describing new genetic variants, our approach provided groundwork for unraveling key candidate genes and mutations.