Analysis of the Genome and Transcriptome of var. Reveals Complex RNA Expression and Microevolution Leading to Virulence Attenuation

1. BrownGD, DenningDW, LevitzSM (2012) Tackling human fungal infections. Science 336: 647.

2. ParkBJ, WannemuehlerKA, MarstonBJ, GovenderN, PappasPG, et al. (2009) Estimation of the current global burden of cryptococcal meningitis among persons living with HIV/AIDS. AIDS 23: 525–530.

3. Kwon-Chung KJ, Boekhout T, Wickes BL, Fell JW, Systematics of the genus Cryptococcus and its type species C. neoformans, in Cryptococcus from human pathogen to model yeast, J. Heitman, et al., Editors. 2011, ASM Press: Washington DC. p.3–15.

4. HeitmanJ, CasadevallA, LodgeJK, PerfectJR (1999) The Cryptococcus neoformans genome sequencing project. Mycopathologia 148: 1–7.

5. LoftusB, FungE, RoncagliaP, RowleyD, AmedeoP, et al. (2005) The genome and transcriptome of Cryptococcus neoformans, a basidiomycetous fungal pathogen of humans. Science 307: 1321–1324.

6. D'SouzaCA, KronstadJW, TaylorG, WarrenR, YuenM, et al. (2011) Genome variation in Cryptococcus gattii, an emerging pathogen of immunocompetent hosts. mBio 2: e00342–10.

7. PerfectJR (2006) Cryptococcus neoformans: the yeast that likes it hot. FEMS Yeast Res 6: 463–468.

8. Casadevall A, Perfect JR, Cryptococcus neoformans. 1998, Washington, D.C.: American Society for Microbiology Press.

9. IdnurmA, BahnYS, NielsenK, LinX, FraserJA, et al. (2005) Deciphering the model pathogenic fungus Cryptococcus neoformans. Nat Rev Microbiol 3: 753–764.

10. SunS, XuJ (2009) Chromosomal rearrangements between serotype A and D strains in Cryptococcus neoformans. PLoS One 4: e5524.

11. KavanaughLA, FraserJA, DietrichFS (2006) Recent evolution of the human pathogen Cryptococcus neoformans by intervarietal transfer of a 14-gene fragment. Mol Biol Evol 23: 1879–1890.

12. LangmeadB, TrapnellC, PopM, SalzbergSL (2009) Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol 10: R25.

13. TrapnellC, PachterL, SalzbergSL (2009) TopHat: discovering splice junctions with RNA-Seq. Bioinformatics 25: 1105–1111.

14. PaoSS, PaulsenIT, SaierMJ (1988) Major facilitator superfamily. Microbiol Mol Biol Rev 62: 1–34.

15. RepM (2005) Small proteins of plant-pathogenic fungi secreted during host colonization. FEMS Microbiol Lett 253: 19–27.

16. MagnyEG, PueyoJI, PearlFMG, CespedesMA, NivenJE, et al. (2013) Conserved regulation of cardiac calcium uptake by peptides encoded in small Open Reading Frames. Science 341: 1116–1120.

17. IngoliaNT, GhaemmaghamiS, NewmanJRS, WeissmanJS (2009) Genome-wide analysis in vivo of translation with nucleotide resolution using ribosome profiling. Science 324: 218–223.

18. ChackoN, LinX (2013) Non-coding RNAs in the development and pathogenesis of eukaryotic microbes. Appl Microbiol Biotechnol 97: 7989–7997.

19. DavidL, HuberW, GranovskaiaM, ToedlingJ, PalmCJ, et al. (2006) A high-resolution map of transcription in the yeast genome. Proc Natl Acad Sci USA 103: 5320–5325.

20. YangL, DuffM, GraveleyB, CarmichaelG, ChenL-L (2011) Genomewide characterization of non-polyadenylated RNAs. Genome Biol 12: R16.

21. CsurosM, RogozinIB, KooninEV (2011) A detailed history of intron-rich eukaryotic ancestors inferred from a global survey of 100 complete genomes. PLoS Comput Biol 9: e1002150.

