Brain Gene Regulatory Networks Coordinate Nest Construction in Birds.

Nest building is a vital behavior exhibited during breeding in birds, and is possibly induced by environmental and social cues. Although such behavioral plasticity has been hypothesized to be controlled by adult neuronal plasticity, empirical evidence, especially at the neurogenomic level, remains limited. Here, we aim to uncover the gene regulatory networks that govern avian nest construction and examine whether they are associated with circuit rewiring. We designed an experiment to dissect this complex behavior into components in response to pair bonding and nest material acquisition by manipulating the presence of mates and nest materials in 30 pairs of zebra finches. Whole-transcriptome analysis of 300 samples from five brain regions linked to avian nesting behaviors revealed nesting-associated gene expression enriched with neural rewiring functions, including neurogenesis and neuron projection. The enriched expression was observed in the motor/sensorimotor and social behavior networks of female finches, and in the dopaminergic reward system of males. Female birds exhibited predominant neurotranscriptomic changes to initiate the nesting stage, while males showed major changes after entering this stage, underscoring sex-specific roles in nesting behavior. Notably, major neurotranscriptomic changes occurred during pair bonding, with minor changes during nest material acquisition, emphasizing social interactions in nest construction. We also revealed gene expression associated with reproductive behaviors and tactile sensing for nesting behavior. This study presents novel neurogenomic evidence supporting the hypothesis of adult neural plasticity underlying avian nest-construction behavior. By uncovering the genetic toolkits involved, we offer novel insights into the evolution of animals' innate ability to construct nests.

A chromosome-level genome assembly of the Asian house martin implies potential genes associated with the feathered-foot trait.

The presence of feathers is a vital characteristic among birds, yet most modern birds had no feather on their feet. The discoveries of feathers on the hind limbs of basal birds and dinosaurs have sparked an interest in the evolutionary origin and genetic mechanism of feathered feet. However, the majority of studies investigating the genes associated with this trait focused on domestic populations. Understanding the genetic mechanism underpinned feathered-foot development in wild birds is still in its infancy. Here, we assembled a chromosome-level genome of the Asian house martin (Delichon dasypus) using the long-read High Fidelity sequencing approach to initiate the search for genes associated with its feathered feet. We employed the whole-genome alignment of D. dasypus with other swallow species to identify high-SNP regions and chromosomal inversions in the D. dasypus genome. After filtering out variations unrelated to D. dasypus evolution, we found six genes related to feather development near the high-SNP regions. We also detected three feather development genes in chromosomal inversions between the Asian house martin and the barn swallow genomes. We discussed their association with the wingless/integrated (WNT), bone morphogenetic protein, and fibroblast growth factor pathways and their potential roles in feathered-foot development. Future studies are encouraged to utilize the D. dasypus genome to explore the evolutionary process of the feathered-foot trait in avian species. This endeavor will shed light on the evolutionary path of feathers in birds.

Genomic architecture underlying morphological and physiological adaptation to high elevation in a songbird.

Organisms often acquire physiological and morphological modifications to conquer ecological challenges when colonizing new environments which lead to their adaptive evolution. However, deciphering the genomic mechanism of ecological adaptation is difficult because ecological environments are often too complex for straightforward interpretation. Thus, we examined the adaptation of a widespread songbird-the rufous-capped babbler (Cyanoderma ruficeps)-to a relatively simple system: distinct environments across elevational gradients on the mountainous island of Taiwan. We focused on the genomic sequences of 43 birds from five populations to show that the Taiwan group split from its sister group in mainland China around 1-2 million years ago (Ma) and colonized the montane habitats of Taiwan at least twice around 0.03-0.22 Ma. The montane and lowland Taiwan populations diverged with gene flow between them, suggesting strong selection associated with different elevations. We found that the montane babblers had smaller beaks than the lowland ones, consistent with Allen's rule, and identified candidate genes-COL9A1 and SOX11-underlying the beak size changes. We also found that altitudinally divergent mutations were mostly located in noncoding regions and tended to accumulate in chromosomal inversions and autosomes. The altitudinally divergent mutations might regulate genes related to haematopoietic, metabolic, immune, auditory and vision functions, as well as cerebrum morphology and plumage development. The results reveal the genomic bases of morphological and physiological adaptation in this species to the low temperature, hypoxia and high UV light environment at high elevation. These findings improve our understanding of how ecological adaptation drives population divergence from the perspective of genomic architecture.