22. GoebelsC, ThonnA, Gonzalez-HilarionS, RollandO, MoyrandF, et al. (2013) Introns regulate gene expression in Cryptococcus neoformans in a Pab2p dependent pathway. PLoS Genet 9: e1003686.

23. StajichJE, DietrichFS, RoySW (2007) Comparative genomic analysis of fungal genomes reveals intron-rich ancestors. Genome Biol 8: R223.

24. BahnYS, KojimaK, CoxGM, HeitmanJ (2006) A unique fungal two-component system regulates stress responses, drug sensitivity, sexual development, and virulence of Cryptococcus neoformans. Mol Biol Cell 17: 3122–3125.

25. SchwartzSH, SilvaJ, BursteinD, PupkoT, EyrasE, et al. (2008) Large-scale comparative analysis of splicing signals and their corresponding splicing factors in eukaryotes. Genome Res 18: 88–103.

26. HonCC, WeberC, SismeiroO, ProuxC, KouteroM, et al. (2013) Quantification of stochastic noise of splicing and polyadenylation in Entamoeba histolytica. Nucleic Acids Res 41: 1936–1952.

27. PauwsE, van KampenAH, van de GraafSA, de VijlderJJ, Ris-StalpersC (2001) Heterogeneity in polyadenylation cleavage sites in mammalian mRNA sequences: implications for SAGE analysis. Nucleic Acids Res 29: 1690–1694.

28. Di GiammartinoDC, NishidaK, ManleyJL (2011) Mechanisms and consequences of alternative polyadenylation. Mol Cell 43: 853–866.

29. WuX, LiuM, DownieB, LiangC, JiG, et al. (2011) Genome-wide landscape of polyadenylation in Arabidopsis provides evidence for extensive alternative polyadenylation. Proc Natl Acad Sci USA 108: 12533–12538.

30. MoqtaderiZ, GeisbergJV, JinY, FanX, StruhlK (2013) Species-specific factors mediate extensive heterogeneity of mRNA 3′ ends in yeasts. Proc Natl Acad Sci USA 110: 11073–11078.

31. JiZ, LeeJY, PanZ, JiangB, TianB (2009) Progressive lengthening of 3′ untranslated regions of mRNAs by alternative polyadenylation during mouse embryonic development. Proc Natl Acad Sci USA 106: 7028–7033.

32. ElkonR, UgaldeAP, AgamiR (2013) Alternative cleavage and polyadenylation: extent, regulation and function. Nat Rev Genet 14: 496–506.

33. BaileyTL (2011) DREME: motif discovery in transcription factor ChIP-seq data. Bioinformatics 27: 1653–1659.

34. ZhaoJ, HymanL, MooreC (1999) Formation of mRNA 3′ ends in eukaryotes: mechanism, regulation, and interrelationships with other steps in mRNA synthesis. Microbiol Mol Biol Rev 403: 405–445.

35. WilkeningS, PelechanoV, JärvelinAI, TekkedilMM, AndersS, et al. (2013) An efficient method for genome-wide polyadenylation site mapping and RNA quantification. Nucleic Acids Res 41: e65.

36. PrescottEM, ProudfootNJ (2002) Transcriptional collision between convergent genes in budding yeast. Proc Natl Acad Sci USA 99: 8796–8801.

37. HobsonDJ, WeiW, SteinmetzLM, SvejstrupJQ (2012) RNA Polymerase II collision interrupts convergent transcription. Mol Cell 48: 365–374.

38. DonaldsonME, SavilleBJ (2012) Natural antisense transcripts in fungi. Mol Microbiol 85: 405–417.

39. van DijkEL, ChenCL, d'Aubenton-CarafaY, GourvennecS, KwapiszM, et al. (2011) XUTs are a class of Xrn1-sensitive antisense regulatory non-coding RNA in yeast. Nature 475: 114–119.

40. JanbonG, MaengSY, D-F , KoY-J, JungK-W, MoyrandF, et al. (2010) Characterizing the role of RNA silencing components in Cryptococcus neoformans. Fungal Genet Biol 47: 1070–1080.

41. WangX, HsuehYP, LiW, FloydA, SkalskyR, et al. (2010) Sex-induced silencing defends the genome of Cryptococcus neoformans via RNAi. Genes Dev 24: 2566–2582.