Microbowls with Controlled Concavity for Accurate Microscale Mass Spectrometry.

Patterned surfaces can enhance the sensitivity of laser desorption ionization mass spectrometry by segregating and concentrating analytes, but their fabrication can be challenging. Here, a simple method to fabricate substrates patterned with micrometer-scale wells that yield more accurate and sensitive mass spectrometry measurements compared to flat surfaces is described. The wells can also concentrate and localize cells and beads for cell-based assays.

Domestication obscures genomic estimates of population history.

Domesticated species are valuable models to examine phenotypic evolution, and knowledge on domestication history is critical for understanding the trajectories of evolutionary changes. Sequentially Markov Coalescent models are often used to infer domestication history. However, domestication practices may obscure the signal left by population history, affecting demographic inference. Here we assembled the genomes of a recently domesticated species-the society finch-and its parent species-the white-rumped munia-to examine its domestication history. We applied genomic analyses to two society finch breeds and white-rumped munias to test whether domestication of the former resulted from inbreeding or hybridization. The society finch showed longer and more runs of homozygosity and lower genomic heterozygosity than the white-rumped munia, supporting an inbreeding origin in the former. Blocks of white-rumped munia and other ancestry in society finch genomes showed similar genetic distance between the two taxa, inconsistent with the hybridization origin hypothesis. We then applied two Sequentially Markov Coalescent models-psmc and smc++-to infer the demographic histories of both. Surprisingly, the two models did not reveal a recent population bottleneck, but instead the psmc model showed a specious, dramatic population increase in the society finch. Subsequently, we used simulated genomes based on an array of demographic scenarios to demonstrate that recent inbreeding, not hybridization, caused the distorted psmc population trajectory. Such analyses could have misled our understanding of the domestication process. Our findings stress caution when interpreting the histories of recently domesticated species inferred by psmc, arguing that these histories require multiple analyses to validate.

To Trim or Not to Trim: Effects of Read Trimming on the De Novo Genome Assembly of a Widespread East Asian Passerine, the Rufous-Capped Babbler (Cyanoderma ruficeps Blyth).

Trimming low quality bases from sequencing reads is considered as routine procedure for genome assembly; however, we know little about its pros and cons. Here, we used empirical data to examine how read trimming affects assembled genome quality and computational time for a widespread East Asian passerine, the rufous-capped babbler (Cyanoderma ruficeps Blyth). We found that scaffolds assembled from raw reads were always longer than those from trimmed ones, whereas computational times for the former were sometimes much longer than the latter. Nevertheless, assembly completeness showed little difference among the trimming strategies. One should determine the optimal trimming strategy based on what the assembled genome will be used for. For example, to identify single nucleotide polymorphisms (SNPs) associated with phenotypic evolution, applying PLATANUS to gently trim reads would yield a reference genome with a slightly shorter scaffold length (N50 = 15.64 vs. 16.89 Mb) than the raw reads, but would save 75% of computational time. We also found that chromosomes Z, W, and 4A of the rufous-capped babbler were poorly assembled, likely due to a recently fused, neo-sex chromosome. The rufous-capped babbler genome with long scaffolds and quality gene annotation can provide a good system to study avian ecological adaptation in East Asia.

Historical demography of four gecko species specializing in boulder cave habitat: Implications in the evolutionary dead end hypothesis and conservation.

Specialization in narrow ecological niches may not only help species to survive in competitive or unique environments but also contribute to their extermination over evolutionary time. Although the "evolutionary dead end" hypothesis has long been debated, empirical evidence from species with detailed information on niche specialization and evolutionary history remains rare. Here we use a group of four closely related Cnemaspis gecko species that depend highly on granite boulder caves in the Mekong Delta to investigate the potential impact of ecological specialization on their evolution and population dynamics. Isolated by unsuitable floodplain habitats, these boulder-dwelling geckos are among the most narrowly distributed Squamata in the world. We applied several coalescence-based approaches combined with the RAD-seq technique to estimate their divergence times, gene flow and demographic fluctuations during the speciation and population differentiation processes. Our results reveal long-term population shrinkage in the four geckos and limited gene flow during their divergence. The results suggest that the erosion and fragmentation of the granite boulder hills have greatly impacted population divergence and declines. The habitat specialization of these geckos has led to fine-scaled speciation in these granite rocky hills; in contrast, specialization might also have pushed these species toward the edge of extinction. Our study also emphasizes the conservation urgency of these vulnerable, cave-dependent geckos.

Whole-genome sequencing and comparative analysis of two plant-associated strains of Rhodopseudomonas palustris (PS3 and YSC3).

Rhodopseudomonas palustris strains PS3 and YSC3 are purple non-sulfur phototrophic bacteria isolated from Taiwanese paddy soils. PS3 has beneficial effects on plant growth and enhances the uptake efficiency of applied fertilizer nutrients. In contrast, YSC3 has no significant effect on plant growth. The genomic structures of PS3 and YSC3 are similar; each contains one circular chromosome that is 5,269,926 or 5,371,816 bp in size, with 4,799 or 4,907 protein-coding genes, respectively. In this study, a large class of genes involved in chemotaxis and motility was identified in both strains, and genes associated with plant growth promotion, such as nitrogen fixation-, IAA synthesis- and ACC deamination-associated genes, were also identified. We noticed that the growth rate, the amount of biofilm formation, and the relative expression levels of several chemotaxis-associated genes were significantly higher for PS3 than for YSC3 upon treatment with root exudates. These results indicate that PS3 responds better to the presence of plant hosts, which may contribute to the successful interactions of PS3 with plant hosts. Moreover, these findings indicate that the existence of gene clusters associated with plant growth promotion is required but not sufficient for a bacterium to exhibit phenotypes associated with plant growth promotion.

The Rho signalling pathway mediates the pathogenicity of AHPND-causing V. parahaemolyticus in shrimp.

An emerging bacterial disease, acute hepatopancreatic necrosis disease (AHPND), is caused by strains of Vibrio parahaemolyticus with an additional AHPND-associated plasmid pVA1 encoding a virulent toxin (Pirvp ) that damages the shrimp's hepatopancreas. Like other species of Vibrio, these virulent strains initially colonise the shrimp's stomach, but it is not yet understood how the bacteria or toxins are subsequently able to cross the epithelial barrier and reach the hepatopancreas. Here, by using transcriptomics and system biology methods, we investigate AHPND-induced changes in the stomach of AHPND-causing V. parahaemolyticus (5HP)-infected shrimp and identify host molecular mechanisms that might explain how the integrity of the stomach barrier is compromised. We found that the expression of 376 unique genes was differentially regulated by AHPND infection. Gene ontology, protein interaction, and gene-to-gene correlation expression interaction analyses indicated that in addition to the immune system, a number of these genes were involved in cytoskeleton regulation by Rho GTPase. The involvement of Rho pathway regulation during AHPND pathogenesis was further supported by experiments showing that while Rho inhibitor pretreatment delayed the infection, pretreatment with Rho activator enhanced the pathogenicity of 5HP, and both the bacteria and toxin were detected sooner in the hepatopancreas. Further, disruption of the stomach epithelial structure was found in both Rho preactivated shrimp and in 5HP-infected shrimp. Taken together, we interpret our results to mean that Rho signalling helps to mediate AHPND pathogenesis in shrimp.

Dissemination of carbapenem-resistant Klebsiella pneumoniae harboring KPC-carrying plasmid pKPC_P16, a pKPC_LK30 variant, in northern Taiwan.

The prevalence of carbapenem-resistant Klebsiella pneumoniae (CRKP) was up to 30% between 2014 and 2016 in the study hospital. Of these 77 CRKP isolates, 22 isolates with sequence type ST11 carried the new pKPC_P16 plasmid, a pKPC_LK30 variant, and were widely disseminated between 2014 and 2015 in northern Taiwan.

Revealing the Saline Adaptation Strategies of the Halophilic Bacterium Halomonas beimenensis through High-throughput Omics and Transposon Mutagenesis Approaches.