42. DumesicPA, NatarajanP, ChenC, DrinnenbergIA, SchillerBJ, et al. (2013) Stalled spliceosomes are a signal for RNAi-mediated genome defense. Cell 152: 957–968.

43. WangX, DarwicheS, HeitmanJ (2013) Sex-induced silencing operates during opposite-sex and unisexual reproduction in Cryptococcus neoformans. Genetics 193: 1163–1174.

44. HouseleyJ, RubbiL, GrunsteinM, TollerveyD, VogelauerM (2008) A ncRNA modulates histone modification and mRNA induction in the yeast GAL gene cluster. Mol Cell 32: 685–695.

45. MoyrandF, LafontaineI, FontaineT, JanbonG (2008) UGE1 and UGE2 regulate the UDP-glucose/UDP-galactose equilibrium in Cryptococcus neoformans. Eukaryot Cell 7: 2069–2077.

46. HenikoffJG, AhmadK, MalikHS (2001) The centromere paradox: stable inheritance with rapidly evolving DNA. Science 293: 1098–1102.

47. RhindN, ChenZ, YassourM, ThompsonDA, HaasBJ, et al. (2011) Comparative functional genomics of the fission yeasts. Science 332: 930–936.

48. RoyB, SanyalK (2011) Diversity in requirement of genetic and epigenetic factors for centromere function in fungi. Eukaryot Cell 10: 1384–1395.

49. PadmanabhanS, ThakurJ, SiddharthanR, SanyalaK (2008) Rapid evolution of Cse4p-rich centromeric DNA sequences in closely related pathogenic yeasts, Candida albicans and Candida dubliniensis. Proc Natl Acad Sci USA 105: 19597–19802.

50. StajichJE, HarrisT, BrunkBP, BrestelliJ, FischerS, et al. (2012) FungiDB: an integrated functional genomics database for fungi. Nucleic Acids Res 40: D675–681.

51. KozubowskiL, YadavV, ChatterjeeG, SridharS, YamaguchiM, et al. (2013) Ordered kinetochore assembly in the human-pathogenic basidiomycetous yeast Cryptococcus neoformans. MBio 4: e00614–13.

52. VollrathD, DavisRW, ConnellyC, HieterP (1988) Physical mapping of large DNA by chromosome fragmentation. Proc Natl Acad Sci USA 85: 6027–6031.

53. IdnurmA (2010) A tetrad analysis of the basidiomycete fungus Cryptococcus neoformans. Genetics 185: 153–163.

54. ShiJ, WolfSE, BurkeJM, PrestingGG, Ross-IbarraJ, et al. (2010) Widespread gene conversion in centromere cores. PLoS Biol 8: e1000327.

55. ThakurJ, SanyalK (2013) Efficient neocentromere formation is suppressed by gene conversion to maintain centromere function at native physical chromosomal loci in Candida albicans. Genome Res 23: 638–652.

56. FolcoHD, PidouxAL, UranoT, AllshireRC (2008) Heterochromatin and RNAi are required to establish CENP-A chromatin at centromeres. Science 319: 94–97.

57. VolpeTA, KidnerC, HallIM, TengG, GrewalSI, et al. (2002) Regulation of heterochromatic silencing and histone H3 lysine-9 methylation by RNAi. Science 297: 1833–1837.

58. NiL, FeretzakiM, LiW, Floyd-AveretteA, MieczkowskiP, et al. (2013) Unisexual and heterosexual meiotic reproduction generates phenotypic and genotypic diversity de novo in Cryptococcus neoformans involving aneuploidy. PLoS Biol 11: e1001653.

59. MondonP, ChangYC, VarmaA, Kwon-ChungKJ (2000) A novel episomal shuttle vector for transformation of Cryptococcus neoformans with the ccdB gene as a positive selection marker in bacteria. FEMS Microbiol Lett 187: 41–45.

60. EdmanJC (1992) Isolation of telomerelike sequences from Cryptococcus neoformans and their use in high efficiency transformation. Mol Cell Biol 12: 2777–2783.