Studies on the halotolerance of bacteria are attractive to the fermentation industry. However, a lack of sufficient genomic information has precluded an investigation of the halotolerance of Halomonas beimenensis. Here, we describe the molecular mechanisms of saline adaptation in H. beimenensis based on high-throughput omics and Tn5 transposon mutagenesis. The H. beimenensis genome is 4.05 Mbp and contains 3,807 genes, which were sequenced using short and long reads obtained via deep sequencing. Sixteen Tn5 mutants with a loss of halotolerance were identified. Orthologs of the mutated genes, such as nqrA, trkA, atpC, nadA, and gdhB, have significant biological functions in sodium efflux, potassium uptake, hydrogen ion transport for energy conversion, and compatible solute synthesis, which are known to control halotolerance. Other genes, such as spoT, prkA, mtnN, rsbV, lon, smpB, rfbC, rfbP, tatB, acrR1, and lacA, function in cellular signaling, quorum sensing, transcription/translation, and cell motility also shown critical functions for promoting a halotolerance. In addition, KCl application increased halotolerance and potassium-dependent cell motility in a high-salinity environment. Our results demonstrated that a combination of omics and mutagenesis could be used to facilitate the mechanistic exploitation of saline adaptation in H. beimenensis, which can be applied for biotechnological purposes.

Clonal spread of blaOXA-72-carrying Acinetobacter baumannii sequence type 512 in Taiwan.

This is the first report to show an insidious outbreak of armA- and blaOXA-72-carrying Acinetobacter baumannii sequence type 512 (ST512) at a study hospital in northern Taiwan. Multilocus sequence typing revealed that this was a ST512 clone. All of the isolates with ST512 carried a novel 12,056-bp repGR2 in combination with a repGR12-type plasmid. This plasmid, designated pAB-ML, had one copy of the blaOXA-72 gene that was flanked by XerC/XerD-like sites and conferred resistance to carbapenems.

Identification of miRNAs and Their Targets in the Liverwort Marchantia polymorpha by Integrating RNA-Seq and Degradome Analyses.

Bryophytes (liverworts, hornworts and mosses) comprise the three earliest diverging lineages of land plants (embryophytes). Marchantia polymorpha, a complex thalloid Marchantiopsida liverwort that has been developed into a model genetic system, occupies a key phylogenetic position. Therefore, M. polymorpha is useful in studies aiming to elucidate the evolution of gene regulation mechanisms in plants. In this study, we used computational, transcriptomic, small RNA and degradome analyses to characterize microRNA (miRNA)-mediated pathways of gene regulation in M. polymorpha. The data have been integrated into the open access ContigViews-miRNA platform for further reference. In addition to core components of the miRNA pathway, 129 unique miRNA sequences, 11 of which could be classified into seven miRNA families that are conserved in embryophytes (miR166a, miR390, miR529c, miR171-3p, miR408a, miR160 and miR319a), were identified. A combination of computational and degradome analyses allowed us to identify and experimentally validate 249 targets. In some cases, the target genes are orthologous to those of other embryophytes, but in other cases, the conserved miRNAs target either paralogs or members of different gene families. In addition, the newly discovered Mpo-miR11707.1 and Mpo-miR11707.2 are generated from a common precursor and target MpARGONAUTE1 (LW1759). Two other newly discovered miRNAs, Mpo-miR11687.1 and Mpo-miR11681.1, target the MADS-box transcription factors MpMADS1 and MpMADS2, respectively. Interestingly, one of the pentatricopeptide repeat (PPR) gene family members, MpPPR_66 (LW9825), the protein products of which are generally involved in various steps of RNA metabolism, has a long stem-loop transcript that can generate Mpo-miR11692.1 to autoregulate MpPPR_66 (LW9825) mRNA. This study provides a foundation for further investigations of the RNA-mediated silencing mechanism in M. polymorpha as well as of the evolution of this gene silencing pathway in embryophytes.

Isolation and Identification of Post-Transcriptional Gene Silencing-Related Micro-RNAs by Functionalized Silicon Nanowire Field-effect Transistor.

Many transcribed RNAs are non-coding RNAs, including microRNAs (miRNAs), which bind to complementary sequences on messenger RNAs to regulate the translation efficacy. Therefore, identifying the miRNAs expressed in cells/organisms aids in understanding genetic control in cells/organisms. In this report, we determined the binding of oligonucleotides to a receptor-modified silicon nanowire field-effect transistor (SiNW-FET) by monitoring the changes in conductance of the SiNW-FET. We first modified a SiNW-FET with a DNA probe to directly and selectively detect the complementary miRNA in cell lysates. This SiNW-FET device has 7-fold higher sensitivity than reverse transcription-quantitative polymerase chain reaction in detecting the corresponding miRNA. Next, we anchored viral p19 proteins, which bind the double-strand small RNAs (ds-sRNAs), on the SiNW-FET. By perfusing the device with synthesized ds-sRNAs of different pairing statuses, the dissociation constants revealed that the nucleotides at the 3'-overhangs and pairings at the terminus are important for the interactions. After perfusing the total RNA mixture extracted from Nicotiana benthamiana across the device, this device could enrich the ds-sRNAs for sequence analysis. Finally, this bionanoelectronic SiNW-FET, which is able to isolate and identify the interacting protein-RNA, adds an additional tool in genomic technology for the future study of direct biomolecular interactions.