61. MesnerLD, CrawfordEL, HamlinJL (2006) Isolating apparently pure libraries of replication origins from complex genomes. Mol Cell 5: 719–726.

62. PoloumienkoA, DershowitzA, DeJ, NewlonCS (2001) Completion of replication map of Saccharomyces cerevisiae chromosome III. Mol Biol Cell 12: 3317–3327.

63. PatelPK, ArcangioliB, BakerSP, BensimonA, RhindN (2006) DNA replication origins fire stochastically in fission yeast. Mol Biol Cell 17: 308–316.

64. WoodsJP, GoldmanWE (1993) Autonomous replication of foreign DNA in Histoplasma capsulatum. J Bacteriol 175: 636–641.

65. LabibK (2011) Building a double hexamer of DNA helicase at eukaryotic replication origins. EMBO J 30: 4853–4855.

66. AvesSJ, LiuY, RichardsTA (2012) Evolutionary diversification of eukaryotic replication machinery. Subcell Biochem 62: 19–35.

67. LiH, StillmanB (2012) The origin recognition complex: a biochemical and structural view. Subcell Biochem 62: 37–58.

68. NielsenK, CoxGM, WangP, ToffalettiDL, PerfectJR, et al. (2003) Sexual cycle of Cryptococcus neoformans var. grubii and virulence of congenic a and alpha isolates. Infect Immun 71: 4831–4841.

69. FranzotSP, MukherjeeJ, CherniakR, ChenL, HamdanJS, et al. (1998) Microevolution of a standard strain of Cryptococcus neoformans resulting in differences in virulence and other phenotypes. Infect Immun 66: 89–97.

70. LiuOW, ChunCD, ChowED, ChenC, MadhaniHD, et al. (2008) Systematic genetic analysis of virulence in the human fungal pathogen Cryptococcus neoformans. Cell 135: 174–188.

71. OrmerodKL, MorrowCA, ChowEWL, LeeIR, ArrasSDM, et al. (2013) Comparative genomics of serial isolates of Cryptococcus neoformans reveals gene associated with carbon utilisation and virulence. G3 3: 675–686.

72. HassBJ, ZengQ, PearsonMD, CuomoCA, WortmanJR (2011) Approaches to fungal genome annotation. Mycology 3: 118–141.

73. FoissacS, GouzyJ, RombautsS, MathéC, AmselemJ, et al. (2008) Genome annotation in plants and fungi: EuGène as a model platform. Current Bioinformatics 3: 87–97.

74. WilhelmBT, MargueratS, WattS, SchubertF, WoodV, et al. (2008) Dynamic repertoire of a eukaryotic transcriptome surveyed at single-nucleotide resolution. Nature 453: 1239–1243.

75. BachmannSP, VandeWalleK, RamageG, PattersonTF, WickesBL, et al. (2002) In vitro activity of caspofungin against Candida albicans biofilms. Antimicrob Agents Chemother 46: 3591–3596.

76. DelmasS, PullanST, GaddipatiS, KokolskiM, MallaS, et al. (2012) Uncovering the genome-wide transcriptional responses of the filamentous fungus Aspergillus niger to lignocellulose using RNA sequencing. PLoS Genet 8: e1002875.

77. MartinezD, LarrondoLF, PutnamN, GelpkeMD, HuangK, et al. (2004) Genome sequence of the lignocellulose degrading fungus Phanerochaete chrysosporium strain RP78. Nat Biotechnol 22: 695–700.

78. XuJ, SaundersCW, HuP, GrantRA, BoekhoutT, et al. (2007) Dandruff-associated Malassezia genomes reveal convergent and divergent virulence traits shared with plant and human fungal pathogens. Proc Natl Acad Sci USA 104: 18730–18735.

79. KämperJ, KahmannR, BölkerM, MaLJ, BrefortT, et al. (2006) Insights from the genome of the biotrophic fungal plant pathogen Ustilago maydis. Nature 444: 97–101.

80. ZuccaroA, LahrmannU, GüldenerG, LangenG, PfiffiS, et al. (2011) Endophytic life strategies decoded by genome and transcriptome analyses of the mutualistic root symbiont Piriformospora indica. PLoS Pathog 7: e1002290.