Dnmt3l-knockout donor cells improve somatic cell nuclear transfer reprogramming efficiency.

Nuclear transfer (NT) is a technique used to investigate the development and reprogramming potential of a single cell. DNA methyltransferase-3-like, which has been characterized as a repressive transcriptional regulator, is expressed in naturally fertilized egg and morula/blastocyst at pre-implantation stages. In this study, we demonstrate that the use of Dnmt3l-knockout (Dnmt3l-KO) donor cells in combination with Trichostatin A treatment improved the developmental efficiency and quality of the cloned embryos. Compared with the WT group, Dnmt3l-KO donor cell-derived cloned embryos exhibited increased cell numbers as well as restricted OCT4 expression in the inner cell mass (ICM) and silencing of transposable elements at the blastocyst stage. In addition, our results indicate that zygotic Dnmt3l is dispensable for cloned embryo development at pre-implantation stages. In Dnmt3l-KO mouse embryonic fibroblasts, we observed reduced nuclear localization of HDAC1, increased levels of the active histone mark H3K27ac and decreased accumulation of the repressive histone marks H3K27me3 and H3K9me3, suggesting that Dnmt3l-KO donor cells may offer a more permissive epigenetic state that is beneficial for NT reprogramming.

Transcriptome profiling in imipenem-selected Acinetobacter baumannii.

BACKGROUND: Carbapenem-resistance in Acinetobacter baumannii has gradually become a global challenge. To identify the genes involved in carbapenem resistance in A. baumannii, the transcriptomic responses of the completely sequenced strain ATCC 17978 selected with 0.5 mg/L (IPM-2 m) and 2 mg/L (IPM-8 m) imipenem were investigated using RNA-sequencing to identify differences in the gene expression patterns. RESULTS: A total of 88 and 68 genes were differentially expressed in response to IPM-2 m and IPM-8 m selection, respectively. Among the expressed genes, 50 genes were highly expressed in IPM-2 m, 30 genes were highly expressed in IPM-8 m, and 38 genes were expressed common in both strains. Six groups of genes were simultaneously expressed in IPM-2 m and IPM-8 m mutants. The three gene groups involved in DNA recombination were up-regulated, including recombinase, transposase and DNA repair, and beta-lactamase OXA-95 and homologous recombination. The remaining gene groups involved in biofilm formation were down-regulated, including quorum sensing, secretion systems, and the csu operon. The antibiotic resistance determinants, including RND efflux transporters and multidrug resistance pumps, were over-expressed in response to IPM-2 m selection, followed by a decrease in response to IPM-8 m selection. Among the genes over-expressed in both strains, blaOXA-95, previously clustered with the blaOXA-51-like family, showed 14-fold (IPM-2 m) to 330-fold (IPM-8 m) over-expression. The expression of blaOXA-95 in IPM-2 m and IPM-8 m cells was positively correlated with the rate of imipenem hydrolysis, as demonstrated through Liquid Chromatography-Mass Spectrometry/Mass Spectrometry, suggesting that blaOXA-95 plays a critical role in conferring carbapenem resistance. In addition, A. baumannii shows an inverse relationship between carbapenem resistance and biofilm production. CONCLUSION: Gene recombination and blaOXA-95 play critical roles in carbapenem resistance in A. baumannii. Taken together, the results of the present study provide a foundation for future studies of the network systems associated with carbapenem resistance.

Insertion sequence transposition determines imipenem resistance in Acinetobacter baumannii.

This study employed genomewide analysis to investigate potential resistance mechanisms in Acinetobacter baumannii following imipenem exposure. Imipenem-selected mutants were generated from the imipenem-susceptible strain ATCC 17978 by multistep selection resistance. Antibiotic susceptibilities were examined, and the selected mutants originated from the ATCC 17978 strain were confirmed by pulsed-field gel electrophoresis. The genomic sequence of a resistant mutant was analyzed using a next-generation sequencing platform, and genetic recombination was further confirmed by PCR. The result showed that phenotypic resistance was observed with carbapenem upon exposure to various concentrations of imipenem. Genomewide analysis showed that ISAba1 transposition was initiated by imipenem exposure at concentrations up to 0.5 mg/L. Transposition of ISAba1 upstream of blaOXA-95 was detected in all the selected mutants. The expression of blaOXA-95 was further analyzed by quantitative PCR, and the results demonstrated that a 200-fold increase in gene expression was required for resistance to imipenem. This study concluded that imipenem exposure at a concentration of 0.5 mg/L mediated the transposition of ISAba1 upstream of the blaOXA-95 gene and resulted in the overexpression of blaOXA-95 gene, which may play a major role in the resistance to imipenem in A. baumannii.

The origin of Acinetobacter baumannii TYTH-1: a comparative genomics study.

There have been increasing reports of bla(OXA-23)-carrying strains of carbapenem-resistant Acinetobacter baumannii (CRAB), which has become a significant public health concern in Taiwan. To determine the origin of these CRAB strains, the prevalence of CRAB and bla(OXA-23)-carrying CRAB in a regional hospital was analysed retrospectively. The genome of A. baumannii TYTH-1 was completely sequenced and annotated. Multiple comparative genomics studies, including phylogenetic analysis, functional comparison via the Clusters of Orthologous Groups (COGs) database, and determination of variance in GC profiles in the whole genome and gene arrangements in resistance islands, were performed using 11 completely sequenced A. baumannii genomes. bla(OXA-23)-carrying CRAB isolates became dominant clones in 2007. A comparative genomics analysis revealed a common strain lineage between Taiwanese strains (TYTH-1 and TCDC-AB0715) and Chinese strains (MDR-TJ and MDR-ZJ06). Phylogenetic studies and GC profiles showed that the genome of TYTH-1 was closest to MDR-ZJ06. However, the resistance island of TYTH-1 (RI(TYTH-1)) was nearly identical to that of RI(MDT-TJ). The functional category for COGs was similar in the tested genomes. The results reveal that dissemination of bla(OXA-23)-carrying CRAB in Taiwan may have been mediated by the transfer of people between Taiwan and China during 2007. The global spread of CRAB is now a worldwide public health problem. In Taiwan, the government needs to focus more attention on the importance of identifying and tracing resistant pathogens and issuing notifications of CRAB infections.

Genome sequence of Acinetobacter baumannii TYTH-1.

Acinetobacter baumannii has emerged recently as a major cause of health care-associated infections due to the extent of its antimicrobial resistance and its propensity to cause large nosocomial outbreaks. Here we report the genome sequence of Acinetobacter baumannii TYTH-1 isolated in Taiwan during 2008.

Sequence of Leptospira santarosai serovar Shermani genome and prediction of virulence-associated genes.

Leptospirosis, a widespread zoonosis, is a re-emerging infectious disease caused by pathogenic Leptospira species. In Taiwan, Leptospira santarosai serovar Shermani is the most frequently isolated serovar, causing both renal and systemic infections. This study aimed to generate a L. santarosai serovar Shermani genome sequence and categorize its hypothetical genes, particularly those associated with virulence. The genome sequence consists of 3,936,333 nucleotides and 4033 predicted genes. Additionally, 2244 coding sequences could be placed into clusters of orthologous groups and the number of genes involving cell wall/membrane/envelope biogenesis and defense mechanisms was higher than that of other Leptospira spp. Comparative genetic analysis based on BLASTX data revealed that about 73% and 68.8% of all coding sequences have matches to pathogenic L. interrogans and L. borgpetersenii, respectively, and about 57.6% to saprophyte L. biflexa. Among the hypothetical proteins, 421 have a transmembrane region, 172 have a signal peptide and 17 possess a lipoprotein signature. According to PFAM prediction, 32 hypothetical proteins have properties of toxins and surface proteins mediated bacterial attachment, suggesting they may have roles associated with virulence. The availability of the genome sequence of L. santarosai serovar Shermani and the bioinformatics re-annotation of leptospiral hypothetical proteins will facilitate further functional genomic studies to elucidate the pathogenesis of leptospirosis and develop leptospiral vaccines.