81. Kwon-ChungKJ, BennettJE, RhodesJC (1982) Taxonomic studies of Filobasidiella species and their anamorphs. Antonie van Leeuwenhoek 48: 25–38.

82. ChaskesS, TyndallRL (1975) Pigment production of Cryptococcus neoformans from para-diphenols and ortho-diphenols: Effect of nitrogen source. J Clin Microbiol 1: 509–514.

83. NielsenK, De ObaldiaAL, HeitmanJ (2007) Cryptococcus neoformans mates on pigeon guano: implications for the realized ecological niche and globalization. Eukaryot Cell 6: 949–959.

84. JaffeDB, ButlerJ, GnerreS, MauceliE, Lindblad-TohK, et al. (2003) Whole-genome sequence assembly for mammalian genomes: Arachne 2. Genome Res 13: 91–96.

85. KentWJ (2002) BLAT–the BLAST-like alignment tool. Genome Res 12: 656–664.

86. AltschulSF, MaddenTL, SchäfferAA, ZhangJ, ZhangZ, et al. (1997) Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 25: 3389–3402.

87. Ter-HovhannisyanV, LomsadzeA, ChernoffYO, BorodovskyM (2008) Gene prediction in novel fungal genomes using an ab initio algorithm with unsupervised training. Genome Res 18: 1979–1990.

88. StankeM, SchoffmannO, MorgensternB, WaackS (2006) Gene prediction in eukaryotes with a generalized hidden Markov model that uses hints from external sources. BMC Bioinformatics 7: 62.

89. KorfI, FlicekP, DuanD, BrentMR (2001) Integrating genomic homology into gene structure prediction. Bioinformatics 7: S140–S148.

90. BlancoE, ParraG, GuigoR (2007) Using geneid to identify genes. Current Protocols in Bioinformatics 4 doi:10.1002/0471250953.bi0403s18

91. ElsikCG, MackeyAJ, ReeseJT, MilshinaNV, RoosDS, et al. (2007) Creating a honey bee consensus gene set. Genome Biol 8: R13.

92. BirneyE, ClampM, DurbinR (2004) GeneWise and Genomewise. Genome Res 14: 988–995.

93. LiM, MaB, KismanD, TrompJ (2004) Patternhunter II: highly sensitive and fast homology search. J Bioinformatics Comput Biol 2: 417–439.

94. LagesenK, HallinPF, RødlandE, StærfeldtHH, RognesT, et al. (2007) RNammer: consistent annotation of rRNA genes in genomic sequences. Nucleic Acids Res 35: 3100–3108.

95. LoweTM, EddySR (1997) tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res 25: 955–964.

96. LiH, DurbinR (2009) Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25: 1754–1760.

97. DePristoMA, BanksE, PoplinR, GarimellaKV, MaguireJR, et al. (2011) A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet 43: 491–498.

98. McKennaA, HannaM, BanksE, SivachenkoA, CibulskisK, et al. (2010) The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res 20: 1297–1303.

99. ChenK, WallisJW, McLellanMD, LarsonDE, KalickiJM, et al. (2009) BreakDancer: an algorithm for high-resolution mapping of genomic structural variation. Nature Meth 6: 677–681.

100. AlbersCA, LunterG, MacArthurDG, McVeanG, OuwehandWH, et al. (2011) Dindel: accurate indel calls from short-read data. Genome Res 21: 961–973.

101. ParkhomchukD, BorodinaT, AmstislavskiyV, BanaruM, HallenL, et al. (2009) Transcriptome analysis by strand-specific sequencing of complementary DNA. Nucleic Acids Res 37: e123.

102. LevinJZ, YassourM, AdiconisX, NusbaumC, ThompsonDA, et al. (2010) Comprehensive comparative analysis of strand-specific RNA sequencing methods. Nature Meth 7: 709–715.

103. CuomoCA, DesjardinsCA, BakowskiMA, GoldbergJ, MaAT, et al. (2012) Microsporidian genome analysis reveals evolutionary strategies for obligate intracellular growth. Genome Res 22: 2478–2488.

104. GrabherrMG, HaasBJ, YassourM, LevinJZ, ThompsonDA, et al. (2011) Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol 29: 644–652.

105. HaasBJ, PapanicolaouA, YassourM, GrabherrM, BloodPD, et al. (2013) De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nature Protocols 8: 1494–1512.

106. LiB, DeweyCN (2011) RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. Bioinformatics 12: 323.

107. RobinsonMD, McCarthyDJ, SmythGK (2010) edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26: 139–140.

108. KadotaK, NishiyamaT, ShimizuK (2012) A normalization strategy for comparing tag count data. Algorithms Mol Biol 7: 1–13.

109. BenjaminiY, HochbergY (1995) Controlling the false discovery rate: a practical and powerfull approach to multiple testing. J R Statis Soc B 57: 289–300.

110. EddySR (2011) Accelerated profile HMM searches. PLoS Computl Biol 7: e1002195.

111. StoreyJD, TibshiraniR (2003) Statistical significance for genomewide studies. Proc Natl Acad Sci USA 100: 9440–9445.

112. LengelerKB, WangP, CoxGM, PerfectJR, HeitmanJ (2000) Identification of the MATa mating-type locus of Cryptococcus neoformans reveals a serotype A MATa strain thought to have been extinct. Proc Natl Acad Sci USA 97: 14455–14460.

113. SouthernE (2006) Southern blotting. Nature Protocols 1: 518–525.

114. GoodwinTJ, PoulterRT (2001) The diversity of retrotransposons in the yeast Cryptococcus neoformans. Yeast 18: 865–880.

115. PitkinJW, PanaccioneDG, WaltonJD (1996) A putative cyclic peptide efflux pump encoded by the TOXA gene of the plant-pathogenic fungus Cochliobolus carbonum. Microbiology 142: 1557–1565.

116. FraserJA, SubaranRL, NicholsCB, HeitmanJ (2003) Recapitulation of the sexual cycle of the primary fungal pathogen Cryptococcus neoformans var. gattii: implications for an outbreak on Vancouver island, Canada. Eukaryot Cell 2: 1036–1045.

117. ToffalettiDL, RudeTH, JohnstonSA, DurackDT, PerfectJR (1993) Gene transfer in Cryptococcus neoformans by use of biolistic delivery of DNA. J Bacteriol 175: 1405–1411.

118. SanyalK, BaumM, CarbonJ (2004) Centromeric DNA sequences in the pathogenic yeast Candida albicans are all different and unique. Proc Natl Acad Sci USA 101: 11374–11379.

119. DubinM, FuchsJ, GrafR, SchubertI, NellenW (2010) Dynamics of a novel centromeric histone variant CenH3 reveals the evolutionary ancestral timing of centromere biogenesis. Nucleic Acids Res 38: 7526–7537.

120. HubermanJA, SpotilaLD, NawotkaKA, El-AssouliSM, DavisLR (1987) The in vivo replication origin of the yeast 2 µm plasmid. Cell 51: 473–481.

121. TheisJF, NewlonCS (2001) Two compound replication origins in Saccharomyces cerevisiae contain redundant origin recognition complex binding sites. Mol Cell Biol 21: 2790–2801.

122. ChristensenWB (1946) Urea decomposition as a means of differentiating Proteus and paracolon cultures from each other and from Salmonella and Shigella types. J Bacteriol 52: 461–466.

123. AlspaughJA, PerfectJR, HeitmanJ (1997) Cryptococcus neoformans mating and virulence are regulated by the G-protein a subunit GPA1 and cAMP. Genes Dev 11: 3206–3217.

124. CoxGM, MukherjeeJ, ColeGT, CasadevallA, PerfectJR (2000) Urease as a virulence factor in experimental cryptococcosis. Infect Immun 68: 443–448.

125. MylonakisE, MorenoE, El KhouryJB, IdnurmA, HeitmanJ, et al. (2005) Galleria mellonella as a model system to study Cryptococcus neoformans pathogenesis. Infect Immun 73: 3842–3850.

126. AndersS, HuberW (2010) Differential expression analysis for sequence count data. Genome Biol 11: R106.

127. LiL, StoeckertCJJ, RoosDS (2003) OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res 13: 2178–2189.

128. EdgarRC (2004) MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res 19: 1792–1797.